variational methods ( vm ) have long been known as powerful tools in theoretical physics , especially in the context of reduced modeling . as opposed to approximating differential equations , approximating the action functional guarantees that the corresponding model too has a variational structure by definition .
this implies that conservation laws are also generally preserved , so a model inherits key features of the original system notwithstanding the reduction @xcite .
one of the areas that particularly benefit from this fact is wave theory @xcite , where even crude approximations to an action functional often yield insightful reduced models @xcite .
the advantages of vm become especially obvious when one deals with media such as plasmas , which can be nonstationary , inhomogeneous , anisotropic , and exhibit both temporal and spatial dispersion simultaneously @xcite .
first - principle differential equations are often unrealistic to handle in this case , whereas vm allow for simple and intuitive modeling .
specifically , a wave theory can be constructed as an axiomatic field theory without even specifying the waves nature ( electromagnetic , acoustic , etc . ) . defining the action functional
explicitly is needed only to connect the resulting general equations with specific quantities of interest @xcite .
for example , this provides a convenient alternative to using maxwell s equations for electromagnetic waves @xcite , which are generally complicated @xcite ( as opposed to the corresponding vm ) even in the geometrical - optics ( go ) limit  @xcite .
but which exact variational principle does one begin with ?
the microscopic least action principle ( lap ) is typically available in the form @xmath0 = 0 $ ] , where @xmath1 describes the wave field , @xmath2 describes the medium , and @xmath3 is a known functional .
however , practical applications in the context described above require a different variational principle , where @xmath1 is the _ only _ independent variable , while the medium s degrees of freedom are `` integrated out '' , i.e.,somehow expressed through @xmath1 .
if the medium response is adiabatic ( e.g.,if @xmath4 $ ] can be approximated with a local function ) , an effective lap for @xmath1 can be formulated simply as @xmath5 = 0 $ ] , as will also be discussed below .
this principle has enjoyed insightful applications , e.g.,in the theory of generalized ponderomotive forces @xcite and nonlinear plasma waves @xcite . however
, when the medium response is not adiabatic , the dynamics of @xmath1 becomes dissipative .
then , the mentioned effective lap becomes inapplicable , and the problem requires a more subtle approach .
this issue has been recently receiving attention in plasma physics @xcite and is also important for advancing variational formulations of plasma and wave kinetics such as in ref .
the purpose of the present paper is to find the appropriate modification of the effective lap that would incorporate dissipation under the above assumptions .
accommodating dissipative effects within the lap has a long history ( e.g.,cf .
@xcite ) both as a formal problem of finding a variational principle for a given equation @xcite and also in modeling physical systems , for example , quantizing radiation in dissipative media @xcite .
it is hardly possible to overview the topic here , so the given references are intended as representative but not exhaustive . in any case
, the subject of our paper is different from those discussed in literature , particularly , in the following aspects :    * we do not address the problem of finding a lap for a system with general dissipation .
instead , we consider only emergent dissipation in nonconservative subsystems of systems that are overall conservative .
* we avoid the standard approach ( known under many different names @xcite ) where additional independent functions are introduced that serve as lagrange multipliers @xcite . in our formulation , only time - independent lagrange multipliers are used .
* we avoid variable transformations that are used ( as in , e.g.,refs .
@xcite ) to compensate the loss of the phase space volume caused by dissipation . *
in contrast to theories discussed in quantum contexts , our theory is not statistical .
instead , we search for a lap that describes deterministic dynamics on a finite time interval for a given initial state of the medium .
thus , we can not describe wave fields using the standard fourier representation , and we also have to specify boundary conditions ( bc ) in time , which are often ignored .
in fact , it is the attention to these bc that helps accommodate dissipation within a lap . * instead of modeling media as collections of stationary oscillators ( as in , e.g.,refs .
@xcite ) , we allow for media that are nonstationary , inhomogeneous , anisotropic , and have both temporal and spatial dispersion simultaneously . moreover , our approach allows treating classical and quantum degrees of freedom on the same footing .
* our theory is formulated within the standard calculus of variations .
we avoid exotic constructs such as complex actions or fractional derivatives @xcite .
feynman integrals , which are used in related quantum theories ( e.g.,see ref .
@xcite ) , are also avoided .
our main results are as follows .
we formulate a variational principle for a dissipative subsystem of an overall - conservative system by introducing constant lagrange multipliers and lagrangians nonlocal in time .
we call it the variational principle for projected dynamics ( vppd ) . although the work is originally motivated by needs of the plasma wave theory , our results are applicable to general dynamical systems just as well . in the first part of the paper ( secs .  [
sec : notation]-[sec : routh ] ) , we introduce the general approach where @xmath4 $ ] is an arbitrary nonlinear functional and derive an exact effective lap for @xmath1 . in the second part ( sec .
[ sec : losc ] ) , we elaborate on applications of this approach to the important special case of linear @xmath4 $ ] , which corresponds to linear media .
our formulation is equally applicable to both classical and quantum media . in the third part ( secs .
[ sec : asym ] and [ sec : cont ] ) , we discuss various asymptotic approximations of the effective lap for oscillations and waves in linear media , including the quasistatic and go approximation . to our knowledge , this is the _ first variational formulation of dissipative go _ in a general linear medium , which is allowed to be nonstationary , inhomogeneous , nonisotropic , and exhibits both temporal and spatial dispersion simultaneously .
also notably , the `` spectral representation '' of the dielectric tensor @xmath6 that enters the go equations is shown to be , strictly speaking , the weyl symbol of the medium dielectric permittivity .
( previously , it has been known that weyl symbols emerge naturally in wave kinetics @xcite . )
this can be considered as an invariant definition of @xmath6 , as opposed to _
definitions used in literature that have been a source of a continuing debate @xcite .
our work is also intended as a stepping stone for extending the variational theory of modulational stability in general wave ensembles @xcite to dissipative systems .
details will be described in a separate publication .
the following notation is used throughout the paper .
the symbol @xmath7 denotes definitions , and @xmath8 denotes transposition . for example , for a column vector @xmath1 [ eq .  ( [ eq : q ] ) ] , one has @xmath9 , which is a row vector ; i.e.,in this case , @xmath8 merely lowers the index .
we assume @xmath10 for any @xmath11 , where @xmath12 is the kronecker delta .
( summation over repeated indices is assumed unless specified otherwise . ) hence , in principle , @xmath13 and @xmath14 could be considered equivalent , but we prefer to retain the distinction because some quantities ( such as coordinates ) are naturally defined with upper indexes , while others ( such as momenta ) are naturally defined with lower indexes . introducing a more fundamental geometrical interpretation will not be needed for our purposes .
note also that some indexes will be used _
just _ to introduce new symbols .
this applies to all scalars ( e.g.,@xmath15 and @xmath16 ) and to some vectors ( e.g.,@xmath17 , @xmath18 , @xmath19 , @xmath20 ) .
we also define @xmath21 .
the symbol @xmath22 will denote @xmath23 .
the notation `` @xmath24 '' will denote that the equation is obtained by extremizing a corresponding action functional with respect to @xmath11 .
the remaining notation is introduced within the main text .
finally , the abbreviations used in the text are summarized as follows : +    [ cols= " > , < " , ]      +
we start by introducing the general formalism of lagrangian mechanics . although this formalism is commonly found in textbooks ( e.g.,see ref .
@xcite ) , it is essential to restate it here , so we can explicitly reference specific definitions and equations later in the text .      consider a lagrangian system @xmath25 characterized by some real coordinate vector @xmath26 the dynamics of @xmath25 on a time interval @xmath27 is governed by the lap .
we call this interval @xmath28 , and we also use the symbol @xmath28 to denote its duration , @xmath29 . we assume the action functional @xmath30 to have the form @xmath31 = \int^{t_2}_{t_1 } l^{(q)}(t , q , \dot{q})\,dt,\end{gathered}\ ] ] where @xmath32 is called a lagrangian
. then , the standard formulation of the lap is as follows : @xmath33 = 0 , \\
\delta q(t_{1,2 } ) = 0.\label{eq : dqzero2}\end{gathered}\ ] ] [ here and further , the notation @xmath34 means `` both @xmath35 and @xmath36 '' ] .
note that eq .
( [ eq : dqzero2 ] ) implies bc @xmath37 where @xmath38 and @xmath39 are some given constants .    from eq .
( [ eq : aq ] ) , one obtains @xmath40\,\delta q\,dt,\end{gathered}\ ] ] where we introduced the `` canonical momentum '' @xmath41 with @xmath42 , or , more compactly , @xmath43 . the first term on the right in eq .
( [ eq : aux1 ] ) vanishes because @xmath34 are fixed .
hence , the lap leads to the following euler - lagrange equations ( ele ) : @xmath44 this is a second - order equation for the @xmath45-dimensional coordinate @xmath1 , so the @xmath46 bc ( [ eq : bcq ] ) ( two per each component of @xmath1 ) provide just the right number of parameters to define a solution of eq .
( [ eq : ele ] ) .
alternatively , one can also approach eq .
( [ eq : ele ] ) as an initial - value problem ; then , the condition on on @xmath47 would need to be replaced with a condition on @xmath48 .
it is straightforward to extend this formulation to lagrangians that depend also on higher - order derivatives , but we will not consider such extensions for the sake of brevity
.      now suppose that @xmath25 is coupled to another lagrangian system @xmath49 , which we call the `` medium '' .
suppose also that @xmath49 is characterized by some @xmath50-dimensional coordinate @xmath2 .
the dynamics of the resulting system @xmath51 is governed by the lap @xmath0 = 0 $ ] , where @xmath3 is the total action . for clarity , we adopt it in the form @xmath52 , where @xmath53 is the lagrangian given by @xmath54 the assumed bc are @xmath55 which introduce @xmath56 bc similar to eq .
( [ eq : bcq ] ) , namely , @xmath57 like in sec .
[ sec : single ] , one then obtains @xmath58   = ( p \cdot \delta q)\,\big|_{t_1}^{t_2 } + ( \eta \cdot \delta \xi ) \,\big|_{t_1}^{t_2}\\ + \int^{t_2}_{t_1 } \big [ ( { \partial}_q l^{(q ) } - \dot{p})\,\delta q +   ( { \partial}_\xi l^{(\xi ) } - \dot{\eta})\,\delta \xi\ , \big]\,dt,\end{gathered}\ ] ] where we introduced @xmath59 the first two terms on the right side of eq .
( [ eq : auxa2 ] ) vanish due to eqs .
( [ eq : dqzero ] ) , so the lap yields the following ele : @xmath60 equations ( [ eq : ele21 ] ) and ( [ eq : ele22 ] ) are second - order equations for @xmath61 and @xmath62 , so the mentioned bc provide just enough parameters to define a solution .
although eqs .
( [ eq : ele21 ] ) and ( [ eq : ele22 ] ) form a self - consistent system , it is also convenient to approach eq .
( [ eq : ele22 ] ) formally as if @xmath1 were as a prescribed function . in that case , one can , in principle , solve for @xmath2 and express it as some functional @xmath63 $ ] , provided that one is given a @xmath64-dimensional integration constant @xmath65 to determine a solution unambiguously .
the resulting system is characterized by @xmath1 alone and is not conservative .
our goal is to derive an effective variational principle for this subsystem , or , in other words , for the dynamics of @xmath51 `` in projection '' on @xmath25 .
first , consider the variation of @xmath3 evaluated on @xmath66 $ ] : @xmath67 = ( p \cdot \delta q)\ , \big|_{t_1}^{t_2 } + ( \hat{\eta } \cdot \delta \hat{\xi } ) \,\big|^{t_2}_{t_1 } + \delta { \mathscr{f}}^{(q ) } + \delta { \mathscr{f}}^{(\xi)},\notag\\ \delta { \mathscr{f}}^{(q ) } \doteq \int^{t_2}_{t_1 } f^{(q)}(t , q , \dot{q } , \hat{\xi } , \dot{\hat{\xi}})\,\delta q\,dt,\notag\\ \delta { \mathscr{f}}^{(\xi ) } \doteq \int^{t_2}_{t_1 } f^{(\xi)}(t , q , \dot{q } , \hat{\xi } , \dot{\hat{\xi}})\,\delta \xi\,dt.\notag\end{gathered}\ ] ] by definition of @xmath68 , one has @xmath69 .
we will also adopt @xmath70 , as usual .
then , @xmath67 = ( \hat{\eta } \cdot \delta \hat{\xi})\,\big|^{t_2}_{t_1 } + \delta { \mathscr{f}}^{(q)}.\label{eq : a10}\end{gathered}\ ] ] if @xmath65 is a local function of @xmath71 , then the first term on the right of eq .
( [ eq : a10 ] ) vanishes . in this case , by requiring that @xmath3 is extremized , one arrives at the correct equation ( [ eq : ele21 ] ) . however , such choice of @xmath65 is inconvenient .
it is more practical to choose @xmath65 such that it characterizes the _ initial _ state of the medium , namely , @xmath72 [ for example , if @xmath1 is electric field , and @xmath2 is linear polarization , then @xmath73 , where @xmath74 describes free oscillations independent of @xmath1 , and @xmath75 is some kernel ( secs .  [
sec : losc]-[sec : cont ] ) . ] in this case , @xmath76 is zero , but @xmath77 is not even though @xmath34 are fixed .
hence , @xmath67 = ( \hat{\eta } \cdot \delta \hat{\xi})\big|^{t_2 } + \delta { \mathscr{f}}^{(q)},\label{eq : a1}\end{gathered}\ ] ] where the first term on the right generally does not vanish .
[ an exception is the special case when @xmath78 is a local function of @xmath1 ; then , fixing @xmath36 implies fixing @xmath79 too . ]
this means that satisfying eq .
( [ eq : ele21 ] ) at fixed @xmath34 is _ not _ sufficient to extremize @xmath80}$ ] . in order to ensure that the action is extremized at fixed @xmath34
, we must constrain the value of @xmath81 separately .    the problem of extremizing @xmath3 under an additional constraint @xmath82
is a standard isoperimetric problem that is solved by introducing a lagrange multiplier ( * ? ? ?
ii.14 ) . in our case , it is an @xmath50-dimensional multiplier , which we will call @xmath83 .
[ the fact that @xmath83 itself needs to be found increases the number of independent variables in the system , which is why @xmath84 can be constrained independently from @xmath34 . ] specifically , the dynamics `` projected '' on @xmath25 is governed by the `` effective '' action @xmath85 \doteq a[q , \hat{\xi } ] + \lambda \cdot \hat{\xi}(t_2),\end{gathered}\ ] ] and the corresponding variational principle , vppd ,  is @xmath86 = 0 , \quad \delta q(t_{1,2 } ) = 0.\end{gathered}\ ] ]    the value of @xmath83 can be readily found in the general case as follows . under the bc ( [ eq : bceff ] ) ,
one has @xmath87 = \delta { \mathscr{f}}^{(q ) } + ( \eta_t + \lambda ) \cdot \delta \hat{\xi}(t_2).\end{gathered}\ ] ] here , @xmath88 is the value of the functional @xmath89 corresponding to the physical trajectory for given @xmath17 , @xmath18 , @xmath19 , and @xmath90 .
the fact that the assumed fixed value of @xmath84 must be consistent also with the assumed @xmath36 implies @xmath91 .
hence , to satisfy eq .
( [ eq : bceff ] ) , we adopt @xmath92 this ensures that @xmath93 = 0 $ ] , so the variational principle ( [ eq : bceff ] ) indeed leads to the correct equation  ( [ eq : ele21 ] ) .    the vppd ( [ eq : bceff ] ) , combined with the solution for @xmath83 given by eq .
( [ eq : aux71 ] ) , is the main conceptual result of this paper . in what follows ,
we discuss various modifications of this principle and its applications to specific types of problems of practical interest .
first , consider the following alternative formulation . we start by introducing the hamiltonian of @xmath94 , @xmath95 as a function of @xmath96 . at fixed @xmath97 ,
one has @xmath98 where we substituted eqs .
( [ eq : petadef ] ) and ( [ eq : ele22 ] ) .
this leads to @xmath99 which are known as hamilton s equations .
they are equivalent to the combination of eqs .
( [ eq : petadef ] ) and ( [ eq : ele22 ] ) , as one can also recheck by direct calculation using eq .
( [ eq : hxi ] ) .
consider rewriting @xmath3 in terms of these new variables : @xmath100\,dt,\end{gathered}\ ] ] where @xmath101 can be recognized as a routhian ( * ? ? ?
41 ) . then , @xmath102\,dt.\notag\end{gathered}\ ] ] let us assume that @xmath34 are kept fixed , as in sec .  [ sec : coupled ] .
then , the requirement @xmath103 leads to correct equations [ eqs .
( [ eq : ele21 ] ) and ( [ eq : he ] ) ] if @xmath104 are treated as @xmath105 independent variables . for this reason
, the lap for our coupled systems can be reformulated as @xmath106 = 0\end{gathered}\ ] ] with bc still given by eqs .
( [ eq : dqzero ] ) .
( note that , although @xmath107 is now an independent variable too , it yet does _ not _ need to be fixed at the end points .
otherwise , the problem would be overdefined . )
accordingly , @xmath108 is defined  as @xmath109 \doteq a[q , \hat{\xi } , \hat{\eta } ] + \lambda \cdot \hat{\xi}(t_2),\end{gathered}\ ] ] where @xmath83 is a constant lagrange multiplier . like
earlier , one can show that the solution for @xmath83 is given by eq .
( [ eq : aux71 ] ) .
let us also introduce a `` symmetrized '' representation @xcite of the above equations that we will need in sec .
[ sec : losc ] .
notice that the action ( [ eq : apq ] ) can be rewritten as @xmath110 = \frac{1}{2}\,\eta \cdot \xi\,\big|_{t_1}^{t_2 } + { \mathscr{a}}[q , \xi , \eta],\\ { \mathscr{a}}[q , \xi , \eta ] \doteq \int^{t_2}_{t_1 } \big[\lambda^{(q ) }   + \frac{1}{2}\,(\eta \cdot \dot{\xi } - \dot{\eta } \cdot \xi ) - h^{(\xi)}\big]\,dt.\label{eq : aux33}\end{gathered}\ ] ] using @xmath111 and eq .
( [ eq : mla ] ) , one  gets @xmath112\big|_{t_1}^{t_2 } + \delta { \mathscr{f}}^{(q ) } + \delta { \mathscr{f}}^{(\xi)}.\label{eq : aux34}\end{gathered}\ ] ] as usual , we require @xmath70 . then , we have @xmath50 bc per end point , and the first term in the square brackets vanishes .
the remaining bc can be defined , e.g. ,  as @xmath113 ( summation over indexes is _ not _ assumed here ; also , the symbol @xmath114 denotes that the equality must be satisfied independently at @xmath115 and @xmath116 . ) . with obvious reservations , eq .
( [ eq : aux102 ] ) is also equivalent to @xmath117 this provides precisely @xmath50 bc per each end point , as needed , and eliminates the second term in the square brackets in eq .
( [ eq : aux34 ] ) .
hence , the requirement of extremal @xmath118 leads to correct equations [ eqs .
( [ eq : ele21 ] ) and ( [ eq : he ] ) ] .
therefore , the lap can be reformulated as @xmath119 = 0\end{gathered}\ ] ] with bc given by eqs .
( [ eq : dqzero2 ] ) and ( [ eq : aux102 ] ) . since the phase space variables @xmath2 and @xmath107 enter @xmath118 on the same footing , we call @xmath118 a phase - space action .    for @xmath118 , the analog of the constrained variational principle that we introduced in sec .
[ sec : constrained ] is @xmath120   = 0 , \quad \delta q(t_{1,2 } ) = 0,\end{gathered}\ ] ] where the action can be taken in the form @xmath121
\doteq   { \mathscr{a}}[q , \hat{\xi } , \hat{\eta } ]   + \sum_m \lambda_m ( \hat{\xi}^m/\hat{\eta}_m)\big|^{t_2}.\end{gathered}\ ] ] it is easy to see that this vppd leads to correct equations for @xmath1 provided that @xmath122 satisfy @xmath123 where @xmath124 is the value of @xmath125 corresponding to the physical trajectory for given @xmath17 , @xmath18 , @xmath19 , and @xmath90 .
let us also introduce an alternative form of @xmath126 $ ] that we will need further on . instead of fixing @xmath127 , one can fix @xmath128 and @xmath79 independently at the expense of doubling the number of unknown lagrange multipliers .
specifically , one can introduce two @xmath50-dimensional lagrange multipliers , @xmath129 and @xmath130 , and define the effective action as follows : @xmath131 \doteq   { \mathscr{a}}[q , \hat{\xi } , \hat{\eta } ]   + \hat{\eta}(t_2 ) \cdot \lambda^{(\eta ) } + \lambda^{(\xi ) } \cdot \hat{\xi}(t_2).\end{gathered}\ ] ] then , one can show that eq .
( [ eq : aux101 ] ) yields correct equations for @xmath1 provided that @xmath132 where @xmath133 and @xmath134 are the values of @xmath125 and @xmath135 corresponding to the physical trajectory for given @xmath17 , @xmath18 , @xmath19 , and @xmath90 .
these equations can also be presented in a more compact form , as discussed in appendix  [ app : psp ] .
let us now consider , as an example of practical interest , the case when @xmath49 is an ensemble of linear oscillators . in this case
, the above phase - space representation facilitates the formulation of a general theory of linear dispersion in a convenient `` quantumlike '' form that is introduced below .
note that truly quantum equations are also subsumed under this formulation . as a side note ,
the special case of _ adiabatic _ oscillators can be described by a simpler machinery discussed in appendix  [ app : osc ] .
when @xmath49 is an ensemble of linear oscillators , there exists @xcite a linear transformation of the variables @xmath136 to some variables @xmath137 ( roman indexes are used to distinguish them from standard italic indexes that denote vector components ) such that @xmath138 = \int_{t_1}^{t_2 } \lambda^{(\xi)}\,dt}$ ] , @xmath139 where @xmath140 and @xmath141 .
( here , we assume that the underlying oscillators have positive energies for simplicity .
negative - energy oscillations @xcite , which are somewhat exotic and typically unstable yet not impossible , could be included by generalizing the notation as explained in ref .
the matrix @xmath142 is hermitian and serves as a hamiltonian . also , @xmath143 is a bilinear form of @xmath144 , and @xmath145 is the complex - conjugate bilinear form of @xmath146 .
the terms @xmath143 and @xmath145 are determined by the rate at which parameters of the medium evolve in time and vanish when those parameters are time - independent .
also  @xcite , @xmath147 so the bc ( [ eq : aux102 ] ) can be expressed as @xmath148\big|^{t_{1,2}},\end{gathered}\ ] ] or , equivalently , as @xmath149_{t_{1,2 } } = 0 $ ] @xcite .
the corresponding ele are as follows : @xmath150 yet @xmath151 and @xmath152 can be expressed as linear combinations of @xmath144 and @xmath146 .
hence , one can also rewrite eqs .
( [ eq : aux64 ] ) in an equivalent complex form @xmath153 where @xmath154 is considered a functional of @xmath144 and @xmath146 . in this sense , @xmath144 and @xmath146
can be formally treated as independent variables .
then , it is easy to see that , under the bc ( [ eq : bcpsi ] ) , one has @xmath155\,dt.\end{gathered}\ ] ] accordingly , the ele can be cast as follows : @xmath156 an additional equation is obtained for @xmath146 , but it is simply the adjoint of eq .
( [ eq : aux65 ] ) .
( this is understood because we have only two independent real functions , @xmath151 and @xmath152 , so there can be only one independent complex equation . )
thus , the equation for @xmath146 will be omitted for brevity .
the terms @xmath143 and @xmath145 describe effects caused by parametric resonances .
( for example , if @xmath144 were the complex probability amplitude of a free quantum particle , they would describe particle annihilation and production . )
we ignore such effects for simplicity , so @xmath143 and @xmath145 are henceforth omitted .
then , the lagrangian of @xmath49 can be represented simply as follows : @xmath157 this leads to a quantumlike equation for @xmath144 , @xmath158 which conserves @xmath159 even though @xmath142 may depend on time .
the quantity @xmath159 is understood as the total number of `` quanta '' , or `` particles '' , in @xmath160 @xcite . [ another name of this quantity is the total action of oscillations in @xmath160 .
we prefer not to use this term here to avoid confusion , since @xmath161 is also called the action of @xmath160 . ]      now , let us introduce the interaction between @xmath49 and @xmath25
. we will describe it by an additional hamiltonian @xmath162 . assuming that the coupling is sufficiently weak
, one can model @xmath162 with a function linear in @xmath2 and @xmath107 .
since @xmath162 is also real , its representation through @xmath144 and @xmath146 must then have the form @xmath163 where @xmath164 is some complex vector and the minus sign is added for convenience .
we will assume @xmath164 to depend on @xmath165 for clarity , but the dependence on derivatives of @xmath1 could also be included . as a side note , this particular model does not conserve @xmath159 . a conservative model , which describes parametric interactions and quantum interactions in particular , is discussed in appendix  [ app : q ] .
since @xmath142 in this model can depend only on time , it is convenient to eliminate it using a variable transformation .
specifically , consider a new variable @xmath166 defined via @xmath167 .
let us choose the operator @xmath168 ( a `` propagator '' ) such that @xmath169 and require that @xmath168 be a unit matrix at @xmath115 .
this ensures that @xmath168 remains unitary at all times , so it can be cast as @xmath170 , where @xmath171 is some hermitian operator .
more specifically , the propagator can be expressed in terms of the following time - ordered exponential @xmath172 : @xmath173.\end{gathered}\ ] ] also importantly , for any @xmath174 , @xmath175 , and @xmath176 , one has @xmath177    in the new variables , the symmetrized action has the form @xmath178 = \int_{t_1}^{t_2 } l\,dt$ ] , where @xmath179 and @xmath180 . the corresponding lap is formulated  as @xmath181 = 0\end{gathered}\ ] ] with bc given by eq .
( [ eq : dqzero2 ] ) and @xmath182\big|^{t_{1,2 } } = 0.\end{gathered}\ ] ] it is easy to see that the corresponding ele are @xmath183      now suppose that one has been able to express @xmath166 through @xmath1 .
we denote the resulting functional @xmath184 and assume that it is parameterized by some @xmath50 complex constants ( i.e.,@xmath64 real constants ) , say , @xmath185 .
hence , we can introduce a vppd for @xmath25 as in sec .
[ sec : sym ] with @xmath186 = \int^{t_2}_{t_1}\big\ { \lambda^{(q ) } ( t , q , \dot{q } )   \\
+ \frac{1}{2}\,[\gamma^\dag(t , q ) \cdot \hat{\psi } + \hat{\psi}^\dag \cdot \gamma(t , q ) ] \big\}\,dt + { \mathcal{b}},\end{gathered}\ ] ] @xmath187.\label{eq : aux334}\end{gathered}\ ] ] here , @xmath188 is a complex lagrange multiplier introduced such that the action remains real .
this representation implies that we seek to fix the values of @xmath189 and @xmath190 independently from @xmath36 by introducing just the right number of lagrange multipliers ( namely , the @xmath64 components of @xmath188 , where the real and imaginary parts are counted separately )
. the values of these multipliers are to be found from the condition that the values of @xmath189 and @xmath190 must be consistent with the value of @xmath36 .
it is easily seen that choosing @xmath191 where @xmath192 is the value of @xmath193 corresponding to the physical trajectory for given @xmath194 , leads to a correct equation for @xmath1 [ eq .
( [ eq : qfq ] ) ] .
hence , @xmath195 is indeed the sought effective action for the projected dynamics of @xmath25 .        by integrating eq .
( [ eq : psiq ] ) , one can express @xmath184 as follows : @xmath196 here , @xmath197 is a linear integral operator , @xmath198\,dt',\end{gathered}\ ] ] which can be understood as the polarizability of @xmath49 with respect to @xmath199 .
( the symbol @xmath200  means `` applied to the expression on the right '' . )
hence , the effective lagrangian can be cast as follows : @xmath201 = \int^{t_2}_{t_1 } \lambda^{(q)}_{0 } \,dt + \xi + { \mathcal{b}},\\ \xi \doteq \text{re } \int^{t_2}_{t_1 } \gamma^\dag \cdot ( \alpha_\gamma \circ \gamma)\,dt.\end{gathered}\ ] ] the difference between @xmath202 and @xmath203 produces an effect similar to that of an external force , while @xmath204 can be described as follows .
suppose a `` small - amplitude approximation '' , namely , that the coordinate @xmath1 is chosen such that @xmath205 , where @xmath11 is a small parameter . assuming that the interaction is weak
, a typical interaction hamiltonian will then be linear in @xmath1 , so we assume @xmath206 , where @xmath207 is an @xmath208 complex matrix that may depend only on @xmath176 , if at all .
( the coefficient @xmath209 is introduced to simplify the interpretation of the final result . )
this gives @xmath210 where @xmath211 is another @xmath208 complex matrix , namely , @xmath212 ( so that , while @xmath213 is typically slow , @xmath211 is oscillatory ) .
then , @xmath214 and @xmath204 can be written as follows : @xmath215 where @xmath216 is an integral operator defined as @xmath217 $ ] .
it is convenient to express this operator in terms of a response ( or green s ) function @xmath216 as @xmath218 here , @xmath219 is the heaviside step function , and @xmath216 is a real matrix function given by @xmath220 \notag \\ & = - \text{im } \left[v^\dag(t ) \cdot v(t')\right ] \notag \\ & = ( i/2)\,\big[v^\dag(t ) \cdot v(t ' ) -
v^{\rm t}(t ) \cdot v^*(t')\big ] , \label{eq : akern}\end{aligned}\ ] ] which satisfies @xmath221 note that `` re '' and `` i m '' above actually mean real and imaginary parts , _ not _ hermitian and anti - hermitian parts , and , in general , @xmath216 is _ not _ hermitian .
as a side remark , the polarizability can be understood as a pseudodifferential operator expressible through @xmath176 and @xmath222 .
the fact that the polarizability operator is not necessarily hermitian can be related to the fact that @xmath222 is not hermitian on a finite time interval , as opposed to the whole time axis . a related discussion for spatial ( rather than temporal ) dispersion can be found in refs .
@xcite .
also note that , using eq .
( [ eq : alphaq ] ) , one can cast @xmath204  as @xmath223 or , equivalently , as @xmath224 where we used eq .
( [ eq : attpsym ] ) and @xmath225 . these formulas will be used below .
if @xmath1 is understood as a vector , then @xmath216 serves as a tensor of rank @xmath226 .
hence , it is convenient to introduce also the related covariant tensor of rank @xmath227 , i.e.,``lower the index '' of the matrix @xmath228 .
we denote such index lowering by underlining . using the corresponding matrix @xmath229 , one can write @xmath230 for any @xmath11 and @xmath231 . accordingly
, we can lower the index also in @xmath216 , i.e.,define @xmath232 this allows one to write @xmath233 ( because transposing a vector is the same as lowering its index ) , and also @xmath204 can be expressed without involving @xmath234 : @xmath235 if one does not need to distinguish upper and lower indexes , the distinction between @xmath216 and @xmath236 can be ignored .
it is instructive to explain at this point how the ele ( [ eq : qfq ] ) that we anticipate from the general theory can be derived also by a brute - force calculation using @xmath238 = \int_{t_1}^{t_2 } \left[{\partial}_q \lambda^{(q ) } - \dot{p }   + \frac{\text{re}\,(v^\dag \cdot \psi_0)}{\sqrt{2}}\right ] \cdot \delta q\,dt \\
+ \delta \xi + \delta{\mathcal{b}}.\end{gathered}\ ] ]    first , using eq .
( [ eq : xittp ] ) and the fact that @xmath239 for any scalar @xmath11 , one gets    @xmath240\,dt\,dt ' \notag\\ & = \frac{1}{2 } \int_{t_1}^{t_2 } \int^{t_2}_{t_1 } \left [ { q^{\rm t}}(t ' ) \cdot { \alpha^{\rm t}}(t , t')\,\theta(t - t ' ) \cdot \delta q(t ) + { q^{\rm t}}(t ) \cdot \alpha(t , t')\,\theta(t - t ' ) \cdot \delta q(t ' ) \right]\,dt\,dt ' \notag\\ & = \frac{1}{2 } \int_{t_1}^{t_2 } \int^{t_2}_{t_1 } \left [ { q^{\rm t}}(t ' ) \cdot { \alpha^{\rm t}}(t , t')\,\theta(t - t ' ) \cdot \delta q(t ) + { q^{\rm t}}(t ' ) \cdot \alpha(t ' , t)\,\theta(t ' - t ) \cdot \delta
q(t ) \right]\,dt\,dt ' \notag\\ & = \frac{1}{2 } \int_{t_1}^{t_2 }   \left [ \int_{t_1}^{t } { q^{\rm t}}(t ' ) \cdot { \alpha^{\rm t}}(t , t')\,dt ' + \int_{t}^{t_2 } { q^{\rm t}}(t ' ) \cdot \alpha(t ' , t)\,dt ' \right ] \cdot \delta q(t ) \ , dt \notag\\ & \equiv \int_{t_1}^{t_2 } \xi^{(q){\rm t}}(t ) \cdot \delta q(t)\,dt.\label{eq : xiq}\end{aligned}\ ] ]    using eq .
( [ eq : attpsym ] ) , we can also rewrite @xmath241 as follows : @xmath242    next , one gets from eq .
( [ eq : aux334 ] ) that @xmath243,\end{gathered}\ ] ] where @xmath244 can be obtained from eq .
( [ eq : fnpsi ] ) ; namely , @xmath245 let us substitute eq .
( [ eq : mupsi ] ) and again use eq .
( [ eq : fnpsi ] ) , now to express @xmath192 .
this gives @xmath246 \\
\cdot \frac{1}{\sqrt{2 } }   \int^{t_2}_{t_1 } v(t)\cdot \delta q(t)\,dt\bigg\},\end{gathered}\ ] ] or , equivalently , @xmath247}{\sqrt{2 } } + { \mathcal{b}}^{(q){\rm t}}\bigg\}\cdot \delta q(t)\,dt,\label{eq : dbq}\\ { \mathcal{b}}^{(q ) } \doteq \frac{1}{2}\int^{t_2}_{t_1 } \alpha(t , t ' ) \cdot q(t')\,dt',\end{gathered}\ ] ] and we also notice that @xmath248    by combining eqs .
( [ eq : xiq ] ) and ( [ eq : dbq ] ) , we obtain @xmath249}{\sqrt{2}}\\   + { [ \xi^{(q)}(t ) + { \mathcal{b}}^{(q)}(t)]^{\rm t}}\bigg\ } \cdot \delta q(t)\,dt,\end{gathered}\ ] ] and eqs .
( [ eq : xiq ] ) and ( [ eq : biq ] ) give @xmath250 .
thus , @xmath251}{\sqrt{2 } }   + ( \underline{\alpha } \circ q)(t ) \bigg\}\cdot \delta q(t)\,dt,\end{gathered}\ ] ] so eq .
( [ eq : aeffxib ] ) becomes @xmath252 = \int_{t_1}^{t_2 } \big[{\partial}_q \lambda^{(q ) } - \dot{p }   + \underline{\alpha } \circ q   \\
+ \frac{v^\dag \cdot \psi_0}{\sqrt{2 } } + \frac{\psi_0^\dag \cdot v}{\sqrt{2 } } \big ]
\cdot \delta q\,dt.\end{gathered}\ ] ] then , the requirement @xmath253 = 0 $ ] leads to @xmath254 considering the assumed expression for @xmath199 [ eq .
( [ eq : gamma ] ) ] , this agrees with the general ele [ eq .  ( [ eq : qfq ] ) ] , as expected .
in this section , we make approximations to the vppd and derive reduced theories for certain limiting cases .
although the calculations below may seem cumbersome , producing equivalent results by a brute - force reduction of eq .
( [ eq : pg ] ) ( as opposed to the underlying functional ) would still have been harder to do in the general case .
consider a representation of the response function in terms of the `` symmetrized '' time coordinates @xmath255 specifically , we define @xmath256 which has the following important property : @xmath257    using eqs .
( [ eq : akern])-([eq : gdeftau ] ) along with eq .
( [ eq : vdef ] ) , one can write @xmath258,\label{eq : avu}\\ { \mathcal{u}}(\bar{t } , \tau ) \doteq u(\bar{t } + \tau/2 , t_1 ) \cdot u^\dag(\bar{t } - \tau/2 , t_1).\end{gathered}\ ] ] since @xmath168 is unitary , one has @xmath259 . using this along with eq .
( [ eq : mcu2 ] ) , we can rewrite @xmath260 as follows : @xmath261 or , in other words , @xmath262,\end{gathered}\ ] ] where @xmath263 is a time - ordered exponential .
first , consider a stationary medium , i.e.,a medium where @xmath213 and @xmath142 are constant .
then , eq .  ( [ eq : mcu3 ] ) gives @xmath264 , and eq .
( [ eq : avu ] ) gives @xmath265 for simplicity , assume the basis in which @xmath142 is diagonal ( `` energy basis '' ) ; i.e.,@xmath266 .
then , matrix elements of @xmath267 can be written as @xmath268 provided that the set of @xmath269 is dense enough , the sum over  @xmath270 can be replaced with an integral [ @xmath271 : @xmath272 where @xmath273 $ ] . since @xmath274 is thereby a fourier image of @xmath275 ( or vice versa ) , it approaches zero at @xmath276 , at least if @xmath275 is continuous . if the characteristic gap @xmath277 between neighboring eigenfrequencies is nonzero ( and @xmath278 ) , the same conclusion applies at @xmath279 , when the spectrum discreteness is inessential .
thus , as long as larger @xmath280 are not of interest ( as will be assumed below ) , there is no need to distinguish continuous and discrete spectra .    the effect @xmath281 is known as phase mixing .
although introduced here for a stationary medium , it extends also to media evolving at large enough time scales @xmath282 , namely , @xmath283 where @xmath284 is the characteristic eigenvalue of @xmath142 .
the former requirement ensures that @xmath213 and @xmath142 be approximately constant on the time scale @xmath285 .
the latter requirement allows approximating eq .
( [ eq : mcu3 ] ) much like in the case of a stationary medium : @xmath286.\end{gathered}\ ] ] our intention is to discuss approximations to the vppd that become possible in this limit . in order to do that ,
let us introduce some definitions first .
let us introduce an auxiliary matrix function @xmath287 where @xmath288 is assumed real , @xmath267 is given by eq .
( [ eq : gdeftau ] ) , and @xmath289 is some constant .
( in later sections , we will treat @xmath289 as a function of @xmath290 , but this distinction is not important here . )
let us also introduce the hermitian and anti - hermitian parts of @xmath291 , namely , @xmath292 such that @xmath293 .
( in this notation , both @xmath294 and @xmath295 are hermitian , so @xmath296 and @xmath297 are anti - hermitian . )
as shown in appendix  [ app : alphas ] , @xmath298    provided efficient phase mixing , the function @xmath291 has a well - defined limit at large @xmath289 .
we call this limit @xmath299 ; i.e.,@xmath300 the function @xmath299 is recognized as the laplace image of @xmath267 evaluated at real @xmath288 and also , more generally , as the _ weyl image of the polarizability operator _ @xcite . like for @xmath291 , the hermitian and anti - hermitian parts of @xmath301 satisfy @xmath302 approximate formulas for @xmath291 and @xmath301 in the limit ( [ eq : vro ] ) can be found in appendix  [ app : alphas ] .
it is also shown there that @xmath303 and @xmath304 satisfy the kramers - kronig relations .
first , let us consider the case when @xmath1 and , possibly , the medium evolve on time scales @xmath282 that satisfy @xmath305 ( to simplify the notation , we will assume all the three small parameters to be comparable to each other , but this is not essential . ) below , we derive an asymptotic expression for @xmath195 using @xmath306 as a small parameter .          to approximate the term @xmath204 as given by eq .
( [ eq : xisgn ] ) , let us first map the integration variables @xmath307 in the double integral to @xmath308 using eqs .
( [ eq : ttau ] ) .
the integration domain is then mapped as shown in figs .
[ fig : domain](a ) and ( b ) , namely , to @xmath309 , where @xmath310 is a piecewise - linear function .
this gives @xmath311 where tr stands for `` trace '' , and @xmath312 is a matrix given by @xmath313 to the lowest ( zeroth ) order in @xmath306 , one has @xmath314 the function @xmath315 [ given by eq .
( [ eq : aproph ] ) at @xmath316 can be approximated with @xmath317 everywhere except for an @xmath318 part of the integration domain . neglecting this @xmath318 correction
, we get @xmath319 \cdot q(\bar{t})\,d\bar{t},\end{gathered}\ ] ] which corresponds to the approximation illustrated in fig .  [ fig : domain](c ) .
we have omitted the index @xmath320 in @xmath317 because , in the zero - frequency limit , the polarizability s weyl symbol ( [ eq : aw ] ) becomes hermitian ( appendix  [ app : alphas ] ) ; i.e.,@xmath321 .
specifically , @xmath322 as seen from eq .
( [ eq : zerow ] ) in the limit when the oscillating term is eliminated by phase mixing .
using eq .
( [ eq : fnpsi ] ) , let us express @xmath184 as @xmath324 the function @xmath171 is such that , at zero @xmath306 , it has a constant derivative equal to @xmath142 .
hence , at small nonvanishing @xmath306 , one can expect @xmath325 to be slow and close to @xmath142 .
then , the integrand in eq .
( [ eq : intpsi0 ] ) is a rapidly oscillating exponent times a slow function @xmath326 .
such integral is entirely determined by values of the integrand at the ends of the integration domain and is of order @xmath327 @xcite .
this gives the following estimate for the term @xmath323 in eq .
( [ eq : aux334 ] ) : @xmath328 assuming the effect of @xmath194 vanishes due to phase mixing . hence , @xmath329 .      by combining the above results and neglecting corrections @xmath318 and smaller , one obtains @xmath121 = \int^{t_2}_{t_1 } [ \lambda^{(q)}_{0 } ( t , q , \dot{q } ) - \phi_0(t , q)]\ , dt,\\ \phi_0(t , q ) \doteq - \frac{1}{2}\,[\underline{\alpha}_0(t ) \cdot q ] \cdot q.\label{eq : udef}\end{gathered}\ ] ] clearly , the corresponding ele is as follows : @xmath330 these results indicate that the interaction of a quasistatic system with a linear medium @xmath49 is equivalent to the interaction with an effective potential @xmath331 given by eq .
( [ eq : udef ] ) .
the symmetric matrix @xmath332 represents the covariant form of  @xmath333 given by eq .
( [ eq : alpha0est ] ) and serves as the polarizability of @xmath49 with respect to @xmath1 . also note that , for positive - definite hamiltonian @xmath142 , the polarizability @xmath333 is positive - definite too , and thus @xmath334 .      here
, we derive an asymptotic variational principle for the case when @xmath1 oscillates rapidly compared to the evolution of the medium .
the method is akin to whitham s averaging @xcite , but the averaging procedure is justified differently .
specifically , this is done as follows .
suppose @xmath203 is bilinear in @xmath335 .
then , eq .  ( [ eq : pg ] ) is linear and allows complex solutions .
let us consider a complex solution @xmath336 of the form @xmath337 , where @xmath338 is a complex vector and @xmath339 is a real phase .
provided that the time scale @xmath282 at which the medium evolves is large enough , in principle it is possible to construct asymptotic @xmath338 and @xmath339 such that @xmath340 is slow compared to @xmath339 at all @xmath176 ; i.e.,@xmath1 is quasimonochromatic . (
the effect of @xmath194 will be assumed negligible , like in sec .
[ sec : qs ] . )
more specifically , harmonics of @xmath288 remain suppressed with exponential accuracy in the small parameter @xmath341 .
( if the system supports more than one mode , we assume that only one mode is excited , and the spectral distance to other modes is @xmath342 . ) here , we assume that @xmath282 is also the time scale of the evolution of @xmath338 and @xmath288 , and @xmath343 .    clearly , if @xmath337 is a solution , then @xmath344 a solution too , and so is @xmath345 .
since the latter is also real , it can be adopted as an _ exponentially accurate asymptotic model _ of a realizable physical process , which we can also express as follows : @xmath346.\end{gathered}\ ] ] let us substitute this model into eq .
( [ eq : lgamma ] ) .
as usual @xcite , the rapidly oscillating terms that couple @xmath336 and @xmath347 can be eliminated , which is also expected from the definition of @xmath336 and @xmath347 .
( such `` @xmath339-averaging '' is denoted below with angular brackets , @xmath348 . )
this implies that the medium responds to @xmath336 and @xmath347 independently , as if there were two independent subsystems @xmath49 that interact with @xmath336 and @xmath347 , respectively .
the corresponding responses will be denoted @xmath349 .
hence , instead of the lagrange multipliers @xmath188 and @xmath350 that we introduced before , now we must introduce twice as much lagrange multipliers .
we denote them @xmath351 and @xmath352 .
hence , like in sec .
[ sec : pdaeff ] , the effective action ( [ eq : lgq ] ) can be represented as follows : @xmath353 @xmath354 [ where we used eq .
( [ eq : xisgn ] ) ] , and @xmath355 , where @xmath356,\end{gathered}\ ] ] assuming the notation @xmath357 and @xmath358 .
also like in sec .
[ sec : pdaeff ] , one finds @xmath359 [ the quantity @xmath351 must not be confused with @xmath188 given by eq .
( [ eq : mupsi ] ) .
in particular , note the extra factor @xmath360 . ]    as usual , the functions @xmath339 , @xmath338 , and @xmath361 are considered independent . then
, @xmath362\,dt.\end{gathered}\ ] ] we seek to approximate @xmath363 to the zeroth order in the small parameter @xmath364 where @xmath285 is the phase - mixing time ( sec .
[ sec : pm ] ) .
( as in sec .
[ sec : qs ] , we adopt for clarity that all the three parameters are comparable to each other . )
we assume @xmath365 and also the standard ordering for @xmath299 ( sec .
[ sec : weyl ] ) , which will be be verified _ a  posteriori _ : @xmath366 then , the individual terms on the right side of eq .
( [ eq : aeffgo1 ] ) can be calculated as follows .
since @xmath203 is a local function of @xmath335 , the function @xmath368 can not depend on @xmath339 other than through @xmath288 .
hence , @xmath369\ , dt,\end{gathered}\ ] ] where we used @xmath370 and integrated by parts with @xmath371 .
since @xmath372 while @xmath373 , @xmath374 , and @xmath375 , all the three terms in the square brackets are @xmath376 .
hence , within the accuracy of our model , it is enough to approximate @xmath367 to the zeroth order in @xmath377 .
since @xmath203 was assumed to be a bilinear function of @xmath335 , the functional @xmath367 then has the following form : @xmath378\,dt,\end{gathered}\ ] ] where the factor @xmath360 is introduced for convenience , and @xmath379 is some matrix that is determined by the local frequency @xmath380 and , in nonstationary systems , also @xmath176 . [ below
, we also refer to this matrix using a shortened notation @xmath381 . ]
the anti - hermitian part of this matrix provides zero contribution to @xmath382}$ ] , so we can assume that @xmath381 is hermitian for all real @xmath288 without loss of generality . then , `` re '' in eq .
( [ eq : ad ] ) can be omitted .
let us now approximate @xmath383 given by eq .
( [ eq : xiqsgn ] ) .
first , note that , due to eq .
( [ eq : ttau ] ) , we have @xmath384 hence , @xmath385 where the second term on the right is estimated as @xmath386 = o(\epsilon_\omega^2 ) \ll 1.\end{gathered}\ ] ] this gives @xmath387 omitting @xmath388 , one can rewrite eq .
( [ eq : xiqsgn ] ) as follows : @xmath389 here , the matrix @xmath143 is given by @xmath390 since the integrand in eq .
( [ eq : xiaux ] ) does not depend on @xmath339 explicitly , the correction @xmath391 can be omitted for the same reason as in sec .
[ sec : aq ] .
hence , we switch from integration in @xmath307 to the integration in @xmath308 and get @xmath392 note that @xmath393 ( sec .
[ sec : gobasic ] ) . this function can be approximated with @xmath303 everywhere except for an @xmath391 part of the integration domain .
since corrections @xmath391 in @xmath383 are deemed negligible , we obtain @xmath394    as a side note , one can also rewrite eq .
( [ eq : xixi ] ) as @xmath395 , so @xmath383 is understood as a contribution of an effective potential @xmath396 that the subsystem @xmath25 experiences due to the interaction with the subsystem @xmath49 .
the potential @xmath397 represents the ponderomotive energy ( or , more precisely , the ponderomotive term in the lagrangian ) of @xmath49 . in the context of particle dynamics with hermitian @xmath301 , the obtained linear relation between @xmath397 and @xmath301
is also known as the @xmath398-@xmath399 theorem @xcite .    since @xmath383 has the same form as @xmath367 given by eq .
( [ eq : ad ] ) , let us combine these two terms into a `` dressed action '' @xmath400 . for that , let us introduce @xmath401 , whose hermitian and anti - hermitian parts are @xmath402 then , @xmath403 this can also be written as @xmath404 , where @xmath405 is understood as the dressed lagrangian of @xmath25 .
accordingly , one gets @xmath406\ , dt,\end{gathered}\ ] ] where @xmath407 is understood as the number of quanta of the oscillation mode , or the action of this mode @xcite [ again , not to be confused with the action such as @xmath195 ] .      a straightforward calculation ( appendix  [ app : db ] ) gives @xmath410.\label{eq : dy}\end{gathered}\ ] ] using @xmath411 , one further gets @xmath412 { \partial}_\omega \big\ } \cdot \widetilde{\alpha}_{\omega , a}.\end{gathered}\ ] ] note that @xmath413 ( sec .  [
sec : gobasic ] ) and @xmath414 [ eq .
( [ eq : dgosc ] ) ] . since we need to calculate @xmath415 to the zeroth order in @xmath377 , the contribution of @xmath416
can then be neglected .
hence , @xmath417\,\delta\theta(\bar{t})\,d\bar{t}.\end{gathered}\ ] ] similarly , one obtains that @xmath418 , so @xmath419 is entirely negligible . in summary , substituting this result into eq .
( [ eq : aux3001 ] ) gives @xmath420 where we omitted corrections @xmath391 and , for the same reason , replaced @xmath295 with @xmath304 .
using @xmath421 along with eqs .
( [ eq : dadr ] ) and ( [ eq : dbgo ] ) and also requiring @xmath422 , we get the following set of ele : @xmath423 plus the conjugate equation for @xmath361 .
note that eq .
( [ eq : it ] ) implies @xmath424 , which is in agreement with the original assumption ( [ eq : alphaord ] ) .
it is also instructive to rewrite these in terms of the unit polarization vector @xmath425 and the scalar amplitude @xmath11 , assuming @xmath426 : @xmath427   = -2({\mathsf{e}}_\omega^\dag\cdot { \mathsf{d}}_{\omega , a}\cdot { \mathsf{e}}_\omega)\,a^2.\label{eq : ampleq}\end{gathered}\ ] ] equation ( [ eq : dhe ] ) can be understood as a local dispersion relation , which determines @xmath428 and the local frequency @xmath429 .
equation ( [ eq : ampleq ] ) can be understood as an amplitude equation . in a dissipationless medium , which we _ define _ as a medium with @xmath430 ,
this equation becomes @xmath431 , which is known as the action conservation theorem @xcite .
for example , consider a one - dimensional stationary system . in this case , eq .
( [ eq : dhe ] ) becomes @xmath432 and the amplitude equation gives @xmath433 , where @xmath434 ( this quantity is not to be confused with @xmath164 that was introduced in sec .  [ sec : co ] . )
one hence obtains that @xmath1 has a well - defined complex frequency ; i.e.,@xmath435 , where @xmath436 . within the accuracy of our theory , this can be understood as a solution of the complex dispersion relation @xmath437 as seen from the fact that @xmath438 combined with eqs .
( [ eq : dh ] ) and ( [ eq : gammainst ] ) .
hence , @xmath439 can be understood as the complex dispersion function .
although the above results were derived for discrete oscillators , they can be readily extended to waves in continuous media .
a continuous medium is understood as a system where the scalar product can be expressed as an integral over some continuous space @xmath440 ( we assume that @xmath440 is euclidean for clarity , but see refs .
@xcite for a more general case ) , specifically as @xmath441 here , @xmath442 , and @xmath443 and @xmath444 are fields of some finite dimension @xmath445 .
for example , they can be electromagnetic fields in physical space ; then @xmath446 .
( we use bold symbols to distinguish @xmath447-dimensional coordinates and @xmath445-dimensional vectors from @xmath45- and @xmath50-dimensional coordinates and vectors that we introduced earlier , since @xmath45 and @xmath50 are infinite in the continuous limit . )
accordingly , the polarizability operator can be rewritten  as @xmath448 the symmetrized kernel is introduced as @xmath449 , where @xmath290 and @xmath280 are defined as usual [ eq .  ( [ eq : ttau ] ) ] , and @xmath450 we also introduce the corresponding weyl symbol ( the term is used with the same reservations as earlier @xcite ) , @xmath451 the modification of the remaining notation is obvious , so it will not be presented in detail .
let us consider a special case of practical interest where the medium is weakly inhomogeneous in both time and space .
( anisotropy and general dispersion , including both temporal and spatial dispersion , are implied . ) unlike in sec .  [
sec : fastglobal ] , where all oscillators were assumed to have the same @xmath288 at given @xmath176 , we now allow the frequency to be spatially inhomogeneous . specifically , we adopt @xmath452\end{gathered}\ ] ] and assume that the frequency and the wave vector , @xmath453 are slow functions of @xmath454 . explicitly , we require @xmath455 where @xmath377 is defined in eq .
( [ eq : epsilonom ] ) , and @xmath456 is the analogous parameter that involves spatial scales instead of temporal scales .
then , one obtains @xmath457 , where @xmath458 , and @xmath459 is the lagrangian density given  by @xmath460 this leads to @xcite @xmath461 + \frac{1}{4}\,[{{\boldsymbol{\rm q}}}^\dag \cdot { { \boldsymbol{\rm { \mathsf{d}}}}}_{(\omega , { { \boldsymbol{\rm k } } } ) , h } ] \cdot \delta { { \boldsymbol{\rm q}}}\bigg\},\end{gathered}\ ] ] where @xmath462 and @xmath463 .
the expression for @xmath464 is obtained like in sec .
[ sec : fastglobal ] and can be written as @xmath465\,\delta \theta.\end{gathered}\ ] ] hence , one obtains the following ele : @xmath466,\label{ref : act}\\ \delta { { \boldsymbol{\rm q}}}^\dag & : \quad { { \boldsymbol{\rm { \mathsf{d}}}}}_{(\omega,{{\boldsymbol{\rm k } } } ) , h}\cdot { { \boldsymbol{\rm q } } } = 0,\label{eq : qdisp}\end{aligned}\ ] ] plus the conjugate equation for @xmath467 .
it is also instructive to rewrite these equations as @xmath468 here , @xmath469 is the unit polarization vector .
we also added an equation for @xmath470 that flows from its definition [ eq .
( [ eq : defwk ] ) ] , introduced @xmath471 , which represents the group velocity , and also introduced @xmath472 which represents the local dissipation rate .
equations ( [ eq : dhe2])-([eq : gammago ] ) form a complete set of equations that describe dissipative go waves in inhomogeneous nonstationary media .
they are in agreement with results of ref .
@xcite , where dissipative effects were described as an addition to the variational formulation rather than as a part of it .
( the variational formulation of _ dissipationless _ wave dynamics is also discussed in refs .
@xcite . )
applications of these general equations to electromagnetic waves are discussed in detail in ref .
@xcite . here
, we only point out that one can choose @xmath473 to be the complex amplitude of the electric field , @xmath474 ; then , @xmath475 the second term in eq .
( [ eq : lem ] ) is understood as @xmath476 , where we used faraday s law @xmath477 to express the complex amplitude of the magnetic field @xmath478 , and @xmath479 is the speed of light . also , @xmath6 is the dielectric tensor in the spectral representation ( or , more specifically , the weyl symbol of the dielectric permittivity operator ) , and the indexes @xmath320 and @xmath3 denote its hermitian and anti - hermitian parts , as usual .
accordingly , eq .  (
[ ref : act ] ) can be understood as a restatement of poynting s theorem @xcite .    as a side note , an even simpler variational
derivation of the results reported in this section is possible by means of the weyl calculus combined with a technique _
  la _  ref .
