Anomalous phase ordering of a quenched ferromagnetic superfluid
L. A. Williamson P. B. Blakie






Anomalous phase ordering of a quenched ferromagnetic superfluid








L. A. Williamson1,2*,
P. B. Blakie1







1 Dodd-Walls Centre for Photonic and Quantum Technologies, Department of Physics, University of Otago, Dunedin 9016, New Zealand


2 Department of Physics, Lancaster University, Lancaster, LA1 4YB, United Kingdom



* l.a.williamson@lancaster.ac.uk





December 30, 2023







§ ABSTRACT


Coarsening dynamics, the canonical theory of phase ordering following a quench across a symmetry breaking phase transition, is thought to be driven by the annihilation of topological defects. Here we show that this understanding is incomplete. We simulate the dynamics of an isolated spin-1 condensate quenched into the easy-plane ferromagnetic phase and find that the mutual annihilation of spin vortices does not take the system to the equilibrium state. A nonequilibrium background of long wavelength spin waves remain at the Berezinskii-Kosterlitz-Thouless temperature, an order of magnitude hotter than the equilibrium temperature. The coarsening continues through a second much slower scale invariant process with a length scale that grows with time as t^1/3. This second regime of coarsening is associated with spin wave energy transport from low to high wavevectors, bringing about the the eventual equilibrium state. Because the relevant spin waves are noninteracting, the transport occurs through a dynamic coupling to other degrees of freedom of the system. The transport displays features of a spin wave energy cascade, providing a potential profitable connection with the emerging field of spin wave turbulence. Strongly coupling the system to a reservoir destroys the second regime of coarsening, allowing the system to thermalise following the annihilation of vortices.








1pt
fancy
1pt






§ INTRODUCTION



Quenching a system across a continuous phase transition from a high to low symmetry phase causes the system to spontaneously break symmetry. Immediately after the quench causally disconnected regions of the system will break symmetry independently, resulting in the formation of domains with independent order parameter orientation. The subsequent growth of these domains toward the global equilibrium state is known as coarsening dynamics. Although the microscopic details of coarsening are usually extremely complicated, at a macroscopic level a much simpler scaling regime can emerge for large average domain size L. Spatial correlations of the order parameter at different times t then collapse onto a single curve when rescaled by L, and the domains grow as L∼ t^1/η with the scaling exponent η determined by the dynamic universality class <cit.>. Such universal dynamics has been explored in a vast variety of systems, ranging from the early universe <cit.> to superfluid formation <cit.> to opinion spreading in sociology <cit.>. When the quench produces topological defects, the decay of these defects has long been thought to provide a unifying framework for understanding the coarsening <cit.>.


Recently, there has been much interest in coarsening dynamics in ultracold atom systems, which are well isolated from their environment and present a pristine system for studying nonequilibrium phase transitions <cit.>. Of particular interest are multicomponent condensates, which support a rich variety of order parameter manifolds and associated topological defects <cit.>. Theoretical studies of coarsening in a variety of cold atom systems <cit.> have culminated in the recent experimental observation of universal dynamics in a quenched quasi-1D scalar Bose gas <cit.> and in a quenched quasi-1D spin-1 condensate <cit.>. Simulations of a homogeneous quasi-2D spin-1 condensate quenched from the polar phase to the easy-plane ferromagnetic phase, see Fig. <ref>(a), identified coarsening dynamics driven by the mutual annihilation of transverse spin vortices with domain size growing as L∼ t/log t <cit.>. A log correction to scaling is familiar from two dimensional systems supporting vortices <cit.>.

In this work we study the easy-plane ferromagnetic ordering of a homogeneous quasi-2D spin-1 condensate after all vortices have annihilated. Remarkably, we find that the annihilation of vortices does not take the system to the equilibrium state. Instead, a nonequilibrium background of spin waves remain at the Berezinskii-Kosterlitz-Thouless (BKT) temperature, an order of magnitude hotter than the eventual equilibrium temperature. The coarsening then continues via spin wave energy transport from low to high wavevectors that displays features of a novel turbulent cascade, relevant to the emerging area of spin turbulence (e.g. see <cit.>). As the transverse spin waves do not interact, their dynamics arises from a dynamic coupling to interacting axial spin degrees of freedom. Order parameter correlations show dynamic scale invariance during the spin wave coarsening, with a length scale that grows as t^1/3. This scaling is distinct from that during the vortex driven coarsening, showing that there are two renormalisation group fixed points affecting the phase ordering of this system. Strongly coupling the system to a reservoir of energy and particles destroys the second scaling regime, allowing the system to thermalise following the annihilation of vortices. Our results give new insights into the phase ordering dynamics of isolated systems and provide a potential profitable connection between phase ordering and wave turbulence.




§ BACKGROUND

A spin-1 condensate can be described by three interacting classical fields ψ_m for condensates in the three spin components with spin projections m=-1,0,1. The quasi-2D Hamiltonian density within a uniform trap <cit.> is <cit.>,

    ℋ=-∑_m=-1^1ψ_m^*p̂^2/2Mψ_m+g_n/2n^2+ℋ_s

where p̂=-iħ∇ is the momentum operator, M is the atom mass, n=∑_m|ψ_m|^2 is the areal density, g_n is the quasi-2D density interaction strength and ℋ_s encompasses the spin dependent terms,

    ℋ_s=g_s/2n^2|𝐅|^2+∑_m=-1^1qm^2|ψ_m|^2.

