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: turbulence is a generic phenomenon observed in the relaxation dynamics of many - body systems far from thermal equilibrium @xcite . it comprises a quasi - stationary flow of energy within certain inertial regimes in momentum space @xcite . correlation functions like the energy spectrum and higher momenta of the velocity distribution exhibit universality and scaling @xcite . this , and quasistationarity are the key characteristics rendering turbulence a nonthermal fixed point of the system s dynamics .    in quantum many - body systems , from the formation of bose - einstein condensates in ultracold gases to quark - gluon plasmas produced in heavy - ion collisions and reheating after early - universe inflation , nonequilibrium dynamics governs many interesting phenomena . in this context , turbulence during thermalization is being studied with increasing effort @xcite . a nonthermal fixed point of the evolution of a many - body system has the potential to strongly affect the equilibration process by forcing the evolution to critically slow down before the system can thermalize . new scaling laws were found by analyzing non - perturbative quantum field dynamic equations @xcite . analogous predictions for a dilute ultracold bose gas were given in @xcite , proposing strong matter - wave turbulence in the regime of long - range excitations . superfluid turbulence , also referred to as quantum turbulence ( qt ) has been the subject of extensive studies in the context of helium @xcite and dilute bose gases  @xcite . in contrast to eddies in classical fluids , vorticity in a superfluid is quantized @xcite , and the creation and annihilation processes of quantized vortices are distinctly different @xcite . the observation of a kolmogorov 5/3-law @xcite in experiments with superfluid helium  @xcite received much attention @xcite . in particular , the role of the normal - fluid as compared to the superfluid component in the turbulent flow is under debate  @xcite . superfluid turbulence plays an important role in the context of the kinetics of condensation and the development of long - range order in a dilute bose gas . this , as well as turbulence in its acoustic excitations has been discussed in refs . a possible observation of qt in ultracold atomic gases presently poses an exciting task for experiments  @xcite . here we emphasize that the experimental study of superfluid turbulence and , more generally , of nonthermal fixed points in ultracold bose gases has strong potential implications for many other areas of physics . besides vortical excitations this also includes other ( quasi-)topological excitations such as solitary waves in one - dimensional gases . a satisfactory ab - initio mathematical description of both quantum and classical turbulence is inherently difficult due to the strong correlations building up within the system . analytical results are known , however , in regimes where kinetic theory applies : in a dilute , degenerate bose gas the normal - fluid component can vary at the expense or gain of the superfluid part . as a consequence , the gas is compressible and so - called weak wave - turbulence phenomena can occur for which scaling laws can be derived by analyzing kinetic equations @xcite .    generically , however , the description in terms of wave kinetic equations such as the quantum boltzmann equation breaks down in the infrared ( ir ) regime of long wavelengths . for a bose gas in this regime , single - particle occupation numbers grow large and the description in terms of , for example , elastic two - to - two collisions becomes unreliable . in the infrared limit , so - called strong wave turbulence is expected to occur . recent developments presented in refs . @xcite allow one to set up a unifying description of scaling , both in the ultraviolet ( uv ) quantum boltzmann kinetic regime and in the infrared limit . in the ir regime , new scaling laws were found for a relativistic scalar field by analyzing non - perturbative kadanoff - baym dynamic equations with respect to nonthermal stationary solutions  @xcite .    in this article , we present the details of our studies of the relaxation of two and three dimensional dilute bose gases through stages of superfluid turbulence and the approach of a nonthermal fixed point , by means of simulations in the classical - wave limit of the underlying quantum field theory . in section  [ sec : dyncritpoints ] , we summarize the quantum - field theoretical predictions of refs . @xcite , in particular for the scaling exponents of the single - particle momentum distribution and compare these with the numerically determined scaling . while we find excellent agreement , our results provide an interpretation of the nonthermal fixed points proposed in refs . @xcite for the case of an ultracold bose gas : the appearance of nonperturbative infrared scaling reflects the presence of statistically independent vortices , as previously reported in brief in ref . @xcite , see also  figs . [ fig:2dphase512]-[fig : modeoccupation3g8 ] . this phenomenon appears to be distinctly different from the weak wave turbulence we observe in the uv part of the spectrum .    in section  [ sec : randomvortexmodel ] , we describe a model of independent as well as pair - correlated vortical excitations @xcite which allows the more refined interpretation of the scaling behavior during the different stages of the evolution presented in section  [ sec : results ] . we show that the stationary scaling is maintained by the presence of energy ( uv ) and particle ( ir ) fluxes , further supporting the analytic theory . in section [ sec : vorticesand ] , we study the power - law distributions of the compressible and incompressible contributions to the flow pattern . this adds to the clear understanding of the bimodal scaling laws in the overall momentum spectrum . ir power - law spectra of subdominant compressible excitations suggest the presence of acoustic turbulence @xcite on the top of the vorticity - bearing quasicondensate . in comparison with recent experiments and analytical predictions , the velocity field probability distribution as well as the velocity statistics of individual vortices are studied . this observable is of great interest , since it has recently been used to experimentally verify a distinction between classical and quantum turbulence @xcite . we finally show that the complete decay of the turbulent scaling is anticipated by the appearance of weaker ir power laws reflecting vortex - antivortex pairing correlations . we begin with a brief summary of the analytical and numerical results on dynamical fixed points and matter - wave turbulence reported in refs . @xcite before we discuss these results in terms of a statistical model of vortex excitations . motivated by the original work presented in ref . @xcite nonthermal fixed points of kadanoff - baym dynamic equations for time - dependent green functions were analysed , for the case of an ultracold bose gas , in ref .  @xcite . at these fixed points of the dynamical evolution of the system , the single - particle momentum distribution @xmath0 becomes ( quasi-)stationary , @xmath1 . it furthermore exhibits a characteristic universal power - law behavior , i.e. , it scales according to @xmath2 in a certain regime of momenta @xmath3 . here @xmath4 is some positive , real number and @xmath5 a universal exponent which was determined from the dynamic equations for the field correlation functions . different exponents resulted in different momentum regimes . exponents valid in the ultraviolet regime of large @xmath6 were found to correspond to well - known fixed points of weak wave turbulence @xcite in the respective systems . in addition to these , new , larger exponents were predicted in the infrared regime of small @xmath6 on the basis of a non - perturbative analysis of the kadanoff - baym dynamic equations derived from the two - particle - irreducible ( 2pi ) effective action @xcite . simulations of the classical field equations for a relativistic @xmath7-symmetric scalar model , for @xmath8 , confirmed the existence of the analytically derived scalings in the infrared regime  @xcite . as reported in ref . @xcite , corresponding simulations were performed for an ultracold , degenerate bose gas in two and three spatial dimensions which demonstrated that for these systems the infrared scalings predicted in @xcite reflect and are caused by the presence of quantized vortical excitations of the superfluid . in this way , the nonthermal fixed points can be related in a clear manner to topological excitations of the interacting coherent matter - wave field which shows that both , weak wave turbulence and macroscopic , topological excitations of the field can be described within a unified field - theoretic approach . differently expressed , both , weak turbulent flow and non - linear solitary bulk excitations are described in a unified manner as representing a nonthermal critical fixed point of the system . in turn , the approach also implies that superfluid turbulence can be studied in a new way , in the frame of a universal quantum field theoretical approach . we remark that the relation between the new infrared exponents and topological excitations of the background ( i.e. , condensate ) field gains further support by the topological pattern formation described in ref . @xcite for a relativistic @xmath9-symmetric scalar model , which serves as a particular model for the reheating period after early - universe inflation . we will discuss in detail vortical excitations in a degenerate bose gas and their relation to the above nonthermal fixed points . before this we summarize briefly the relevant results of refs . @xcite . to be specific , throughout this article we consider an ultracold bose gas of atoms in @xmath10 , @xmath11 space dimensions interacting through @xmath4-wave collisions which is described by the hamiltonian @xmath12 \,,\ ] ] where the time and space dependent fields @xmath13 satisfy bose commutation relations , and where the coupling @xmath14 in three dimensions is defined in terms of the @xmath4-wave scattering length @xmath15 . here and in the following we set @xmath16 , see appendix [ app : classfields ] for more details . suppose the generic case that for sufficiently large momenta @xmath17 occupation numbers @xmath18 are small enough such that for a given coupling @xmath19 the quantum boltzmann equation ( qbe ) @xmath20 ,    \label{eq : kinscattint}\end{aligned}\ ] ] describes the evolution of @xmath21 under the effects of collisions . in our case , the transition matrix element squared @xmath22 is a numerical constant proportional to @xmath23 and thus independent of momenta . zeroes of the scattering integral @xmath24 correspond to fixed points of the time evolution within the regime of applicability of the qbe @xcite . most prominent amongst these are the thermal fixed point corresponding to the system in thermal equilibrium and the trivial fixed point where the occupation number is independent of @xmath3 . at both fixed points the scattering integral vanishes and @xmath0 becomes independent of @xmath25 . note that both , the trivial distribution and the thermal bose - einstein distribution in the rayleigh - jeans regime , for @xmath26 , show a power - law behavior of the form ( [ eq : scale ] ) with @xmath27 and @xmath28 , respectively . the theory of weak wave turbulence @xcite allows to analytically derive further , nonthermal fixed points at which the occupation number @xmath29 obeys a scaling law of the form ( [ eq : scale ] ) and , in general , @xmath30 . as in classical turbulence of an incompressible fluid one assumes that universal scaling appears within a certain regime of momenta , the inertial range . according to this picture , outside the scaling regime kinetic energy enters the system from an external source and/or is dissipated into heat , whereas there are no sources and sinks within the inertial interval . instead , kinetic energy is transported from momentum shell to momentum shell without loss or gain . to a good approximation , this process is described by a continuity equation in momentum space , with a momentum - independent , radially oriented current vector . a central aspect of weak - wave - turbulence theory is that the qbe can be cast into different such equations  @xcite , for the radial number density @xmath31 and the energy density @xmath32 , @xmath33 , @xmath34 depending on whether the radial particle current @xmath35 or energy current @xmath36 is taken to be independent of @xmath37 , one derives different scaling exponents . the resulting exponents are @xmath38 these exponents can be obtained by simple power counting : combining eqs . ( [ eq : qkineq ] ) and ( [ eq : baleqq ] ) gives the radial relation @xmath39 which implies that stationarity requires @xmath40 to become @xmath37-independent , i.e. scale like @xmath41 . counting all powers of @xmath37 in @xmath42 , eq . ( [ eq : kinscattint ] ) , in the wave - kinetic regime where the terms of third order in the occupation numbers dominate the scattering integral , this requires @xmath43 . analogously one infers the exponent @xmath44 from the balance equation ( [ eq : baleqp ] ) for the energy density @xmath45 . despite this simple procedure , the existence of the respective scaling solutions has to and can be derived rigorously from the qbe by means of zakharov conformal integral transforms  @xcite . we note that , similar to classical turbulence , the case @xmath10 is special , where the scaling exponent @xmath46 equals that of a thermal distribution in the rayleigh - jeans regime . given a positive scaling exponent @xmath5 momentum occupation numbers @xmath47 grow large in the ir regime of small @xmath37 . as a consequence , for a given coupling @xmath19 , the qbe fails in this regime , where contributions to the scattering integral @xmath42 which are of higher order than @xmath48 become important . to find scaling solutions in the ir , an approach beyond kinetic theory is required . this is available through quantum - field dynamic equations derived from the two - particle irreducible ( 2pi ) effective action or @xmath49-functional @xcite beyond the 2-loop order in the expansion of the self - energy . see ref . @xcite for details of the procedure summarized in the following . the 2pi equations include the dyson equation for the time - ordered green function @xmath50 ( here , we use four - vector notation @xmath51 ) , from which a time evolution equation ( [ eq : qkineq ] ) for @xmath29 is derived . as before one considers zeros of the scattering integral which in the dynamic theory reads @xmath52 . \label{eq : stationarity}\ ] ] here , @xmath53 , and @xmath54 and @xmath55 are the spectral and statistical components of @xmath56 , respectively , defined in coordinate space by @xmath57 , @xmath58\rangle$ ] , @xmath59 . the corresponding contributions to the self energy @xmath60 are defined in terms of @xmath56 through a loop expansion of the 2pi effective action , see fig . [ fig:2pi ] . resumming an infinite set of such loop diagrams contributing to the 2pi effective action @xcite leads to a non - perturbative , effectively renormalized coupling in the dynamic equations @xcite .    to derive the scaling solutions of the dynamic equations for green functions @xmath61 one assumes separate scaling of its spectral and statistical components according to @xmath62 , @xmath63 , @xmath64 . here , @xmath65 is the dynamical scaling exponent accounting for a different scaling in @xmath66 as compared to @xmath3 . the scaling exponent @xmath5 is related to @xmath67 by @xmath68 . @xmath67 is derived from the condition that the scattering integral ( [ eq : stationarity ] ) vanishes , making use of zakharov transformations , while the anomalous scaling exponent @xmath69 remains undetermined by this and , for the first , can be assumed to vanish for the @xmath10 and @xmath11 dimensional cases considered in this paper . the ir scaling exponents for radial quasiparticle flow ( @xmath70 ) and radial energy flow ( @xmath71 ) in @xmath72 dimensions were predicted in ref . @xcite to be @xmath73 in situations where a quasiparticle picture applies the dynamical exponent @xmath65 corresponds to the homogeneity index of @xmath66 : @xmath74 . in this case the nonperturbatively resummed effective coupling can be related to the diagonal elements of an effective many - body @xmath75-matrix , @xmath76 in the kinetic boltzmann formulation . in the scaling regimes this scales as @xmath77|$ ] , @xmath78 , where @xmath79 is some constant which fine - tunes the position of the transition from uv to ir scaling : for small @xmath80 and @xmath81 one recovers the uv case discussed in the previous section , i.e. , @xmath82 is a constant independent of @xmath37 . for large @xmath80 , the second term in the denominator dominates which , assuming scaling of @xmath83 , implies a power - law behavior @xmath84 and , in turn , the modified scaling ( [ eq : kappair ] ) of @xmath85 in the infrared regime of small wave numbers as compared to the uv regime discussed before . physically , the renormalized @xmath75-matrix implies a reduction of the effective interaction strength in the ir regime of strongly occupied modes @xcite . as a consequence , single - particle occupation numbers rise , towards smaller wave numbers , in an even steeper way than in the weak - turbulence regime .     at the time shown in fig . [ fig:2dphase512 ] . @xmath37 is in lattice units , see app . [ app : classfields ] . note the double - logarithmic scale . the development of a ir scaling with @xmath86 coincides with the presence of superfluid vortices , with additional wave - turbulent or thermal background short - wavelength fluctuations in the uv . dashed lines indicate the filling of the initially occupied modes with @xmath87.,scaledwidth=40.0% ]     at the time shown in fig . [ fig:2dphase512 ] . @xmath37 is in lattice units , see app . [ app : classfields ] . note the double - logarithmic scale . the development of a ir scaling with @xmath86 coincides with the presence of superfluid vortices , with additional wave - turbulent or thermal background short - wavelength fluctuations in the uv . dashed lines indicate the filling of the initially occupied modes with @xmath87.,scaledwidth=45.0% ]     for the snapshot in fig . [ fig : phase512 ] . @xmath37 is in lattice units . note the double - logarithmic scale . the development of a ir scaling @xmath88 coincides with the presence of superfluid vortex lines . dashed lines indicate the filling of the initially occupied modes with @xmath87.,scaledwidth=41.0% ]     for the snapshot in fig . [ fig : phase512 ] . @xmath37 is in lattice units . note the double - logarithmic scale . the development of a ir scaling @xmath88 coincides with the presence of superfluid vortex lines . dashed lines indicate the filling of the initially occupied modes with @xmath87.,scaledwidth=45.0% ]      we now briefly review the relation between wave - turbulent scaling and the appearance of vortical excitations in an ultracold degenerate bose gas . in @xcite the results of semiclassical simulations of the gas dynamics were reported , obtained by solving the classical field equation @xmath89 \phi(\mathbf{x},t)\ ] ] in a box with periodic boundary conditions . see appendix [ app : classfields ] for details on the simulations and on lattice units . the initial field @xmath90 was prepared by macroscopically populating a few of the lowest momentum modes in the computation such that the resulting condensate density in configuration space varied between zero and some maximum value . quantum noise is taken into account by adding a small random contribution to each field mode . vortical excitations were created in large numbers , within shock - waves forming during the non - linear evolution of the coherent matter - wave field . [ fig:2dphase512 ] shows an example of the condensate phase distribution over a numerical lattice in @xmath10 dimensions , with vortices and antivortices indicated by white crosses and circles , respectively . [ fig : phase512 ] depicts the location of vortex - ring cores in a @xmath91-dimensional system . for videos of the evolution see  @xcite . both figures show a snapshot of the system at some time after the vortical excitations have formed . at these times , in @xmath10 , part of the vortex - antivortex pairs have re - annihilated with each other already , and in @xmath91 some of the rings have undergone reconnections and also disappeared by shrinking to zero size . runs were repeated many times ( @xmath92 times in @xmath10 and @xmath93 times in @xmath91 ) for an ensemble of initial configurations differing by statistical noise . the number of runs was chosen such that the statistical error arising from run - to - run fluctuations was reduced to a value on the order of the size of the symbols used in the figures . in the regime of momenta with large occupation numbers where quantum statistical fluctuations play little role , correlation functions like the spectrum @xmath80 could be computed in a quasi exact way by averaging at a given time over the ensemble . in this way , the flow was analysed in terms of the ensemble- and angle - averaged single - particle spectrum @xmath94 as a function of radial momentum @xmath78 . [ fig:2dmodeoccupation ] and [ fig : modeoccupation3g8 ] show this spectrum for @xmath10 and @xmath11 , at the times chosen in figs . [ fig:2dphase512 ] and [ fig : phase512 ] , respectively . most importantly , these spectra show scaling , eq . ( [ eq : scale ] ) , within a range of momenta between about the maximum of the initially occupied momentum modes indicated by the vertical lines and @xmath95 in lattice units . the exponent @xmath5 is in agreement with the field theoretical prediction @xmath96 given in eq . ( [ eq : kappauv ] ) . it was found that the annihilation processes which eventually destroy the vortical structure while coherence builds up in the system are very slow and thus stabilize the scaling solution over a long time . in particular , only after the last vortex - antivortex pair has annihilated and the last ring shrunk to zero the power law @xmath97 breaks down and a thermal distribution of particles is left over . characteristic times are @xmath98 for vortex formation , @xmath99 for stabilization of the scaling solution and @xmath100 for the last vortices to annihilate . as we have seen above , the scaling of the momentum occupation numbers in the ir regime of large wave numbers , predicted within a non - perturbative analysis of strong wave - turbulence , correspond , for @xmath101- and @xmath11-dimensional degenerate bose gases , to the appearance of macroscopic vortical excitations . in the following , we analyse the observed scaling spectra by comparing them to the single - particle momentum distributions for a set of randomly positioned point vortices ( vortex lines ) in @xmath10 ( @xmath91 ) dimensions . we will show that uncorrelated vortices ( vortex lines ) are sufficient to yield the infrared scaling with exponent @xmath102 , see eq . ( [ eq : kappair ] ) . we will also show that beyond this , pair correlations between vortices and antivortices as well as configurations with small rings well separated from each other can give rise to a further scaling exponent deviating from @xmath102 ( see fig . [ fig : vortexsystemrandom ] ) .    the point ( line ) vortex model employed here was introduced by onsager in 1949 @xcite . it describes the complex flow pattern in terms of the statistical mechanics of interacting classical point objects . this model has been constructed as a discrete vorticity approximation of classical fluid turbulence , but it is even more suitable to describe superfluid turbulence consisting of quantized vortices . further applications include plasma physics or stellar dynamics @xcite . in the following , we study the spectrum of a set of vortices , finding different scaling regimes dependent on vortex correlation functions . with this approach , nonthermal fixed points in an ultracold bose gas can be related to the statistics of vortices . numerical calculations , sampling field configurations of static randomly positioned vortices , confirm the analytical result . finally , the velocity probability distribution of a vortex - dominated flow  @xcite is derived and numerically confirmed .           in two dimensions , an isolated , singly quantized vortex is described by the complex field @xmath103 . as the @xmath104-dependence of the density @xmath105 only becomes important at small scales on the order of the healing length @xmath106 where our simulations are dominated by thermal excitations , we will omit it in the following and assume @xmath107 to be uniform . a system of @xmath108 vortices in @xmath10 dimensions can be defined as @xmath109 , where @xmath110 is the single vortex field centered around @xmath111 . we derive the corresponding particle spectrum by considering the hydrodynamic velocity field @xmath112 . denoting the velocity field of a single vortex as @xmath113 ( see appendix [ app : singlevortex ] ) we can express the mean classical kinetic energy density of the velocity field as @xmath114 where @xmath115 defines the spatial distribution of vortices with winding number @xmath116 . here and in the following , @xmath117 denotes an ensemble average over different realizations of the classical field @xmath118 . we derive the low-@xmath37 scaling of @xmath119 from the kinetic - energy spectrum @xmath120 , given by the angle - averaged fourier transform of @xmath121 , taking into account that at low @xmath37 , the single - particle spectrum is dominated by the superfluid velocity field @xmath122 , i.e. , @xmath123 one has , from eq . ( [ eq : energyspec ] ) @xmath124 with @xmath125 ( lattice units ) , the modulus of the velocity field of a single vortex scales as @xmath126 and is radially symmetric in momentum space . hence , the angle - averaged single - particle spectrum scales like @xmath127 here , @xmath128 denotes the zeroth - order bessel function and @xmath129 is the distance between vortices @xmath130 and @xmath131 .    assuming that the positions @xmath132 of the vortices are uncorrelated one can take the average over relative positions @xmath133 within the area @xmath134 , @xmath135 and , for fixed @xmath37 , the limit @xmath136 . hence , the second term in brackets in eq . ( [ eq : ve ] ) vanishes and one finally obtains the scaling  @xcite @xmath137      in a two - dimensional superfluid containing vortices and antivortices an effectively attractive force between the two species can lead to pair correlations . we study a signature of this feature in the single - particle momentum spectrum by applying the point vortex model introduced above to the case of vortex - antivortex pairs .    as a first step we calculate the velocity field @xmath138 for a vortex - antivortex pair with the vortex situated at @xmath139 and the antivortex at @xmath140 , @xmath141 the squared velocity field far away from the center of the pair can be obtained via a dipole approximation @xmath142 which yields the scaling @xmath143 . hence , in fourier space , the pair velocity field scales as @xmath144 for low momenta . in this regime , the vortex - antivortex pair can again be treated as a point - like object with modified velocity scaling . to obtain the infrared scaling of a set of random vortex - antivortex pairs , we define a spatial pair distribution @xmath145 , with @xmath111 denoting the center of the @xmath130-th vortex - antivortex pair . then , the analysis performed for random vortices above can be adopted . therefore , in the case of independently distributed pairs , the approach predicts the infrared scaling of the occupation number to be the same as for a single pair , i.e. @xmath146 .      the considerations from sect . [ subsec : indpairs ] are only valid in the far ir ( or equivalently for pair size going to zero ) . in our numerical simulations we found , however , that for large evolution times , a finite minimum distance emerged between vortices and antivortices , and also between vortices with the same circulation , see sect . [ subsec : pairingspectrum ] . to take into account this observation and to analyse the full spectrum , we go back to writing the distribution @xmath147 as the sum of distributions @xmath148 of @xmath108 vortices and @xmath149 of @xmath108 antivortices . hence , @xmath150 with @xmath151 . this allows for a derivation of the kinetic - energy distribution in terms of correlation functions of vortex positions . we model pairing by the density - density correlation functions @xmath152 where @xmath153 is the area in which the theta function equals one . the contributions @xmath154 take into account that , besides the pairing , vortices and antivortices avoid each other in the dilute gas , keeping a minimum distance @xmath155 . this is due to vortex - vortex repulsion and fast vortex - antivortex annihilation on small distances . the functions @xmath156 cancel out in eq . ( [ eq : cf ] ) . apply for vortices and antivortices , the terms do not cancel , but the remaining term does not alter the results for pair scaling derived here . ]    inserting eqs . ( [ eq : rhovrhov ] ) , ( [ eq : rhovrhoa ] ) into eq . ( [ eq : cf ] ) , evaluating the fourier transform and angular averaging gives @xmath157 the expansion of the integral for @xmath158 yields the leading - order result @xmath159 and , from eqs . ( [ eq : nkfromevk ] ) and ( [ eq : evk ] ) , the same occupation - number scaling as for independent pairs , @xmath160 at momenta @xmath161 the independent - vortex scaling @xmath86 is restored . we note that the infrared result ( [ eq : ps ] ) can also be achieved by choosing more general pair correlations of the form @xmath162^{-1}|\mathbf{x}-\mathbf{x'}|^{\epsilon-2}\theta(\lambda_{\mathrm{max } } - |\mathbf{x}-\mathbf{x'}| ) \theta (   |\mathbf{x}-\mathbf{x'}| -\lambda_{\mathrm{min } } ) +   p_{\mathbf{x},\mathbf{x}'}$ ] , with @xmath163 and @xmath164 , which is closer to our numerical observations and includes the case of a fixed `` binding length '' @xmath165 . as before , pair scaling @xmath166 is found for @xmath167 , critical scaling @xmath168 for @xmath169 , irrespective of @xmath170 . bound vortex - antivortex pairs , on a @xmath171 grid . note the double - logarithmic scale . the vortices are positioned according to the probability distribution ( [ eq : pd ] ) . at low momenta the power law is consistent with the existence of random vortex pairs @xmath172 . above @xmath173 , the distribution exhibits the scaling of an ensemble of independent vortices , @xmath174 , while at momenta larger than the healing - length scale @xmath175 one observes the vortex - core scaling @xmath176 . , scaledwidth=45.0% ]    we supplement our discussion of scaling in the point vortex model with numerical data obtained by averaging over an ensemble of field configurations in @xmath10 dimensions . these configurations were constructed by multiplying uncorrelated single - vortex fields centered at positions according to the probability distribution @xmath177 with @xmath178 , where we neglect that unpaired vortices avoid each other , i.e. , choose @xmath179 . they do not represent stable solutions of the classical field equation nor do they contain sound - wave and related excitations yet which would build up through the interactions of the field and the vortices .    the resulting momentum spectrum is presented in fig .  [ fig : randomvortex ] . three scaling regimes can be observed . at low momenta the power law is consistent with the existence of random vortex pairs @xmath172 . above @xmath180 , choosing @xmath181 ( lattice units ) , the distribution exhibits the scaling of an ensemble of independent vortices , @xmath174 , while at momenta larger than the healing - length scale @xmath175 one observes the vortex - core scaling @xmath176 . the above results reflect that in a vortex dominated flow , particles with low momenta are found far away from vortex cores . in the case of pairing , the flow field far away from the cores is given by the field of a vortex pair , and the low - momentum scaling follows the pair - field scaling . particles closer to the vortex cores pick up a higher momentum . above @xmath173 the field will be dominated by the field of a single vortex . note that fig . [ fig : randomvortex ] shows the result of a numerical calculation in which we have sampled field configurations of random static vortices . dynamical simulations of the classical field equation eq . ( [ eq : gpe ] ) will be shown in sect . [ sec : results ] . we close with recalling the picture onsager developed in ref . @xcite of thermodynamic equilibrium states of a fixed number of vortices and antivortices in two dimensions . he used the hamiltonian of vortical flow in two dimensions @xcite , @xmath182 to describe the dynamics of a system of @xmath108 vortices in a superfluid which hence behave like a coulomb gas . here , the position of the @xmath130-th vortex is denoted as @xmath183 . due to the fact that the @xmath184 and @xmath185 coordinates of each vortex are canonical conjugates , phase space is identical with configuration space of the vortex positions . hence , for vortices moving in a volume @xmath186 the total phase space is given by @xmath187 . the hamiltonian ( [ eq : hamiltonian ] ) implies that low - energy configurations feature vortices of opposite sign close to each other , whereas high - energy configurations require vortices of equal sign to group . due to these constraints , the number of configurations @xmath188 available for the system at a given energy @xmath189 decreases towards high and low energies , with a maximum at some intermediate @xmath190 . according to boltzmann , the entropy is @xmath191 and the inverse temperature @xmath192 is positive for @xmath193 and negative for @xmath194 . it follows that positive - temperature states are characterized by vortex - antivortex pairing , while negative - temperature states feature vortices of the same circulation to cluster . at the point of maximum entropy @xmath195 and infinite temperature , onsager expected a state of uncorrelated vortices and antivortices . this understanding of the nonthermal fixed point as a quasi - equilibrium state of vortices corroborates our result that the ir scaling behavior at the fixed point corresponds to the appearance of topological excitations .       grid . note the double - logarithmic scale . the major and minor radius scales @xmath196 , for @xmath197 , and @xmath198 , for @xmath199 , respectively , are indicated . these scales separate pair scaling @xmath146 as for a near - circular vortex ring , scaling @xmath200 for two anti - circulating vortex lines , and @xmath201 , as for a vortex ring or of a pair of straight vortex lines , corresponding to the scaling exponent @xmath202 at the nonthermal fixed point . see sect . [ subsec : vortexloops ] for more details.,scaledwidth=45.0% ]    we now consider the three - dimensional case of vortex lines and loops . a formulation similar to the onsager point vortex model in sect . [ subsec : independentvortices ] is possible  @xcite . see also refs . @xcite for an extension including the dynamics of tangles of vortex lines . we write the classical kinetic energy spectrum of the velocity field in terms of the vorticity density @xmath203 as @xmath204 the vorticity is vanishing everywhere but on the vortex lines , i.e. , assuming @xmath108 individual vortex loops , @xmath205 in this expression , the vortex filaments are represented by the connected curves @xmath206 , parametrized by the one - dimensional coordinate @xmath207 , @xmath208 being the arc length of filament @xmath130 . the loops are closed , @xmath209 , possibly also accross the walls of the 3-dimensional volume in accordance with periodic boundary conditions . @xmath210 is the tangent vector along the filament at @xmath206 , of unit length @xmath211 . we parametrize the @xmath130-th vortex loop in terms of a single center coordinate and a relative curve , @xmath212 and write the vorticity as @xmath213 hence , in fourier space , @xmath214 where @xmath215 is the fourier transform of the vorticity of the @xmath130-th vortex loop , given by @xmath216 the ensemble - averaged vorticity @xmath217 becomes @xmath218 assuming that the shapes of the individual vortex loops are statistically independent of their position , it follows that @xmath219 . if the loops are also uncorrelated among themselves , then @xmath220 . hence , for statistically identical loops , @xmath221 which means that the vorticity spectrum scales in the same way as the average vorticity of a vortex loop centered at the origin . finally , the scaling of the momentum spectrum @xmath222 follows from that of the angle - averaged vorticity , @xmath223 for the case of two straight parallel vortex lines of opposite circulation eq . ( [ eq : energyspeccorrfinal ] ) is evaluated in appendix [ app : vortexline ] , for an ensemble of such paired lines in appendix [ app : vortexlines ] . making use of the procedure developed there , the case of circular vortex rings of radius @xmath104 is discussed in appendix [ app : vortexring ] . for the latter , the resulting angle - averaged momentum spectrum scales like @xmath224 for momenta @xmath225 and @xmath88 for momenta @xmath226 . to include effects from squeezed vortex loops , we consider elliptical filaments . the ellipse is defined by a major radius @xmath227 and minor radius @xmath228 . in fig . [ fig : vortexellipse ] , the angle - averaged momentum spectrum of an elliptical vortex loop is shown . three scaling regimes can be distinguished . for the lowest momenta , one has @xmath146 , which equals the infrared scaling for a vortex ring . for momenta @xmath229 , one finds @xmath200 , which coincides with the infrared scaling of two anti - circulating vortex lines ( see app .  [ app : vortexline ] ) . for the ellipse , above @xmath198 , the momentum distribution scales like @xmath201 , which is the same as the high - momentum scaling of a vortex ring ( see app .  [ app : vortexring ] ) or of a pair of straight vortex lines ( app . [ app : vortexlines ] ) .      ) as obtained from a random distribution of @xmath230 bound vortex - antivortex pairs , on a @xmath171 grid , as in fig . [ fig : randomvortex ] . note the double - logarithmic scale . at low momenta the velocity field follows the scaling of vortex - antivortex pairs ( @xmath231 ) , whereas at higher momenta the distribution reflects the presence of vortices ( @xmath232 ) . note , that there is a deviation from @xmath231 in the regime of low velocities . from simulations with higher resolution , we find the best agreement for @xmath233 . we remark , that in the limit @xmath234 the velocity distribution is expected to be gaussian @xcite.,scaledwidth=45.0% ]    the vortex model can provide insight into another observable accessible in our numerical simulations , the velocity distribution . @xcite for discussions of the velocity distribution and ref . @xcite for recent experimental results . as discussed above , the velocity field scales , far away from any core , as @xmath235 , with @xmath232 for a single vortex and @xmath231 for a vortex - antivortex pair in two dimensions . the velocity probability distribution @xmath236 is calculated as @xmath237 , with spatial vortex distribution function @xmath238 . for a uniform distribution @xmath239 , it follows after angular averaging that @xmath240 it is numerically convenient to calculate the probability density of a single component of the field , e.g. @xmath241 in fig . [ fig : vortexvelocity ] , the corresponding velocity distribution @xmath242 is shown . in accordance with the analytical predictions two scaling regimes are observed . at low momenta the velocity field follows the scaling of vortex - antivortex pairs , whereas at higher momenta the distribution reflects a random distribution of vortices . in this section we return to the process of the formation of vortical excitations and of the bose gas approaching the nonthermal fixed point characterized by the wave - turbulent scaling solutions discussed in sect . [ sec : dyncritpoints ] and refs . thereby we first focus on the evolution of the single - particle momentum distributions and the emergence of nonthermal power laws . we relate these to structure formation in the form of vortical excitations . in order to identify in a clearer way the relevant processes in approaching the fixed point we compute the fluxes in momentum space . remarkably , we find that the appearance of the particle scaling exponent @xmath243 in the ir and of the energy exponent @xmath244 in the uv are compatible with predictions based on general arguments of irreversibility as well as energy and number conservation in a collision process  @xcite .      , for four different times of a run in @xmath10 dimensions . parameters are : @xmath245 , @xmath246 , @xmath247 . note the double - logarithmic scale . an early development of a scaling @xmath248 is followed by a quasi - stationary period of bimodal scaling with @xmath86 in the ir , due to the presence of vortices , and @xmath224 in the uv , corresponding to weak wave turbulence or thermal equilibrium.,scaledwidth=48.0% ]    , for the four different times of a run in @xmath91 dimensions . parameters are : @xmath249 , @xmath250 , @xmath251 . note the double - logarithmic scale . a period of bimodal scaling @xmath88 ( vortex lines , ir ) and @xmath248 ( weak wave turbulence , uv ) is followed by trimodal scaling , which also exhibits pairing , i.e. , @xmath248 in the far ir , induced by a set of far separated small vortex rings , in addition to @xmath88 ( ir ) , and thermal scaling @xmath224 in the uv.,scaledwidth=48.0% ]    as discussed in sect . [ subsec : summaryvortexnumerics ] , vortical excitations can be created in large numbers , within shock waves forming during the non - linear evolution of the coherent matter - wave field . [ fig : modeoccupationsequence2 ] and [ fig : modeoccupationsequence3 ] show the angle- and ensemble - averaged radial momentum spectra , eq . ( [ eq : spectrumdef ] ) , for a bose gas in a box with periodic boundary conditions , in @xmath10 and @xmath91 dimensions , respectively . four snapshots are shown , taken at the dimensionless times @xmath252 as indicated in each panel . see appendix [ app : classfields ] for details on the simulations and on lattice units . the initial field configurations were prepared by macroscopically populating a few of the lowest momentum modes in the computation such that the resulting condensate density in configuration space varied between zero and some maximum value . in figs . [ fig : modeoccupationsequence2 ] and [ fig : modeoccupationsequence3 ] , the panels representing the earliest times show the system after a brief initial evolution during which momentum is rapidly transported from the few initially occupied modes near zero over a comparatively large range of wave numbers .    shortly after this , vortical excitations are created which immediately causes the spectrum in @xmath10 dimensions to exhibit a power - law behavior @xmath253 within a range of momenta @xmath254 $ ] , see the upper right panel of figs . [ fig : modeoccupationsequence2 ] . subsequently , the evolution slows down and a quasi - stationary period is entered . during an intermediate stage ( bottom left panel of fig . [ fig : modeoccupationsequence2 ] and upper right panel of fig . [ fig : modeoccupationsequence3 ] ) of the vortex - bearing phase two distinct power laws develop which are in excellent agreement with the analytical prediction in eqs . ( [ eq : kappauv ] ) and ( [ eq : kappair ] ) . while in the ultraviolet the exponent @xmath255 exhibits weak wave turbulence , eq . ( [ eq : kappauv ] ) , in the infrared , the exponent confirms the field theory prediction @xmath256 , cf . ( [ eq : kappair ] ) . more specifically , in two dimensions , at @xmath257 , one observes scaling exponents @xmath258 within a range of momenta @xmath259 $ ] and @xmath260 within a range of momenta @xmath261 $ ] , see the lower left panel of figs .  [ fig : modeoccupationsequence2 ] . at @xmath262 , @xmath263 within a range of momenta @xmath259 $ ] , see the lower right panel of figs .  [ fig : modeoccupationsequence2 ] . in @xmath91 dimensions , at @xmath264 , one observes @xmath265 within a range of momenta @xmath266 $ ] and @xmath267 within a range of momenta @xmath268 $ ] , see the upper right panel of figs . [ fig : modeoccupationsequence3 ] . at @xmath269 , @xmath270 within a range of momenta @xmath271 $ ] , see the lower left panel of figs . [ fig : modeoccupationsequence3 ] . the appearance of the bimodal power laws corroborates results for a relativistic @xmath7 model reported in refs . @xcite .    during the ensuing evolution , the weak - wave - turbulence scaling decays towards @xmath28 , reflecting a thermal uv tail . note that in @xmath10 , the weak - turbulence exponent @xmath272 is identical to that in thermal equilibrium in the rayleigh - jeans regime , @xmath273 @xcite . in @xmath91 we observe , at late times , a change of the infrared scaling behavior from @xmath274 to @xmath275 , pointing to the development of pairing correlations , cf .  sects . [ subsec : correlatedvortices ] , [ subsec : vortexloops ] , and sect . [ subsec : pairingspectrum ] below . more specifically , at @xmath276 , one observes a scaling exponent @xmath277 within a range of momenta @xmath278 $ ] , see the lower right panel of figs .  [ fig : modeoccupationsequence3 ] . at late times , after the last vortical excitations have disappeared , we observe the entire spectrum to become thermal , i.e. , exhibit rayleigh - jeans scaling with @xmath28 ( not shown ) . we emphasize that thermal scaling of the single - particle occupation number with @xmath28 applies despite the fact that quasiparticles with a linear dispersion are expected to thermalize in the regime of wave numbers smaller than the inverse healing length @xmath279 ( @xmath10 ) and @xmath280 ( @xmath10 ) . in the bogoliubov approximation this is seen by taking into account the power - law dependence of the coefficients @xmath281 which contribute to @xmath282 @xcite . we are going to study more details of the time evolution of the system in a forthcoming publication and focus in the remainder of this paper on properties of the stationary scaling distributions . an important question in this context is , why the system selects the particular exponents @xmath255 and @xmath256 from the set of four possible exponents given in eqs . ( [ eq : kappauv ] ) and ( [ eq : kappair ] ) . for this , the fluxes underlying the stationary but non - equilibrium distributions are relevant . ( upper set ) and @xmath91 ( lower set ) , for the four different snapshots of figs . [ fig : modeoccupationsequence2 ] and [ fig : modeoccupationsequence3 ] , respectively . note the logarithmic @xmath37-axes . the appearance of the bimodal scaling coincides with a positive kinetic energy flux in the uv , and a negative particle flux in the ir . flux units : @xmath283=[gnq]=(4m^2a_s^{d+4})^{-1}$ ] , cf . [ app : classfields ] . , title="fig:",scaledwidth=45.0% ] +   ( upper set ) and @xmath91 ( lower set ) , for the four different snapshots of figs . [ fig : modeoccupationsequence2 ] and [ fig : modeoccupationsequence3 ] , respectively . note the logarithmic @xmath37-axes . the appearance of the bimodal scaling coincides with a positive kinetic energy flux in the uv , and a negative particle flux in the ir . flux units : @xmath283=[gnq]=(4m^2a_s^{d+4})^{-1}$ ] , cf . [ app : classfields ] . , title="fig:",scaledwidth=45.0% ]    the timeline of distributions shown in figs . [ fig : modeoccupationsequence2 ] and [ fig : modeoccupationsequence3 ] suggests that the evolution of the gas involves a transport process from intermediate momenta around @xmath284 both towards lower and higher wave numbers , building up a bimodal power - law distribution . to describe the character of these transport processes we plot , in fig . [ fig : pflux2d3d ] for @xmath10 ( upper set of four panels ) and @xmath91 ( lower set ) , the radial particle and kinetic - energy flux distributions @xmath285 and @xmath286 , respectively . note that the radial particle flux density @xmath285 is multiplied by @xmath287 to have the same units as the energy flux density @xmath288 . these flux densities are defined through the balance equations ( [ eq : baleqq ] ) and ( [ eq : baleqp ] ) , respectively , with kinetic energy density @xmath289 . the initial stage is governed by an infrared particle transport towards larger @xmath37 , also causing a flux @xmath290 of interaction energy . at intermediate times , @xmath291 changes sign . this is accompanied by a positive kinetic energy transport in the uv , as observed in two - dimensional simulations in ref . @xcite . we emphasize that the negative particle flux in the ir and the positive kinetic - energy flux in the uv coincide with the appearance of the bimodal momentum distributions in figs . [ fig : modeoccupationsequence2 ] and [ fig : modeoccupationsequence3 ] ( bottom left panels ) . although the derivation of the ir exponents requires the full dynamical theory with non - perturbatively resummed self - energies , the signs of the fluxes correspond to the respective scaling exponents , i.e. , @xmath202 in the ir , and @xmath292 in the uv , cf . [ subsec : wwt ] . moreover , at late times , the kinetic - energy flux @xmath71 vanishes due to a thermalized uv momentum distribution , but @xmath70 still reshuffles particles and therefore interaction energy , keeping the system out of equilibrium close to the nonthermal fixed point . we finally remark that , as discussed in detail in ref . @xcite , one can show on the general grounds of the boltzmann @xmath293-theorem that a necessary condition for a non - equilibrium stationary distribution in our systems is energy damping in the region of large @xmath37 . moreover , energy and particle number conservation in the interaction of different momentum modes can be shown to imply the existence of at least one more sink , i.e. , a region where the right - hand sides of eq . ( [ eq : baleqq ] ) effectively has an additional damping term @xmath294 , with negative @xmath295 . in between these sinks , a source region supplies the input to the bidirectional flux pattern towards the uv and ir . it was shown , moreover , in ref . @xcite , within wave - kinetic theory , that generally , a positive , @xmath37-independent flux transports energy , @xmath296 while a negative flux transfers particles , @xmath297 . the only exception is the case of @xmath10 where the thermal and the weak - wave - turbulence exponent @xmath292 can not be distinguished . remarkably , this pattern remains valid in our case , beyond the uv weak - wave - turbulence regime , in the ir region where the exponent emerges from a fixed point of the full dynamic equations . as already pointed out in ref . @xcite , however , the derivation of the ir exponent @xmath256 requires the existence of a sufficiently well defined quasi - particle dispersion relation , suggesting a treatment in terms of the quantum boltzmann equation with a momentum dependent scattering matrix element to be applicable .    from this point of view , the negative flux @xmath70 and scaling in the ir and the positive flux @xmath71 and weak wave turbulence in the uv , as observed in the numerics , emerge as a necessary consequence of conservation laws and transport processes described by wave - kinetic transport equations . in this last section we analyse the single - particle momentum spectra obtained in our numerical simulations in view of their interpretation in terms of vortical excitations and wave - turbulence as implied by the point- and line - vortex models introduced in sect .  [ sec : randomvortexmodel ] . for this we first decompose the flow pattern of the system which has closely approached the nonthermal fixed point , into transverse ( incompressible ) and longitudinal ( compressible ) contributions . in this way we can show that the ir scaling is dominated by the incompressible part while in the uv the particles mainly belong to the compressible as well as a quantum pressure components . we find a further scaling exponent @xmath298 for the subdominant compressible component in the ir which is interpreted as a signature of acoustic turbulence .    moreover , going beyond the results presented in @xcite we identify signs of pair formation in the final stage . the scaling caused by this goes beyond the predictions from dynamical quantum field theory summarized in section  [ sec : dyncritpoints ] , and we interpret it as a signature for the system leaving the fixed point again for final thermalization . , as defined in eqs . ( [ eq : ev])([eq : decomposition ] ) : total single - particle occupation number @xmath119 ( black dots ) , incompressible ( solenoidal - flow ) component @xmath299 ( red circles ) , compressible ( rotationless ) component @xmath300 ( blue filled squares ) , quantum - pressure component @xmath301 ( grey open squares ) . parameters are the same as in fig . [ fig : modeoccupationsequence2 ] ( lower right panel ) , for the run in @xmath10 dimensions at the time @xmath262 . note the double - logarithmic scale . a scaling with @xmath302 corresponds to a power - law exponent @xmath303 for the kinetic energy in @xmath10 . see text for more details on the scaling exponents . inset : phase angle @xmath304 as in fig . [ fig:2dphase512].,scaledwidth=45.0% ]    , as defined in eqs . ( [ eq : ev])([eq : decomposition ] ) ( see caption of fig .  [ fig : energyspectrum ] for more information ) . parameters are the same as in fig . [ fig : modeoccupationsequence3 ] ( lower left panel ) , for the run in @xmath91 dimensions at the time @xmath269 . note the double - logarithmic scale . a scaling with @xmath305 corresponds to a power - law exponent @xmath303 for the kinetic energy in @xmath91 . inset : the black spots mark points around the vortex line cores where the density falls below 5% of the mean density @xmath107.,scaledwidth=45.0% ]    to exhibit vortical flow and define the decomposition we use the polar representation @xmath306 of the field in terms of the density @xmath307 and a phase angle @xmath308 . this allows to express the particle current @xmath309 in terms of the velocity field @xmath310 .    with this , we decompose the kinetic - energy spectrum following ref . @xcite , splitting the total kinetic energy @xmath311 as @xmath312 into a ` classical ' part @xmath313 and a ` quantum - pressure ' component @xmath314 the radial energy spectra for these fractions involve the fourier transform of the generalised velocities @xmath315 and @xmath316 , @xmath317 following ref . @xcite , the velocity @xmath318 , which due to the multiplication of @xmath122 with the density @xmath107 becomes regularized and can be fourier transformed , is furthermore decomposed into ` incompressible ' ( divergence free ) and ` compressible ' ( solenoidal ) parts , @xmath319 , with @xmath320 , @xmath321 , to distinguish vortical superfluid and rotationless motion of the fluid . for comparison of the kinetic - energy spectrum with the single - particle spectra @xmath119 , we determine occupation numbers corresponding to the different energy fractions as @xmath322 the resulting spectra @xmath323 , @xmath324 , and @xmath325 add up to @xmath326 , which agrees with the single - particle spectrum up to small corrections , see appendix [ app : hydro ] .    in figs . [ fig : energyspectrum ] and [ fig : modeoccupationboltzmann3 ] , we depict the momentum distributions of the occupation numbers @xmath323 , @xmath324 , and @xmath325 , together with the previously shown total single - particle spectrum @xmath119 , each at a late time when the system is close to the fixed point , i.e. , shows the predicted scaling both in the ir and the uv . red circles denote @xmath327 , filled blue squares @xmath328 , and open grey squares @xmath329 . qualitatively , the results are similar for @xmath10 and @xmath91 .    in the range of large wavenumbers , the spectrum is dominated by the compressible and quantum - pressure components . the scaling of these excitations exhibits the weak - wave - turbulence exponent @xmath255 . for smaller wave numbers the scaling changes to @xmath330 . the decomposition into the various components now shows that this switching to a different scaling in the ir is clearly due to the take - over of a different character of the excitations with a modified flow pattern accounted for by the incompressible ( solenoidal ) component of @xmath318 . the fact that in this regime contributions from vortical flow @xmath327 dominate is in accordance with our interpretation of the strong ir scaling by a model of an ensemble of vortical excitations . as we have discussed in the previous section the analytically predicted infrared power laws @xmath330 are consistent with a finite density of independent vortices and antivortices ( @xmath10 ) or vortex lines ( @xmath91 ) . moreover , we find that the compressible component in the ir represents the second strongest contribution to the flow and develops a nonthermal scaling as @xmath331 . we will discuss this result in more detail in sect . [ subsec : acousticturbulence ] below .    for momenta larger than about @xmath332 for @xmath10 and @xmath333 for @xmath91 the compressible and quantum - pressure components dominate the momentum distributions . in the regime of intermediate momenta , above the scale @xmath334 with @xmath335 the mean inter - vortex spacing , where the incompressible flow and the rest are of roughly equal strength one finds a scaling of approximately @xmath336 , corresponding to @xmath337 for both the two- and the three - dimensional case . scaling with this exponent , which is the same as in classical kolmogorov turbulence in an incompressible fluid , has been found in simulations of superfluid turbulence in @xmath338 before @xcite . in some of these cases , turbulent flow appeared in the time evolution of a system starting from a taylor - green configuration of vortex line tangles . we emphasize that our configuration after the creation of vortical excitations resembles more a state of the kind usually termed chaotic turbulence . it is unclear whether chaotic turbulence can be related to classical kolmogorov turbulence as this does not bear near - classical bundles of equal - orientation vortex lines . nonetheless we observe such @xmath303 scaling in the range of intermediate momenta where the interplay between the incompressible and compressible components is most pronounced .      ) ([eq : decomposition ] ) are shown as functions of the radial momentum @xmath37 , from an average over @xmath339 single runs in @xmath10 dimensions ( upper panel ) and @xmath91 ( lower panel ) , at times @xmath340 ( @xmath10 ) and @xmath341 ( @xmath91 ) , triggered on the decay of the last vortical excitation . see caption of fig . [ fig : energyspectrum ] for more information . parameters are : @xmath10 : @xmath245 , @xmath342 , @xmath247 ; @xmath91 : @xmath249 , @xmath343 , @xmath344 . note the double - logarithmic scale . the figure shows that , shortly after the last vortex ring has disappeared ( incompressible component breaks down ) , compressible excitations exhibiting acoustic turbulence scaling @xmath345 remain present.,title="fig:",scaledwidth=45.0% ]   + ) ([eq : decomposition ] ) are shown as functions of the radial momentum @xmath37 , from an average over @xmath339 single runs in @xmath10 dimensions ( upper panel ) and @xmath91 ( lower panel ) , at times @xmath340 ( @xmath10 ) and @xmath341 ( @xmath91 ) , triggered on the decay of the last vortical excitation . see caption of fig . [ fig : energyspectrum ] for more information . parameters are : @xmath10 : @xmath245 , @xmath342 , @xmath247 ; @xmath91 : @xmath249 , @xmath343 , @xmath344 . note the double - logarithmic scale . the figure shows that , shortly after the last vortex ring has disappeared ( incompressible component breaks down ) , compressible excitations exhibiting acoustic turbulence scaling @xmath345 remain present.,title="fig:",scaledwidth=45.0% ]    the ir scaling @xmath331 of the compressible component ( blue filled squares ) in figs . [ fig : energyspectrum ] , [ fig : modeoccupationboltzmann3 ] , and also fig . [ fig : pairingspectrum ] below suggests an interpretation in terms of acoustic turbulence @xcite and corroborates the numerical findings reported in ref . the scaling is persistent until late times , see , e.g. , fig . [ fig : pairingspectrum ] , but decays after the vortices have disappeared . to check that the nonthermal scaling of the incompressible component is not an artefact of the decomposition of the vortical flow into compressible and incompressible parts , we show the momentum spectrum for an average over selected runs in @xmath91 , see fig . [ fig : acoustic ] ( lower panel ) , where the snapshot is taken shortly after the last vortex ring has disappeared . one can observe that acoustic turbulence can survive for a limited period , @xmath346 , beyond the time when vortical excitations are negligible . in the case of @xmath10 , the effect is weaker but still present , see fig . [ fig : acoustic ] ( upper panel ) , where we show the averaged spectrum shortly after the last vortex - antivortex annihilation . our interpretation of this finding is that acoustic turbulence coexists with the dominant vortical flow . it is driven by vortex motion and especially vortex annihilation processes , which are known to produce compressible excitations . in @xmath10 , signatures of acoustic turbulence are less pronounced , which we attribute to reduced vortex dynamics or , equivalently , weaker driving of compressible excitations . , for the run in @xmath91 dimensions , at the time @xmath269 . note the double - logarithmic scale . the black line is a gaussian fit to the data which shows that in the low - velocity regime , the distribution is dominated by gaussian fluctuations . a high - velocity scaling with @xmath347 reflects the presence of uncorrelated vortices , cf . ( [ eq : pvelx]).,scaledwidth=45.0% ]     dimensions at @xmath348 . at this time an average of @xmath349 vortices are present in the system . parameters : @xmath350 @xmath351 , @xmath352 . note the double - logarithmic scale . a high - velocity scaling with @xmath347 reflects the presence of uncorrelated vortices , cf . ( [ eq : pvelx ] ) . gaussian fluctuations are suppressed in the distribution of velocities of the individual vortices , measured through the motion tracking of the vortex cores.,scaledwidth=45.0% ]    we furthermore investigated the velocity - field probability distribution in the turbulent regime to compare our data with the expectations summarized in sect . [ subsec : spectrumvelo ] . in @xmath91 dimensions , see fig . [ fig : velocitydistribution3 ] , the results corroborate theoretical and experimental results reported in refs .  @xcite . the low - momentum distribution is characterized by gaussian statistics , whereas the uv regime clearly shows scaling @xmath353 predicted by the model of randomly distributed , uncorrelated vortex lines , see eq . ( [ eq : pvelx ] ) . note that even at late times , when vortex rings shrink in size , no signs of vortex pairing , @xmath354 , can be seen in the velocity distribution . this is due to dominating gaussian fluctuations in the low momentum region . in @xmath10 dimensions , see fig . [ fig : velocitydistribution2 ] , we compare the velocity - field distribution with the velocity distribution of the positions of the individual vortices . similar to the three - dimensional case the ir part of the velocity - field distribution is characterized by gaussian statistics , whereas the uv regime shows single - vortex scaling  @xcite . the vortex velocity distribution is obtained from the statistical analysis of a vortex - tracking algorithm which is designed to detect regions of low density accompanied by a winding number equal to one . the method is analogous to the experimental set - up of ref . @xcite , where particles are trapped inside vortex cores to study vortex velocities in @xmath91 . in agreement with the experiment , fig . [ fig : velocitydistribution2 ] shows that gaussian random motion is suppressed for vortices , but @xmath355 scaling clearly persists . in this way , quantum turbulence can be distinguished from classical turbulence , where a continuous distribution of vorticity favors a gaussian velocity distribution .      . parameters are the same as in fig . [ fig : modeoccupationsequence3 ] , for the run in @xmath91 dimensions , at the time @xmath276 ( lower right panel ) . note the double - logarithmic scale . the scaling @xmath356 in the far ir reflects pair correlations present in far - separated small vortex loops . ( see caption of fig . [ fig : energyspectrum ] for more information).,scaledwidth=45.0% ]     in @xmath10 dimensions , at the time @xmath262 : occupation number obtained from averaging over all runs ( black dots ) ; occupation number obtained from averaging over selected runs featuring tightly bound vortex antivortex pairs with maximum pair correlation length @xmath357 ( red circles ) . parameters : @xmath358 , @xmath359 , @xmath251 . note the double - logarithmic scale . a scaling with @xmath166 reflects the presence of correlated vortex - antivortex pairs in @xmath10.,scaledwidth=45.0% ]    in fig . [ fig : pairingspectrum ] we show the decomposition of the single - particle spectrum into the previously defined incompressible , compressible , and quantum - pressure components , see eq . ( [ eq : decomposition ] ) , for the 3-dimensional system at @xmath360 , as in the lower right panel of fig . [ fig : modeoccupationsequence3 ] . at this time , only a small density of vortex rings has survived , which decays under the influence of the noise field , cf . also refs . @xcite . during the period when only a few small vortex rings remain one can observe a decrease of the infrared scaling exponents of the occupation number and its incompressible part to the value @xmath275 . since vortical excitations are still dominating the spectrum , we can interpret this observation in terms of the statistical point vortex model introduced in sect . [ sec : randomvortexmodel ] , as vortex - antivortex pair correlations , see sect . [ subsec : vortexloops ] . this is consistent with the snapshots of individual runs showing small vortex ellipses ( see inset of fig .  [ fig : pairingspectrum ] ) . the scaling transition can be identified with the scale @xmath361 , with estimated minor radius @xmath362 . [ fig : pairingspectrum2 ] shows the average occupation number spectrum at @xmath262 in @xmath10 dimensions , as in the lower right panel of fig . [ fig : modeoccupationsequence2 ] . in the average over all generated runs pairing effects can hardly be observed . however , for selected runs with small maximum pair correlation length @xmath357 , as introduced in sect . [ subsec : correlatedvortices ] , the pair scaling appears . we find that in the selected runs ( @xmath363 of total number of runs ) , at the chosen point of time one or at most a few vortex - antivortex pairs with pairing length smaller than the distance between pairs are present . this constitutes the final period of the evolution , shortly before the last pair has disappeared through mutual annihilation and the system fully thermalizes . the generic configuration in @xmath10 dimensions during the preceding time interval of critical slowing down close to the nonthermal fixed point is characterized by randomly positioned vortices and antivortices , with pairing correlations nevertheless present . pairing is seen in the density - density correlation function between vortices and antivortices shown in fig . [ fig : paircorr2d ] , from which we read off a largest scale @xmath364 , cf .  sect . [ subsec : correlatedvortices ] . this implies pair scaling below @xmath365 , not visible in the few ( black ) full circles in fig . [ fig : pairingspectrum2 ] .    . parameters as in fig . [ fig : modeoccupationsequence2 ] , at the time @xmath366 . red crosses : pair correlations of vortices and antivortices within the same circulation . blue stars : pair correlations between vortices and antivortices . the rise above @xmath367 indicates an effective `` binding '' . black crosses : resulting correlation function @xmath368 introduced in eq . ( [ eq : cf ] ) . , scaledwidth=45.0% ]    we observe that , in @xmath10 , scattering between the vortices of equal and opposite circulation can lead to the increase of the pairing length between correlated vortices and antivortices . this is in contrast to the situation for @xmath91 where the scattering between vortex rings rarely leads to an increase in the size of the rings . it is energetically advantageous for the rings to shrink in size , and thus , at a particular late time of the evolution the presence of separated small rings is generic and  pairing  is seen in the spectrum .    in summary , our analysis gives a picture of the necessary conditions and the character of the nonthermal fixed point . it also shows that pairing effects are a first signature for the system leaving this fixed point again for the final move to thermal equilibrium . we have studied in detail superfluid turbulence in two- and three - dimensional dilute bose gases by means of simulations in the classical - wave limit of the underlying quantum field theory . a focus is set on the identification and characterization of stationary scaling solutions , in particular for single - particle momentum distributions . characteristic exponents @xmath5 of @xmath119 corroborate the analytical predictions made in ref . our findings suggest that local field expectation values and short- to intermediate - range coherence , including topological excitations , are at the basis of the infrared power laws predicted within a full nonperturbative dynamical field theory @xcite employing two - particle - irreducible effective action techniques . we have shown that the stationary scaling is maintained by the presence of particle ( ir ) and energy ( uv ) fluxes , originating at intermediate scales and directed towards the low- and high - frequency limits , respectively . the respective fluxes were found to be consistent with the particular scaling exponents of the momentum distribution in the respective regimes . our main result is the identification of the nonthermal fixed - point with the appearance of topological excitations in the system . we could successfully employ statistical models of point vortices in two dimensions and vortex rings in three dimensions to interpret the ir scaling exponent of the spectra . moreover , we observed in our simulations , during the final thermalization stage , decay of the turbulence scaling into a scaling derived from the assumption of vortex - antivortex pairing . to characterize further the bimodal power spectra developed close to the nonthermal fixed point we have analysed the decomposition of the overall flow pattern into compressible and incompressible components . we found that the ir power spectra of underlying compressible excitations suggest an understanding in terms of acoustic turbulence on top of the vorticity - bearing quasicondensate . finally , to compare with the results of recent experiments and analytical predictions , the velocity - field probability distribution as well as the velocity statistics of individual vortices has been studied . this observable is of great interest , as it can be used to experimentally distinguish between classical and quantum turbulence . the connection of qt phenomena with ab - initio dynamical field theoretic methods points a way to unified analytical studies of turbulence . moreover , it provides hints of how the proposed nonthermal fixed points in relativistic systems @xcite are realized in nature . a manifestation of the approach of a nonthermal fixed - point in terms of quasi - solitary pattern formation and charge separation in a relativistic scalar model has recently been reported in ref . @xcite . experimental studies of universal phenomena in nonequilibrium dynamics of ultracold atoms have great potential since universal effects do not depend significantly on initial conditions and details of the system . following this idea , our numerical protocol was chosen such that an experimental verification of our findings is within reach of present - day cold - atom experiments . the study of turbulence in ultracold gases may have great impact on many other fields of physics . prominent examples are strongly correlated nuclear matter produced in heavy - ion collisions and early - universe cosmology . t.  g. and d.  s. thank j. berges and c. scheppach for collaboration on related work . the authors thank j.  berges , n.  berloff , e.  bodenschatz , m.  j.  davis , g.  falkovich , s .- c . gou , h.  horner , r.  kerr , g.  krstulovic , m.  k.  oberthaler , j.  m. pawlowski , b.  shivamoggi , b.  svistunov , and m.  tsubota for useful discussions . they acknowledge support by the deutsche forschungsgemeinschaft ( ga 677/7,8 ) , by the university of heidelberg ( frontier , excellence initiative , center for quantum dynamics ) , by the helmholtz association ( ha216/emmi ) , and by the university of leipzig ( grawp - cluster ) . the authors furthermore thank kitp , santa barbara , for its hospitality . this research was supported in part by the national science foundation under grant no . phy05 - 51164 . [ appendix ] a singly quantized vortex with circulation @xmath369 in two dimensions is described , in polar coordinates , by the solution @xmath370 of the gross - pitaevskii equation ( [ eq : gpe ] ) , where @xmath371 is real and approaches the square root of the bulk density @xmath372 for large distances @xmath104 from the vortex core where @xmath373 . a vortex is a stationary solution of eq . ( [ eq : gpe ] ) , evolving as @xmath374 with @xmath375 . the velocity field of the vortex is @xmath376 the fourier transform of this field can be conveniently evaluated by the help of the first - order bessel function , @xmath377 , as @xmath378      based on eq . ( [ eq : vorticity ] ) , we calculate the vorticity of two straight vortex lines in three dimensions , with circulation @xmath379 , @xmath380 , placed at positions @xmath381 and @xmath382 in the @xmath184-@xmath185 plane and relate it to the momentum spectrum , see eq . ( [ eq : spectrumfromvorticity ] ) . the lines are parametrized by @xmath383 with @xmath384 $ ] . in fourier space , the vorticity density , defined in eq . ( [ eq : vorticity ] ) , reads @xmath385 the vorticity spectrum follows from eq . ( [ eq : vorticity ] ) , @xmath386 ) @xmath387 ) one obtains @xmath388 we regularize the delta distribution , giving it a finite width @xmath389 . the angle average of eq . ( [ eq : lines1 ] ) then reads @xmath390,\end{aligned}\ ] ] where @xmath391 is the length of the system in @xmath65-direction and @xmath392 . one obtains the scaling behaviour @xmath393 which is the scaling of a single line . with eq . ( [ eq : spectrumfromvorticity ] ) , we conclude @xmath88 . following the same reasoning , the angle average of eq . ( [ eq : lines2 ] ) is given as @xmath394\,.\end{aligned}\ ] ] for @xmath395 @xmath396 which gives @xmath248 . for @xmath397 @xmath398 which gives @xmath88 . these scalings are confirmed in fig . [ fig : straightvortexline ] .     on a @xmath399 lattice . the distance scale @xmath400 is indicated . note the double - logarithmic scale.,scaledwidth=45.0% ]     on a @xmath399 lattice . the radius scale @xmath401 is indicated . note the double - logarithmic scale.,scaledwidth=45.0% ]      in sects . [ subsec : indpairs ] and [ subsec : correlatedvortices ] , we have discussed the scaling of the momentum distribution for ensembles of independent and pair correlated vortices in two dimensions , respectively . in this appendix we generalize the scaling derived above for straight vortex lines to the case of many such lines in @xmath91 . we derive an expression for the momentum spectrum of a system of @xmath108 straight vortex lines oriented along the @xmath65 direction . the average vorticity squared is given by eq . ( [ eq : energyspeccorr3 ] ) where @xmath402 are vectors in the two - dimensional @xmath184-@xmath185 plane pointing to the @xmath130-th vortex line with circulation @xmath403 . the vorticity of the single line is @xmath404 hence , eq .  ( [ eq : energyspeccorr3 ] ) has the same form as eq . ( [ eq : vortexdensity ] ) for @xmath10 , @xmath405 up to a @xmath406 term , which arises from the infinite extent of the vortex lines in @xmath65-direction . as in app . [ app : vortexline ] , we regularize the delta distributions , and averaging ( [ eq : energyspeccorr3line ] ) over solid angles yields @xmath407 which is analogous to eq . ( [ eq : ve ] ) in two dimensions . as the position of the @xmath130-th vortex line is determined by @xmath402 , statistical averaging of ( [ eq : energyspeccorr3linefinal ] ) can be done in the same way as in @xmath10 . a vortex ring with radius @xmath104 lying in the @xmath184-@xmath185 plane can be parametrized as @xmath408 with @xmath409 $ ] . in fourier space , similar steps as above lead to the angle - averaged vorticity density @xmath410 for @xmath411 and @xmath412 for @xmath413 , corresponding to @xmath224 and @xmath88 . this is shown in fig . [ fig : vortexring ] . both scalings are also found in the momentum spectrum of a vortex ellipse , discussed in sect . [ subsec : vortexloops ] . in this paper , a dilute superfluid gas of bosons of mass @xmath414 is considered , with contact interactions quantified by a coupling constant @xmath19 . its dynamics is described , in the classical wave limit , by the gross - pitaevskii equation ( gpe ) ( [ eq : gpe ] ) . we use units where @xmath16 . we consider the gas to be contained in a box of size @xmath415 , @xmath416 , with periodic boundary conditions . in @xmath91 , the coupling constant @xmath19 is related to the @xmath4-wave scattering length @xmath15 by @xmath417 . in @xmath10 , one has @xmath418^{-1}$ ] where @xmath419 is the chemical potential and @xmath420 is a scattering length in two dimensions which , for a gas of hard - spheres of radius @xmath15 is given by @xmath421 , with the euler - mascheroni constant @xmath422 @xcite . for a two - dimensional gas created by trapping a three - dimensional one tightly in one dimension , with harmonic - oscillator length @xmath423 , the effective 2d scattering length is given by @xmath424 , where @xmath425 @xcite . diluteness implies that @xmath426 , the interparticle spacing @xmath427 being determined by the mean density @xmath428 .    the initial values for the real and imaginary parts of the field @xmath429 are randomly chosen from a gaussian distribution with width @xmath430 , centered around @xmath431 , where @xmath432 is the occupation number at time @xmath25 and @xmath433 is a random phase angle . correlation functions including @xmath0 are obtained by averaging over many trajectories . to induce transport from small to large wave numbers , only a few modes near @xmath434 are chosen to be macroscopically occupied at the initial time , @xmath435 . such an initial state can be prepared , e.g. , by bragg scattering of photons from a bose - einstein condensate . our numerical simulations are performed on space - time lattices with side lengths @xmath436 , with lattice spacing @xmath437 and in total @xmath438 gridpoints . we use periodic boundary conditions and @xmath439 for @xmath10 and @xmath440 for @xmath91 . ( [ eq : gpe ] ) is written in terms of the dimensionless variables @xmath441 , @xmath442 and @xmath443 . all quantities in the figures are either dimensionless or shown in lattice units , with the length unit given by @xmath444 . the dimensionless lattice momenta are @xmath445^{1/2 } , \mathbf{k}=2\pi\mathbf{n}/n_s , \ , ,n_d),\ ,   n_i= -n_s/2, ... ,n_s/2 $ ] .    in order to relate our simulations to a typical situation in experiment , we give parameters for rb-87 . we estimate the total energy of the gas in equilibrium to be given by the interaction energy which , by equipartition is related to the temperature , @xmath446 . at the transition to degeneracy the thermal de broglie wave length @xmath447 is of the order of the interparticle spacing and hence @xmath448 . inserting this into the expression for the interaction energy allows to express the lattice spacing in terms of the density and the scattering length as @xmath449 inserting parameters for a typical rb-87 experiment in @xmath91 , viz . , @xmath450 nm , @xmath451 kg , @xmath452m@xmath453 we obtain @xmath454@xmath419 m . hence , time scale as @xmath455s . our lattice time typically runs until @xmath456 , which therefore represents a realistic observation timescale in experiments . typical parameter choices in our simulations are @xmath343 , @xmath344 , @xmath249 for @xmath91 and @xmath342 , @xmath247 , @xmath245 for @xmath10 . as a function of the radial momentum @xmath37 at different times , for @xmath10 ( left panel ) and @xmath91 ( right panel ) . parameters are the same as in fig . [ fig : modeoccupationsequence2 ] , for the run in @xmath10 dimensions and in fig . [ fig : modeoccupationsequence3 ] , for @xmath91 . note the logarithmic scale . , scaledwidth=48.0% ]    in sect . [ subsec : superfluidturbulence ] , we have defined occupation numbers corresponding to the incompressible , compressible and quantum pressure components of the kinetic energy . numerically , the resulting spectra @xmath323 , @xmath324 , and @xmath325 add up to @xmath458 which is equal to the full single - particle spectrum @xmath119 up to small corrections . in particular , the scaling behavior emerges as the same for both , @xmath458 and @xmath119 . to make this more explicit , we calculate @xmath459 in terms of generalized velocities @xmath460v , q@xmath461 , @xmath462 @xmath463 denotes the real part and @xmath464 the fourier transform . the middle term vanishes since the expectation value is purely imaginary due to isotropy . it follows that @xmath465 since phase and phase velocity are expected not to be correlated , the main contribution to the 4-point correlation functions is assumed to arise from terms of the form @xmath466 , @xmath467v , q@xmath461 . in the superfluid regime , the phase of the field is a slowly varying function in position space . therefore , @xmath468-terms are strongly peaked at zero momentum , effectively acting as regularized delta functions under the convolution . hence , @xmath469 this approximative identity is supported by fig . [ fig : differencespectrum2d3d ] .
Abstract:

Answer: nonthermal fixed points of the dynamics of a dilute degenerate bose gas far from thermal equilibrium are analysed in two and three spatial dimensions . universal power - law distributions , previously found within a nonperturbative quantum - field theoretical approach and recently shown to be related to vortical dynamics and superfluid turbulence [ phys . rev . b * 84 * , 020506(r ) ( 2011 ) ] , are studied in detail . the results imply an interpretation of the scaling behavior in terms of independent vortex excitations of the superfluid and show that the statistics of topological excitations can be described in the framework of wave turbulence . the particular scaling exponents observed in the single - particle momentum distributions are found to be consistent with irreversibility as well as conservation laws obeyed by the wave interactions . moreover , long - wavelength acoustic excitations of the vortex - bearing condensate , driven by vortex annihilations , are found to follow a nonthermal power law . considering vortex correlations in a statistical model , the long - time departure from the nonthermal fixed point is related to vortex - antivortex pairing . the studied nonthermal fixed points are accessible in cold - gas experiments . the results shed light on fundamental aspects of superfluid turbulence and have strong potential implications for related phenomena , e.g. , in early - universe inflation or quark - gluon plasma dynamics .