Evolution of Vortex Strings after a Thermal Quench in a Holographic Superfluid

The formation of topological defects during continuous phase transitions exhibits nonequilibrium universality. While the Kibble-Zurek mechanism (KZM) predicts universal scaling of point-like defect numbers under slow driving, the statistical properties of extended defects remain largely unexplored across both slow and fast protocols. We investigate vortex string formation in a three-dimensional holographic superfluid. For slow quenches, the vortex string number follows KZM scaling, while for rapid quenches, it exhibits complementary universal scaling governed by the final temperature. Beyond the vortex string number, the loop-length distribution reveals a richer structure: individual loops follow the first-return statistics of three-dimensional random walks, $P(\ell) \sim \ell^{-5/2}$. While the total vortex length distribution remains Gaussian, its cumulants obey universal scaling laws with varying power-law exponents, and thus differ markedly from those observed in point-defect systems, indicating distinct statistical features of extended topological defects.

💡 Research Summary

In this work the authors study the non‑equilibrium formation of vortex strings in a three‑dimensional strongly coupled superfluid using holographic duality. The bulk theory consists of a charged scalar field in a five‑dimensional AdS‑Schwarzschild black‑hole background; the Hawking temperature of the black hole sets the temperature of the dual 3+1‑dimensional field theory. By linearly lowering the temperature across the critical point (T(t)=T_c(1−t/τ_Q) until a final temperature T_f is reached) they drive a continuous U(1) symmetry‑breaking transition. Near the critical point the equilibrium correlation length ξ∝|ε|^{-ν} and relaxation time τ∝|ε|^{-zν} obey mean‑field exponents ν=1/2 and z=2. The Kibble‑Zurek mechanism (KZM) predicts a freeze‑out time ˆt∝τ_Q^{zν/(1+zν)} and a corresponding freeze‑out length ˆξ∝τ_Q^{ν/(1+zν)} that set the size of causally disconnected domains.

Within each domain vortex loops are nucleated independently. The authors model the number of loops N as a binomial process with success probability p determined by the closure condition of a random walk. Consequently all cumulants κ_n(N) scale with the mean, leading to universal power‑law behavior both for slow quenches (τ_Q≫100) and fast quenches (τ_Q≲100):

• Slow regime: κ_1, κ_2 ∝ τ_Q^{-3/4}, in agreement with the generalized KZM prediction for line‑like defects (D=3, d=1).
• Fast regime: κ_1, κ_2 ∝ ε_f^{3/2}, where ε_f=(T_f−T_c)/T_c is the quench depth. This reflects the breakdown of the usual KZM scaling and the dominance of the final temperature in setting the freeze‑out length.

Beyond the loop count, the distribution of individual loop lengths ℓ is examined. By invoking the first‑return statistics of a three‑dimensional random walk, the authors predict a power‑law tail \