we leave details to future publications .
in summary , we have formulated a variational principle for a dissipative subsystem of an overall - conservative system by introducing constant lagrange multipliers and lagrangians nonlocal in time .
we call it the variational principle for projected dynamics , or vppd .
we have also elaborated on applications of the vppd to the special case of linear systems , particularly in the context of wave propagation in general linear media .
the focus was on how the variational formulation helps in deriving reduced models , such as the quasistatic and go models .
in particular , we have proposed a variational formulation of dissipative go in a linear medium that is inhomogeneous , nonstationary , nonisotropic , and exhibits both temporal and spatial dispersion simultaneously .
the `` spectral representation '' of the dielectric tensor @xmath6 that enters go equations is shown to be , strictly speaking , the weyl symbol of the medium dielectric permittivity .
this can be considered as an invariant definition of @xmath6 , as opposed to _
ad  hoc _ definitions that are commonly adopted in literature and have been a source of a continuing debate .
our work is also intended as a stepping stone for extending the variational theory of modulational stability in general wave ensembles @xcite to dissipative systems .
details will be described in a separate publication .
the authors thank d.  bnisti for stimulating discussions .
the work was supported by the nnsa ssaa program through doe research grant no .
de - na0002948 , by the u.s .
doe through contract no .
de - ac02 - 09ch11466 , and by the u.s .
dod ndseg fellowship through contract no .
as a side note , equations of sec .  [
sec : sym ] can be expressed in a more compact form in terms of the phase space coordinate @xmath480 and the poincar two - form  @xmath481 , @xmath482 \eta \end{array } \right ) , \quad \omega_{\alpha\beta } \doteq   \left ( \begin{array}{c @{\quad } c } 0 & - \mathbb{i}\\[3pt ]   \mathbb{i } & 0 \end{array } \right),\end{gathered}\ ] ] where @xmath483 is a @xmath484 unit matrix
. in particular , one can cast eq .
( [ eq : aux33 ] ) as @xmath485 = \int^{t_2}_{t_1 } \big[\lambda^{(q ) }   + \frac{1}{2}\,\omega_{\alpha\beta } z^\alpha \dot{z}^\beta - h^{(z)}\big]\,dt\end{gathered}\ ] ] ( greek indexes span from 1 to @xmath64 ) and eq .
( [ eq : aux34 ] ) as @xmath486 = \big[p \cdot \delta q + \frac{1}{2}\,\omega_{\alpha\beta}z^\alpha\,\delta z^\beta\big]\,\big|_{t_1}^{t_2 } + \delta { \mathscr{f}}^{(q ) } + \delta { \mathscr{f}}^{(z ) } , \notag\\ \delta { \mathscr{f}}^{(z ) } \doteq \int^{t_2}_{t_1}(\omega_{\alpha\beta}\dot{z}^\beta - { \partial}_\alpha h^{(z)})\,\delta z^\alpha\,dt,\end{gathered}\ ] ] where @xmath487 and @xmath488 . accordingly
, hamilton s equations ( [ eq : he ] ) are cast as follows : @xmath489 also , eq .  ( [ eq : aeffsr ] ) can be represented as follows : @xmath490   = { \mathscr{a}}[q , \hat{z } ]   + \frac{1}{2}\,\omega_{\alpha\beta } \hat{z}^\alpha(t_2 ) z^\beta_t,\end{gathered}\ ] ] where @xmath491 \eta_t \end{array } \right).\end{gathered}\ ] ] on physical trajectories one has @xmath492 . since @xmath481 is antisymmetric
, this implies that the second term on the right side of eq .
( [ eq : auxx ] ) vanishes on such trajectories .
here , we discuss an instructive special case , namely , the situation when the hamiltonian @xmath493 in eq .
( [ eq : apq ] ) does not depend on @xmath2 explicitly . in this case ,
@xmath107 is conserved and can be treated as a constant parameter , which we denote @xmath494 .
hence , one can write @xmath495 = \int_{t_1}^{t_2 } [ \lambda^{(q)}(t , q , \dot{q } ) - h^{(\xi)}(t , q , \dot{q } ; j)]\,dt.\end{gathered}\ ] ] an alternative derivation of the same result , also known as the routh reduction can be found , e.g.,in ref .
@xcite .    in particular , eq .
( [ eq : aeffl ] ) can be applied @xcite when @xmath49 is an _
adiabatic oscillator _
in which @xmath136 are the angle - action variables .
( then , @xmath494 is called an adiabatic invariant @xcite . ) in the special case when @xmath49 is a _ linear _ oscillator , one can further write @xmath496 , where @xmath497 is the corresponding canonical frequency ( i.e.,such that @xmath498 ) , which is independent of @xmath494 but may depend on @xmath97 .
then , @xmath499 = \int_{t_1}^{t_2 } [ \lambda^{(q)}(t , q , \dot{q } ) - j \cdot \omega(t , q , \dot{q})]\,dt.\end{gathered}\ ] ]    for example , one can use this to describe a charged particle in a weakly inhomogeneous field @xmath500 . in this case , @xmath497 is the gyrofrequency , and @xmath494 is proportional to the particle magnetic moment @xmath501 , such that @xmath502 @xcite .
then , the last term in eq .
( [ eq : leff ] ) is nothing but the effective potential responsible for the diamagnetic force .
in addition to the linear coupling discussed in sec .  [ sec : co ] , it is instructive to consider parametric interactions , particularly because they subsume quantum interactions and thus can be understood as most general .
specifically , let us consider @xmath503 where @xmath504 is some hermitian operator .
( we assumed that @xmath162 is independent of derivatives of @xmath1 for brevity , but the corresponding generalization is straightforward . )
the total lagrangian in this case can be written as @xmath505 and the corresponding ele are @xmath506 notably , eq .  ( [ eq : psiq2 ] )
conserves @xmath507 even in the presence of coupling , in contrast with eq .
( [ eq : psiq ] ) .    for projected dynamics ,
the effective action @xmath195 is introduced just like in sec .
[ sec : pdaeff ] .
when the solution for @xmath166 and @xmath508 are substituted into eq .
( [ eq : aux103 ] ) , the second term and the third term cancel each other due to eq .
( [ eq : psiq2 ] ) , so one ends up with @xmath121 = \int^{t_2}_{t_1 } \lambda^{(q ) } ( t , q , \dot{q})\,dt
+ { \mathcal{b } } , \\
{ \mathcal{b } } \doteq \frac{i}{2}\,\big[\nu^\dag \cdot \hat{\psi}(t_2 ) - \hat{\psi}^\dag(t_2 ) \cdot \nu\big].\end{gathered}\ ] ]      let us show how a correct ele for @xmath1 is obtained from here by a straightforward calculation .
first , @xmath252 = \int_{t_1}^{t_2 } [ { \partial}_q \lambda^{(q ) } - \dot{p}]\,\delta q\,dt + \delta{\mathcal{b}},\label{eq : daeffq}\\ \delta { \mathcal{b } } = \frac{i}{2}\,\big[\nu^\dag \cdot \delta \hat{\psi}(t_2 ) - \delta \hat{\psi}^\dag(t_2 ) \cdot \nu\big].\end{gathered}\ ] ] in order to calculate @xmath244 , let us use the fact that @xmath509 , where @xmath510 is a unitary propagator given by the following time - ordered exponential : @xmath511.\end{gathered}\ ] ] if the variation @xmath512 is localized to an infinitesimal time interval @xmath513 following some time @xmath280 , one can readily write @xmath514 where one can further substitute @xmath515 . then , for a general variation @xmath512 , one gets @xmath516\,\delta q(\tau)\,d\tau.\notag\end{gathered}\ ] ] for @xmath517 at any given time @xmath176 , this gives @xmath518.\end{gathered}\ ] ] then , after substituting eq .
( [ eq : mupsi ] ) , one gets @xmath519 \\ = - \text{re}\,\big[\hat{\psi}^\dag(t ) \cdot { \partial}_q h(t , q(t ) ) \cdot \hat{\psi}(t)\big].\label{eq : aux112}\end{gathered}\ ] ] here , we also used @xmath520 and the fact that @xmath521^\dag \cdot u_h(t_2 , t)$ ] is a unit operator , because @xmath510 is unitary . since @xmath522 is hermitian , one can omit `` re '' , so eq .  ( [ eq : aux112 ] ) becomes @xmath523 after substituting this into eq .
( [ eq : daeffq ] ) and requiring @xmath524 = 0 $ ] , one obtains @xmath525 this is equivalent to eq .
( [ eq : qleffq ] ) , as expected .
suppose that the coordinate @xmath1 is chosen such that @xmath205 , where @xmath11 is a small parameter .
assuming that the interaction is weak , a typical interaction hamiltonian is then linear in @xmath1 . like @xmath142 ,
any zeroth - order term can be removed by a variable transformation , so we assume @xmath526 , where @xmath527 may depend only on @xmath176 , if at all .
we now seek to derive an approximation of @xmath53 accurate up to @xmath528 ; i.e.,terms @xmath529 will be neglected .    to construct such an approximation , we adopt an asymptotic representation of @xmath166 as a power series in @xmath11 ; i.e.,@xmath530 , where @xmath531 . according to eq .
( [ eq : psiq2 ] ) , @xmath532 , so truncating this series as @xmath533 is sufficient .
since @xmath194 is constant , @xmath49 is then fully described by @xmath534 , which serves as the new independent variable instead of the original @xmath166 .
hence , after omitting complete time derivatives ( which do not affect ele ) , eq .
( [ eq : aux103 ] ) becomes similar to eq .
( [ eq : lgamma ] ) , namely , @xmath535 where the new @xmath536 and @xmath199 are given by @xmath537 ( one may also recognize this as the leading - order born approximation of the original problem . ) this result shows that theory of linear dispersion for interaction hamiltonians of the form ( [ eq : hint2 ] ) can be constructed identically to that for interaction hamiltonians of the form ( [ eq : hint ] ) .
in particular , this means that general predictions of linear dispersion theory are independent of whether oscillators comprising a medium are classical or quantum .
here , we discuss some properties of the functions @xmath291 and @xmath301 that are introduced in sec .  [ sec : weyl ] .
first , eq .
( [ eq : aproph ] ) is proved as follows : @xmath541 where we used that @xmath267 is real by definition .
equation ( [ eq : apropa ] ) is proved similarly .
second , let us explicitly derive an approximation for @xmath291 in the limit ( [ eq : vro ] ) . using eq .
( [ eq : avu ] ) for @xmath267 together with eq .
( [ eq : mcuasym ] ) for @xmath260 , we obtain @xmath542.\notag\end{gathered}\ ] ] then , using @xmath543 , one gets @xmath544}{h_0 - \omega }   \\ + \frac{1- \cos[(h_0 + \omega)\tau_m]}{h_0 + \omega }   \bigg\ } \cdot v,\notag\end{gathered}\ ] ] @xmath545}{h_0 - \omega }   - \frac{\sin[(h_0 + \omega)\tau_m]}{h_0 + \omega }   \right\ } \cdot v,\notag\end{gathered}\ ] ] or , in the continuous - spectrum limit ( sec .
[ sec : pm ] ) , @xmath546}{2(\omega - \omega ) } \,[f(\bar{t } , \omega ) - f(\bar{t } , -\omega)]\,d\omega,\label{eq : ath}\\ \widetilde{\alpha}_{\omega , a } = \int \frac{\sin[(\omega - \omega)\tau_m]}{2(\omega - \omega ) } \,[f(\bar{t } , \omega ) - f(\bar{t } , -\omega)]\,d\omega.\label{eq : ata}\end{gathered}\ ] ] notably , the anti - hermitian part of @xmath291 vanishes in the zero-@xmath288 limit ; i.e.,@xmath547 , so @xmath548 , and @xmath549 \cdot v \\
\to \int \omega^{-1}[1 - \cos(\omega \tau_m)]\,f(\bar{t } , \omega)\,d\omega . \end{gathered}\ ] ]      asymptotic expressions for @xmath301 under the condition ( [ eq : vro ] ) can be obtained by taking the large-@xmath289 limit of eqs .
( [ eq : ath ] ) and ( [ eq : ata ] ) .
the oscillating terms vanish for @xmath550 . for @xmath551 ,
one has @xmath552 \equiv 0 $ ] , so , instead of averaging the corresponding cosine to zero , we entirely exclude the point @xmath551 from the integration domain .
also , @xmath553/(\omega - \omega ) \to \pi \delta(\omega - \omega)$ ] .
hence , one gets @xmath554,\end{gathered}\ ] ] where @xmath555 denotes the cauchy principal value .
equation ( [ eq : ahaa ] ) can be understood also as @xmath556 , where @xmath557 is the hilbert transform . using that @xmath558 ,
we then obtain @xmath559 .
thus , in summary , @xmath560 which can be recognized as the kramers - kronig relations .
to calculate @xmath562 used in sec .
[ sec : fastglobal ] , we proceed as follows .
first , using eqs .
( [ eq : auxb ] ) and ( [ eq : nupm ] ) , we get @xmath563    \right\}.\end{gathered}\ ] ] using the notation @xmath564 , we further rewrite this as follows : @xmath565      \right\ }       \notag\\ & = \frac{1}{16 } \bigg\{\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , \tau)\,q^{\dag}(t ) \cdot v^{\dag}(t ) \cdot v(t ' ) \cdot q(t')\,\delta \theta ( t ' )    + \int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      [ e(\bar{t } , \tau)\,q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t')\cdot q(t')]^{\dag}\,\delta \theta ( t ' )      \notag \\ &   \mbox { } \quad        - i\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , \tau)\,q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t')\cdot \delta q(t ' )      + i\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      [ e(\bar{t } , \tau)\,q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t')\cdot \delta q(t')]^{\dag }      \bigg\ }      \notag \\ & = \frac{1}{16}\bigg\ {      \int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , \tau)\,q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t')\cdot q(t')\,\delta \theta ( t ' )      + \int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , -\tau)\,q^{\dag}(t ' ) \cdot v^{\dag}(t ' ) \cdot v(t)\cdot q(t)\,\delta \theta ( t ' )      \notag \\ &   \mbox { } \quad        - i\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , \tau)\,q^{\dag}(t ) \cdot v^{\dag}(t)\cdot v(t')\cdot \delta q(t ' )      + i\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , -\tau)\,\delta q^{\dag}(t')\cdot v^{\dag}(t')\cdot v(t)\cdot q(t )      \bigg\ }      \notag \\ & = \frac{1}{16}\bigg\ {       \int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,       e(\bar{t } , \tau)\,q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t ' ) \cdot q(t')\,\delta \theta ( t ' )       + \int_{t_1}^{t_2 } dt ' \int_{t_1}^{t_2 } dt\ ,       e(\bar{t } , \tau)\,q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t')\cdot q(t')\,\delta \theta ( t )       \notag \\ &   \mbox { } \quad         - i\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,       e(\bar{t } , \tau)\,q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t')\cdot \delta q(t ' )       + i\int_{t_1}^{t_2}dt \int_{t_1}^{t_2 } dt'\ ,       e(\bar{t } , \tau)\,\delta q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t')\cdot q(t ' )      \bigg\ }      \notag \\ & = \frac{1}{8}\bigg\ {      \int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , \tau)\,q^{\dag}(t ) \cdot v^{\dag}(t)\cdot v(t')\cdot q(t')\,\frac{\delta \theta ( t')+\delta \theta(t)}{2 }      \notag \\ &   \mbox { } \quad        -\frac{i}{2}\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , \tau)\,q^{\dag}(t ) \cdot v^{\dag}(t ) \cdot v(t')\cdot \delta q(t ' )      + \frac{i}{2}\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ , e(\bar{t } , \tau)\,\delta q^{\dag}(t)\cdot v^{\dag}(t)\cdot v(t')\cdot q(t ' )      \bigg\ }      \notag \\ & = \frac{1}{8}\,\text{tr } \int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,       e(\bar{t } , \tau ) \left\{w(\bar{t } , \tau)\ ,       \frac{\delta \theta ( t')+\delta \theta ( t)}{2 }       + \frac{i}{2 } \left[q(t')\cdot \delta q^\dag ( t)-\delta q(t')\cdot q^\dag ( t)\right ]       \right\}\cdot v^{\dag}(t)\cdot v(t').\notag\end{aligned}\ ] ] next , we notice that @xmath566 = \delta { \mathcal{b}}_+[q^ * , -\theta ] = ( \delta { \mathcal{b}}_+[q^ * , -\theta])^*$ ] .
hence , @xmath567\right\ }    \cdot v^{\dag}(t)\cdot v(t')\right]^*\\ = -\frac{1}{8}\,\text{tr}\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,    e(\bar{t } , \tau)\,\left[w(\bar{t } , \tau)\,\frac{\delta \theta ( t')+\delta \theta ( t)}{2 } + \delta y(\bar{t } , \tau)\right ]    \cdot v^{\rm t}(t)\cdot v^*(t'),\notag\end{gathered}\ ] ] where @xmath143 is given by eq .
( [ eq : ww ] ) and @xmath568 is given by eq .
( [ eq : dy ] ) . for @xmath569 ,
this leads to @xmath570        \cdot \left[v^{\dag}(t)\cdot v(t')-v^{\rm t}(t)\cdot v^*(t')\right]\notag\\ & = \frac{1}{8}\,\text{tr}\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , \tau ) \left[w(\bar{t } , \tau)\,\frac{\delta \theta ( t')+\delta \theta(t)}{2 } + \delta y(\bar{t } , \tau)\right ]      \cdot \left[-2i \alpha(t , t')\right]\notag\\ & = - \frac{i}{4}\,\text{tr } \int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,      e(\bar{t } , \tau ) \left[w(\bar{t } , \tau)\,\frac{\delta \theta ( t')+\delta \theta(t)}{2 } + \delta y(\bar{t } , \tau)\right ]      \cdot \alpha(t , t')\notag\\ & = - \frac{i}{4}\,\text{tr}\int_{t_1}^{t_2 } dt \int_{t_1}^{t_2 } dt'\ ,
e^{i \omega \left(\bar{t}\right ) \tau }      \left[w(\bar{t } , \tau)\,\delta \theta \left(\bar{t}\right)+\delta y(\bar{t } , \tau)\right ]      \cdot \bar{\alpha } ( \bar{t } , \tau ) + o\left(\epsilon_\omega^2\right)\notag\\ & \approx \delta { \mathcal{b}}_{\theta}+\delta { \mathcal{b}}_q,\end{aligned}\ ] ] where @xmath571 and @xmath572 are given by eqs .
( [ eq : dbtheta ] ) and ( [ eq : dbq ] ) , respectively .
the only exception are media with extremely simple dispersion , such as those described by a local refraction index .
although often presented as characteristic models of dielectric media ,
they , in fact , do not provide a truly representative picture of general wave dynamics .
d.  bnisti , _ envelope equation for the linear and nonlinear propagation of an electron plasma wave , including the effects of landau damping , trapping , plasma inhomogeneity , and the change in the state of wave _ , phys .
plasmas * 23 * , 102105 ( 2016 ) .                      c.  m. bender , m.  gianfreda , n.  hassanpour , and h.  f. jones , _ comment on `` on the lagrangian and hamiltonian description of the damped linear harmonic oscillator '' [ j. math .
48 , 032701 ( 2007 ) ] _ , j. math
. phys . * 57 * , 084101 ( 2016 ) .
r. l. dewar , r. f. abdullatif , and g. g. sangeetha , _ complex nonlinear lagrangian for the hasegawa - mima equation _ , in proceedings of the 20th iaea fusion energy conference ( 2004 ) , th / p6 - 1 , http://www-naweb.iaea.org/napc/physics/fec/fec2004/datasets/th_p6-1.html .
this implies @xmath50 _ real _ bc per each end point , as in the original problem ( sec .
[ sec : coupled ] ) .
this is in contrast with ref .
@xcite , where @xmath50 _ complex _ bc were assumed per end point in the quantumlike representation [ namely , @xmath578 , so the problem was overdefined .
a brief overview of the weyl calculus can be found , for instance , in refs .
@xcite . here
, we only notice that , for an operator @xmath3 given by @xmath579 , its weyl image is @xmath580 .
the polarizability operator can be brought to the form assumed for @xmath581 by making two replacements in eq .
( [ eq : alphaq ] ) : ( i ) @xmath582 , which is justified provided we consider the operator on functions that are identically zero at @xmath583 ; and ( ii ) @xmath584 , which is justified by the presence of the heaviside function in the integrand .    as asymptotic approximation of @xmath585 \doteq \int_a^b    g(t)\,e^{i\vartheta(t ) } \,dt}$ ] , where @xmath586 evolves at some rate @xmath587 , can be found via integration by parts .
this gives @xmath588 = \smash{\int_a^b } [ ( d / dt)(ig    e^{i\vartheta}/\omega )
- i g'(t)/\omega]\,dt$ ] or , equivalently , @xmath588 =    \smash{(ig e^{i\vartheta}/\omega)|_a^b } - \smash{{\mathscr{j}}[i g'(t)/\omega]}$ ] , where the functional on the right is smaller than the functional on the left by a factor @xmath589 . the calculation can also be iterated . then
, @xmath588 $ ] can be represented as an asymptotic series of terms that are determined by @xmath590 , @xmath339 and their derivatives evaluated locally at the ends of the integration domain .          in linear systems , @xmath591 , where @xmath592 is the energy .
thus , the action conservation law that we presented is also known as the conservation of an adiabatic invariant @xmath593 . for details ,
see , e.g. , ref .
@xcite .
h.  qin , j.liu , j.  xiao , r.  zhang , y.  he , y.  wang , y.  sun , j.  w. burby , l  ellison , and y  zhou , _
canonical symplectic particle - in - cell method for long - term large - scale simulations of the vlasov - maxwell equations _ , nucl .
fusion * 56 * , 014001 ( 2015 ) .
a. stern , y. tong , m. desbrun , and j.  e. marsden , _ geometric computational electrodynamics with variational integrators and discrete differential forms _ , in _ geometry , mechanics , and dynamics : the legacy of jerry marsden _ , edited by d. e. chang , d. d. holm , g. patrick , and t. ratiu ( springer , new york , 2015 ) , p. 437 .