Fundamental Connection Between Temperature-Quenched 2D Superfluids and 2D Quantum Turbulence

There is a fundamental connection between temperature-quenched 2D superfluids and 2D quantum turbulence: the mechanism responsible for the decay of the vorticity after the quench is the enstrophy cascade of 2D turbulence. The range of the cascade is shown to increase with time after the quench as t1/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t^{1/2}$$\end{document}, being equal to the dynamic scaling length characterizing the quench. These results may apply to other phase-ordering transitions involving the decay of topological objects.


Introduction
The static equilibrium properties of the two-dimensional (2D) superfluid phase transition are now well understood in terms of the Kosterlitz-Thouless theory [1] involving thermally excited vortex-antivortex pairs. However, the nonequilibrium dynamics of such a system is less well characterized. An example of this is a sudden temperature quench of the superfluid that leaves a nonequilibrium high density of vortices at the new low temperature, and the only way they can decay is by vortexantivortex annihilation, but this takes time due to the friction on the normal vortex cores. Early studies of the vortex dynamics by Minnhagen and co-workers [2] in 2D XY model simulations found puzzling results, that there seemed to be be two different dynamic exponents, depending on the boundary conditions used in the simulations. One dynamic exponent found only with fluctuating-twist boundary conditions took the form z scale = 2 K − 2 , where K = ℏ 2 s (T)∕(m 2 k B T) with s the macroscopic 2D superfluid areal density and m the atomic mass. This exponent varies with the temperature T, taking a value z scale = 2 at the critical temperature T KT where K = 2∕ , and then increasing at temperatures below T KT . A second dynamic exponent found with all boundary conditions was the constant value z ≈ 2 . The presence of two different dynamic exponents was unusual, and it was not entirely clear how these came into the vortex dynamics. An exact solution [3] for temperature quenches starting from T KT and below shed light on some of these questions, where the vortex density was found to decay with time t after the quench as t −z scale ∕z . Here, z scale is evaluated at the initial temperature before the quench, and z = 2 is exact.
Another nonequilibrium situation in 2D superfluids is the case of quantum turbulence, which is the case when extra vortices are injected into the superfluid. By analogy with turbulence results in 2D classical fluids [4,5], it was expected that the turbulence in superfluids would self-organize into constant-flux cascades of both energy and vorticity (enstrophy). This was indeed observed in computer simulations, where a cascade with a constant flux of energy from small to large scales was found with a Kolmogorov k −5∕3 energy spectrum [6] . A constant-flux enstrophy cascade from large scales to small was also observed in a different simulation, with the expected k −3 energy spectrum [7], as had been predicted much earlier [8]. Analytic solutions have been found for the enstrophy cascade [9], which is a forward cascade where vortex pairs of large separation are injected into the superfluid at a constant rate, and due to the friction on the vortex cores they diffuse at a constant flux to smaller separations, finally annihilating at the vortex core scale at the same rate they are being injected.
We have found that these seemingly unrelated nonequilibrium states have a fundamental connection: the mechanism of the decay of the vorticity after the temperature quench is the enstrophy cascade. There is a constant flux of vortices from the large separation scales present at the initial temperature, moving to the small scale of the vortex core size where they annihilate, to finally get to the very low vortex density of the final temperature. This is precisely the definition of the enstrophy cascade, and the range of the cascade is found to grow as the dynamic length, varying as t 1∕z .