The first term in ℋ_s is the spin interaction energy, with spin density 𝐅=∑_m m^'ψ_m^* 𝐟_m m^'ψ_m^'/n^2 for spin-1 matrices (f_x,f_y,f_z)≡𝐟, and quasi-2D spin interaction strength g_s. The sign of g_s determines whether the interactions are ferromagnetic (g_s<0), which occurs in ^87Rb <cit.>, or antiferromagnetic (g_s>0), which occurs in ^23Na <cit.>. Here we consider the ferromagnetic case. The second term in ℋ_s is a quadratic Zeeman splitting of the spin components, which can be induced using either DC magnetic fields or AC microwave stark shifts <cit.>. A linear Zeeman term pnF_z can also be included, but conservation of nF_z means this term does not affect the system dynamics and can be removed via the unitary transformation e^-ipmt/ħψ_m→ψ_m. The quasi-2D regime is obtained from a 3D system by tightly confining the system in one direction and integrating over the resulting spatial profile along that direction <cit.>.


The relative strength of the two terms  in ℋ_s produces a rich phase diagram, from which a variety of quenches can be explored. The zero temperature mean field phase diagram for ferromagnetic interactions and with q>0 is shown in Figure <ref>(a). A quantum critical point at q=q_0≡ 2|g_s|n_0 (n_0 is the mean condensate density) separates the unmagnetised polar phase (all atoms in the m=0 condensate) from the easy-plane ferromagnetic phase with spin order parameter 𝐅_⊥≡ (F_x,F_y) (for quantization along F_z). The order parameter manifold of 𝐅_⊥ is SO(2) with transverse spin vortices as topological defects. These vortices consist of a positive or negative phase winding of the transverse spin angle θ (tanθ=F_y/F_x), and can only decay via the mutual annihilation of two vortices of opposite sign. Vortices with negative phase winding are also termed antivortices. The energy scale q_0 defines a time scale t_s≡ħ/q_0 and the spin healing length ξ_s≡ħ/√(Mq_0).



§ RESULTS



 §.§ Quench dynamics: anomalous phase ordering

We simulate the condensate dynamics following an instantaneous quench of the quadratic Zeeman energy from deep in the polar phase to q=0.3q_0 in the easy-plane ferromagnetic phase, see Fig. <ref>(a),(b). Symmetry breaking and the production of transverse spin vortices following such a quench have been observed in experiments with ^87Rb <cit.>. Conservative dynamics of our system is simulated by numerically integrating the three coupled Gross-Pitaevskii equations (GPEs) obtained from Eq. (<ref>) <cit.>,

    iħ∂ψ_m/∂ t=(p̂^2/2M+qm^2+g_nn)ψ_m+g_sn∑_m^'𝐅·𝐟_mm^'ψ_m^'.

Further numerical details are described in Appendix <ref>. A homogeneous system can be realised in experiments using a flat bottomed trap <cit.>.






We quantify order in the system by spatial correlations of 𝐅_⊥,

    G(r,t)≡<𝐅_⊥(𝐫,t)·𝐅_⊥(0,t)>,

where angular brackets denote an ensemble average (see Appendix <ref>). Figure <ref>(c) shows the evolving correlation function Eq. (<ref>). For times 10^2t_s≲ t≲ 10^3t_s, the growth of order is scale invariant and driven by the mutual annihilation of transverse spin vortices of opposite sign, see Fig. <ref>(b), with correlations decaying to zero at a length scale on the order of the intervortex spacing. This vortex driven coarsening has been described in previous work <cit.>. We find that all vortices have annihilated by a time t≈ 2.8× 10^3t_s, after which correlations extend to the boundary. The correlations can then be compared to the equilibrium (thermalised) prediction <cit.>,


    G_eq(r)∼ r^-ν,ν=T_eq/4T_BKT.

Here T_BKT=π K/2k_B is the BKT temperature associated with the unbinding of transverse spin vortices <cit.>, with K=ħ^2n_0(1-q/q_0)/2M the spin wave stiffness and k_B Boltzmann's constant. The equilibrium temperature T_eq of our microcanonical system is calculated by equipartitioning the energy liberated by the quench amongst all collective modes of the system <cit.>, which gives ν≈ 0.011.  This equilibrium prediction is shown in Fig. <ref>(c). Surprisingly, even after very long simulation times t=10^5 t_s, correlations of transverse spin only agree with the equilibrium prediction for length scales r≲ 5ξ_s. For larger length scales the correlations decay more rapidly. This absence of equilibrium following the annihilation of topological defects is not predicted by the current theory of coarsening dynamics <cit.>.




 §.§ Spin wave energy transport driving phase ordering




To identify the origin of the unexpectedly slow ordering displayed in Fig. <ref>(c) we look at the distribution of energy in the gradient of the transverse spin angle ∇θ (this vector field is proportional to currents of F_z magnetization <cit.>). We firstly perform a Helmholtz decomposition ∇θ=𝐯_i+𝐯_c with ∇·𝐯_i=0 and ∇×𝐯_c=0. The first contribution 𝐯_i, known as the incompressible field, arises from vortex excitations while the second contribution 𝐯_c, known as the compressible field, arises from transverse spin wave excitations. The spectral energies of the incompressible and compressible fields are given by,

    ϵ_μ(k,t)=K/2<|𝐯̃_μ(𝐤,t)|^2>,μ=i, c

where 𝐯̃_μ(𝐤)=l^-1∫ d^2𝐫 𝐯_μ(𝐫)e^-i𝐤·𝐫 is the Fourier transform of 𝐯_μ(𝐫) and angular brackets denote an ensemble average (see Appendix <ref>).

