Formation of holographic vortex in a rotating shell-shaped superfluid

Technical Summary: Formation of Holographic Vortex in a Rotating Shell-Shaped Superfluid

Problem Statement
This work investigates the dynamics of a superfluid confined to a spherical shell geometry under external rotation. While superfluidity in planar geometries and Bose-Einstein condensates (BECs) in various shapes (disk, cigar, toroidal) are well-studied, shell-shaped systems have recently become experimentally accessible via microgravity platforms (e.g., the Cold Atom Lab on the ISS). Theoretical understanding of these systems faces limitations: the Gross-Pitaevskii (GP) equation is restricted to ground-state dynamics and requires phenomenological additions for finite-temperature effects, while conventional statistical mechanics is limited to equilibrium properties. Furthermore, standard hydrodynamics often fails to describe far-from-equilibrium processes involving topological defects. To address these gaps, the authors employ the AdS/CFT correspondence (holographic duality) to model finite-temperature superfluid dynamics in a rotating shell, mapping the problem to a gravitational system in an Anti-de Sitter (AdS) black hole spacetime.

Methodology
The study utilizes a holographic superfluid model based on the Abelian-Higgs model coupled to Einstein's gravity in an asymptotically AdS$_4$ spacetime.

• Setup: The authors consider a spherically symmetric Schwarzschild-AdS$_4$ background. They work in the probe limit ($q \to \infty$), neglecting the backreaction of matter fields on the metric. The boundary system lives on a sphere with radius $L$.
• Static Solutions: To model the rotating shell, they impose an ansatz for the bulk scalar field $\Psi$ and gauge field $A_\mu$ that depends on the radial coordinate $z$ and polar angle $\theta$. The rotation is introduced via boundary conditions where the gauge field component $A_\phi$ at the boundary corresponds to an external angular velocity $\Omega$. The system is analyzed in the canonical ensemble with a fixed average charge density $\bar{\rho}$.
• Linear Stability Analysis: The authors perform a linear perturbation analysis around the static rotating solutions. They calculate Quasi-Normal Modes (QNMs) using ingoing Eddington-Finkelstein coordinates to ensure regularity at the black hole horizon. The perturbations are decomposed into modes characterized by an azimuthal quantum number $p$ and parity (even or odd under reflection across the equator).
• Nonlinear Dynamics: Fully nonlinear numerical simulations are conducted to evolve the system from static configurations perturbed by the dominant unstable modes. The simulations track the time evolution of the condensate amplitude and gauge fields to observe the formation and migration of topological defects.

Key Results

1. Critical Temperature Enhancement: The study confirms that the critical temperature for the onset of superfluidity is enhanced on a spherical geometry compared to a planar one. This enhancement persists and becomes more pronounced as the external angular velocity increases.
2. Linear Instability and Phase Diagram: The linear stability analysis reveals that the rotating superfluid becomes dynamically unstable beyond specific critical angular velocities. The instability is driven by modes with positive imaginary frequencies.
   • The unstable modes are classified into even and odd sectors based on their symmetry with respect to the equatorial plane.
   • Even modes are found to be dominant over odd modes.
   • Three critical angular velocities ($\Omega_a, \Omega_b, \Omega_c$) are identified, marking transitions where the dominant unstable mode shifts from $p=1$ to $p=2$ and then to $p=3$.
   • A dynamical phase diagram in the $T-\Omega$ plane is constructed, delineating regions where the system's evolution is controlled by specific even unstable modes ($p=1, 2, 3$).
3. Nonlinear Dynamical Evolution:
   • Vortex Nucleation: When the system is perturbed by unstable modes, vortex-antivortex pairs nucleate at the equator, where the superfluid velocity is maximal and the condensate amplitude is minimal (Landau instability).
   • Mode Dependence: The number of nucleated pairs is determined by the symmetry of the triggering mode.
     • For even modes with azimuthal number $p$, $p$ vortex-antivortex pairs nucleate.
     • For odd modes with azimuthal number $p$, $2p$ vortex-antivortex pairs nucleate.
   • Migration and Accumulation: Following nucleation, the vortices and antivortices separate and spiral along the direction of rotation toward the poles. For $p > 1$, this process results in the accumulation of vortices at the poles, forming a "giant vortex."
   • Winding Numbers: The final winding number $w$ of the vortex at the north pole is directly linked to the mode number $p$ and the rotation direction. For $\Omega > 0$, $w = -p$ for even modes and $w = 2p$ for odd modes. The sign reverses for $\Omega < 0$.
4. Density Constraints: The authors identify a critical average number density $\bar{\rho}^c_{\min} = 1.891$ below which no superfluid solution exists. A higher threshold, $\bar{\rho}^s_{\min} = 2.1998$, is required to support a stable rotating shell superfluid; below this density, the system is unstable even at low angular velocities.

Significance and Claims
The paper claims to extend holographic superfluid studies from uniform configurations to inhomogeneous, rotating shell geometries. By combining linear instability analysis with fully nonlinear simulations, the work provides a first-principles framework for understanding how rotation induces topological defects in finite-temperature superfluids on curved manifolds. The authors assert that their results offer a foundation for exploring other inhomogeneous structures (such as solitons, skyrmions, and turbulence) within the AdS/CFT correspondence. They note that the phase transition phenomena and dynamical behaviors described await future experimental verification in shell-shaped BEC systems. The study highlights the universality of the spherical geometry's effect on critical temperatures and the specific role of mode symmetry in determining the topology of the resulting vortex configurations.