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 ] .