Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: theoretical study of the space - time evolution of matter in high energy nuclear collisions has a long history since the early pioneering works of fermi and landau employing thermodynamics and fluid mechanics  @xcite . more recent works  @xcite have been motivated by the prospect of studying a new form of matter experimentally by means of very energetic collisions of heavy nuclei as those currently underway at brookhaven with relativistic heavy - ion collider ( rhic ) and planned at cern with large hadron collider ( lhc ) . the data taken from the rhic experiments has shown existence of a strong anisotropic ( elliptic ) collective flow of hadrons in non - central collisions  @xcite , indicating early local equilibration of dense matter produced by the collision  @xcite . this supports the hydrodynamic picture of matter evolution and gives us a hope to learn information of the equation of state of dense matter from systematic analyses of the data . much attention has been paid to the early stage of matter evolution where a dense plasma of unconfined quarks and gluons has been expected to be formed through some complex non - equilibrium processes  @xcite . one anticipates also a breakdown of the hydrodynamic behavior of the system when it is diluted sufficiently and the collision time exceeds the characteristic time scale of expansion . this stage of the matter evolution is usually refered to as the freezeout stage . the aim of this work is to present a kinetic theory which is designed to describe the freezeout stage of the expanding hadron gas . the freezeout of the expanding hadronic matter would proceed in several steps ; chemical freezeout of the relative abandunce of hadron spieces may occur before the kinetic freezeout of the momentum distribution of hadrons . observed relative abandunce of hadrons fit very nicely to a simple picture that the chemical freezeout occurs at a certain temperature and baryon chemical potential  @xcite . some of the other observables , such as two particle momentum correlations , the pion analogue of the hanbury brown and twiss two photon intensity interferometry , may be very sensitive to the dynamics of the final kinetic freezeout as indicated by the  hbt puzzle " found in the recent data analysis  @xcite . if the initial state of the dense matter formed in the collision is a plasma of deconfined quarks and gluons , freezeout of the color degrees of freedom should proceed these processes when the quark - gluon plasma hadronizes . chiral symmmetry  @xcite , an approximate global symmetry of qcd which becomes exact in the limit of vanishing quark masses and is considered to be broken spontaneously in the qcd vacuum , may also play an important role in the freezeout dynamics .    as a dense hadronic matter formed by ultrarelativistic nuclear collision is diluted by the expansion , one expects that the system undergoes a phase transition associated with the spontaneous breakdown of the chiral symmetry which is restored temporalily after the collision by the formation of a quark - gluon plasma . as the quark - gluon plasma hadronizes and the system turns into the confining phase , the system would gradually develop a vacuum chiral condensate and remaining excitaions would expand in the influence of the growing chiral condensate . this physical picture has been elaborated in terms of a classical equation of motion for the chiral condensate ; effects of excitations were described in terms of statistical fluctuations in the classical fields  @xcite . the fluctuating condensate described by classical pion field has been termed disoriented chiral condensate ( dcc ) , emphasizing the symmetry aspect of the problem  @xcite . what was missing in these classical treatments of the meson fields is the existence of particle excitations in addition to the condensate . inclusion of particle excitations requires the quantization of the fields . the effect of the quantum fluctuations of meson fields has been studied by tsue , vautherin and one of the present authors ( tm ) @xcite by the functional schrdinger picture formalism @xcite . the physical picture we have described is very similar to what happens when the dilute gas of magnetically trapped atoms of alkali metals cools down by evapolation @xcite . some of the atoms condense into the lowest single particle level in the trapping external potential forming a bose - einstein condensate . the dynamics of such a system may be described by the coupled equations of motion for the condensate , or the gross - pitaevskii equation @xcite , and the kinetic equation , or the boltzmann - vlasov equations , for the phase space distribution of the excitations in the presence of the condensate  @xcite .    in this paper we will show that a similar set of equations can be derived for a system of interacting mesons described by the relativistic quantum field theory by the mean field approximation . this approximation corresponds to a neglect of all correlations in the system  @xcite . we will make no attempt in this work to justify this approximation and leave it as an open problem for further study .    in this context we note that field theoretical derivation of a boltzmann - type kinetic equation has been given by many others  @xcite . most of these works focuses on the derivation of the collision terms . the present work is distinguished from these works in the emphasis of the role of the mean field in the evolution of the system coupled with the quasi - particle excitations in the same spirit as in @xcite . we omit the effect of the quasi - particle collisions in this work . this procedure may be reasonable , at least as a first step , for describing dynamical aspect of the freeze - out process : even in the absence of collisions , interactions between the quasi - particles and the evolving condensate would affect the final particle distribution . but the present approach is not adequate , however , for the early thermalization problem , where the collision terms play essential role  @xcite . we adopt the standard heisenberg picture instead of the functional schrdinger picture used in @xcite since we found it more straightforward to see the connection to familiar semi - classical kinetic equation . in this first of a series of papers we shall use one - component real scalar field model interacting via a @xmath0 self - interaction term in order to concentrate on the presentation of the basic features of the theory . analysis with multi - component scalar model with @xmath1 symmetry will be deferred for the forthcoming paper  @xcite .    in the next section we formulate a quantum kinetic theory for quantum scalar meson field , starting from the heisenberg equations for quantum field of one - component real scalar . the mean field approximation replaces the products of quantum fields by the products of classical fields and the statistical average of bilinear forms of the quantum fluctuations . the latter is expressed in terms of the wigner functions which is reduced to the single particle distribution function in the classical boltzmann equation . our theory contains another forms of the wigner functions which have no classical counter parts and arises due to the coherent pair creation and annihilation processes in non - uniform systems . only in uniform systems , these _ off - diagonal _ components of the fluctuation can be eliminated by suitable redefinition of the particle mass . appearance of the off - diagonal wigner functions is reminiscent of the anomalous propagators in the microscopic theory of superconductivity  @xcite . similar structure also appears in the theory of bose - einstein condensate @xcite . some details of the mean field calculation is given in appendix . in section 3 , we apply our method to uniform systems and show that each mode characterized by the particle momentum obeys non - linear forced oscillatory motion . in section 4 we show that in equilibrium these kinetic equations are reduced to a gap equation  @xcite which determine the equilibrium amplitude of the condensate and the mass parameter as a function of the temperature . in this paper , we ignore the effect of the divergent vacuum polarizations which requires a subtle renormalization procedure in the mean field approximation  @xcite . the solution exhibits characteristic features of the first order phase transition .    in section 5 , we study slowly varying non - uniform systems . taking long wavelength approximation , the equations of motion of the diagonal components of the wigner functions are reduced to a vlasov equation in a form generalized by landau for quasi - particles excitations in quantum fermi liquid  @xcite with the quasi - particle energy given in the mean field approximation . non - vanishing off - diagonal components of the wigner functions generate extra terms in the vlasov equation which may be interpreted as particle source and sink terms .    in section 6 , we compute the dispersion relations of excitations of the system near equilibrium . solving the coupled kinetic equations by linearizing the equations with respect to small oscillatory deviations from equilibrium solution we obtain dispersion relation of the excitation modes in the system near equilibrium . we find the continuum of quasi - particle excitations in the entire space - like energy - momentum region in addition to the continuum in time - like region due to the ( thermally induced ) pair creation . we found that in the low temperature phase the meson pole shifts due to the coupling to the quasi - particle continua . the effective meson mass vanishes at the edge of the spinodal instability line of the first order transition . meson mode has been studied first by hatsuda and kunihiro as a precursory phenomenon of the chiral phase transition using the nambu - jona - lasinio model with quark fields.@xcite ] we also found that the coupling of the off - diagonal components of the wigner function plays an important role to prevent appearance of undamped tachyonic sound mode which propagates with a velocity greater than that of light . a short summary of the paper is given in section 7 with remarks on the remaining problems . in this section we derive quantum kinetic equations which describes the time evolution of the meson condensate coupled with mesonic quasi - particle excitations . we use the natural unit @xmath2 throughout this paper . we first take a simple model of a self - interacting real scalar field in the heisenberg picture . the hamiltonian is given by @xmath3 \right]\ ] ] where @xmath4 = \frac{1}{2 } m^2 { { \hat \phi}}^2 + \frac{\lambda}{4 ! } { { \hat \phi}}^4 .\ ] ] here the scalar field @xmath5 and its canonical conjugate momentum field @xmath6 are quantized by the equal - time commutation relations : @xmath7   & = &    i \delta ( { { \mathbf r}}- { { \mathbf r } } ' ) \\ \left [ { { \hat \phi } } ( { { \mathbf r } } , t ) , { { \hat \phi}}({{\mathbf r } } ' , t ) \right ] & = & [ { { \hat \pi } } ( { { \mathbf r } } , t ) , { { \hat \pi } } ( { { \mathbf r } } ' , t ) ] = 0   \ ] ] the heisenberg equation of motion of the quantum field @xmath8 is given by @xmath9 =   { { \hat \pi } } ( { { \mathbf r } } , t ) \ ] ] while the equation of motion of the canonical conjugate field @xmath10 becomes @xmath11 = ( \nabla^2 - m^2 ) { { \hat \phi}}({{\mathbf r } } , t )   - \frac{1}{3 ! } \lambda { { \hat \phi}}^3 ( { { \mathbf r } } , t ) \ ] ] elimination of the field momentum @xmath6 from these equations yields a modified klein - gordon equation for the quantum scalar field @xmath5 : @xmath12 where @xmath13 .      we are interested in the time evolution of the system described by the density operator @xmath14 where @xmath15 are a set of normalized wave functions and @xmath16 is the probability distribution for a mixed state described by this density matrix so that it satisfies @xmath17 this density matrix can be also expressed in terms of some complete set @xmath18 of the wave functions of our hilbert - fock space as @xmath19 where @xmath20    all physical information of the system to be described are contained in a specific form of the density matrix . in thermodynamic equilibrium , the density matrix is given by @xmath21 with @xmath22 = e^{\beta f ( t ) } \ ] ] where @xmath23 gives the helmholtz free energy of the system at temperature @xmath24 . in the heisenberg picture the density matrix is time - independent since the wave functions are time - independent ; all time dependence arises from the time dependence of an operator : @xmath25\ ] ] for example , we define the classical condensate fields by the statistical average of the quantum fields , @xmath26 in the following we take a gaussian ansatz for the density matrix : @xmath27 and @xmath28 where @xmath29 and @xmath30 are shifted field operators defined by @xmath31 the shifted field operators obey the same equal - time commutation relations as the original fields : @xmath32   & = &    i \delta ( { { \mathbf r}}- { { \mathbf r } } ' ) \label{q1 } \\ \left [ { { \tilde \phi } } ( { { \mathbf r } } , t ) , { { \tilde \phi}}({{\mathbf r } } ' , t ) \right ] & = & [ { { \tilde \pi } } ( { { \mathbf r } } , t ) , { { \tilde \pi } } ( { { \mathbf r } } ' , t ) ] = 0   \label{q2 } \end{aligned}\ ] ]    we will show that this choice of the density matrix is a non - equilibrium generalization of the hartree approximation in equilibrium . it has been alternatively introduced in the functional schrdinger representation  @xcite      we introduce the one - particle wigner function by @xmath33 here the particle creation and annihilation operators may be defined in terms of the fourier transforms of the shifted fields @xmath34 as @xmath35 \\ a^\dagger_{{\mathbf p}}(t ) & = & \frac{1}{\sqrt{2 { { \omega_{{\mathbf p } } } } } }     \left [ { { \omega_{{\mathbf p}}}}{{\tilde \phi}}_{{\mathbf p}}(t )   -   i { { \tilde \pi}}_{- { { \mathbf p } } } ( t ) \right ] \end{aligned}\ ] ] with @xmath36 these relations are rewritten as : @xmath37 \\ { { \tilde \pi}}({{\mathbf r } } , t ) & = & i \sum_{{\mathbf p}}e^{i { { \mathbf p}}\cdot { { \mathbf r } } } \sqrt{\frac{{{\omega_{{\mathbf p}}}}}{2 } }   \left [ a_{{\mathbf p}}^\dagger ( t ) - a_{-{{\mathbf p } } }   ( t ) \right]\end{aligned}\ ] ] the quantization rules ( [ q1 ] , [ q2 ] ) are transcribed to @xmath38 & = & \delta_{{{\mathbf p } } , { { \mathbf p } } ' } \\ \left [ a_{{\mathbf p}}(t ) , a_{{{\mathbf p } } ' } \right ( t ) ] & = &   [ a_{{\mathbf p}}^\dagger ( t ) , a_{{{\mathbf p}}'}^\dagger ( t ) ] = 0 .\end{aligned}\ ] ] with which we may interpret @xmath39 ( @xmath40 ) as annihilation ( creation ) operator of  particle excitation " with momentum @xmath41 . we also introduce the following other forms of the wigner functions : @xmath42 these wigner functions are not independent but are related to each other . the complex conjugate of the wigner functions are given by @xmath43 where the asterisk ( @xmath44 ) stands for the complex conjugate . the commutation relations imply also that @xmath45 and @xmath46 are even functions of @xmath41 @xmath47 and @xmath48 the four wigner functions may be grouped together to form a matrix form of the wigner function @xmath49    the appearance of the  off - diagonal " components of the wigner functions is reminiscent of the anomalous propagators in the bcs theory of superconductivity which arises due to the presence of fermion pair condensate  @xcite . note that our definition of the matrix components of the wigner function is slightly different from those for the propagators . we write the fourier transforms of the wigner functions as @xmath50 the ( [ fstar ] ) and ( [ fbarstar ] ) imply that @xmath51 and @xmath52 are real functions and are related to each other by @xmath53 while ( [ gstar ] ) implies that @xmath54 and @xmath55 are complex conjugate to each other : @xmath56    here we have chosen the particle  mass " @xmath57 to be different from the mass parameter @xmath58 in the original hamiltonian . physical particle mass for interacting fields is generally different from the mass parameter in the hamiltonian or lagrangian due to the effect of interaction , e.g. renormalization with or without the spontaneous symmetry breaking . it may also depend on the physical conditions described by the statistical average with the density matrix @xmath59 . since we are interested in non - equilibrium time - evolution of the system where the physical particle mass may not have a definite meaning in the intermediate states , we consider here @xmath57 just as a parameter to be chosen at our discretion , for an appropriate choice of the initial conditions specified by the gaussian density matrix . a different choice of this mass parameter would give different definition for  particle excitations " ; the wigner functions thus depend on the particular choice of the mass parameter . suppose we take a different particle mass @xmath60 to define the particle creation and annihilation operators @xmath61 \\ b^\dagger_{{\mathbf p}}(t ) & = & \frac{1}{\sqrt{2 { { \omega_{{\mathbf p } } } } ' } }     \left [ { { \omega_{{\mathbf p } } } } ' { { \tilde \phi}}_{{\mathbf p}}(t )   -   i { { \tilde \pi}}_{- { { \mathbf p } } } ( t ) \right ] \end{aligned}\ ] ] with @xmath62 these new particle creation and annihilation operators should also obey the commutation relations , @xmath63 & = & \delta_{{{\mathbf p } } , { { \mathbf p } } ' } \\ \left [ b_{{\mathbf p}}(t ) , b_{{{\mathbf p } } ' } \right ( t ) ] & = &   [ b_{{\mathbf p}}^\dagger ( t ) , b_{{{\mathbf p}}'}^\dagger ( t ) ] = 0 .\end{aligned}\ ] ] so that they are related to the original ones by the bogoliubov transformation : @xmath64 where the real parameter @xmath65 is determined by requiring that they describe the same fields : @xmath66 = \sum_{{\mathbf p}}\frac{e^{i { { \mathbf p}}\cdot { { \mathbf r } } } } { \sqrt{2 \omega'_{{\mathbf p } } } } \cdot \left [ b_{{\mathbf p}}(t ) + b_{-{{\mathbf p}}}^\dagger   ( t ) \right ] \nonumber \\\end{aligned}\ ] ] and this gives @xmath67 or @xmath68    the new wigner functions defined by replacing the creation and annihilation by the new ones are related to the original wigner functions by : @xmath69 where @xmath70 with @xmath71 for a small change of the mass parameter @xmath72 the wigner function will change by @xmath73 in particular , @xmath74    in uniform equilibrium system the mass parameter may be chosen to `` diagonalize '' the one - body mean field hamiltonian . as will be shown later , this procedure will lead to the well - known gap equation , a self - consistency condition to determine @xmath57 . if the system is slowly changing in time , one may still use such procedures adjusted to slowly varying quasi - equiliblium conditions , introducing a time dependent effective mass as a dynamical parameter to describe such adiabatic process . for such calculations , the relation ( [ df ] ) may be used to describe the adiabatic change of the wigner functions by the change of the mass parameter .    in more general non - equilibrium situations as we expect to encounter at the freeze - out stage of expanding matter , however , there may be no such appropriate condition to determine the mass parameter . the situation could even be worse : if the system goes through unstable state with respect to small fluctuation of the field , then the adiabatically determined mass parameter would become pure imaginary , reflecting the extremum , instead of the minimum , of the effective potential . in such case we may keep the value of @xmath57 at a certain real value reflecting initial conditions . the instability would then show up as appearance of a growing solution to our kinetic equations ; which would eventually be stabilized by the non - linear interaction . to extract the physical information , such as the particle distribution in the final asymptotic state , we should use the wigner function defined with the physical particle mass in the vacuum . but these asymptotic physical wigner functions may be calculated from the wigner functions with a different choice of the mass parameter by the relation ( [ rotation ] ) . the equation of motion of the classical mean field @xmath75 is obtained by taking the quantum statistical average of the field equation ( [ kg ] ) . with the gaussian ansatz for the density matrix , we find @xmath76\ ] ] this equation corresponds to the non - linear schdinger equation ( also called the gross - pitaevskii equation @xcite ) in the theory of bose - einstein condensates . so we may call this equation _ non - linear klein - gordon equation_. the non - linearity arises due to the self - interaction of the classical field @xmath75 ( condensate ) and also due to the interaction with fluctuations @xmath77 which also depends on @xmath75 implicitly . the latter may be interpreted as due to  particle excitations " , since the fluctuation can be expressed by the wigner functions as @xmath78 \label{phi2}\end{aligned}\ ] ] the time - evolution of the classical mean field @xmath79 is thus coupled with the time - evolution of the wigner functions .    to derive the equation of motion of the wigner functions , we need to compute the time - derivative of the bilinear forms of the operators @xmath80 and @xmath81 which in turn requires computation of the commutators of these operators with the hamiltonian . we decompose the original hamiltonian as @xmath82 where @xmath83 is the classical hamiltonian obtained from @xmath84 by replacing the quantum fields by their classical expectation values and @xmath85 contain the @xmath86-thrth power of the quantum fluctuation @xmath5 ( or @xmath39 and @xmath87 ) . a straightforward calculation yields @xmath88 \\ h_2 & = & \int d { { \mathbf r}}\left [ \frac{1}{2 } { { \tilde \pi}}^2 + \frac{1}{2 } ( \nabla { { \tilde \phi}})^2   + \frac{1}{2 } m^2 { { \tilde \phi}}^2 ( { { \mathbf r } } , t ) + \frac{\lambda}{4 } \phi_c^2 ( { { \mathbf r } } , t ) { { \tilde \phi}}^2 ( { { \mathbf r } } , t ) \right ] \\ h_3 & = & \frac{\lambda}{3 ! } \int d { { \mathbf r}}\phi_c ( { { \mathbf r } } , t ) { { \tilde \phi}}^3 ( { { \mathbf r } } , t ) \\ h_4 & = & \frac{\lambda}{4 ! } \int d { { \mathbf r}}{{\tilde \phi}}^4 ( { { \mathbf r } } , t ) \end{aligned}\ ] ]    the commutators of bilinear forms of @xmath80 and @xmath81 with @xmath89 vanish and the commutators with @xmath90 would give either a linear term or the third power of the fluctuation , both of which may vanish when taking the average with the gaussian density matrix . what remain to be computed are then the commutators with @xmath91 and with @xmath92 . they will give either the bilinear form of @xmath80 and @xmath87 or the fourth power of the fluctuations . the gaussian average of the resultant equations of motion of the bilinear field operators would give the desired equations of motion of the wigner functions . details of this computation is given in appendix a.    the resultant equation of motion of the wigner functions may be obtained more easily by introducing the mean field hamiltonian defined by @xmath93   \nonumber \\ & = & \int d { { \mathbf r}}\left [ \frac{1}{2 } { { \tilde \pi}}^2 + \frac{1}{2 } ( \nabla { { \tilde \phi}})^2   + \frac{1}{2 } \mu^2 { { \tilde \phi}}^2 + \frac{1}{2 } { { \delta \pi } } ( { { \mathbf r } } , t ) { { \tilde \phi}}^2 \right ]   \end{aligned}\ ] ] where @xmath94 and @xmath95 in the momentum representation this mean field hamiltonian may be written as @xmath96 where @xmath97    the commutator of a bilinear operator product of @xmath39 and @xmath87 with this mean - field hamiltonian is given by @xmath98    & = & - ( \omega_{{{\mathbf p}}_1 } - \omega_{{{\mathbf p}}_2 } )   a^\dagger_{{{\mathbf p}}_1 }   a_{{{\mathbf p}}_2 }   - \sum_{{\mathbf q}}{{\delta \pi}}_{{\mathbf q}}\cdot \frac { ( a_{-{{\mathbf p}}_1-{{\mathbf q } } } +   a_{{{\mathbf p}}_1 + { { \mathbf q}}}^\dagger ) a_{{{\mathbf p}}_2 } } { \sqrt{2 \omega_{{{\mathbf p}}_1 + { { \mathbf q } } } } \sqrt{2\omega_{{{\mathbf p}}_2 } } } \nonumber \\ & & \qquad   + \sum_{{\mathbf q}}{{\delta \pi}}_{{\mathbf q}}\cdot \frac{a_{{{\mathbf p}}_1}^\dagger ( a_{{{\mathbf p}}_2 - { { \mathbf q } } }    + a_{- { { \mathbf p}}_2 + { { \mathbf q}}}^\dagger ) }    { \sqrt{2 \omega_{{{\mathbf p}}_1 } } \sqrt{2\omega_{{{\mathbf p}}_2 - { { \mathbf q } } } } } \ ] ] we show in appendix that the quantum statistical average of this commutator with the gaussian density matrix gives precisely the same result for the same statistical average of the commutator with the original hamiltonian : @xmath99   { \rangle}= { \langle } [ a^\dagger_{{{\mathbf p}}_1 } a_{{{\mathbf p}}_2 } , h ]   { \rangle}\ ] ] therefore one can compute the equations of motion of the wigner functions using this effective hamiltonian , @xmath100 { \rangle}= { \langle } [ a^\dagger_{{{\mathbf p}}+ { { \mathbf k}}/2 } a_{{{\mathbf p}}- { { \mathbf k}}/2 } , h_{\rm mf } ] { \rangle}\ ] ] using this we find , @xmath101 equations of motion of other three wigner functions @xmath102 , @xmath103 , @xmath104 can be also computed from the commution relations of the product operators @xmath105 , @xmath106 , @xmath107 , with the mean field hamiltonian @xmath108 , respectively : we obtain @xmath109 and the equation of motion of @xmath110 can be obtained from ( [ eomf ] ) by the substitution @xmath111 : @xmath112 these equations form a closed system of coupled differential equations with the non - linear klein - gordon equation ( [ nlkg ] ) which may be rewritten as @xmath113    these four equations of motion of the wigner function may be combined into a single matrix form as @xmath114 for a uniform system , we expect that the classical mean field and the self - energy become functions only of time : @xmath115 thus the non - linear klein - gordon equation ( [ nlkg ] ) becomes @xmath116 and the mean - field hamiltonian is reduced to @xmath117    in this case the wigner functions contain non - vanishing components only for the diagonal elements ( @xmath118 ) so that they may be written as @xmath119 then the equations of motion of the wigner functions become @xmath120    these coupled equations have a time - independent solution of the form @xmath121 with @xmath122 thus if we define the particle mass @xmath57 so as to satisfy @xmath123 , then two components of the wigner functions @xmath124 and @xmath125 vanish for all @xmath41 . in this case , @xmath126 may be interpreted as the momentum distribution of physical particle excitations .    to find time - dependent solutions , we write @xmath127 and rewrite the equations as @xmath128 these equations are to be solved together with the non - linear klein - gordon equation ( [ nlkg0 ] ) with @xmath129\ ] ]    from these equations we can derive second order differential equation : @xmath130 with @xmath131\ ] ] we observe that each of these coupled equations looks like an equation of a forced oscillator . for example the frequency of the oscillator @xmath132 is given by @xmath133 while the frequency of the oscillator @xmath134 is given by @xmath135 so that @xmath134 oscillates much rapidly than @xmath132 . the coupling of these oscillators may produce an interesting interesting effects which may be studied by numerical integration . in statistical equilibrium , we have the bose distribution @xmath136 together with @xmath137 where @xmath138 is the inverse temperature . the condition of vanishing @xmath139 implies from ( [ self ] ) @xmath140 where the thermal fluctuation of the quantum field is given by the relation ( [ phi2 ] ) as @xmath141 on the other hand , the condensate amplitude @xmath142 in equilibrium should also satisfy the static non - linear klein - gordon equation , @xmath143\ ] ] from these two conditions , we find , for @xmath144 , @xmath145 and @xmath146 we note that , since the thermal distribution @xmath147 depends on the mass ( gap ) parameter @xmath57 through @xmath148 , the equation ( [ gap ] ) determines the mass parameter @xmath57 self - consistently as a function of the temperature . this equation is called the _ gap equation _  @xcite . as the temperature increases the mass gap @xmath57 and the condensate amplitude @xmath142 decreases and vanishes at high temperature .    to find the behavior of the mass parameter as a function of the temperature we can use the following formula  @xcite for the integral of the bose distribution function : @xmath149 where the dimensionless function @xmath150 is given as @xmath151 ( @xmath152 is euler s number ) . the gap equation ( [ gap ] ) has a solution @xmath153 at @xmath154 determined by @xmath155 however the solution exhibits the behavior of the first order transition due to the non - analytic behavior of the function @xmath156 @xcite . this is a generic feature of the mean field approximation @xcite which may be a theoretical artifact and may not survive the inclusion of correlations missing in the mean field approximation . an improvement of the mean field approximation has been proposed in @xcite and later it has been shown to lead to the second order transition in two - loop approximations @xcite .     for @xmath157 ( dashed line ) and @xmath158 ( solid line ) ]    in fig . 1 we plot @xmath142 as determined by ( [ phi0 ] ) as a function of temperature . it shows a behavior characteristic of the first order phase transition : there is a region @xmath159 where there are three solutions of ( [ phi0 ] ) , one at @xmath160 and other two at @xmath161 as well as @xmath162 . one expects that two solutions @xmath163 and @xmath163 correspond to two local minima of the effective potential @xmath164 while the other one @xmath165 corresponds to the local maximum of @xmath166 . indeed we can show this explicitly by constructing the effective potential @xmath167 by demanding that the local extremal condition @xmath168 coincides with the condition ( [ phi0 ] ) . @xmath169 the integration of the first two terms are trivial . beside irrelevant integration constant , which we may choose to be zero , we obtain @xmath170 to carry out the remaining integral over @xmath171 we demand that the mass parameter @xmath57 depends on @xmath171 through the relation @xmath172 which appears to be equivalent to the condition ( [ mu20 ] ) , but we assume that this relation holds not only for the equilibrium value of @xmath171 , namely at @xmath173 , but also for any value of @xmath171 . noting @xmath174 which results from ( [ mu2 ] ) , we transform the integration variable from @xmath171 to @xmath175 : @xmath176 remaining integral over @xmath175 can be carried out as @xmath177     = - \sum_{{\mathbf p}}\frac{{{\mathbf p}}^2}{3 \omega_{{\mathbf p } } } f_{\rm eq . } ( { { \mathbf p}})\ ] ] where in deriving the final expression we performed integration by part in @xmath178 and omitted again the irrelevant integration constant .    to sum up , we find the following expression for the effective potential : @xmath179 where @xmath171 dependence of the last two terms are implicitly given from the @xmath171-dependence of @xmath57 . we note that this result coincides , in the neglect of the renormalization effects due to the meson mass shift , with the effective potential obtained by amelino - camelia and pi  @xcite using the cornwall - jackiw - tombolis ( cjt ) composite operator effective potential formalism  @xcite , we note that our mass parameter @xmath57 corresponds to the variational mass parameter @xmath180 of the cjt potential . the third term may be interpreted as the pressure of an ideal gas of the quasi - particles obeying the dispersion @xmath181 . the effective potential given by the formula ( [ veff ] ) takes in general a complex value since @xmath175 becomes negative at low temperature and at small value of @xmath171 . at low temperatures the potential develops a second local minimum at non - zero value of @xmath171 , which becomes the absolute minimum below certain temperature @xmath182 where the first order transition takes place in equilibrium and the order parameter @xmath142 jumps discontinuously from @xmath183 to @xmath142 as the temperature is lowered in the presence of inhomogeneity the wigner functions acquire non - vanishing elements with @xmath184 . if we assume that this inhomogeneity are due to the long - wavelength fluctuations in the system , @xmath185 in this case we can obtain a familiar form of the vlasov equation @xcite from the equation of motion of the wigner function by the procedure usually referred to as the gradient expansion @xcite . to show this we make the following approximations : @xmath186 and @xmath187 we also make use of the taylor expansion of the wigner functions : @xmath188    then the equations of motion ( [ eomf ] ) of the wigner function @xmath189 becomes @xmath190 on the right hand side , the first term and the second term are disguised forms of the drift term and the vlasov term respectively .    to obtain more familiar form we make the fourier transforms ( 42 ) and ( 43 ) of the wigner functions . then we find @xmath191 now it is clear that the first term is the drift term which describes the change of the particle position by the drift with the velocity @xmath192 . the second term can be interpreted as the vlasov term which represents the change of particle momentum due to the continuous acceleration by the velocity dependent equivalent potential , @xmath193 acting on the particle . the third term appears to be the correction to the drift term due to the local change of particle mass @xmath194 which causes change in particle velocity . other three terms are associated with the other components of the wigner function and has no counter parts in non - relativistic vlasov equation .    noting that @xmath195 and @xmath196 the above equation may be rewritten in a more compact form : @xmath197 where @xmath198 and @xmath199 the quantity @xmath200 plays the same role as the quasi - particle energy which appears in the kinetic equation of landau s fermi - liquid theory @xcite . we note that because of the relation ( [ gstar - gbar ] ) , @xmath201 is a pure imaginary function , while @xmath202 , @xmath203 , @xmath204 , and @xmath200 are all real functions .    in the long wavelength approximation , the equation of motion of @xmath45 becomes @xmath205 which may be rewritten for the fourier transforms in a compact form as @xmath206 where @xmath207 . thus the kinetic equation for @xmath54 does not look like a vlasov equation ; it takes a form of a simple oscillator equation with the oscillator frequency @xmath208 exposed to the external perturbation created by the particle distribution @xmath209 . similar equation is derived for @xmath210 which is the complex conjugate of @xmath54 . taking the complex conjugate of the eq . ( [ g - eq ] ) , we obtain the kinetic equation for @xmath211 : @xmath212 adding or subtracting ( [ g - eq ] ) and ( [ gbar - eq ] ) , we find for the real function @xmath213 and the pure imaginary function @xmath214 @xmath215    in the long wavelength approximation , we also have @xmath216   \nonumber \\ & = & \sum_{{\mathbf p}}\frac{1}{2\omega_{{\mathbf p } } } \left [ f_+ ( { { \mathbf p } } , { { \mathbf r } } , t ) +   g_+ ( { { \mathbf p } } , { { \mathbf r } } , t ) \right ] \end{aligned}\ ] ] which appears in the non - linear klein - gordon equation ( [ nlkg ] ) for the condensate and in the self energy term ( [ pi ] ) in the vlasov equation . we now apply our coupled kinetic equations , consisting of the non - linear klein - gordon equation and the vlasov equations , to find a dispersion relation of the excitations in the system near equilibrium . the dispersion relations of the excitations of the finite temperature system may be computed using the propagator ( two - time ) formalism , either with real - times ( schwinger - keldysh ) formalism @xcite or with imaginary - times ( matsubara ) formalism with subsequent analytic continuation @xcite . although such calculations have been performed by various authors@xcitenumerical results have been presented focusing on the time - like region of the spectral functions . we study with our kinetic theory the entire range of the @xmath217 plane including the space - like excitations of the system . as has been shown by one of the present authors @xcite the calculation of the zero sound mode in the degenerate fermi liquid using the relativistic landau kinetic equation@xcite , similar to ours , reproduces the same result as the propagator theory in the long wavelength limits . we thus expect that our result may also reproduce the long wavelength behaviors of the excitations obtained from the propagator theory .    for this purpose we assume that the distribution functions consist of uniform equilibrium term and a small deviation from it : @xmath218 where @xmath219 we also assume the mean meson field consists of uniform equilibrium term and a deviation from it @xmath220 where @xmath142 is determined by the static non - linear klein - gordon equation ( [ phi0 ] ) together with the solution of the gap equation ( [ gap ] ) for the mass gap @xmath57 . linearization of the non - linear klein - gordon equation with respect to @xmath221 yields @xmath222    \nonumber \\\end{aligned}\ ] ] where the thermal fluctuation of the quantum field in equilibrium is given by @xmath141 and the deviation @xmath223 is expressed in term of the small deviation of the distribution function from equilibrium value . @xmath224   \nonumber \\   & = &    \sum_{{\mathbf p}}\frac{1}{2\omega_{{\mathbf p } } }   \left [ \delta f ( { { \mathbf p } } , { { \mathbf r } } , t ) + \delta f ( - { { \mathbf p } } , { { \mathbf r } } , t )   + \delta g ( { { \mathbf p } } , { { \mathbf r } } , t ) + \delta g^ * ( { { \mathbf p } } , { { \mathbf r } } , t ) \right ]   \nonumber \\\end{aligned}\ ] ] where in deriving the last line we have used @xmath225 which follow from ( [ fbar ] ) and ( [ gbar ] ) respectively .    by the linearization with respect to the small increments , @xmath226 , @xmath227 , the vlasov equations become @xmath228 the right hand side of the vlasov equation ( [ vlasov1 ] ) vanishes because @xmath229 vanishes in equilibrium . the fluctuation in the mean field potential @xmath230 which appears in these equations is related to the fluctuations of the condensate amplitude and the distribution functions @xmath231 the linearized klein - gordon equation ( [ linkg ] ) and the linearized vlasov equations ( [ lin - vlasov1 ] ) and ( [ lin - vlasov2 ] ) , supplimented by the  constitutive relations " ( [ deltau ] ) and ( [ deltaphi2 ] ) , form a closed set of equations to determine a small fluctuation propagating the system in equilibrium .    noting that @xmath232 and @xmath233 are real functions , we seek a solution in the form : @xmath234 where we have introduced the landau prescription  @xcite , @xmath235 with a positive infinitesimally small constant @xmath236 , to set an adiabatic switching - on of the fluctuation , namely @xmath237 , @xmath238 all vanish slowly as @xmath239 . the constitutive relations for the fluctuations ( [ deltaphi2 ] ) and the mean field potential ( [ deltau ] ) now read @xmath240 and @xmath241 e^ { i ( { { \mathbf k}}\cdot { { \mathbf r}}- \omega_+ t ) }   \nonumber \\ & & \quad + \frac{\lambda}{2\omega_{{\mathbf p } } }   \left [ 2 \phi^*_0 \delta \phi^ * + \sum_{{\mathbf p}}\frac{1}{2\omega_{{\mathbf p } } }   \left ( \delta f^*_{{\mathbf p}}+ \delta f^ * _ { - { { \mathbf p } } } + \delta g^*_{{\mathbf p}}+ \delta g_{-{{\mathbf p } } } \right ) \right ] e^ { - i ( { { \mathbf k}}\cdot { { \mathbf r}}- \omega_- t ) }   \nonumber \\\end{aligned}\ ] ] respectively .    using these relations , we find from the linearized klein - gordon equation ( [ linkg ] ) , @xmath242   \delta \phi    =   - \frac{1}{2 } \lambda \phi_0 \sum_{{\mathbf p}}\frac{1}{2\omega_{{\mathbf p } } }   \left [ \delta f_{{\mathbf p}}+ \delta f_{-{{\mathbf p } } } + \delta g_{{\mathbf p}}+ \delta g^*_{-{{\mathbf p } } } \right ] , \ ] ] and its complex conjugate relation , where we have used the equilibrium relation @xmath243 also , from the linearized vlasov equations we obtain @xmath244   \nonumber \\ & & \qquad \qquad \times \left ( 1 + f_{\rm eq . } ( { { \mathbf p } } ) \right ) f_{\rm eq . } ( { { \mathbf p } } )   \nonumber \\ \label{deltaf}\end{aligned}\ ] ] and @xmath245 f_{\rm eq . } ( { { \mathbf p } } )   \nonumber \\ \label{deltag}\end{aligned}\ ] ] where we have used @xmath246 using ( [ deltaphi ] ) , we eliminate @xmath247 in ( [ deltaf ] ) and ( [ deltag ] ) and find @xmath248 where we have introduced the notation @xmath249 for the deviation of the fluctuation . solving ( [ df ] ) and ( [ dg ] ) for @xmath250 and @xmath251 respectively , and inserting the results into ( [ f ] ) we obtain the relation @xmath252 where @xmath253   ~.   \nonumber \\\end{aligned}\ ] ]    applying the prescription , @xmath254 with @xmath255 implying that the principal part should be taken in integration , we find that @xmath256 consists of the real and imaginary parts : @xmath257 the real part is given by @xmath258 with @xmath259   ~ ,   \nonumber   \\\end{aligned}\ ] ] while the imaginary part is given by @xmath260 with @xmath261   \nonumber   \\   \end{aligned}\ ] ] where the integration can be carried out analytically , yielding @xmath262 for @xmath263 . the values of this function for @xmath264 can be found by noting that it is an odd function of @xmath265 . the first term in ( [ phi2 ] ) corresponds to the space - like ( scattering ) continuum with @xmath266 . while the second term corresponds to the continuum of thermally induced pair creation / annihilation of mesons with energy @xmath267 , the function @xmath268 ( and @xmath269 ) has non - vanishing supports in the regions on @xmath217 plane as indicated by the shaded areas in fig . 2 . we plotted in fig . 3 @xmath270 and @xmath268 as functions of @xmath271 at @xmath272 for two different values of @xmath273 .     and @xmath268 have non - vanishing value is shown by shaded areas ( two - quasi - particle continua ) . the solid hyperbola corresponding to the meson poles ( @xmath274 ) on which @xmath269 has a @xmath275-function singularity in the low temperature phase . the shaded area in the time - like region below the curve @xmath276 is only an artifact of the long wavelength approximation . ]     and @xmath277 as a function of @xmath271 at @xmath278 and @xmath279 ( dashed line ) , @xmath280 ( solid line ) . both @xmath281 and @xmath282 vanish at @xmath283 . @xmath282 also vanishes at the pair creation threshold @xmath284 ( in the long wavelength approximation ) while @xmath281 has a cusp at this point.,title="fig : " ]   and @xmath277 as a function of @xmath271 at @xmath278 and @xmath279 ( dashed line ) , @xmath280 ( solid line ) . both @xmath281 and @xmath282 vanish at @xmath283 . @xmath282 also vanishes at the pair creation threshold @xmath284 ( in the long wavelength approximation ) while @xmath281 has a cusp at this point.,title="fig : " ]    the condition @xmath285 which follows from ( [ omegaf ] ) determines the dispersion relation of a possible long wavelength collective excitation of the system .    before examining the solutions of the dispersion relation ( [ dispersion ] ) . we note that at zero temperature , namely in the absence of quasi - particle excitations , the only solution of our coupled kinetic equations is the one which satisfies the linearized klein - gordon equation ( [ linkg ] ) or ( [ deltaphi ] ) with @xmath286 and @xmath287 . this gives a simple meson pole @xmath288 which appears in the time - like region ( @xmath289 ) . there is no space - like mode of excitations in the absence of quasi - particle excitations due to the mass gap ( @xmath290 ) for the excitations of the vacuum . the situation is different in the case of the ordinary bose - einstein condensate ( bec ) which possesses , even at zero temperature , a low - lying acoustic mode ( bogoliubov phonon ) in the mean - field approximation to repulsive two - body contact interaction@xcite . the ordinary bec contains space - like quasi - particle excitations even at absolute zero temperature which corresponds to the excitation of one of the zero energy particles forming the condensate . the bogoliubov phonon is a collective excitation of such space - like quasi - particle excitation modes modified by the repulsive mean - field interaction .    at finite non - zero temperature , the system contains continuum of the space - like excitations for all @xmath291 satisfying @xmath292 as signified by the non - vanishing value of @xmath293 in addition to the continuum in the time - like region @xmath294 . ( see fig . 2 ) any solution of ( [ dispersion ] ) in the space - like region therefore is subject to the collisionless dissipation known as the landau damping@xcite . we plot the function @xmath295 in fig . 4 as a function of @xmath265 at fixed value of @xmath272 . in the low temperature phase , it contains a singularity at the position of the mesonic pole at @xmath296 while this singularity disappears in the high temperature phase . in the presence of the condensate , a disturbance in the quasi - particle distribution may be absorbed into an excitation of the condensate which then propagates with the mesonic dispersion relation and is converted back again to the quasi - particle excitations . this coupling between the condensate and the quasiparticle excitations generates a collective mode in the low temperature phase.he due to the excitation of the condensate coupled with quasi - particle excitations . ] in the high temperature phase , the quasi - particle excitations couple each other only through their direct interaction , hence no mixing with single mesonic pole . the high temperature behavior of @xmath297 is the same as that of the function @xmath281 : it is negative at @xmath298 and increases with @xmath265 and reaches a positive maximum slightly below @xmath283 . in the low temperature phase , @xmath295 reverts the sign for @xmath265 below this singularity . one may interpret that effective coupling strength @xmath299 changes its sign below the meson pole . note that here we have used the relation @xmath300 ( 111 ) for the condensate amplitude . we note that @xmath301 always vanishes at @xmath283 since @xmath302 . the cusp at @xmath284 appears at the threshold of the two quasiparticle creation . the condition ( [ dispersion ] ) is fulfilled only near the meson pole in the low temperature phase and it gives a shift of the mesonic excitation spectrum . we could not find any additional solution satisfying ( [ dispersion ] ) . we like to note here that without the coupling to the fluctuations to pair creation or annihilation , we would not have gotten the first term in the function @xmath270 and then @xmath303 would have become 1 at @xmath304 creating a undamped tachyonic sound mode . hence the off - diaginal wigner functions plays an important role in making our framework consistent with causality .     as a function of @xmath271 at fixed @xmath305,title="fig : " ]   as a function of @xmath271 at fixed @xmath305,title="fig : " ]    we introduce the `` response function '' defined by @xmath306 =   - \frac { \omega_2 ( \omega , \mu ) } { \left[1 - \omega_1(\omega , \mu ) \right]^2   + \left [ \omega_2 ( \omega , \mu ) \right]^2}\ ] ] it is plotted in fig . 5 as a function of @xmath265 and compared with the `` bare '' response function given by @xmath307 the response function in the space - like momentum region gives the dynamic form factor of the cross section of the scattering of particles coupled to the excitations of the system @xcite , while its time - like component may give the rate of pair annihilation of the quasi - particles . we note that in the low temperature phase the response function changes its sign below the meson pole due to the sign change of the effective coupling ( [ lambda ] ) . we observe some enhancement ( depletion ) of the strength in the time - like pair annihilation near threshold and small depletion ( enhancement ) in the space - like region in the low ( high ) temperature phase . ( [ response ] ) ( solid lines ) at @xmath278 . for comparison the bare response function @xmath308 is also shown ( dashed lines ) . note that in the low temperature phase @xmath309 changes sign in the space - like region and at the meson pole ( @xmath310 ) due to the sign change of the effective coupling strength ( [ lambda]).,title="fig : " ]   ( [ response ] ) ( solid lines ) at @xmath278 . for comparison the bare response function @xmath308 is also shown ( dashed lines ) . note that in the low temperature phase @xmath309 changes sign in the space - like region and at the meson pole ( @xmath310 ) due to the sign change of the effective coupling strength ( [ lambda]).,title="fig : " ]    as we have noted , in the low temperature phase with non - vanishing condensate amplitude @xmath142 , the mesonic excitations couple with the quasi - particle excitations and this give the additional shift of the meson mass from the one obtained by the gap equation . the shift of the meson pole from the bare spectrum @xmath311 may be computed by the condition @xmath312 which in the long wavelength limit ( @xmath313 ) gives a solution at @xmath314 . the shift of the meson mass @xmath315 is determined by @xmath316 where we have used @xmath300 . in the case of @xmath317 this may be solved approximately @xmath318    the equation ( [ mass - shift ] ) may be solved by iteration starting from the first approximation ( [ mass - shift0 ] ) by inserting it for @xmath319 in the integrand on the right hand side . we plot in fig . 6 the meson mass shift determined by this method . the meson mass becomes smaller as @xmath320 increases and eventually becomes zero at the temperature which satisfies the condition @xmath321 the vanishing of the effective mass of mesonic excitation may be identified as the onset of instability in the metastable low temperature phase by small fluctuation ( spinodal decomposition ) . the same instability occurs when one approaches to @xmath322 from high temperature phase . we note that the temperature @xmath323 at the spinodal point given by ( [ spinodal ] ) does not coincide with the `` backbending '' temperature @xmath324 where @xmath325 diverges . the latter temperature is determined by the condition @xmath326 which can be obtained from the gap equation ( [ gap ] ) . we found that @xmath324 is slightly above @xmath323 and the unstable region appears associated with the lower solutions @xmath57 of the gap equation . these unstable solutions thus appear only in the region which is not easily accessible and may well be considered as another pathology of the mean field approximation . in the case of the second order transition of the landau - type . such behavior has been obtained in @xcite by a temperature - dependent loop expansion method formulated in @xcite . ]    ) is the mass of the excitation determined by ( [ mass - shift ] ) and the dashed line is the solution of the gap equation evaluated at @xmath32710 . the slopes of two curves diverge at the temperature @xmath328 while @xmath60 vanishes at @xmath329 . ] in this paper , we have developed a kinetic theory for a system of interacting quantum fields in the mean field approximation taking into account the existence of quasi - particle excitations . we have obtained a set of coupled equations of motion , one for the meson condensate in the form of non - linear klein - gordon equation which contains extra non - linearlity due to the particle excitations . the equations of motion for quasi - particle excitations are described in terms of the wigner functions , which reduce to a semiclassical vlasov equation for one - particle distribution function with a modification due to the coherent pair - creation and pair - annihilation expressed by the off - diagonal components in the wigner function . these off - diagonal components may be eliminated by a suitable bogoliubov transformation of the particle creation and annihilation operators for a uniform , time - independent system . however , they remain non - vanishing in general non - uniform , time - dependent systems .    we have shown that in equilibrium these equations are reduced to a gap equation in the hartree approximation . this implies that our kinetic equations are natural extension of the hartree approximation to the non - equilibrium situation . it is well known , however , that in this approximation the phase transition becomes first order . this is a generic feature of the mean field approximation  @xcite which persists also for multi - components scalar fields with continuous @xmath1 symmetry for arbitrary finite integer @xmath330 .    in the present work , we studied also the excitations spectrum in the system near equilibrium . in the high temperature phase there is no collective excitation mode in the system besides the two continua of the quasi - particle excitations in the entire space - like energy - momentum region and the time - like region above the pair creation threshold . in the low temperature phase we found that coupling of the meson pole to the quasi - particle continua give rise to the shift of the meson mass which becomes zero at the edge of the spinodal instability line of the first order transition . the spinodal point appears deep inside the meta - stable region of the first order transition which may well be an artifact of the mean field approximation . we found that the coupling to the off - diagonal components of the wigner functions should be properly included to avoid the appearance of the undamped tachyonic mode .    in this work we used a single component real scalar field model which possesses only discrete symmetry . it is straightforward to extend the present analysis to models with continuous symmetry such as the sigma model with @xmath1 symmetry . basic features of the present analysis is preserved in such extension . it is known however that the goldstone theorem is apparently violated in this approximation for a system with continuous symmetry . we will show in the forthcoming paper  @xcite , that the missing nambu - goldstone mode may be retrieved in the collective excitations of the system . our coupled kinetic equations may be solved for an arbitrary initial conditions . we plan to study the freezeout dynamics with these equations with more realistic interactions . it would be interesting to see in particular how much flow is generated by the acceleration by the mean field as the vacuum condensate is restored . 0.5 cm    * acknowledgment *    we thank gordon baym , hirotsugu fujii , tetsuo hatsuda , tetsufumi hirano , osamu morimatsu and koichi ohta for helpful discussions . we are indebted to gordon baym for correcting a sign error in the original manuscript and to tetsuo hatsuda and osama morimatsu for calling our attention to some works on the spectral functions at finite temperature which prompted us to extend our analyses of the excitation spectrum in section 6 . we are also grateful to francois gelis , kazunori itakura , yusuke kato , and tetsuro nikuni for their interests in this work . this work has been supported by the grants - in - aid of the japanese ministry of education , culture , sports , science and technology no . 13440067 and 19540269 . here we present some details of the computation of the equation of motion of the wigner functions and show that with the gaussian ansatz for the density matrix the result is equivalent to what we obtain from the mean - field hamiltonian ( [ hmf ] ) . we like to compute the time - derivative of the operator product @xmath331 which appears in the definition of @xmath332 with @xmath333 and @xmath334 . @xmath335 a_{{{\mathbf p}}_2 }   +   a^\dagger_{{{\mathbf p}}_1 } [ a_{{{\mathbf p}}_2 } , h ] \end{aligned}\ ] ] the commutators @xmath336 $ ] is decomposed into a sum of the commutators with @xmath85 ( @xmath337 ) among which the commutators with @xmath89 and @xmath90 do not survive the average with the gaussian density matrix since they only contain the odd power of the field operators . we only need to compute the commutators with @xmath91 and @xmath92 .    to compute the commutator with @xmath91 it is convenient to rewrite @xmath91 as @xmath338 we then find @xmath339 & = & \omega_{{{\mathbf p } } } a_{{\mathbf p}}(t ) +   \frac{\lambda}{4 } \int d { { \mathbf r}}\phi_c^2 ( { { \mathbf r } } , t )   [ a_{{\mathbf p}}(t ) , { { \tilde \phi}}^2 ( { { \mathbf r } } , t ) ]   \nonumber \\ & = & \omega_{{{\mathbf p } } } a_{{\mathbf p}}(t ) + \frac{\lambda}{2 } \int d { { \mathbf r}}\phi_c^2 ( { { \mathbf r } } , t )   \frac{e^{- i { { \mathbf p}}\cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p } } } } } { { \tilde \phi}}({{\mathbf r } } , t ) \\ \left [ a_{{\mathbf p}}^\dagger ( t ) , h_2 \right ] & = & - \omega_{{{\mathbf p } } } a_{{\mathbf p}}^\dagger ( t ) +    \frac{\lambda}{4 } \int d { { \mathbf r}}\phi_c^2 ( { { \mathbf r } } , t )   [ a_{{\mathbf p}}^\dagger ( t ) , { { \tilde \phi}}^2 ( { { \mathbf r } } , t ) ]   \nonumber \\ & = & - \omega_{{{\mathbf p } } } a_{{\mathbf p}}(t ) - \frac{\lambda}{2 } \int d { { \mathbf r}}\phi_c^2 ( { { \mathbf r } } , t )   \frac{e^ { i { { \mathbf p}}\cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p } } } } } { { \tilde \phi}}({{\mathbf r } } , t ) \end{aligned}\ ] ] where we have used the following formulae : @xmath340   & = & n { { \tilde \phi}}^{n-1 }   ( { { \mathbf r } } , t ) [ a_{{\mathbf p}}^\dagger , { { \tilde \phi}}({{\mathbf r}},t ) ]   = -   n { { \tilde \phi}}^{n-1 }   ( { { \mathbf r } } , t ) \frac{e^ { i { { \mathbf p}}\cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p } } } } } \\   \left [ a_{{\mathbf p } } , { { \tilde \phi}}^n ( { { \mathbf r } } , t ) \right ]   & = & n \phi^{n-1 }   ( { { \mathbf r}},t ) [ a_{{\mathbf p } } , { { \tilde \phi}}({{\mathbf r}},t ) ]   = n \phi^{n-1 }   ( { { \mathbf r}},t ) \frac{e^{- i { { \mathbf p}}\cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p } } } } } \end{aligned}\ ] ] with these results we obtain @xmath341   & = & - ( \omega_{{{\mathbf p}}_1 } - \omega_{{{\mathbf p}}_2 } ) a^\dagger_{{{\mathbf p}}_1 }   a_{{{\mathbf p}}_2 }   \nonumber \\ & & \quad - \frac{\lambda}{2 } \int d { { \mathbf r}}\phi_c^2 ( { { \mathbf r } } , t ) { { \tilde \phi}}({{\mathbf r } } , t )   \frac{e^ { i { { \mathbf p}}_1 \cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p}}_1 } } }   a_{{{\mathbf p}}_2 }    \nonumber \\ & & \quad \qquad   + \frac{\lambda}{2 } \int d { { \mathbf r}}\phi_c^2 ( { { \mathbf r } } , t ) { { \tilde \phi}}({{\mathbf r } } , t ) \frac{e^{- i { { \mathbf p}}_2 \cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p}}_2 } } } a_{{{\mathbf p}}_1}^\dagger \end{aligned}\ ] ] the gaussian average of this equation gives @xmath342 { \rangle } & = & - ( \omega_{{{\mathbf p}}_1 } - \omega_{{{\mathbf p}}_2 } ) { \langle}a^\dagger_{{{\mathbf p}}_1 }   a_{{{\mathbf p}}_2 } { \rangle}\nonumber \\ & & \quad - \frac{\lambda}{2 } \int d { { \mathbf r}}\phi_c^2 ( { { \mathbf r } } , t )    \frac{e^ { i { { \mathbf p}}_1 \cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p}}_1 } } } { \langle}{{\tilde \phi}}({{\mathbf r } } , t ) a_{{{\mathbf p}}_2 } { \rangle}\nonumber \\ & & \quad \qquad   + \frac{\lambda}{2 } \int d { { \mathbf r}}\phi_c^2 ( { { \mathbf r } } , t )   \frac{e^{- i { { \mathbf p}}_2 \cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p}}_2 } } } { \langle}{{\tilde \phi}}({{\mathbf r } } , t ) a_{{{\mathbf p}}_1}^\dagger { \rangle}\end{aligned}\ ] ]    next we compute the commutators with @xmath92 . we first compute @xmath343 & = & \frac{\lambda}{3 ! } \int d { { \mathbf r}}{{\tilde \phi}}^3 ( { { \mathbf r } } , t )   \frac{e^ { i { { \mathbf p}}\cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p } } } } } a_{{\mathbf p}}(t ) \\ \left [ a_{{\mathbf p}}^\dagger ( t ) , h_4 \right ] & = & - \frac{\lambda}{3 ! } \int d { { \mathbf r}}{{\tilde \phi}}^3 ( { { \mathbf r } } , t )   \frac{e^ { i { { \mathbf p}}\cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p } } } } } a_{{\mathbf p}}(t ) \end{aligned}\ ] ] we then use these results to obtain @xmath344   & = & - \frac{\lambda}{3 ! } \int d { { \mathbf r}}{{\tilde \phi}}^3 ( { { \mathbf r } } , t ) \frac{e^ { i { { \mathbf p}}_1 \cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p}}_1 } } }   a_{{{\mathbf p}}_2 }    \nonumber \\ & & \qquad + \frac{\lambda}{3 ! } \int d { { \mathbf r}}{{\tilde \phi}}^3 ( { { \mathbf r } } , t ) \frac{e^{- i { { \mathbf p}}_2 \cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p}}_2 } } } a_{{{\mathbf p}}_1}^\dagger \end{aligned}\ ] ] taking the average with the gaussian density matrix , we find @xmath345 { \rangle } & = & - \frac{\lambda}{2 } \int d { { \mathbf r}}{\langle}{{\tilde \phi}}^2 ( { { \mathbf r } } , t ) { \rangle}\frac{e^ { i { { \mathbf p}}_1 \cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p}}_1 } } }   { \langle}{{\tilde \phi}}({{\mathbf r } } , t ) a_{{{\mathbf p}}_2 }   { \rangle}\nonumber \\ & & \qquad + \frac{\lambda}{2 } \int d { { \mathbf r}}{\langle}{{\tilde \phi}}^2 ( { { \mathbf r } } , t ) { \rangle}\frac{e^{- i { { \mathbf p}}_2 \cdot { { \mathbf r}}}}{\sqrt{2 \omega_{{{\mathbf p}}_2 } } } { \langle}{{\tilde \phi}}({{\mathbf r } } , t ) a_{{{\mathbf p}}_1}^\dagger { \rangle}\end{aligned}\ ] ]    we observe the similarity between the last two terms of the commutator with @xmath91 and the commutator with @xmath92 after gaussian average . these terms can be combined by introducing the self - energy function , @xmath346 integration over the space coordinate @xmath347 give a delta function @xmath348 and @xmath349 . performing the integral over @xmath41 yields @xmath350 { \rangle } & = & - ( \omega_{{{\mathbf p}}_1 } - \omega_{{{\mathbf p}}_2 } ) { \langle}a^\dagger_{{{\mathbf p}}_1 }   a_{{{\mathbf p}}_2 } { \rangle}\nonumber \\ & & \quad   - \frac{\lambda}{4 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt{\omega_{{{\mathbf p}}_1 + { { \mathbf q}}}\omega_{{{\mathbf p}}_2 } } }    ( { \langle}a_{-{{\mathbf p}}_1-{{\mathbf q } } } a_{{{\mathbf p}}_2 } { \rangle}+ { \langle}a_{{{\mathbf p}}_1 + { { \mathbf q}}}^\dagger a_{{{\mathbf p}}_2 } { \rangle } ) \nonumber \\ & & \qquad   + \frac{\lambda}{4 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt { \omega_{{{\mathbf p}}_1 } \omega_{p_2 - { { \mathbf q } } } } }    ( { \langle}a_{{{\mathbf p}}_2 - { { \mathbf q } } }   a_{{{\mathbf p}}_1}^\dagger { \rangle}+    { \langle}a_{- { { \mathbf p}}_2 + { { \mathbf q}}}^\dagger a_{{{\mathbf p}}_1}^\dagger { \rangle } )   \nonumber \\\end{aligned}\ ] ] this result coincides with the average of the commutator with the mean field hamiltonian @xmath108 defined by ( [ hmf ] ) . @xmath351 { \rangle}= { \langle } [ a^\dagger_{{{\mathbf p}}_1 } a_{{{\mathbf p}}_2 } , h_{\rm mf } ] { \rangle}\ ] ]    commutators of the four bilinear operator products with the mean - field hamiltonian are listed below : @xmath98    & = & - ( \omega_{{{\mathbf p}}_1 } - \omega_{{{\mathbf p}}_2 } )   a^\dagger_{{{\mathbf p}}_1 }   a_{{{\mathbf p}}_2 }   \nonumber \\ & & \quad   - \frac{1}{2 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt{\omega_{{{\mathbf p}}_1 + { { \mathbf q}}}\omega_{{{\mathbf p}}_2 } } }    (   a_{-{{\mathbf p}}_1-{{\mathbf q } } } a_{{{\mathbf p}}_2 }   +   a_{{{\mathbf p}}_1 + { { \mathbf q}}}^\dagger a_{{{\mathbf p}}_2 }   ) \nonumber \\ & & \qquad   + \frac{1}{2 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt { \omega_{{{\mathbf p}}_1 } \omega_{{{\mathbf p}}_2 - { { \mathbf q } } } } }    ( a_{{{\mathbf p}}_1}^\dagger   a_{{{\mathbf p}}_2 - { { \mathbf q } } }    +     a_{{{\mathbf p}}_1}^\dagger a_{- { { \mathbf p}}_2 + { { \mathbf q}}}^\dagger   )    \nonumber \\    \\ \left [ a_{{{\mathbf p}}_1 } a_{{{\mathbf p}}_2 } , h_{\rm mf } \right ]   & = & ( \omega_{{{\mathbf p}}_1 } + \omega_{{{\mathbf p}}_2 } )   a_{{{\mathbf p}}_1 }   a_{{{\mathbf p}}_2 }   \nonumber \\ & & \quad   + \frac{1}{2 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt{\omega_{{{\mathbf p}}_1 + { { \mathbf q}}}\omega_{{{\mathbf p}}_2 } } }    (   a_{{{\mathbf p}}_1-{{\mathbf q } } }   +   a_{-{{\mathbf p}}_1 + { { \mathbf q}}}^\dagger )   a_{{{\mathbf p}}_2 }    \nonumber \\ & & \qquad   + \frac{1}{2 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt { \omega_{{{\mathbf p}}_1 } \omega_{{{\mathbf p}}_2 - { { \mathbf q } } } } }    a_{{{\mathbf p}}_1 } (   a_{{{\mathbf p}}_2 - { { \mathbf q } } }    + a_{- { { \mathbf p}}_2 + { { \mathbf q}}}^\dagger   )   \nonumber \\   \\ \left [ a^\dagger_{{{\mathbf p}}_1 } a^\dagger_{{{\mathbf p}}_2 } , h_{\rm mf } \right ]    & = & - ( \omega_{{{\mathbf p}}_1 } + \omega_{{{\mathbf p}}_2 } )   a^\dagger_{{{\mathbf p}}_1 }   a_{{{\mathbf p}}_2 }   \nonumber \\ & & \quad   - \frac{1}{2 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt{\omega_{{{\mathbf p}}_1 + { { \mathbf q}}}\omega_{{{\mathbf p}}_2 } } }    (   a_{-{{\mathbf p}}_1-{{\mathbf q } } } a_{{{\mathbf p}}_2 }   +   a_{{{\mathbf p}}_1 + { { \mathbf q}}}^\dagger a_{{{\mathbf p}}_2 }   ) \nonumber \\ & & \qquad   - \frac{1}{2 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt { \omega_{{{\mathbf p}}_1 } \omega_{{{\mathbf p}}_2 - { { \mathbf q } } } } }    (   a_{{{\mathbf p}}_2 - { { \mathbf q } } }   a_{{{\mathbf p}}_1}^\dagger   +     a_{- { { \mathbf p}}_2 + { { \mathbf q}}}^\dagger a_{{{\mathbf p}}_1}^\dagger )   \nonumber \\   \\   \left [ a_{{{\mathbf p}}_1 } a^\dagger_{{{\mathbf p}}_2 } , h_{\rm mf } \right ]    & = & ( \omega_{{{\mathbf p}}_1 } - \omega_{{{\mathbf p}}_2 } )   a_{{{\mathbf p}}_1 }   a^\dagger_{{{\mathbf p}}_2 }   \nonumber \\ & & \quad   + \frac{1}{2 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt{\omega_{{{\mathbf p}}_1 + { { \mathbf q}}}\omega_{{{\mathbf p}}_2 } } }    (   a_{{{\mathbf p}}_1-{{\mathbf q } } } +   a_{- { { \mathbf p}}_1 + { { \mathbf q}}}^\dagger ) a^\dagger_{{{\mathbf p}}_2 }    \nonumber \\ & & \qquad   - \frac{1}{2 } \sum_{{\mathbf q}}\frac{\pi_{{\mathbf q}}}{\sqrt { \omega_{{{\mathbf p}}_1 } \omega_{{{\mathbf p}}_2 - { { \mathbf q } } } } }    a_{{{\mathbf p}}_1}^\dagger    (   a_{{{\mathbf p}}_2 - { { \mathbf q } } }   + a_{- { { \mathbf p}}_2 + { { \mathbf q}}}^\dagger )   \nonumber \\\end{aligned}\ ] ]    # 1ann . phys . * # 1 * # 1phys . rep . * # 1 * # 1act . phys . pol . * # 1 * # 1prog . phys . * # 1 * # 1prog . phys . supplement*#1 * # 1phys . rev . * # 1 * # 1phys . rev . * d#1 * # 1phys . rev . * c#1 * # 1phys . rev . * a#1 * # 1phys . * # 1 * # 1phys . # 1 * # 1rev . phys.*#1 * # 1nucl . phys . * # 1 * # 1nuovo cimento * # 1 * # 1sov . # 1 * # 1eur . j. * a#1 * # 1coll . a#1 *      j. d. bjorken , _ phys . rev . _ * d27 * , 865 ( 1983 ) ; g. baym , b. l. friman , j .- blaizot , m. soyeur , w. czy , _ nucl . _ * a407 * , 541 ( 1983 ) ; k. kajantie and l. mclerran , , 261 ( 1983 ) ; m. gyulassy , t. matsui , _ phys . rev . _ * d29 * , 419 ( 1984 )
Abstract:

Answer: we formulate a kinetic theory of self - interacting meson fields with an aim to describe the freezeout stage of the space - time evolution of matter in ultrarelativistic nuclear collisions . kinetic equations are obtained from the heisenberg equation of motion for a single component real scalar quantum field taking the mean field approximation for the non - linear interaction . the mesonic mean field obeys the classical non - linear klein - gordon equation with a modification due to the coupling to mesonic quasi - particle excitations which are expressed in terms of the wigner functions of the quantum fluctuations of the meson field , namely the statistical average of the bilinear forms of the meson creation and annihilation operators . in the long wavelength limit , the equations of motion of the diagonal components of the wigner functions take a form of vlasov equation with a particle source and sink which arises due to the non - vanishing off - diagonal components of the wigner function expressing coherent pair - creation and pair - annihilation process in the presence of non - uniform condensate . we show that in the static homogeneous system , these kinetic equations reduce to the well - known gap equation in the hartree approximation , and hence they may be considered as a generalization of the hartree approximation method to non - equilibrium systems . as an application of these kinetic equations , we compute the dispersion relations of the collective mesonic excitations in the system near equilibrium . meson condensate ; vlasov equation