The evolving spectral energies ϵ_μ(k,t) are shown in Fig. <ref>(a),(b). The incompressible spectral energy, Fig. <ref>(a), shows a k^-2 decay when vortices are present, in agreement with the infrared (ξ_s k<1) scaling of a distribution of quantum vortices <cit.>. Once all vortices have annihilated the spectral energy drops abruptly. In comparison, the compressible spectral energy, Fig. <ref>(b), shows nonequilibrium features across the duration of the simulation. The initial condition of our simulation results in a flat high energy distribution ϵ_c(k,0)≈ 200k_BT_eq. For times t≳ 10^3t_s, the compressible spectral energy shows three approximate regimes,

    ϵ_c(k,t)={[                  ϵ_lw(k,t),                  k<k_lw(t),; (ϵ_lw(k_lw,t)/k_lw^-α)k^-α,         k_lw(t)≤ k<k_eq(t),;                  k_BT_eq/2,                 k_eq(t)≤ k. ].

We have introduced the evolving wavevectors k_lw(t) and k_eq(t) to signify the boundaries between the three regimes of ϵ_c(k,t). The spectral energy ϵ_lw(k,t) is the long wavelength portion of ϵ_c(k,t), with energy per mode ϵ_lw(k,t)≈ 10 k_BT_eq≈ k_BT_BKT/2 in the wavevector window considered. This nonequilibrium temperature, being approximately at the BKT temperature, corresponds to the typical energy of a single transverse spin vortex <cit.>, and may be a remnant of interactions between spin waves and vortices during the vortex driven coarsening. For k>k_lw ϵ_c(k,t) decays steeply as k^-α with an exponent α≈ 4 until the equilibrium distribution ϵ_c(k,t)=k_BT_eq/2 is reached at a wavevector k_eq(t). The structure of ϵ_c(k,t) is suggestive of a turbulent cascade, with a high temperature long wavelength energy source cascading to a short wavelength thermal field. We provide further evidence of this shortly. With no vortices present, the persistent nonequilibrium features of ϵ_c(k,t) must be responsible for the anomalously slow ordering that we observe in Fig. <ref>(c).

The observed dynamics of ϵ_c(k,t) necessarily involves nonlinear interactions, whereas the transverse spin waves in our system do not interact at any order in the Hamiltonian (which is independent of the phase variable θ). However, the field conjugate to the transverse spin phase θ, i.e. the generator of rotations of θ, is nF_z, leading to dynamic coupling between transverse spin waves and the axial spin waves of F_z <cit.>. Axial spin waves do interact, both between themselves and with other excitations, and must therefore mediate the transverse spin wave interactions. Expanding the system Hamiltonian to quadratic order in F_z and n <cit.> gives the spectral energy of axial spin fluctuations,

    ϵ_F_z(k,t)=n_0ħ^2 k^2/2M+q/2(1-q/q_0)<|F̃_z(𝐤,t)|^2>

where F̃_z(𝐤)≡ l^-1∫ d^2𝐫 F_z(𝐫)e^i𝐤·𝐫 and angular brackets denote an ensemble average (see Appendix <ref>). For times t≳ 400t_s the spectral energy ϵ_F_z(k,t) closely follows ϵ_c(k,t), see Fig. <ref>(b), indicating that the dynamics of the two spectra are coupled and in equilibrium with each other. The nonlinear interactions of axial spin waves allow the redistribution of energy in ϵ_F_z(k,t) and then dynamic coupling to transverse spin waves actuates the same effect in ϵ_c(k,t).

To provide further evidence for the presence of an energy cascade in Fig. <ref>(b) we decompose the total spin wave energy

    E_sw(t)=∑_k 2π k(ϵ_c(k,t)+ϵ_F_z(k,t))

into a low wavevector portion

    E_low(t)=∑_k<k_mid 2π k(ϵ_c(k,t)+ϵ_F_z(k,t))

and a high wavevector portion

    E_high(t)=∑_k≥ k_mid 2π k(ϵ_c(k,t)+ϵ_F_z(k,t)),

where we choose k_mid=0.5ξ_s^-1. In Fig. <ref>(c) we plot the energy changes Δ E(t)≡ E(t)-E(10^5 t_s) of these three quantities for times after all vortices have annihilated. The energy E_low decreases in time while E_high increases, consistent with an energy cascade from k<k_mid to k≥ k_mid. There is also a net decrease in the total spin wave energy E_sw, showing that energy is also lost from the spin wave excitations, either to other quadratic excitations <cit.> or to excitations beyond quadratic order. In principle, one could solve for the dynamics of these additional excitations to obtain effective spin wave dynamics that would transport energy from low to high wavevectors. Figure <ref>(b) shows that the spin wave energy transport is associated with an approximate k^-4 scaling of ϵ_c(k,t) and ϵ_F_z(k,t). (Note the spectral energies in most studies of turbulence include a k phase space factor so that the k^-4 scaling observed here would normally be described as k^-3 scaling.) There are currently no predictions for such a cascade within weak wave turbulence theory <cit.>. To confirm that the energy transport shown in Figs. <ref>(b),(c) is a turbulent cascade would require showing that the energy transport is local in wavevector space.




 §.§ A second regime of scale invariance




The robust shape of ϵ_c(k,t) for times t≳ 10^3t_s (see Fig. <ref>(b)) suggests a regime of scale invariance driven by spin waves, beyond the scale invariant coarsening dynamics driven by vortex annihilation. To explore this we consider the late time dynamics of correlations of transverse spin, Eq. (<ref>), which in a scale invariant regime will evolve as <cit.>,

    G(r,t)=r^-νf(r/L(t))

for some universal function f and growing length scale L(t). The r^-ν correction factor ensures G(r,∞)∼ r^-ν, consistent with equilibrium. Since ν≈ 0.011≪ 1, the correction is only significant when G(r) is close to ordered.

The evolving correlation function for times after all vortices have annihilated is shown in the inset to Fig. <ref>(a). The correlations exhibit a short wavelength ordered portion that grows slowly in time and a nonequilibrium long wavelength portion. The correlation functions collapse onto a single curve after rescaling according to Eq. (<ref>), see Fig. <ref>(a). We define the rescaling factor L_sw(t) by G(L_sw,t)=0.8G(0,t), which follows the boundary between the ordered portion of the correlation function and the nonequilibrium portion. This length scale is governed by spin waves and grows as a power law L_sw∼ t^1/3 for times t≳ 10^3 t_s, i.e. times after all vortices have annihilated, see Fig. <ref>(b). The length scale 2π k_eq(t)^-1, where k_eq(t) is introduced in Eq. (<ref>) and defined more precisely in Appendix <ref>, follows the growth of L_sw(t). For comparison, the scale invariance during vortex driven coarsening is associated with the more rapidly growing average intervortex spacing L_v(t) (defined in Appendix <ref>) <cit.>. The nonequilibrium portions of the correlation functions in Fig. <ref>(a) clearly exhibit an additional algebraic decay G(r)∼ r^-0.21-ν≈ r^-0.22. The value of the decay exponent corresponds to a temperature of T≈ 0.9T_BKT, see Eq. (<ref>), and is consistent with the nonequilbrium temperature of ϵ_lw(k,t) from Eq. (<ref>).




 §.§ Comparison with open system dynamics




Our analysis so far has considered isolated, energy conserving dynamics. It is of interest to compare our results with open system quench dynamics, where the condensate is coupled to a reservoir of energy and particles. Using a stochastic Gross-Pitaevskii theory (see Appendix <ref>), we model a spin-1 condensate strongly coupled to a reservoir with fixed temperature and chemical potential, which we choose such that the equilibriated energy and particle number matches those of the conservative dynamics. We then simulate the same quench as for the isolated system dynamics. Figure <ref>(a) shows the evolution of transverse spin correlations, Eq. (<ref>), for the open system dynamics. The vortex driven coarsening is comparable to the isolated system case, with correlations showing scale invariant growth. For times after t≈ 2× 10^3 t_s correlations in the open system dynamics show excellent agreement with the equilibrium prediction Eq. (<ref>). For comparison, all vortices have annihilated by a time t≈ 1.8× 10^3 t_s. The results in Fig. <ref>(a) are in stark contrast to the results in Fig. <ref>(c) for the isolated system. Indeed, differences in the two cases are apparent from the evolving spin domains, Fig. <ref>(b) and Fig. <ref>(b), with the open system being more ordered in the spaces between vortices. For the large reservoir coupling strength we have used here, spin waves in the open system are able to rapidly thermalise directly with the reservoir rather than via interactions with other spin waves. However, we emphasise that microscopically derived reservoir coupling strengths are much smaller than the value we use here <cit.>, and therefore the isolated system dynamics are a realistic approximation to experiments.

The growing length scales L_v and L_sw for the open system dynamics, defined as for the conservative dynamics, are shown in Fig. <ref>(c). The growth of L_v in the open system is very similar to the isolated system growth (denoted by L_v,ISO in this figure). In the open system, however, there is no second growing length scale, and L_sw follows the growth of L_v.

The decay of transverse spin correlations for open system dynamics following quenches to different values of q show good agreement with Eq. (<ref>) once all vortices have annihilated, see Fig. <ref>(c) inset. (The temperature and chemical potential for these quenches have been adjusted as a function of q; see Appendix <ref>.) The small deviation at the smallest q value may be caused by axial spin fluctuations, which become stronger as q→ 0 due to a diminishing energy gap. Indeed, we expect that the physics will be modified in the limit q→ 0, since the ground state manifold changes from SO(2)×U(1) to SO(3), resulting in changes in collective mode excitations <cit.> and vortex topology <cit.>.




§ CONCLUSION

We have shown that vortex driven coarsening of an isolated easy-plane ferromagnetic spin-1 condensate does not take the system to equilibrium. Instead, a second regime of scale invariant coarsening associated with transport of spin wave energy scales more slowly as t^1/3. Strongly coupling the system to a reservoir of energy and particles destroys this second coarsening process and equilibrium is reached after the vortex driven coarsening.

The presence of two dynamic scaling regimes in the isolated dynamics shows that there are two renormalisation group fixed points affecting the phase ordering. The first, associated with vortices, has been ascribed to the model E dynamic universality class <cit.>. The second slower scaling, L_sw∼ t^1/3, matches that of the scalar model B dynamic universality class <cit.>, however the order parameter does not: the scalar model B universality class describes a one component conserved order parameter, whereas the order parameter in our system has two components and is not conserved. There is, however, a second important field in our system that does satisfy the properties of the scalar model B universality class: the conserved nF_z field. It could be possible that the dynamics of the nF_z field belongs to the scalar model B dynamic universality class, even though this is not the order parameter of the system, and that the dynamic coupling between nF_z and θ leads to model B scaling emerging in the correlations of transverse spin. However, model B is also a dissipative model whereas we have shown that strong dissipation destroys the second scaling regime. Hence the second scaling regime might alternatively belong to a currently unidentified dynamic universality class unique to isolated systems. We have shown that this second scaling regime displays features of a spin wave energy cascade, thus identifying a potential connection between the fields of wave turbulence and phase ordering dynamics.

The nonequilibrium background of spin waves that remain after the vortices have annihilated is at a temperature very close to T_BKT. These spin waves may have thermalised with the vortex field during vortex driven coarsening, either via scattering off of vortices or via spin wave production after vortex annihilation (see <cit.>). The absence of interactions once all vortices have annihilated would then leave behind high temperature spin waves, reminiscent of photons decoupling from matter in the early universe to produce the cosmic microwave background. This intriguing process may be ubiquitous in phase ordering systems involving topological defects interacting with collective mode excitations.



§ ACKNOWLEDGEMENTS

We thank Dan Stamper-Kurn, Kazuya Fujimoto, Matthew Reeves and Ashton Bradley for valuable discussions.








  
Funding information
We acknowledge support from the Marsden Fund of the Royal Society of New Zealand. LW acknowledges support from the University of Otago Postgraduate Publishing Bursary.






§ NUMERICAL DETAILS




 §.§ GPE simulations

The GPE simulations are conducted on a 2D square grid with dimensions l× l=400ξ_s× 400ξ_s covered by an N× N=512× 512 grid of equally spaced points. In experiments in ^87Rb, g_n/|g_s|∼ 100 <cit.>. We use a more modest ratio g_n/|g_s|=10, which is sufficient to suppress density fluctuations at the energy scale we are interested in. The mean condensate density is taken to be n_0=10^4ξ_s^-2. We evolve our system using a recently developed fourth order symplectic integrator <cit.> to ensure that energy, atom number and nF_z magnetization are conserved effectively. We find that total energy and atom number are conserved to within a factor of 10^-9 across the full simulation time. The total axial magnetization ∫ d^2𝐫 n(𝐫)F_z(𝐫) remains below 10^-6n_0l^2.  We use a time step of 0.02t_s for each integration step. The kinetic energy time step is evaluated spectrally using fast Fourier transforms, and we employ periodic boundary conditions. Our initial state is the polar state (ψ_1,ψ_0,ψ_-1)=√(n_0)(0,1,0)+δ, where δ is noise added to Bogoliubov modes on top of the ground state at q=∞, as in <cit.>, which seeds the symmetry breaking evolution. Noise added this way corresponds to adding on average half a particle per mode according to the truncated Wigner prescription <cit.>. We then evolve our system using Eq. (<ref>) at a quadratic Zeeman energy q=0.3q_0, so that the quench is effectively instantaneous at t=0.



 §.§ Ensemble averaging

Correlations Eq. (<ref>) and spectral energies Eq. (<ref>) and Eq. (<ref>) are computed using an ensemble average of the form,

    g̅(𝐮)=<g(𝐮)>

where 𝐮=𝐫, 𝐤 and g(𝐮) denotes the result of a single simulation trajectory. In the GPE simulations, the ensemble average is over 30 simulation trajectories conducted with independent initial noise. In the open system simulations, the ensemble average is over 10 simulation trajectories. We also average g̅(𝐮) over azimuthal angles of the coordinate 𝐮, such that g̅(𝐮)→g̅(u) for u≡ |𝐮|. Correlation functions are additionally averaged over space, i.e. we replace 𝐅_⊥(0)·𝐅_⊥(𝐫) by 𝐅_⊥(𝐫^')·𝐅_⊥(𝐫^'+𝐫) in Eq. (<ref>) and average over the spatial coordinate 𝐫^'.



 §.§ Definition of 𝐤_eq

The wavevector k_eq introduced in Eq. (<ref>) is obtained as follows. We firstly skew the spectrum ϵ_c(k) by multiplying by k. We then define k_eq as the position of the local minimum that appears in kϵ_c(k) at the start of the equilibrium portion of the spectral energy. To improve resolution, we firstly interpolate the numerical values for kϵ_c around its minimum and then find the minimum point of the more highly resolved interpolated data.



 §.§ Vortex detection and averaging

We detect vortices by evaluating the phase winding of the transverse spin angle around plaquettes of our simulation grid. The average intervortex spacings in Fig. <ref>(b) and Fig. <ref>(c) are defined as L_v(t)≡<√(l^2/N_v(t))> where N_v(t) is the vortex number for a single simulation trajectory at time t and angular brackets denote an ensemble average over the 30 simulation trajectories for the GPE results and the 10 simulation trajectories for the open system results. The results for L_v in Fig. <ref>(b) are for the single trajectory displayed.



 §.§ Open system simulations

To model open system evolution we couple our condensate to a reservoir of energy and particles with fixed temperature T and chemical potential μ. The dynamics is simulated using the simple growth stochastic Gross-Pitaevskii equations (SGPEs) <cit.>,

    iħ dψ_m=(1-iγ)(ℒ_m[ψ_m]-μψ_m) dt+dW(𝐫,t).

Here

    ℒ_m[ψ_m]=(p̂^2/2M+qm^2+g_nn)ψ_m+g_sn∑_m^'𝐅·𝐟_mm^'ψ_m^'

is the conservative evolution operator from Eq. (<ref>), μ is the chemical potential and γ is a dimensionless damping. The precise value chosen for γ will not affect equilibrium properties, but will affect the rate that equilibrium is approached. The term dW(𝐫,t) is Gaussian distributed complex noise with delta correlations,

    <dW^*(𝐫,t)dW(𝐫^',t)>=2γ k_B T/ħδ(𝐫-𝐫^')dt.

The SGPEs Eq. (<ref>) take the form of Langevin equations.

The temperature (as a function of q) is chosen to be that obtained by equiparitioning the energy liberated by the quench amongst the 3N^2 numerical modes, as was done in the calculation of T_eq for the microcanonical system. The energy liberated is the energy of the polar state evaluated at the final quadratic Zeemen energy q <cit.>. The temperature is then,

    k_BT=q_0/12N^2(1-q/q_0)^2n_0l^2.

The chemical potential (as a function of q) is chosen to be that of a zero temperature spin-1 condensate in the easy-plane phase, which is obtained by solving ℒ_mψ_m=μψ_m. This gives <cit.>,

    μ=g_nn_0+g_sn_0+q/2.

These choices of T and μ give steady state energy and atom number within 1% of the conservative GPE results. We use γ=10^-2, which gives |dW|≲ q_0|ψ_m|dt, and therefore reservoir scattering events occur within the time scale of spin interactions. The microscopically derived value for γ will be considerably smaller than this <cit.>, resulting in the GPE dynamics overwhelming the reservoir interactions.

We evolve Eq. (<ref>) using an interaction picture fourth order Runge-Kutta integrator with periodic boundary conditions and kinetic energy evaluated to spectral accuracy. The noise is added in a single step following the Runge-Kutta integration of the (1-iγ)(ℒ_m[ψ_m]-μψ_m) term <cit.>. Numerical parameters and initial condition sampling are the same as for the GPE simulations.




























10


urlstyle

bray1994
A. Bray,
Theory of phase-ordering kinetics,
Adv. Phys. 43(3), 357 (1994).

boyanovsky2000
D. Boyanovsky, H. De Vega and R. Holman,
Non-equilibrium phase transitions in condensed matter and
  cosmology: Spinodal decomposition, condensates and defects,
In Topological Defects and the Non-Equilibrium Dynamics of
  Symmetry Breaking Phase Transitions, pp. 139–169. Springer (2000).

damle1996
K. Damle, S. N. Majumdar and S. Sachdev,
Phase ordering kinetics of the Bose gas,
Phys. Rev. A 54, 5037 (1996),
10.1103/PhysRevA.54.5037.

castellano2009
C. Castellano, S. Fortunato and V. Loreto,
Statistical physics of social dynamics,
Rev. Mod. Phys. 81, 591 (2009),
10.1103/RevModPhys.81.591.

sadler2006
L. Sadler, J. Higbie, S. Leslie, M. Vengalattore and D. Stamper-Kurn,
Spontaneous symmetry breaking in a quenched ferromagnetic
  spinor Bose–Einstein condensate,
Nature 443(7109), 312 (2006).

erne2018
S. Erne, R. Büecker, T. Gasenzer, J. Berges and J. Schmiedmayer,
Universal dynamics in an isolated one-dimensional Bose gas
  far from equilibrium,
Nature 563(7730), 225 (2018).

prufer2018
M. Prüfer, P. Kunkel, H. Strobel, S. Lannig, D. Linnemann, C.-M. Schmied,
  J. Berges, T. Gasenzer and M. K. Oberthaler,
Observation of universal dynamics in a spinor Bose gas far
  from equilibrium,
Nature 563(7730), 217 (2018).

weiler2008
C. N. Weiler, T. W. Neely, D. R. Scherer, A. S. Bradley, M. J. Davis and B. P.
  Anderson,
Spontaneous vortices in the formation of Bose-Einstein
  condensates,
Nature 455(7215), 948 (2008),
10.1038/nature07334.

lamporesi2013
G. Lamporesi, S. Donadello, S. Serafini, F. Dalfovo and G. Ferrari,
Spontaneous creation of Kibble-Zurek solitons in a
  Bose-Einstein condensate,
Nat. Phys. 9(10), 656 (2013),
10.1038/nphys2734.

navon2015
N. Navon, A. L. Gaunt, R. P. Smith and Z. Hadzibabic,
Critical dynamics of spontaneous symmetry breaking in a
  homogeneous Bose gas,
Science 347(6218), 167 (2015),
10.1126/science.1258676.

chomaz2015
L. Chomaz, L. Corman, T. Bienaimé, R. Desbuquois, C. Weitenberg,
  S. Nascimbène, J. Beugnon and J. Dalibard,
Emergence of coherence via transverse condensation in a uniform
  quasi-two-dimensional Bose gas,
Nat. Commun. 6, 6162 (2015),
10.1038/ncomms7162.

anquez2016
M. Anquez, B. A. Robbins, H. M. Bharath, M. Boguslawski, T. M. Hoang and M. S.
  Chapman,
Quantum kibble-zurek mechanism in a spin-1 Bose-Einstein
  condensate,
Phys. Rev. Lett. 116, 155301 (2016),
10.1103/PhysRevLett.116.155301.

kang2017
S. Kang, S. W. Seo, J. H. Kim and Y. Shin,
Emergence and scaling of spin turbulence in quenched
  antiferromagnetic spinor Bose-Einstein condensates,
Phys. Rev. A 95, 053638 (2017),
10.1103/PhysRevA.95.053638.

kawaguchi2012
Y. Kawaguchi and M. Ueda,
Spinor Bose–Einstein condensates,
Phys. Rep. 520(5), 253 (2012).

stamperkurn2013
D. M. Stamper-Kurn and M. Ueda,
Spinor Bose gases: Symmetries, magnetism, and quantum
  dynamics,
Rev. Mod. Phys. 85, 1191 (2013),
10.1103/RevModPhys.85.1191.

hofmann2014
J. Hofmann, S. S. Natu and S. Das Sarma,
Coarsening dynamics of binary Bose condensates,
Phys. Rev. Lett. 113, 095702 (2014),
10.1103/PhysRevLett.113.095702.

kudo2013
K. Kudo and Y. Kawaguchi,
Magnetic domain growth in a ferromagnetic Bose-Einstein
  condensate: Effects of current,
Phys. Rev. A 88, 013630 (2013),
10.1103/PhysRevA.88.013630.

kudo2015
K. Kudo and Y. Kawaguchi,
Coarsening dynamics driven by vortex-antivortex annihilation in
  ferromagnetic Bose-Einstein condensates,
Phys. Rev. A 91, 053609 (2015),
10.1103/PhysRevA.91.053609.

williamson2016a
L. A. Williamson and P. B. Blakie,
Universal coarsening dynamics of a quenched ferromagnetic
  spin-1 condensate,
Phys. Rev. Lett. 116, 025301 (2016),
10.1103/PhysRevLett.116.025301.

williamson2016b
L. A. Williamson and P. B. Blakie,
Coarsening and thermalization properties of a quenched
  ferromagnetic spin-1 condensate,
Phys. Rev. A 94, 023608 (2016),
10.1103/PhysRevA.94.023608.

williamson2017
L. A. Williamson and P. B. Blakie,
Coarsening dynamics of an isotropic ferromagnetic superfluid,
Phys. Rev. Lett. 119, 255301 (2017),
10.1103/PhysRevLett.119.255301.

bourges2017
A. Bourges and P. B. Blakie,
Different growth rates for spin and superfluid order in a
  quenched spinor condensate,
Phys. Rev. A 95, 023616 (2017),
10.1103/PhysRevA.95.023616.

symes2017
L. M. Symes and P. B. Blakie,
Nematic ordering dynamics of an antiferromagnetic spin-1
  condensate,
Phys. Rev. A 96, 013602 (2017),
10.1103/PhysRevA.96.013602.

schmied2019
C.-M. Schmied, M. Prüfer, M. K. Oberthaler and T. Gasenzer,
Bidirectional universal dynamics in a spinor Bose gas close
  to a nonthermal fixed point,
Phys. Rev. A 99, 033611 (2019),
10.1103/PhysRevA.99.033611.

schmied2019b
C.-M. Schmied, T. Gasenzer and P. B. Blakie,
Violation of single-length scaling dynamics via spin vortices
  in an isolated spin-1 Bose gas,
arXiv:1904.13222  (2019).

fujimoto2018
K. Fujimoto, R. Hamazaki and M. Ueda,
Unconventional universality class of one-dimensional isolated
  coarsening dynamics in a spinor Bose gas,
Phys. Rev. Lett. 120, 073002 (2018),
10.1103/PhysRevLett.120.073002.

fujimoto2019
K. Fujimoto, R. Hamazaki and M. Ueda,
Flemish strings of magnetic solitons and a nonthermal fixed
  point in a one-dimensional antiferromagnetic spin-1 Bose gas,
Phys. Rev. Lett. 122, 173001 (2019),
10.1103/PhysRevLett.122.173001.

kulczykowski2017
M. Kulczykowski and M. Matuszewski,
Phase ordering kinetics of a nonequilibrium exciton-polariton
  condensate,
Phys. Rev. B 95, 075306 (2017),
10.1103/PhysRevB.95.075306.

Fujimoto2012a
K. Fujimoto and M. Tsubota,
Spin turbulence in a trapped spin-1 spinor Bose-Einstein
  condensate,
Phys. Rev. A 85, 053641 (2012),
10.1103/PhysRevA.85.053641.

Fujimoto2014a
K. Fujimoto and M. Tsubota,
Spin-superflow turbulence in spin-1 ferromagnetic spinor
  Bose-Einstein condensates,
Phys. Rev. A 90, 013629 (2014),
10.1103/PhysRevA.90.013629.

fujimoto2016
K. Fujimoto and M. Tsubota,
Direct and inverse cascades of spin-wave turbulence in spin-1
  ferromagnetic spinor bose-einstein condensates,
Phys. Rev. A 93, 033620 (2016),
10.1103/PhysRevA.93.033620.

Simula2015a
T. Mawson, G. Ruben and T. Simula,
Route to non-abelian quantum turbulence in spinor
  Bose-Einstein condensates,
Phys. Rev. A 91, 063630 (2015),
10.1103/PhysRevA.91.063630.

ho1998
T.-L. Ho,
Spinor Bose condensates in optical traps,
Phys. Rev. Lett. 81, 742 (1998),
10.1103/PhysRevLett.81.742.

ohmi1998
T. Ohmi and K. Machida,
Bose-Einstein condensation with internal degrees of freedom
  in alkali atom gases,
J. Phys. Soc. Jpn. 67(6), 1822 (1998).

stenger1998
J. Stenger, S. Inouye, D. M. Stamper-Kurn, H.-J. Miesner, A. P. Chikkatur and
  W. Ketterle,
Spin domains in ground-state Bose-Einstein condensates,
Nature 396(6709), 345 (1998),
10.1038/24567.

barnett2011
R. Barnett, A. Polkovnikov and M. Vengalattore,
Prethermalization in quenched spinor condensates,
Phys. Rev. A 84, 023606 (2011),
10.1103/PhysRevA.84.023606.

vankempen2002
E. G. M. van Kempen, S. J. J. M. F. Kokkelmans, D. J. Heinzen and B. J.
  Verhaar,
Interisotope determination of ultracold rubidium interactions
  from three high-precision experiments,
Phys. Rev. Lett. 88, 093201 (2002),
10.1103/PhysRevLett.88.093201.

gerbier2006
F. Gerbier, A. Widera, S. Fölling, O. Mandel and I. Bloch,
Resonant control of spin dynamics in ultracold quantum gases by
  microwave dressing,
Phys. Rev. A 73, 041602 (2006),
10.1103/PhysRevA.73.041602.

gaunt2013
A. L. Gaunt, T. F. Schmidutz, I. Gotlibovych, R. P. Smith and Z. Hadzibabic,
Bose-Einstein condensation of atoms in a uniform
  potential,
Phys. Rev. Lett. 110, 200406 (2013),
10.1103/PhysRevLett.110.200406.

chaikin1995
P. M. Chaikin and T. C. Lubensky,
Principles of condensed matter physics,
Cambridge University Press (1995).

kosterlitz1972
J. M. Kosterlitz and D. J. Thouless,
Long range order and metastability in two dimensional solids
  and superfluids. (application of dislocation theory),
J. Phys. C: Solid State Phys. 5(11), L124 (1972).

yukawa2012
E. Yukawa and M. Ueda,
Hydrodynamic description of spin-1 Bose-Einstein
  condensates,
Phys. Rev. A 86, 063614 (2012),
10.1103/PhysRevA.86.063614.

bradley2012
A. S. Bradley and B. P. Anderson,
Energy spectra of vortex distributions in two-dimensional
  quantum turbulence,
Phys. Rev. X 2, 041001 (2012),
10.1103/PhysRevX.2.041001.

billam2015
T. P. Billam, M. T. Reeves and A. S. Bradley,
Spectral energy transport in two-dimensional quantum vortex
  dynamics,
Phys. Rev. A 91, 023615 (2015),
10.1103/PhysRevA.91.023615.

kobayashi2019
M. Kobayashi,
Berezinskii-Kosterlitz-Thouless transition of spin-1
  spinor Bose gases in the presence of the quadratic Zeeman effect,
arxiv:1902.00210  (2019).

murata2007
K. Murata, H. Saito and M. Ueda,
Broken-axisymmetry phase of a spin-1 ferromagnetic
  Bose-Einstein condensate,
Phys. Rev. A 75, 013607 (2007),
10.1103/PhysRevA.75.013607.

uchino2010
S. Uchino, M. Kobayashi and M. Ueda,
Bogoliubov theory and Lee-Huang-Yang corrections in
  spin-1 and spin-2 Bose-Einstein condensates in the presence of the
  quadratic Zeeman effect,
Phys. Rev. A 81, 063632 (2010),
10.1103/PhysRevA.81.063632.

symes2014
L. M. Symes, D. Baillie and P. B. Blakie,
Static structure factors for a spin-1 Bose-Einstein
  condensate,
Phys. Rev. A 89, 053628 (2014),
10.1103/PhysRevA.89.053628.

zakharov1992
V. E. Zakharov, V. S. L'vov and G. Falkovich,
Kolmogorov spectra of turbulence I,
Springer,
10.1007/978-3-642-50052-7 (1992).

nazarenko2011
S. Nazarenko,
Wave turbulence,
Springer,
10.1007/978-3-642-15942-8 (2011).

lee1995
J.-R. Lee, S. J. Lee and B. Kim,
Ordering kinetics of two-dimensional O(2) models: Scaling and
  temperature dependence,
Phys. Rev. E 52, 1550 (1995),
10.1103/PhysRevE.52.1550.

bray2000
A. J. Bray, A. J. Briant and D. K. Jervis,
Breakdown of scaling in the nonequilibrium critical dynamics of
  the two-dimensional 𝑋𝑌 model,
Phys. Rev. Lett. 84, 1503 (2000),
10.1103/PhysRevLett.84.1503.

bradley2014
A. S. Bradley and P. B. Blakie,
Stochastic projected gross-pitaevskii equation for spinor and
  multicomponent condensates,
Phys. Rev. A 90, 023631 (2014),
10.1103/PhysRevA.90.023631.

hohenberg1977
P. C. Hohenberg and B. I. Halperin,
Theory of dynamic critical phenomena,
Rev. Mod. Phys. 49, 435 (1977),
10.1103/RevModPhys.49.435.

williamson2016c
L. A. Williamson and P. B. Blakie,
Dynamics of polar-core spin vortices in a ferromagnetic spin-1
  Bose-Einstein condensate,
Phys. Rev. A 94, 063615 (2016),
10.1103/PhysRevA.94.063615.

symes2016
L. M. Symes, R. I. McLachlan and P. B. Blakie,
Efficient and accurate methods for solving the time-dependent
  spin-1 Gross-Pitaevskii equation,
Phys. Rev. E 93, 053309 (2016),
10.1103/PhysRevE.93.053309.

blakie2008
P. Blakie, A. Bradley, M. Davis, R. Ballagh and C. Gardiner,
Dynamics and statistical mechanics of ultra-cold Bose gases
  using c-field techniques,
Adv. Phys. 57(5), 363 (2008).

gardiner2002
C. Gardiner, J. Anglin and T. Fudge,
The stochastic Gross-Pitaevskii equation,
J. Phys. B: At. Mol. Opt. Phys. 35(6), 1555 (2002),
https://doi.org/10.1088/0953-4075/35/6/310.

gardiner2003
C. Gardiner and M. Davis,
The stochastic Gross-Pitaevskii equation: II,
J. Phys. B: At. Mol. Opt. Phys. 36(23), 4731 (2003).

milstein2004
G. N. Milstein and M. V. Tretyakov,
Stochastic Numerics for Mathematical Physics,
Springer (2004).