Temperature Quenches
To characterize an instantaneous temperature quench from an initial temperature T i to a final low temperature T f , it is necessary to solve a Fokker-Planck equation [10] for the distribution Γ(r, t) of vortex-antivortex pairs of separation r at time t, where this is made dimensionless by taking r in units of the vortex core size a 0 , Γ in units a −4 0 , and t in units of the diffusion time 0 = a 2 0 ∕2D with D the vortex diffusion coefficient at the temperature T f set by the thermal bath (which e.g. would be the substrate of a thin superfluid helium film). The vortex motion is diffusive even at low temperatures, as shown by experiments on helium films below 1 K [11][12][13].
1 is determined by simultaneously solving the Kosterlitz recursion relation [14] for the superfluid density For quenches starting from below T KT , we can make the approximation that K i is effectively a constant, since the spatial renormalization from Eq. (2) is rapid, changing from the initial K 0i to the renormalized K i over a length scale that can be less than a core radius. In that limit the initial distribution at t = 0 is given by Γ � 0 r −2 K i with Γ � 0 a constant depending on the vortex core energy. The solution for the pair distribution after the quench where now K = K f ≃ const is given by where 1 F 1 is the confluent hypergeometric function of the first kind, z = 2 exactly, and is a constant depending on K i and K f . Figure 1 shows a plot of Γ(r, t) for a quench from 0.9 T KT to 0.1 T KT , where K i = 0.814 and K f = 7.479 . Just after the quench only the smallest-separation vortices start to disappear, since they get to the core-size annihilation scale r = 1 the quickest, while the distribution at larger scales remains at the equilibrium variation r −2 K i = r −5.11 . Note in particular at later times that the distribution develops a very flat region at shorter separations, and that the length scale of the flat regions increases with increasing time.
The density of vortices can then be found by integrating over r, where ′ is a constant, and z scale = 2 K i − 2 . At times longer than a few diffusion times the hypergeometric function goes to a constant, and the vortex decay is then t −z scale ∕z = t −1.557 . The exponent z scale appears in the dynamics due to the second term on the right-hand side of Eq. (1), which reflects the attractive interaction between the vortices of the pair, while the factor of z = 2 appears because of the second derivative of the first term in Eq. 1, which gives rise to the diffusive motion of the vortices. This confirms the result of Minnhagen et al. [2] that two dynamic exponents control the dynamics of the superfluid, but it is still not entirely clear why the boundary conditions play a role in this. The factor 2 K i appears in z scale because the initial distribution function Γ i is proportional to r −2 K i , and this gives a very simple interpretation of why the vortex decay becomes more rapid at lower temperatures where K i increases: the steeper distribution means most of the pairs are closer together, and hence annihilate more rapidly. This has nothing to do with any Kibble-Zurek "creation" of vortices in the quench [15,16].

Enstrophy Cascade
Analytic solutions for the enstrophy cascade can be found [9] using the Fokker-Planck Eq. (1), but with a 2D Dirac delta function added to the right-hand side to inject vortex pairs of large separation R. The injection rate multiplying the delta function in dimensionless form is = a 2 0Q 0 where Q is the number of vortex pairs of separation R injected per unit area per time. The thermal bath is held at a very low temperature (i.e. 0.1 T KT ) so that no thermal vorticity is present. Due to the contact with the thermal bath, negative-temperature Onsager states of the vortices [17] cannot form, and hence the energy cascade will not be present. For the comparison with quenched superfluids, we only consider the case where is small enough that the resulting vortex density is well below that at the Kosteritz-Thouless transition. In that limit the superfluid density is unaffected by the vortices, so K = K 0 = const , the "bare" superfluid density, and the solution for the distribution function characterizing the enstrophy cascade is For r > R , the solution is a quasi-thermal distribution extending from R, which arises from injected pairs initially at separation R getting a thermal kick to higher separation. Note for comparison with the quench case that the pair distribution of the cascade is completely constant for the cascade range r < R . This is a solution with a constant flux of vortices from the large scale R to the small core size scale, since the flux is r Γ∕ r + 2 KΓ = = const. The energy spectrum corresponding to Eq. (5) varies as k −3 [9], similar to the classical enstrophy spectrum [4,5], though it is not quite the same. It is proportional to 0 where is the k-space enstrophy injection rate and 0 the diffusion time. The classical spectrum is proportional to 2∕3 since the vorticity there has all different diffusion characteristics, whereas the quantum vortices all have the same 0 .
For further comparison with the quench case, we can compute the flux from the quench distributions of Fig. 1, shown in Fig. 2. The flux is seen to be nearly constant over separations that increase with time, before falling off at large scales. We can arbitrarily define a dynamic length as the length where the flux has decreased by 5% from its value at r = 1 , marked by the arrows in Fig. 2. Figure 3 shows the dynamic lengths as a function of time, and is completely well fit by (t) = 0 t 1∕2 with 0 = 1.38 , the known dynamic length associated with z = 2 . This is strong proof that the mechanism of the vortex decay after the quench is an enstrophy cascade, with the range of the constant-flux cascade given by the increasing dynamic length.

Conclusion
We think it is likely that the connection between turbulent cascades and phaseordering dynamics may apply more generally. There are numerous systems with phase-ordering transitions that involve topological excitations [18], and it would be interesting to look for constant-flux cascades in these systems. One wellknown example is quenches of liquid crystals from the isotropic to nematic state [19,20], where disclinations in the director field form annihilating topological objects.
A further note is that the 2D superfluid system considered here is actually not in the same dynamic universality class as the 2D XY model, even though they are in the same 2D O(2) static class. Dissipation in a superfluid vortex is entirely localized to the normal vortex core, while dissipation in XY vortices occurs over finite distances from the core [21]. This leads to a scale-dependent mobility of the XY vortices [22] and gives rise to a log correction term in the dynamic scaling [23].