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: the formation of a quark - gluon plasma ( qgp ) is achieved in ultra - relativistic heavy - ion collision experiments at relativistic heavy - ion collider ( rhic ) in bnl and at the large hadron collider ( lhc ) in cern . more and more information on properties of strongly interacting quark - gluon matter in quantum chromodynamics ( qcd ) is being accumulated . however , despite a number of interesting discoveries and new notions established in experiments , the thermalization mechanism that leads to the qgp still lacks a deep understanding based on qcd first principles . there have been many attempts to estimate the thermalization time scale , that is phenomenologically known to be @xmath0 from analysis in hydrodynamic models  @xcite . it was argued  @xcite that perturbative qcd processes lead to a thermalization time of order @xmath1 where the saturation momentum @xmath2 is a characteristic hard scale in the high - energy reactions ; @xmath3-@xmath4 in the typical rhic kinematics . if the prefactor is of order unity , this estimate would yield a thermalization time around several @xmath5 , which is too slow to account for the rhic observations . it has also been suggested that qcd perturbation theory , no matter at which order one truncates it , may be unable to describe thermalization at all  @xcite . plasma instabilities , particularly the weibel instability  @xcite , could be promising candidates for the early thermalization mechanism  @xcite and have been well investigated in the hard - loop formalism  @xcite . there are also various other approaches  @xcite including a framework based on the ads / cft correspondence  @xcite , in which thermalization appears as the formation of a black hole horizon in the dual theory . parallel to the analysis on plasma instabilities , an important observation was made in a related theoretical framework . the occupation number of gluons at high energy is so huge that color field amplitudes @xmath6 can be realized , and then non - linearities must be taken into account to all orders . instead of a plain perturbative expansion , then , an approximation in terms of classical field equations becomes more suitable because it effectively makes a resummation of large - amplitude fields @xmath7 . this non - linear regime , whose natural time scale is of order @xmath8 at rhic , is expected to be prevalent during the early stages in heavy - ion collisions . the qcd framework known as the color glass condensate ( cgc )  @xcite rearranges the perturbative expansion in terms of classical fields ( that correspond to the resummation of infinite sets of graphs ) , and perturbations on top of these classical fields are of higher order in @xmath9 . the cgc was first developed in order to deal with the saturation effect in the gluon distribution function  @xcite . because we are interested in the early time ( @xmath10 ) dynamics before thermalization , the cgc description seems to be the most suited for our purpose . the initial condition just above the forward light - cone ( i.e.  at @xmath11 ) for the classical field that appears at lowest order in the cgc description can be nicely formulated in the bjorken ( expanding ) coordinates  @xcite . because in the regime of large occupation numbers , classical field theory and kinetic theory essentially describe the same dynamics  @xcite , it is conceivable that the cgc can describe the isotropization and thermalization processes qualitatively . however , for this program to have a chance to work , it is necessary to go beyond the lowest order description , since at zeroth order it is well known that the pressure tensor never isotropizes . it is unfortunately impossible to obtain an analytical solution of the classical yang - mills equations of motion for the heavy - ion collision ( two - source ) problem , and therefore one has to resort to numerical simulations . in refs . @xcite , heavy - ion collisions were simulated in the mclerran - venugopalan ( mv ) model ( i.e.  the cgc at lowest order , with a local gaussian distribution of color sources ) , formulated gauge - invariantly in terms of link variables . these numerical calculations were later confirmed by an independent simulation in ref . the peculiar feature in the cgc description at leading order is that the initial condition for the classical fields is independent of the space - time rapidity variable @xmath12 , and the equations of motion maintain @xmath12-independence during the time evolution . therefore , the space - time evolution is entirely independent of @xmath12 , that is , invariant under boosts in the direction of the collision axis . intuitively , one can expect a longitudinally extended structure of the initial fields , which is indeed a characteristic property of the so - called `` glasma '' state  @xcite , which could be a source for the so - called ridge structure seen in the two - particle correlations at rhic and lhc  @xcite .    within the cgc framework , an instability has been found when the boost invariant glasma fields are disturbed  @xcite . this instability has some similarities to the weibel instability in anisotropic plasmas , but the origin of the instability is not yet fully clarified , though several interpretations have been proposed  @xcite . also , it is still an issue how to properly formulate disturbances to the boost invariant cgc background . however , it is clear that such perturbations arise when one considers higher order corrections to the leading cgc picture  @xcite , and therefore they should not be ignored since they can alter phenomena such as particle production  @xcite , etc . moreover , without including these perturbations , one can not quantify the thermalization processes . one of the reasons why the origin of the glasma instability lacks a clear understanding is that the available numerical data from earlier works  @xcite is rather limited . this motivates us to carry out a more systematic survey of the glasma instability , exploiting a wider range of simulation parameters . in this work we will specifically investigate the dependence of the glasma instability on the following parameters :    * [ i ] * _ transverse grid size @xmath13 : _   it is known that the initial energy density at @xmath11 is sensitive to the ultraviolet and infrared regulators in the transverse direction , and this sensitivity becomes milder at finite @xmath14  @xcite . it is therefore important to check whether continuum - limit results can be obtained reliably by looking systematically at the dependence on the transverse lattice spacing @xmath15 , where @xmath16 is the length of the box and @xmath13 the number of lattice points in the transverse direction .    * [ ii ] * _ longitudinal grid size @xmath17 : _   physical results should also not be contaminated by ultraviolet singularities in the longitudinal direction . to avoid this , we will use initial fluctuations with sufficiently small longitudinal wave - number @xmath18 ( @xmath18 is the fourier conjugate to the rapidity @xmath12 ) . we will confirm that there is no dependence on the longitudinal grid size @xmath17 or the longitudinal lattice spacing @xmath19 , at least until the longitudinal spectrum extends into the ultraviolet regions at very late time . * [ iii ] * _ magnitude @xmath20 of the seeds : _   in general the strength of the instability increases with the magnitude @xmath20 of the initial @xmath12-dependent perturbations . this dependence is found to be very simple until the unstable modes become large enough to be self - interacting : at moderate times , the finite-@xmath18 modes obey the linearized dynamics and their amplitudes grow linearly as @xmath20 increases . * [ iv ] * _ initial wave - number @xmath21 : _   the time evolution depends on the initial condition in the longitudinal direction , and in particular on the wave - number @xmath21 of the seeds . we will show that the instability produced by the smallest wave - number is defined to be an integer between @xmath22 and @xmath23 . the lowest non - zero wave - number is thus @xmath24 . ] ( @xmath25 ) is the fastest and strongest in a wide range of times . we will also confirm the dominance of linearity and the appreciable presence of residual non - linear effects by verifying that the result from an initial condition with multiple @xmath21 s can be approximated as the superposition of the results from initial conditions with single wave - numbers @xmath21 s .    in this paper , we report on numerical results about these dependences for the growing behavior at the late stage of the glasma evolution . these numerical results provide a pool of observations that we hope will help to understand the microscopic mechanisms that drive the glasma instabilities . in the section  [ sec : formulas ] , we recall the most important formulas for a lattice study of this problem . in the section  [ sec : bi ] , we briefly expose some results in the purely boost - invariant case . these are well - known results already , and merely serve as a check of our implementation . in the section  [ sec : notbi ] , we describe what happens when the glasma fields are perturbed by a rapidity - dependent seed , and we investigate thoroughly how the result depends on the parameters listed above . we explore the time evolution of the @xmath18-spectrum in the section  [ sec : lspectrum ] . finally , the section  [ sec : summary ] is devoted to a summary and outlook . in this section , we summarize our equations , with emphasis on the lattice specificities . the main goal of this section is to avoid confusions due to the existence of different conventions in the literature for some quantities . the link variables are defined as @xmath26 with the gauge potential @xmath27 that connects the neighboring sites from @xmath28 and @xmath29 . this expression for the link variable is the most important one in order to establish the correspondence between the continuum and the lattice ( link - variable ) formulations . under a gauge transformation by @xmath30 the link variables transform as @xmath31 and the chromo - electric fields as @xmath32 the plaquette variable is then defined as @xmath33 \ ; , \label{eq : approxf}\ ] ] where the last form is an approximation which gives rise to the field strength @xmath34 $ ] of the continuum formulation . this is a very useful expression to study the continuum limit . the lattice spacing should be the longitudinal one , i.e.  @xmath35 , which is distinct from the transverse one @xmath36 , if the longitudinal direction is involved in the displacements @xmath37 or @xmath38 in the above formula . the canonical momenta in terms of link variables are expressed as @xmath39 here we note that in the above we already use variables made dimensionless by the transverse lattice spacing @xmath36 , and thus @xmath36 will never appear explicitly in our formula . because @xmath12 is a dimensionless number , @xmath35 does not bring any mass dimension unlike @xmath36 . we can also discretize the time in the above equations in order to solve the time evolution numerically for @xmath40 and @xmath41 ; @xmath42 u_i(\tau ) \;,\\   u_\eta(\tau '' ) & = \exp\bigl [ -2{\delta\tau}\cdot { \mathrm{i}}g a_\eta    \tau ' e^\eta(\tau ' ) \bigr ] u_\eta(\tau ) \;,\end{aligned}\ ] ] where @xmath43 and @xmath44 . note that we have not written explicitly the position arguments in the quantities that appear in these equations . here , the exponentiation of @xmath45 is important ; otherwise the updated @xmath46 would not be a unitary matrix . the equations of motion are discretized in the same manner as @xmath47_\tau \\   & \qquad + 2{\delta\tau}\frac{{\mathrm{i}}\tau}{2g}\ ,    \sum_{j\neq i } \bigl [ u_{ji}(x)+u_{-ji}(x ) - \text{(h.c . ) } \bigr]_\tau \ ; , \end{split}\ ] ] for the transverse components and @xmath48_\tau \;,\ ] ] for the longitudinal component . we recall again that all variables including @xmath14 and @xmath49 above are dimensionless , and expressed in unit of the transverse spacing @xmath36 . one can make sure that these equations are equivalent to the ordinary equations of motion in the continuum limit , by using the approximate form  . the initial conditions are derived from the requirement of avoiding the singularity at the collision point . it is written in terms of @xmath50 , the classical solution of the mv model for a single color source ( the label @xmath51 indicates which of the two nuclei produces the corresponding field ) . this single - source classical solution is a transverse pure - gauge , @xmath52 where the gauge rotation matrix @xmath53 can be written as @xmath54 \;. \label{eq : wilson}\ ] ] here , @xmath55 is the solution of the poisson equation , @xmath56 . in the mv model , the random color source @xmath57 is gaussian distributed , with a two - point correlation @xmath58 higher - point correlation functions are obtained as products of two - point contractions , via wick s theorem . in the mv model , @xmath59 is the only dimensionful parameter and it is related to the saturation momentum @xmath2 . then , the initial conditions are  @xcite ; @xmath60 \ ; , \label{eq : initial_e}\end{aligned}\ ] ] written here in the slightly simpler case of the su(2 ) color group . these discretized equations of motion and initial conditions completely define the numerical mv model . our conventions are mostly consistent with those of previous works , except for the fact that the field @xmath61 is not treated as an adjoint scalar - field variable but instead we also describe it in terms of a link variable @xmath62 . we note that our formulation is reduced to the conventional one with the adjoint scalar - field variable if we take @xmath35 to be sufficiently small . in addition , since our interest lies primarily in the qualitative behavior of the glasma instability , we disregard the color neutralization that has been imposed onto the mv model color sources in other works  @xcite . it is now straightforward to implement the numerical calculation for gauge - invariant observables such as the components of the energy - momentum tensor . they are in the continuum convention written down as @xmath63\bigr\rangle \ ; , \label{eq : e } \\   & p_{_t } = \frac{1}{2}\bigl\langle t^{xx } + t^{yy } \bigr\rangle    = \bigl\langle { \mathrm{tr}}\bigl [ e_{_l}^2 + b_{_l}^2 \bigr ] \bigr\rangle \ ; , \label{eq : p_{_t } } \\   & p_{_l } = \bigl\langle \tau^2 t^{\eta\eta } \bigr\rangle    = \bigl\langle{\mathrm{tr}}\bigl [ e_{_t}^2 + b_{_t}^2 - e_{_l}^2 - b_{_l}^2      \bigr]\bigr\rangle \ ; , \label{eq : p_{_l}}\end{aligned}\ ] ] where we define @xmath64 these formulas are implicitly summed over the color index @xmath36 , for @xmath65 and @xmath66 to be gauge invariant . here , the chromo - magnetic field squared should be expressed in a way consistent with eq .  , that is , @xmath67 these definitions are , however , not very convenient when we decompose the energy into its fourier components , as we will see later . we confirm that the energy - momentum tensor is traceless ; @xmath68 . it should also be mentioned that our definition of @xmath69 is different from that in refs . @xcite by a factor @xmath70 so that @xmath71 for an isotropic system . we first check the consistency with the known results in the boost - invariant case , i.e.  in @xmath72 dimensions . this is useful also for later discussions on the spectral decomposition when we consider @xmath12-dependent fluctuations . in this section , we focus on @xmath12-independent classical solutions , with initial conditions specified in eqs .   and that are themselves independent of @xmath12 . eventually , however , our goal is to superimpose @xmath12-dependent fluctuations to these boost invariant initial conditions in order to incorporate quantum effects  @xcite . we note here that our numerical calculations are limited to the su(2 ) color group , so that we can get more samples and thus better statistics for the same computational cost . as discussed in ref . @xcite , the su(2 ) calculations capture the essential properties of the glasma instability , and we expect that the results we report here for the su(2 ) case would be qualitatively the same for su(3 ) . let us first take a look at a typical fixed initial configuration , without taking the ensemble average over the color sources , in order to develop an intuitive picture for the glasma state . in the mv model , the initial color distribution in the transverse plane is akin to white noise , having no correlation between different transverse sites . ] . since the solution of the poisson equation @xmath73 involves a convolution with a massless ( i.e.  long ranged ) transverse propagator , @xmath74 has a smooth structure as shown in the left panel of fig . [ fig : lam ] .    . the mv - model parameters are chosen as @xmath75 ( number of transverse sites ) and @xmath76 . left : the solution @xmath74 of the poisson equation . only the first color component @xmath77 is plotted . right : corresponding gauge field @xmath78 calculated from @xmath79 $ ] . only the first color component @xmath80 is plotted.,title="fig:",scaledwidth=49.0% ] . the mv - model parameters are chosen as @xmath75 ( number of transverse sites ) and @xmath76 .  left : the solution @xmath74 of the poisson equation . only the first color component @xmath77 is plotted . right : corresponding gauge field @xmath78 calculated from @xmath79 $ ] . only the first color component @xmath80 is plotted.,title="fig:",scaledwidth=49.0% ]    although @xmath74 is fairly smooth due to the regularizing effect of this convolution , the corresponding gauge field @xmath81 , which is a phase of @xmath82 defined in eq . , is far from smooth . this is actually a consequence of the pure - gauge form ; contributions from globally smooth parts in @xmath74 cancel each other and the resulting @xmath81 fluctuates strongly . in other words , @xmath81 has a rough structure because a derivative of @xmath74 undoes the regularization provided by the convolution . the @xmath81 shown in the figure [ fig : lam ] is the gauge field from a single source only . in a collision , the initial fields such as @xmath83 and @xmath84 involve a product of two @xmath81 s from two independent sources , and thus have an even rougher structure . it is worth mentioning that the `` color - flux tube '' structures , whose transverse size is expected to be of order @xmath85 , do not emerge manifestly from the mv model since the source distribution has no correlation length in the transverse direction . this means that the nielsen - olesen instability  @xcite is unlikely to occur at least within the mv model . nevertheless , we expect that color - flux tube structures would appear after the jimwlk evolution from the mv - model distribution . in the mean - field approximation , in fact , the transverse correlation function of color source has been already studied in ref . @xcite , which should be taken into account in the future works . also , to make a clear comparison with the weibel instability and the hard - loop estimates ( especially early works  @xcite ) , it would be useful to simplify the mv - model setup in ( 1 + 1 ) dimensions without transverse dependences in such a way not to lose the qualitative features of the glasma instability that we will see in the section  [ sec : notbi ] . in this sense the results reported in the present paper should serve as a foundation for future investigations in these directions .      . an ensemble average is taken over @xmath86 random configurations of @xmath87 for @xmath88 and @xmath89 , and over @xmath90 configurations for @xmath91 , @xmath92 , @xmath93 and @xmath94 . the mv model parameter is chosen as @xmath76 . right : chromo - electric and chromo - magnetic fields   and for the same mv model parameters and @xmath88.,title="fig:",scaledwidth=49.0% ] . an ensemble average is taken over @xmath86 random configurations of @xmath87 for @xmath88 and @xmath89 , and over @xmath90 configurations for @xmath91 , @xmath92 , @xmath93 and @xmath94 . the mv model parameter is chosen as @xmath76 . right : chromo - electric and chromo - magnetic fields   and for the same mv model parameters and @xmath88.,title="fig:",scaledwidth=49.0% ]     defined in eq .   and transverse pressure @xmath69 defined in eq .  , multiplied by @xmath95.,scaledwidth=49.0% ]    figure  [ fig : energy ] shows the time evolution of the energy density   and . we set @xmath76 which means that @xmath96 with @xmath97 according to the convention for the rhic physics in the literature  @xcite . this choice corresponds to @xmath98 [ @xmath99 ( i.e. the radius of the au atom ) ] and @xmath100 ( where @xmath101 as usually chosen , meaning that @xmath102 ) . in most of the numerical results presented in this work , we keep using these parameters . we note that the strength of the glasma instability as seen in the next section is set by the instability seed @xmath20 relative to @xmath103 which characterizes the cgc background . in this sense , for our purpose in later discussions , the precise value of @xmath104 is of no qualitative importance .    in fig . [ fig : energy ] , we investigate the dependence on the transverse grid size @xmath13 : the left plot shows how the energy density increases with increasing @xmath13 , varying from @xmath105 to @xmath88 from the bottom curve to the topmost one . it is already known in the literature  @xcite that the @xmath13-dependence ( or the uv - cutoff dependence ) is significant at @xmath106 , but becomes almost irrelevant for @xmath107 ( except for small lattice sizes @xmath108 ) .    as is well - known from earlier works  @xcite , the energy density decreases as @xmath109 , which is the scaling behavior of a free - streaming expansion . in fact the energy density multiplied by the time , @xmath110 , saturates at a time @xmath111 ( which is sometimes called the formation time , i.e. the time at which partons are freed from the wave - function of the nuclei ) where the free - streaming expansion starts . this is natural on physical grounds , because the classical fields keep interacting strongly until @xmath111 , and become quasi - free at @xmath112 . we see that the longitudinal chromo - electric and chromo - magnetic fields behave in the same way , which is the characteristic feature of the glasma initial state  @xcite . in our convention @xmath65 , @xmath113 , @xmath66 , and @xmath114 are all approaching a common value for @xmath115 , and the system still remains anisotropic ( isotropy requires that @xmath116 and @xmath117 ) . this becomes clearer if we consider the transverse pressure @xmath69 and the longitudinal one @xmath118 defined in eqs .   and , respectively , as shown in fig . [ fig : p ] . then , we can find that @xmath119 and @xmath120 both approach an asymptotic value and @xmath120 is much smaller than @xmath119 . if @xmath118 is exactly zero , this means that the system reaches the limit of the free - streaming expansion leading to @xmath121 . if the system is completely thermalized , on the other hand , the pressure must be isotropic and thus @xmath71 . because an expansion with positive @xmath118 produces work against the expansion of the matter , the energy density decreases faster than in the free - streaming case , e.g. @xmath122 if @xmath123 and if the expansion is purely longitudinal ( i.e. bjorken expansion ) . therefore , the ( leading order ) glasma state alone can not reach isotropization ; a candidate for causing the longitudinal pressure to grow is the instability due to the @xmath12-dependent fluctuations that occur in higher order corrections . let us now consider the time evolution of the spectral composition of the energy content . that is , the energy is written in terms of fourier components as @xmath124 in terms of the continuum variables . if written in terms of the lattice variables , the transverse momenta take discrete values @xmath125 and the integration is replaced by a summation over @xmath126 . the energy density in the mode @xmath127 is given by @xmath128\bigr\rangle , \label{eq : ee_k}\\   \varepsilon_{_b}(k_\perp ) & \equiv \bigl\langle { \mathrm{tr}}\bigl [    b^{\eta a}(-{\boldsymbol{k}_\perp})b^{\eta a}({\boldsymbol{k}_\perp } ) + \tau^{-2}\bigl (    b^{ia}(-{\boldsymbol{k}_\perp})b^{ia}({\boldsymbol{k}_\perp } ) \bigr ) \bigr]\bigr\rangle\ ; , \label{eq : eb_k}\end{aligned}\ ] ] where @xmath129 , @xmath130 , @xmath131 , and @xmath132 are the fourier transforms of the chromo - electric and chromo - magnetic fields . the above expressions are given in terms of the continuum variables , but it is non - trivial how to re - express them in terms of the link variables : equation   for the chromo - electric component is not changed , but we can not use eq .   as it is and @xmath114 in eq . are perfectly fine to calculate the chromo - magnetic energy  the problem is that they are not simply a sum of squares . ] . it is of course possible to introduce the chromo - magnetic field using the relation   in such a way that @xmath133    \notag\\   b_{x}^a & = f_{\eta y}^a \approx    \frac{2}{{\mathrm{i}}g a_\eta}{\mathrm{tr}}\bigl [ t^a \bigl ( 1-u_{\eta y } \bigr )    \bigr]\ ; , \quad   b_{y}^a = f_{\eta x}^a \approx    \frac{2}{{\mathrm{i}}g a_\eta}{\mathrm{tr}}\bigl [ t^a \bigl ( 1-u_{\eta x } \bigr )    \bigr]\ ; . \label{eq : approxb}\end{aligned}\ ] ] this approximation works for sufficiently large number of the transverse sites , i.e.  @xmath134 . then , we in principle have two options here : we abandon the information on @xmath135 and just evaluate @xmath136 to see the spectral pattern , or , we adopt a sufficiently large value of @xmath13 so that we can use the approximation  . in this work we take the former option , because simulations with @xmath134 would be too time - consuming in the forthcoming studies that include @xmath12-dependent fluctuations ( and where one thus needs to solve ( @xmath137 ) dimensional equations ) .    at @xmath138 , the initial spectrum is calculable analytically in the mv model , leading to the expression  @xcite , @xmath139 , \label{eq : ee_pert}\ ] ] with an ir cutoff @xmath140 . this formula means that @xmath136 behaves as @xmath141 at large @xmath127 , which is the expected perturbative tail .     at @xmath142 , @xmath143 , @xmath144 , and @xmath145 . the transverse size is chosen as @xmath146 and the mv model parameter as @xmath147 . the spectrum at @xmath142 shows the perturbative scaling @xmath148 . an ensemble average is taken over @xmath86 configurations.,scaledwidth=49.0% ]    at non - zero @xmath14 , we have evaluated @xmath136 numerically , with results shown in fig . [ fig : time - trans ] . to draw the figure , we have calculated @xmath149 first as a function of @xmath150 and @xmath151 and then have taken an average over the results that lie in a bin , @xmath152 , to obtain @xmath136 with @xmath153 . we use logarithmic horizontal and vertical axes in the plot of fig . [ fig : time - trans ] , to facilitate the identification of the initial - time power - law scaling in the spectrum . we can notice that the energy spectrum at the initial time @xmath142 ( where we started the simulation ) obeys the perturbative scaling @xmath154 up to the uv - cutoff , as indicated by a line for eye guide in fig . [ fig : time - trans ]  in agreement with the perturbative expression in eq .  . the spectral shape quickly approaches the asymptotic one within @xmath155 and it hardly changes after that time . in fact the shape at @xmath156 is already close to that at @xmath157 ( where we stopped the simulation ) . thanks to this , we do not show the results at any intermediate time between @xmath156 and @xmath157 in order to keep the plot legible .    in later discussions in the section  [ sec : lspectrum ] we will perform a similar analysis in the longitudinal direction to find non - trivial asymptotic behavior of the spectral shapes . after this summary of the boost - invariant results , including the time evolution of the energy spectrum , let us proceed with the analysis of the instability with respect to @xmath12-dependent fluctuations . we will first confirm the existence of the instability in a similar setup as the one used in previous studies . for this purpose , we choose the initial fluctuations according to the prescription used in refs . @xcite and we fix the mv model parameter @xmath103 in the same way as in the boost - invariant situation . we take the extent in the @xmath12 direction to be @xmath158 units of rapidity and we use periodic boundary conditions for @xmath12 as well as for transverse coordinates . we note that this extent covers roughly the mid - rapidity region where the boost - invariant plateau is expected , and therefore it is legitimate to have a boost - invariant background . that is , our parameter choice reads ; @xmath159 and we perform calculations at various values of the transverse grid size @xmath13 and the longitudinal grid size @xmath17 later .    introducing an arbitrary function @xmath160 , we can write the @xmath12-dependent chromo - electric fields in such a way that gauss s law is satisfied by construction , @xmath161\ ,    \xi^i({\boldsymbol{x}_\perp } ) \ ; , \label{eq : deltaei}\\   \delta e^\eta(x ) & = -f(\eta)\sum_{i = x , y }    \bigl [ u_i^\dagger(x-\hat{i})\xi^i({\boldsymbol{x}_\perp}-\hat{i})u_i(x-\hat{i } )    -\xi^i({\boldsymbol{x}_\perp } ) \bigr]\;. \label{eq : deltaeeta}\end{aligned}\ ] ] before applying eq . we can also add @xmath12-dependent fluctuations in @xmath162 . for the moment we will discuss the case with @xmath163 and @xmath164 only . here the choice of @xmath165 is arbitrary and we choose it as a random variable in the transverse plane , i.e.@xmath166 it is important to note that @xmath167 has the dimension of a momentum , which is obvious from the above equation . if all the dimensional quantities are scaled with the transverse spacing @xmath36 as is the case in the numerical implementation , @xmath167 scales as @xmath168 , which is artificial . the strength of the longitudinal disturbance should be given independently . to cancel this artificial @xmath36 dependence , we use a dimensionless @xmath169 satisfying @xmath170 ( on the lattice ) and make @xmath160 scale as @xmath171 , so that the fluctuations are proportional overall to @xmath172 which is a constant .    in the literature @xcite @xmath160 is chosen as a random variable too containing all infrared and ultraviolet modes in the spectrum . here , in order to study the instability behavior in a well - controlled situation , let us perturb the system with a single @xmath12-mode , and consider superposition of multiple modes only later . hence , we will adopt the following simple form for the moment , the frequency @xmath21 is an integer between @xmath22 and @xmath23 . ] , @xmath173 where @xmath174 is chosen to be @xmath70 as we explained before and @xmath20 should scale as @xmath175 to make the fluctuations insensitive to the way we discretize the transverse plane . we will later fix @xmath176 and then vary @xmath13 to study the ( unphysical ) dependence on the transverse grid size .      here we fix @xmath177 in eq .   because , as we will confirm later , this lowest non - zero mode leads to the fastest and strongest growth in the longitudinal pressure . figure  [ fig : size_t ] shows how the instability occurs in the longitudinal pressure @xmath118 in the long - time run . the left - panel is the pressure multiplied by time ( @xmath178 ) in the unit of @xmath179 . we note that the plotted quantity is for @xmath180 which maximizes the fourier transform of @xmath120 ( except for the zero mode @xmath181 ) . that is , we define , @xmath182 which is complex valued in general in our definition ranges from @xmath22 to @xmath183 and the corresponding wave - number is @xmath184 . ] . we thus plot its modulus , @xmath185 , in the left panel of fig . [ fig : size_t ] corresponding to an initial disturbance at @xmath177 and @xmath186 in eqs .   and ( i.e.  we choose @xmath187 as a reference point and then @xmath188 for @xmath187 ) . the horizontal axis represents @xmath189 , which is a natural variable in a longitudinally expanding geometry . note that the fourier transform @xmath190 at finite @xmath18 is not a measurable quantity in experiments . its real and imaginary parts both fluctuate from negative to positive values depending on the initial conditions and they would be vanishing if we compute an average without taking the modulus . only the @xmath181 mode has a physical meaning as the longitudinal pressure . for our purpose of studying the instability , we can nevertheless consider @xmath191 ( i.e. we compute the modulus first , and we perform the average over initial conditions next ) in order to confirm some known results reported previously  @xcite . in later discussions , we will also decompose the energy density into @xmath18 components , which is a more well - defined quantity to see the instability .    , where @xmath177 is chosen for the initial condition  . the results are for @xmath105 , @xmath93 , @xmath92 , @xmath192 , @xmath193 in the transverse direction ( from the top to the bottom ) and @xmath194 fixed in the longitudinal direction ( except for @xmath195 for @xmath75 and @xmath196 for @xmath197 ) . the seed magnitude is @xmath188 for @xmath187 , i.e.  @xmath186 in eq .  . an ensemble average is taken over @xmath198 configurations . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199 ( @xmath200).,title="fig:",scaledwidth=49.0% ] , where @xmath177 is chosen for the initial condition  . the results are for @xmath105 , @xmath93 , @xmath92 , @xmath192 , @xmath193 in the transverse direction ( from the top to the bottom ) and @xmath194 fixed in the longitudinal direction ( except for @xmath195 for @xmath75 and @xmath196 for @xmath197 ) . the seed magnitude is @xmath188 for @xmath187 , i.e.  @xmath186 in eq .  . an ensemble average is taken over @xmath198 configurations . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199 ( @xmath200).,title="fig:",scaledwidth=49.0% ]    from the left panel of fig . [ fig : size_t ] , the exponential growth with the square root of proper time is obvious . the onset of the instability depends on the transverse size @xmath13 but the instability slope is rather insensitive to @xmath13 . with increasing @xmath13 , the instability occurs earlier , and eventually the shift in this onset position saturates for @xmath201 . this means that the instability does not start earlier by increasing @xmath13 further . so , we stop increasing @xmath13 and conclude that simulations at @xmath187 are reasonably close to the continuum limit . in the right panel , we see how the fields , initially localized in the modes @xmath202 ( background field itself ) and @xmath203 ( perturbation ) , populate the higher @xmath18 modes due to the non - linearities in the yang - mills equations . for a small perturbation , the amplitude in the higher modes decreases very fast ( exponentially with @xmath18 ) , because in order to produce a field at a given @xmath18 , one needs to multiply @xmath18 seeds at @xmath177 .      . the results are for @xmath188 and @xmath187 fixed in the transverse direction and @xmath194 and @xmath204 in the longitudinal direction . an ensemble average is taken over @xmath198 configurations . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199.,title="fig:",scaledwidth=49.0% ] . the results are for @xmath188 and @xmath187 fixed in the transverse direction and @xmath194 and @xmath204 in the longitudinal direction . an ensemble average is taken over @xmath198 configurations . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199.,title="fig:",scaledwidth=49.0% ]    physical results must also not depend on how we discretize the longitudinal coordinate . in order to see the dependence on the longitudinal grid size , we make a plot in fig . [ fig : size_l ] in the same way as in fig . [ fig : size_t ] , where we compare two longitudinal grid sizes that differ by a factor @xmath70 . as is clear from fig . [ fig : size_l ] we can not find any sizable dependence on @xmath17 , which means that our numerical calculations are safely free from lattice artifacts in the longitudinal direction .    the influence of @xmath17 can be seen more clearly in the spectrum of @xmath205 shown in the right panel of fig . [ fig : size_l ] . a different choice of @xmath17 with fixed @xmath174 changes the uv cutoff ( the spectrum ranges between @xmath22 and @xmath23 ) : in the right panel of fig .  [ fig : size_l ] , @xmath18 varies from @xmath206 to @xmath207 for the @xmath194 calculation and from @xmath208 to @xmath209 for the @xmath204 calculation . since the seeds are at @xmath180 in our initial condition , we naturally anticipate that the uv cutoff at @xmath210 does not affect physical properties around @xmath211 . one may think that the results for a larger seed amplitude @xmath212 in fig . [ fig : size_l ] could be sensitive to the choice of @xmath17 because all the amplitudes in the spectrum are substantially larger and thus a change in @xmath17 leading to a change in the uv cutoff may have an influence on the spectrum . we have also checked this and found that it is not the case . as is clear from fig . [ fig : size_l ] , the difference between the results for @xmath213 with @xmath212 and that for @xmath214 with @xmath212 is almost invisible . from all these results we can conclude that a value of @xmath194 is sufficiently large as long as @xmath21 is set to be of order unity . if @xmath21 takes a larger value , as we will see soon below , the spectral shape has a wider distribution and the uv cutoff at @xmath210 could influence the results .      . the results are for @xmath215 , @xmath216 , @xmath217 and a fixed lattice size @xmath213 . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199 , which has a uv cutoff at @xmath218 for the @xmath194 calculation.,title="fig:",scaledwidth=49.0% ] . the results are for @xmath215 , @xmath216 , @xmath217 and a fixed lattice size @xmath213 . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199 , which has a uv cutoff at @xmath218 for the @xmath194 calculation.,title="fig:",scaledwidth=49.0% ]    in this subsection , we study the dependence on the seed magnitude @xmath20 of the longitudinal fluctuations in eq .  . in fig .  [ fig : magnitude ] , we show some of numerical results for different values of @xmath20 . it is interesting that the results obviously scale as @xmath219 up to the point where the instability growth eventually saturates . this @xmath20 dependence is almost trivial , which suggests that the non - linear effects in terms of finite-@xmath18 modes are small unless the unstable modes substantially develop at large @xmath95 . one might have thought that any physical observables should be expanded in even powers of @xmath20 and thus @xmath220 should be proportional @xmath221 rather than @xmath20 . this is not the case because we took the modulus of complex @xmath222 before performing the average over initial conditions , thus preventing the linear term in @xmath20 from vanishing . we can also notice from the right panel of fig . [ fig : magnitude ] that , once the instability at @xmath223 gets saturated , the mode at @xmath223 stops growing but the spectrum spreads quickly to higher @xmath18-modes . once this happens , the scaling property with @xmath20 is lost . actually the amplitude of the @xmath223 mode reaches one third of the zero - mode amplitude ( for @xmath212 results at @xmath199 ) before this saturation occurs . then the non - linearities or self - interaction effects can not be neglected any longer . now we investigate the @xmath21-dependence of the instability behavior . first , we have chosen three different @xmath21 s as @xmath177 , @xmath224 , and @xmath225 to uncover the general trend . then we study a more general case , that consists in a superposition of these initial wave - numbers . , @xmath224 , and @xmath225 with the initial condition given as a superposition  . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199 . the solid line represents the results from the non - linear evolution with the initial condition given as a superposition , while the dotted one is a simple linear superposition of the results shown in fig . [ fig : nu].,title="fig:",scaledwidth=49.0% ] , @xmath224 , and @xmath225 with the initial condition given as a superposition  . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199 . the solid line represents the results from the non - linear evolution with the initial condition given as a superposition , while the dotted one is a simple linear superposition of the results shown in fig . [ fig : nu].,title="fig:",scaledwidth=49.0% ]    , @xmath224 , and @xmath225 with the initial condition given as a superposition  . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199 . the solid line represents the results from the non - linear evolution with the initial condition given as a superposition , while the dotted one is a simple linear superposition of the results shown in fig . [ fig : nu].,title="fig:",scaledwidth=49.0% ] , @xmath224 , and @xmath225 with the initial condition given as a superposition  . right : spectrum of the longitudinal `` pressure '' as a function of @xmath18 at @xmath199 . the solid line represents the results from the non - linear evolution with the initial condition given as a superposition , while the dotted one is a simple linear superposition of the results shown in fig . [ fig : nu].,title="fig:",scaledwidth=49.0% ]    the figure  [ fig : nu ] shows the numerical results . it is clear from the figure that the onset of the instability is delayed for larger @xmath21 . this is naturally understandable from the kinetic term for finite-@xmath18 modes . the ( positive ) kinetic term generally tends to stabilize the system , and so the ( negative ) term responsible for the instability is overwhelmed if @xmath18 is large enough . for dimensional reasons , however , this kinetic term for the finite-@xmath18 modes is of the form @xmath226 , and therefore it decreases with time . in contrast , the negative term that generates the instability has a constant coefficient of order unity  therefore , it dominates over the kinetic term at sufficiently large time . in other words , the kinetic term in @xmath226 only delays the start of the instability by a time proportional to @xmath18 . we can confirm qualitatively this interpretation from fig . [ fig : nu ] . the growth appears to start around @xmath227 for @xmath228 , respectively , and the starting times are indeed ordered proportionally to @xmath21 . the spectral shape shown in the right panel of fig . [ fig : nu ] looks discontinuous for @xmath229 and @xmath225 , but this is a trivial artifact of the fact that our initial seed contains a single @xmath21 mode . indeed , the non - linear couplings in the yang - mills equations can only produce linear combinations with integer coefficients of the frequencies present in the initial condition . for instance , with @xmath229 , only modes with @xmath18 a multiple of @xmath224 can be produced by the non - linear evolution . the fact that the amplitude if these higher harmonics is suppressed compared to the amplitude of the base frequency @xmath21 is a good indication that we are still in a regime where the non - linearities are rather small . we can also see that the non - linear effects are certainly present but not very large by repeating the above calculation with an initial condition given as a superposition of three modes , namely , @xmath230 with @xmath231 , @xmath232 , @xmath233 , instead of the single - mode fluctuation of eq .   with @xmath177 , @xmath229 , @xmath234 individually . if the non - linearity is significant , the results should not be a simple superposition of the individual calculation . in view of fig . [ fig : nu_mix ] , however , we can see only small deviations from the individual results shown in fig .  [ fig : nu ] at @xmath199 . this is an important observation because the initial fluctuations in reality should have some continuous spectral distribution ( see ref . @xcite for instance ) and we simplified our analysis by considering the special case of a single - mode fluctuation . as long as the non - linearity is a minor effect , we can just add up the single - mode outputs in order to get approximate results for general initial spectra . although they are several orders of magnitude smaller at @xmath199 , we see that higher harmonics mixing the three @xmath21 s are enhanced by non - linear effects in the right panel of fig . [ fig : nu_mix ] . for instance , the amplitude of the mode @xmath235 is much larger than the sum of the amplitudes at @xmath235 for the three seeds considered in isolation . indeed , when the seed @xmath177 is alone for example , the mode @xmath235 is produced only after 20 interactions , which explains why it is totally negligible . furthermore , the seeds @xmath229 or @xmath225 can never lead to @xmath235 by themselves . when the three seeds are superposed in the initial condition , 21 can be reached more efficiently as @xmath236 , i.e. with 2 interactions only . apart from small non - linear effects , it is certainly true that the smallness of the field fluctuations we have considered gives a justification for our simplified , single mode , analysis . however , at the same time , thermalization is hardly expected to occur in these circumstances , precisely because in this regime the interaction effects are small . in realistic conditions , the parameter @xmath20 that controls the fluctuations is related to the strong coupling @xmath9 , and therefore it could be larger than @xmath237 that we have considered here . as a consequence , we expect that the non - linearities would start much sooner and be much more important in applications of this framework to actual heavy - ion collisions . the study with such large @xmath20 is beyond our current scope because a large @xmath20 requires full information on quantum fluctuations and needs appropriate renormalization procedures , which is technically involved  @xcite . in spite of the instability the zero - mode is still dominant among all the modes in @xmath190 . with an initial condition that contains a single mode as in eq .  , the @xmath180 mode is the second largest , which is obvious from the spectrum shown in the right panel of fig . [ fig : size_t ] . moreover , the spectrum is spreading into higher wave - numbers as the time elapses . in this section , we shall discuss the scaling property of the energy spectrum with respect to the frequency @xmath18 , and show that it possibly displays also a kolmogorov s turbulent spectrum . a turbulent spectrum has already been reported in the context of the qcd plasma instability  @xcite , but these results with the power @xmath238 instead of @xmath239 are in conflict with kolmogorov s wave turbulence . on the other hand , the study performed in ref . @xcite has shown that a kolmogorov spectrum appears due to instabilities in yang - mills theory , for a fixed - volume system . therefore , it is still worth addressing the turbulent spectrum in the present glasma simulation , that has longitudinal expansion . this is because the glasma system is expanding along the beam axis , which tends to tame the turbulences and to delay the isotropization or thermalization in general . it is thus non - trivial whether the universal scaling still holds or not in the expanding system . besides , our description of the glasma evolution is purely classical and the precise relationship between the glasma instability and the qcd plasma instability is not fully understood yet .    . the lattice size is @xmath213 and the seed is @xmath188 . left : results at early times up to @xmath240 . note that the horizontal axis is not logarithmic and also that the vertical axis is for not @xmath241 but @xmath242 ( see the text for details ) . right : results at later times up to @xmath243.,title="fig:",scaledwidth=49.0% ] . the lattice size is @xmath213 and the seed is @xmath188 . left : results at early times up to @xmath240 . note that the horizontal axis is not logarithmic and also that the vertical axis is for not @xmath241 but @xmath242 ( see the text for details ) . right : results at later times up to @xmath243.,title="fig:",scaledwidth=49.0% ]    two plots in fig . [ fig : time - scale - long ] show the longitudinal energy spectrum ( for the chromo - electric part only ) , that is , the energy density decomposed in terms of the longitudinal wave - number @xmath18 as @xmath244 with @xmath245\bigr\rangle .\ ] ] ( the transverse momentum @xmath246 of the fields is integrated over . ) the question is then what kind of scaling spectrum we can expect for @xmath247 at later times . for this purpose , we consider three characteristic quantities ; the wave - number , the fourier decomposed energy , and the rate of energy flow in @xmath18-space ( denoted by @xmath248 below ) . the wave - number @xmath18 itself is , however , dimensionless and it seems difficult to cope with the dimensional matching argument by kolmogorov . we propose that one should use a scaled variable @xmath249 instead of @xmath18 , in the case of an expanding geometry . this is because @xmath250 is approximately equal to the longitudinal coordinate @xmath251 if @xmath12 is small enough , and @xmath249 is approximately equal to the momentum @xmath252 . then , the dimension of the scaled wave - number is @xmath253=l^{-1}$ ] . we note that the integration variable in the decomposition   should also be @xmath254 for this scaling analysis , which means that the relevant energy density per mode should be @xmath255 . this in addition should be multiplied by the expanding length @xmath256 . therefore , the dimensions of the three relevant quantities are @xmath257 = l^{-1 } , \qquad   [ v_\perp(c\tau)^2\varepsilon_{_e}(\nu ) ] = l^3 t^{-2 } , \qquad   [ \psi ] = l^2 t^{-3 } , \ ] ] which are identical to the standard analysis of kolmogorov s spectrum . provided that @xmath258 is expressed in a form of @xmath259 , we can fix @xmath260 and @xmath261 uniquely from the dimensions and thus expect the kolmogorov scaling , that is , @xmath262 . in fact we can confirm this scaling property asymptotically from the right panel of fig . [ fig : time - scale - long ] . one might be wondering what the fate of the power - law spectrum should be at even later times . we understand from fig . [ fig : stage ] that the glasma time evolution in the present mv - model simulation has three stages ; the linear regime , the non - linear regime , and the uv - cutoff regime . in the linear regime the unstable modes grow up and the non - linear effects are still minor . when the amplitudes of @xmath263 modes become so large that the non - linearities can affect the time evolution , the instability stops , as shown in the right panel of fig . [ fig : stage ] , and the non - abelianization occurs @xcite . because non - linearities are no longer negligible in this regime , the energy cascade due to non - linearities is naturally expected , which would lead to the kolmogorov - type spectrum . eventually , as the turbulence decays into higher uv modes , the uv - cutoff effect flattens the whole spectrum as plotted in fig . [ fig : stage ] up to @xmath264 .     up to @xmath264 .  left : results including later time spectra . the horizontal axis is logarithmic to make it evident how the power - law spectrum spreads . right : longitudinal `` pressure '' with three regimes identified schematically.,scaledwidth=49.0% ]    to the best of our knowledge the present work is the first to exhibit a kolmogorov spectrum from the glasma time evolution . our results imply that the framework of the glasma , in terms of classical variables together with @xmath12-dependent fluctuations , may encompass a rich physics contents including the qcd plasma instability .    here , as a final remark in this section , we draw attention to related works in analogous systems . in the context of an ultra - cold gas , the kolmogorov spectrum has been obtained as a solution of the gross - pitaevskii equation ( gpe ) , which is a classical equation of motion like the glasma equations  @xcite . the calculational procedure by means of the gpe is similar to the glasma simulation , but in that case the microscopic mechanism of the energy flow is understood in terms of the vortex dynamics . in future studies on the glasma dynamics , it will be important to clarify the structure of the energy flow , and it would be very interesting to see whether the so - called non - thermal fixed point  @xcite appears in the strongly correlated yang - mills theory and glasma . indeed , a possible deviation in the power could be expected in yang - mills theory by non - perturbative resummations  @xcite . in this paper we investigated the time evolution of the classical equations of motion in the su(2 ) pure yang - mills theory , especially in the framework of the glasma or the mclerran - venugopalan model in ( 3 + 1 ) dimensions . our emphasis has been put on the systematic study of the instability behavior that occurs in fluctuations with respect to the space - time rapidity @xmath12 when the initial conditions break boost invariance . computing the longitudinal pressure component at non - zero wave - number @xmath18 , that is a conjugate of @xmath12 , we have numerically confirmed the glasma instability that had been found in preceding works  @xcite . we carefully checked the system size dependence of the instability to conclude that the glasma instability is a robust consequence from @xmath12-dependent disturbances and is not sensitive to how the transverse and longitudinal coordinates are discretized . hence , the glasma instability is certainly a physically well - defined phenomenon . at the same time , however , this means that taking the continuum limit does not help with an earlier onset of the instability . the typical time scales we found for the glasma instability is of order of hundreds in terms of @xmath95 , which is arguably too slow for being useful in heavy ion collisions . one should however keep in mind that the seeds we have used in this numerical study are very small compared to the glasma fields , while in a realistic situation they are suppressed at most by a power of the strong coupling constant . there may also be a chance that the non - linear effects could more or less accelerate the instability . we compared two results for the spectral distribution , one from the superposition of individual calculations with the initial condition that is given by a different wave - number @xmath21 , and the other from a single calculation with the initial condition that contains all @xmath21 s . we could observe the non - linear effects in our numerical results , which would not affect the instability behavior substantially , however , as long as the instability seed @xmath20 is reasonably small . we have decomposed the energy density into fourier modes to see how the low wave - number modes cascade toward higher wave - number modes due to non - linearity in the equations of motion . our motivation to do so is to examine whether the so - called kolmogorov s scaling can appear in the turbulent system described by the yang - mills theory put in an expanding geometry . in contrast to the boost - invariant case , a kolmogorov - type scaling behavior in the energy spectrum a function of the longitudinal wave - number @xmath18 may appear when the boost invariance is broken by the initial fluctuations . in this case , the violation of boost invariance increases with time due to instability , but the @xmath181 mode remains very large compared to the other @xmath263 modes , which enables us to regard the @xmath181 mode as the source for the energy injection . then , we could clearly observe that the energy spectra approach an asymptotic form that shows the kolmogorov s scaling with the power @xmath265 . as far as we know , our work is the first example that exhibits the kolmogorov - type scaling law in the expanding systems . although it is difficult to find any indication that the system indeed becomes isotropic and thermalized in our simulation , the confirmation of the kolmogorov - type spectra is a promising signal for the general tendency heading for thermalization .    there remain very interesting questions to be investigated in the future . in order to unveil the microscopic mechanism for the instability , it would be very useful to make approximations on the classical equations of motion in such a way that the analytical treatment could be feasible . for instance , as we have confirmed , the instability mostly lies in the linear regime as long as its magnitude specified by @xmath20 is small enough . therefore the linearization of the equations of motion in terms of fluctuating fields on top of the boost - invariant cgc background should correctly describe the instability in early times . in this work we did not exploit such a comparison between the numerical output and the analytical result from the linearized equations . this is because it is not quite efficient to linearize the equations of motion given in terms of the lattice variables . the linearization is of course possible , but it would eventually be reduced to similar calculations in the continuum formulation . then , it would make much more sense to make a comparison between the numerical and analytical results that are both formulated in the continuum variables . we have already executed test programs to make sure if the numerical results from the glasma simulations are consistent with each other , given they are written in terms of lattice and continuum variables . we will report on this analysis elsewhere . it is also an intriguing question to look for a link between kolmogorov s turbulent spectrum and the chaotic behavior generally seen in the solution of the classical yang - mills equations . this problem may have some relevance for a deeper understanding of the microscopic description of the entropy generating process and how the thermalization is achieved at all . the time evolution of the glasma and its physics contents deserve much more investigations in the future . we would like to thank kevin dusling , tuomas lappi , and raju venugopalan for discussions . this work has been supported in part by jsps japan - france integrated action program ( sakura ) . thanks the yukawa institute for theoretical physics , where part of this work has been done , for its hospitality and support . p.  f. kolb and r.  rapp , `` transverse flow and hadrochemistry in au+au collisions at @xmath266 , '' http://dx.doi.org/10.1103/physrevc.67.044903[_phys . rev . _ * c67 * ( 2003 ) 044903 ] , http://arxiv.org/abs/hep-ph/0210222 [ arxiv : hep - ph/0210222 [ hep - ph ] ] .      r.  baier , a.  h. mueller , d.  schiff , and d.  t. son , `` bottom - up thermalization in heavy ion collisions , '' http://dx.doi.org/10.1016/s0370-2693(01)00191-5[_phys . lett . _ * b502 * ( 2001 ) 5158 ] , http://arxiv.org/abs/hep-ph/0009237[arxiv:hep-ph/0009237 ] . y.  v. kovchegov , `` can thermalization in heavy ion collisions be described by qcd diagrams ? , '' http://dx.doi.org/10.1016/j.nuclphysa.2005.08.009[_nucl . _ * a762 * ( 2005 ) 298325 ] , http://arxiv.org/abs/hep-ph/0503038[arxiv:hep-ph/0503038 ] . b. arnold , j.  lenaghan , and g.  d. moore , `` qcd plasma instabilities and bottom up thermalization , '' _ jhep _ * 0308 * ( 2003 ) 002 , http://arxiv.org/abs/hep-ph/0307325[arxiv:hep-ph/0307325 [ hep - ph ] ] . erratum added online , sep/29/2004 . b. arnold , j.  lenaghan , g.  d. moore , and l.  g. yaffe , `` apparent thermalization due to plasma instabilities in quark - gluon plasma , '' http://dx.doi.org/10.1103/physrevlett.94.072302[_phys . * 94 * ( 2005 ) 072302 ] , http://arxiv.org/abs/nucl-th/0409068 [ arxiv : nucl - th/0409068 [ nucl - th ] ] . p.  b. arnold , g.  d. moore , and l.  g. yaffe , `` the fate of non - abelian plasma instabilities in 3 + 1 dimensions , '' http://dx.doi.org/10.1103/physrevd.72.054003[_phys . * d72 * ( 2005 ) 054003 ] , http://arxiv.org/abs/hep-ph/0505212 [ arxiv : hep - ph/0505212 [ hep - ph ] ] .    p.  b. arnold and g.  d. moore , `` qcd plasma instabilities : the nonabelian cascade , '' http://dx.doi.org/10.1103/physrevd.73.025006[_phys . rev . _ * d73 * ( 2006 ) 025006 ] , http://arxiv.org/abs/hep-ph/0509206[arxiv:hep-ph/0509206 [ hep - ph ] ] . j.  randrup and s.  mrowczynski , `` chromodynamic weibel instabilities in relativistic nuclear collisions , '' http://dx.doi.org/10.1103/physrevc.68.034909[_phys . rev . _ * c68 * ( 2003 ) 034909 ] , http://arxiv.org/abs/nucl-th/0303021[arxiv:nucl-th/0303021 ] . s.  mrowczynski , `` instabilities driven equilibration of the quark - gluon plasma , '' _ acta phys . * b37 * ( 2006 ) 427454 , http://arxiv.org/abs/hep-ph/0511052[arxiv:hep-ph/0511052 ] . a.  dumitru , y.  nara , and m.  strickland , `` ultraviolet avalanche in anisotropic non - abelian plasmas , '' http://dx.doi.org/10.1103/physrevd.75.025016[_phys . rev . _ * d75 * ( 2007 ) 025016 ] , http://arxiv.org/abs/hep-ph/0604149 [ arxiv : hep - ph/0604149 [ hep - ph ] ] . j.  berges , d.  gelfand , s.  scheffler , and d.  sexty , `` simulating plasma instabilities in su(3 ) gauge theory , '' http://dx.doi.org/10.1016/j.physletb.2009.05.008[_phys . * b677 * ( 2009 ) 210213 ] , http://arxiv.org/abs/0812.3859[arxiv:0812.3859 [ hep - ph ] ] . d.  bodeker , `` the impact of qcd plasma instabilities on bottom - up thermalization , '' http://dx.doi.org/10.1088/1126-6708/2005/10/092[_jhep_ * 0510 * ( 2005 ) 092 ] , http://arxiv.org/abs/hep-ph/0508223 [ arxiv : hep - ph/0508223 [ hep - ph ] ] . a.  rebhan , p.  romatschke , and m.  strickland , `` hard - loop dynamics of non - abelian plasma instabilities , '' http://dx.doi.org/10.1103/physrevlett.94.102303[_phys . lett . _ * 94 * ( 2005 ) 102303 ] , http://arxiv.org/abs/hep-ph/0412016[arxiv:hep-ph/0412016 ] . a.  rebhan , p.  romatschke , and m.  strickland , `` dynamics of quark - gluon plasma instabilities in discretized hard - loop approximation , '' http://dx.doi.org/10.1088/1126-6708/2005/09/041[_jhep_ * 09 * ( 2005 ) 041 ] , http://arxiv.org/abs/hep-ph/0505261[arxiv:hep-ph/0505261 ] . p.  romatschke and a.  rebhan , `` plasma instabilities in an anisotropically expanding geometry , '' http://dx.doi.org/10.1103/physrevlett.97.252301[_phys . _ * 97 * ( 2006 ) 252301 ] , http://arxiv.org/abs/hep-ph/0605064[arxiv:hep-ph/0605064 ] . a.  rebhan , m.  strickland , and m.  attems , `` instabilities of an anisotropically expanding non - abelian plasma : 1d+3v discretized hard - loop simulations , '' http://dx.doi.org/10.1103/physrevd.78.045023[_phys . rev . _ * d78 * ( 2008 ) 045023 ] , http://arxiv.org/abs/0802.1714 [ arxiv:0802.1714 [ hep - ph ] ] .    a.  ipp , a.  rebhan , and m.  strickland , `` non - abelian plasma instabilities : su(3 ) vs. su(2 ) , '' http://arxiv.org/abs/1012.0298[arxiv:1012.0298 [ hep - ph ] ] . b.  muller and a.  trayanov , `` deterministic chaos in nonabelian lattice gauge theory , '' http://dx.doi.org/10.1103/physrevlett.68.3387[_phys * 68 * ( 1992 ) 33873390 ] . t.  s. biro , c.  gong , and b.  muller , `` lyapunov exponent and plasmon damping rate in nonabelian gauge theories , '' http://dx.doi.org/10.1103/physrevd.52.1260[_phys . * d52 * ( 1995 ) 12601266 ] , http://arxiv.org/abs/hep-ph/9409392[arxiv:hep-ph/9409392 ] . t.  kunihiro _ et al . _ , `` chaotic behavior in classical yang - mills dynamics , '' http://dx.doi.org/10.1103/physrevd.82.114015[_phys . rev . _ * d82 * ( 2010 ) 114015 ] , http://arxiv.org/abs/1008.1156[arxiv:1008.1156 [ hep - ph ] ] . a.  nishiyama and a.  ohnishi , `` entropy production in gluodynamics in temporal axial gauge in 2 + 1 dimensions , '' _ prog . * 125 * ( 2011 ) 775795 , http://arxiv.org/abs/1011.4750[arxiv:1011.4750 [ nucl - th ] ] .    p.  m. chesler and l.  g. yaffe , `` horizon formation and far - from - equilibrium isotropization in supersymmetric yang - mills plasma , '' http://dx.doi.org/10.1103/physrevlett.102.211601[_phys . * 102 * ( 2009 ) 211601 ] , http://arxiv.org/abs/0812.2053 [ arxiv:0812.2053 [ hep - th ] ] .    p.  m. chesler and l.  g. yaffe , `` boost invariant flow , black hole formation , and far - from - equilibrium dynamics in n = 4 supersymmetric yang - mills theory , '' http://dx.doi.org/10.1103/physrevd.82.026006[_phys . rev . _ * d82 * ( 2010 ) 026006 ] , http://arxiv.org/abs/0906.4426 [ arxiv:0906.4426 [ hep - th ] ] .    y.  v. kovchegov and s.  lin , `` toward thermalization in heavy ion collisions at strong coupling , '' http://dx.doi.org/10.1007/jhep03(2010)057 [ _ jhep _ * 03 * ( 2010 ) 057 ] , http://arxiv.org/abs/0911.4707[arxiv:0911.4707 [ hep - th ] ] . v.  balasubramanian _ et al . _ , `` thermalization of strongly coupled field theories , '' http://arxiv.org/abs/1012.4753[arxiv:1012.4753 [ hep - th ] ] . f.  gelis , t.  lappi , and r.  venugopalan , `` high energy scattering in quantum chromodynamics , '' http://dx.doi.org/10.1142/s0218301307008331 [ _ int.j.mod.phys . _ * e16 * ( 2007 ) 25952637 ] , http://arxiv.org/abs/0708.0047[arxiv:0708.0047 [ hep - ph ] ] . l.  d. mclerran and r.  venugopalan , `` computing quark and gluon distribution functions for very large nuclei , '' http://dx.doi.org/10.1103/physrevd.49.2233[_phys . rev . _ * d49 * ( 1994 ) 22332241 ] , http://arxiv.org/abs/hep-ph/9309289[arxiv:hep-ph/9309289 ] . l.  d. mclerran and r.  venugopalan , `` gluon distribution functions for very large nuclei at small transverse momentum , '' http://dx.doi.org/10.1103/physrevd.49.3352[_phys * d49 * ( 1994 ) 33523355 ] , http://arxiv.org/abs/hep-ph/9311205[arxiv:hep-ph/9311205 ] . l.  d. mclerran and r.  venugopalan , `` fock space distributions , structure functions , higher twists and small x , '' http://dx.doi.org/10.1103/physrevd.59.094002[_phys . rev . _ * d59 * ( 1999 ) 094002 ] , http://arxiv.org/abs/hep-ph/9809427[arxiv:hep-ph/9809427 ] . a.  kovner , l.  d. mclerran , and h.  weigert , `` gluon production from nonabelian weizsacker - williams fields in nucleus - nucleus collisions , '' http://dx.doi.org/10.1103/physrevd.52.6231[_phys . rev . _ * d52 * ( 1995 ) 62316237 ] , http://arxiv.org/abs/hep-ph/9502289[arxiv:hep-ph/9502289 ] . a.  kovner , l.  d. mclerran , and h.  weigert , `` gluon production at high transverse momentum in the mclerran - venugopalan model of nuclear structure functions , '' http://dx.doi.org/10.1103/physrevd.52.3809[_phys . _ * d52 * ( 1995 ) 38093814 ] , http://arxiv.org/abs/hep-ph/9505320[arxiv:hep-ph/9505320 ] . d.  son , `` semiclassical approach for multiparticle production in scalar theories , '' http://dx.doi.org/10.1016/0550-3213(96)00386-0[_nucl . phys . _ * b477 * ( 1996 ) 378406 ] , http://arxiv.org/abs/hep-ph/9505338[arxiv:hep-ph/9505338 [ hep - ph ] ] . a.  mueller and d.  son , `` on the equivalence between the boltzmann equation and classical field theory at large occupation numbers , '' http://dx.doi.org/10.1016/j.physletb.2003.12.047[_phys . _ * b582 * ( 2004 ) 279287 ] , http://arxiv.org/abs/hep-ph/0212198 [ arxiv : hep - ph/0212198 [ hep - ph ] ] . a.  krasnitz and r.  venugopalan , `` non - perturbative computation of gluon mini - jet production in nuclear collisions at very high energies , '' http://dx.doi.org/10.1016/s0550-3213(99)00366-1[_nucl . _ * b557 * ( 1999 ) 237 ] , http://arxiv.org/abs/hep-ph/9809433[arxiv:hep-ph/9809433 ] . a.  krasnitz and r.  venugopalan , `` the initial energy density of gluons produced in very high energy nuclear collisions , '' http://dx.doi.org/10.1103/physrevlett.84.4309[_phys . * 84 * ( 2000 ) 43094312 ] , http://arxiv.org/abs/hep-ph/9909203[arxiv:hep-ph/9909203 ] . a.  krasnitz and r.  venugopalan , `` the initial gluon multiplicity in heavy ion collisions , '' http://dx.doi.org/10.1103/physrevlett.86.1717 [ _ phys . * 86 * ( 2001 ) 17171720 ] , http://arxiv.org/abs/hep-ph/0007108[arxiv:hep-ph/0007108 ] . t.  lappi , `` production of gluons in the classical field model for heavy ion collisions , '' http://dx.doi.org/10.1103/physrevc.67.054903[_phys . * c67 * ( 2003 ) 054903 ] , http://arxiv.org/abs/hep-ph/0303076[arxiv:hep-ph/0303076 ] . t.  lappi and l.  mclerran , `` some features of the glasma , '' http://dx.doi.org/10.1016/j.nuclphysa.2006.04.001[_nucl . phys . _ * a772 * ( 2006 ) 200212 ] , http://arxiv.org/abs/hep-ph/0602189[arxiv:hep-ph/0602189 ] . a.  dumitru , f.  gelis , l.  mclerran , and r.  venugopalan , `` glasma flux tubes and the near side ridge phenomenon at rhic , '' http://dx.doi.org/10.1016/j.nuclphysa.2008.06.012[_nucl . phys . _ * a810 * ( 2008 ) 91108 ] , http://arxiv.org/abs/0804.3858 [ arxiv:0804.3858 [ hep - ph ] ] . k.  dusling , f.  gelis , t.  lappi , and r.  venugopalan , `` long range two - particle rapidity correlations in a+a collisions from high energy qcd evolution , '' http://dx.doi.org/10.1016/j.nuclphysa.2009.12.044 ,    10.1016/j.nuclphysa.2009.12.044[_nucl . phys . _ * a836 * ( 2010 ) 159182 ] , http://arxiv.org/abs/0911.2720[arxiv:0911.2720 [ hep - ph ] ] . p.  romatschke and r.  venugopalan , `` collective non - abelian instabilities in a melting color glass condensate , '' http://dx.doi.org/10.1103/physrevlett.96.062302[_phys . _ * 96 * ( 2006 ) 062302 ] , http://arxiv.org/abs/hep-ph/0510121[arxiv:hep-ph/0510121 ] . p.  romatschke and r.  venugopalan , `` the unstable glasma , '' http://dx.doi.org/10.1103/physrevd.74.045011[_phys . rev . _ * d74 * ( 2006 ) 045011 ] , http://arxiv.org/abs/hep-ph/0605045[arxiv:hep-ph/0605045 ] . k.  fukushima , `` initial fields and instability in the classical model of the heavy - ion collision , '' http://dx.doi.org/10.1103/physrevc.76.021902[_phys . _ * c76 * ( 2007 ) 021902 ] , http://arxiv.org/abs/0704.3625[arxiv:0704.3625 [ hep - ph ] ] . k.  fukushima , `` erratum : initial fields and instability in the classical model of the heavy - ion collision , '' http://arxiv.org/abs/0711.2634[arxiv:0711.2634 [ hep - ph ] ] . h.  fujii and k.  itakura , `` expanding color flux tubes and instabilities , '' http://dx.doi.org/10.1016/j.nuclphysa.2008.05.016[_nucl . phys . _ * a809 * ( 2008 ) 88109 ] , http://arxiv.org/abs/0803.0410[arxiv:0803.0410 [ hep - ph ] ] . h.  fujii , k.  itakura , and a.  iwazaki , `` instabilities in non - expanding glasma , '' http://dx.doi.org/10.1016/j.nuclphysa.2009.06.020 [ _ nucl . phys . _ * a828 * ( 2009 ) 178190 ] , http://arxiv.org/abs/0903.2930[arxiv:0903.2930 [ hep - ph ] ] . k.  dusling , t.  epelbaum , f.  gelis , and r.  venugopalan , `` role of quantum fluctuations in a system with strong fields : onset of hydrodynamical flow , '' http://dx.doi.org/10.1016/j.nuclphysa.2010.11.009[_nucl . _ * a850 * ( 2011 ) 69109 ] , http://arxiv.org/abs/1009.4363[arxiv:1009.4363 [ hep - ph ] ] . k.  fukushima , f.  gelis , and l.  mclerran , `` initial singularity of the little bang , '' http://dx.doi.org/10.1016/j.nuclphysa.2007.01.086[_nucl . _ * a786 * ( 2007 ) 107130 ] , http://arxiv.org/abs/hep-ph/0610416[arxiv:hep-ph/0610416 [ hep - ph ] ] .    f.  gelis and r.  venugopalan , `` particle production in field theories coupled to strong external sources , '' http://dx.doi.org/10.1016/j.nuclphysa.2006.07.020[_nucl . * a776 * ( 2006 ) 135171 ] , http://arxiv.org/abs/hep-ph/0601209 [ arxiv : hep - ph/0601209 [ hep - ph ] ] . f.  gelis and r.  venugopalan , `` particle production in field theories coupled to strong external sources . ii . generating functions , '' http://dx.doi.org/10.1016/j.nuclphysa.2006.08.015[_nucl . phys . _ * a779 * ( 2006 ) 177196 ] , http://arxiv.org/abs/hep-ph/0605246 [ arxiv : hep - ph/0605246 [ hep - ph ] ] . f.  gelis , s.  jeon , and r.  venugopalan , `` how particles emerge from decaying classical fields in heavy ion collisions : towards a kinetic description of the glasma , '' http://dx.doi.org/10.1016/j.nuclphysa.2008.11.011 [ _ nucl . phys . _ * a817 * ( 2009 ) 6189 ] , http://arxiv.org/abs/0706.3775[arxiv:0706.3775 [ hep - ph ] ] . t.  lappi , `` energy density of the glasma , '' http://dx.doi.org/10.1016/j.physletb.2006.10.017[_phys . lett . _ * b643 * ( 2006 ) 1116 ] , http://arxiv.org/abs/hep-ph/0606207[arxiv:hep-ph/0606207 ] . e.  iancu , k.  itakura , and l.  mclerran , `` a gaussian effective theory for gluon saturation , '' http://dx.doi.org/10.1016/s0375-9474(03)01477-5 [ _ nucl . * a724 * ( 2003 ) 181222 ] , http://arxiv.org/abs/hep-ph/0212123[arxiv:hep-ph/0212123 ] . a.  krasnitz , y.  nara , and r.  venugopalan , `` gluon production in the color glass condensate model of collisions of ultrarelativistic finite nuclei , '' http://dx.doi.org/10.1016/s0375-9474(03)00636-5[_nucl.phys._ * a717 * ( 2003 ) 268290 ] , http://arxiv.org/abs/hep-ph/0209269 [ arxiv : hep - ph/0209269 [ hep - ph ] ] . a.  krasnitz , y.  nara , and r.  venugopalan , `` classical gluodynamics of high - energy nuclear collisions : an erratumn and an update , '' http://dx.doi.org/10.1016/j.nuclphysa.2003.08.004[_nucl.phys._ * a727 * ( 2003 ) 427436 ] , http://arxiv.org/abs/hep-ph/0305112 [ arxiv : hep - ph/0305112 [ hep - ph ] ] . h.  fujii , k.  fukushima , and y.  hidaka , `` initial energy density and gluon distribution from the glasma in heavy - ion collisions , '' http://dx.doi.org/10.1103/physrevc.79.024909[_phys . rev . _ * c79 * ( 2009 ) 024909 ] , http://arxiv.org/abs/0811.0437 [ arxiv:0811.0437 [ hep - ph ] ] . j.  berges , s.  scheffler , and d.  sexty , `` turbulence in nonabelian gauge theory , '' http://dx.doi.org/10.1016/j.physletb.2009.10.032[_phys . _ * b681 * ( 2009 ) 362366 ] , http://arxiv.org/abs/0811.4293[arxiv:0811.4293 [ hep - ph ] ] . j.  berges , a.  rothkopf , and j.  schmidt , `` non - thermal fixed points : effective weak - coupling for strongly correlated systems far from equilibrium , '' http://dx.doi.org/10.1103/physrevlett.101.041603[_phys . lett . _ * 101 * ( 2008 ) 041603 ] , http://arxiv.org/abs/0803.0131 [ arxiv:0803.0131 [ hep - ph ] ] .
Abstract:

Answer: we extensively study the growing behavior of the energy and the pressure components depending on the space - time rapidity in the framework of the glasma , which describes the early - time dynamics in the ultra - relativistic heavy - ion collisions . we simulate the glasma solving the classical equations of motion in the su(2 ) yang - mills theory and systematically investigate the dependence of the glasma instability on the model parameters . we have checked that the transverse and longitudinal grid sizes in our simulation are large enough to handle cutoff effects under control . by comparing the numerical results from several initial conditions with different magnitudes of instability seed and also those with different wave - numbers for rapidity fluctuations , we clearly see that unstable modes dominantly grow up in the linear regime and we also confirm non - linear effects in the time evolution . to extract more detailed information on the evolving glasma , we decompose the energy into the components in terms of rapidity wave - numbers . we observe an energy flow from low wave - number modes into higher wave - number modes due to non - linearity in the equations of motion . we find that the energy spectrum approaches an asymptotic scaling that is consistent with kolmogorov s power - law form even in the expanding system of the glasma .    and    color glass condensate , glasma , relativistic heavy - ion collision , instability , kolmogorov spectrum