Vortex Retention Mediated Turbulent Transitions in Self-Gravitating Bosonic and Axionic Condensates

This study demonstrates that while both self-gravitating bosonic and axionic condensates undergo turbulent transitions from Kolmogorov to Vinen scaling during rapid spin-down, axionic condensates exhibit enhanced vortex retention that increasingly disrupts Kolmogorov scaling as interaction strength grows, driven primarily by quantum pressure during vortex detachment.

The Gist

The Big Picture: Spinning Ice Cream and Dark Matter

Imagine you have two giant, magical bowls of ice cream. These aren't normal ice cream; they are made of quantum fluid (specifically, Bose-Einstein Condensates). In this state, all the particles act like a single, giant super-particle. They are so cold and calm that they flow without any friction, like a ghost.

Now, imagine these bowls are spinning very fast, like a figure skater. But then, someone suddenly tells them to slow down.

This paper is a computer simulation of what happens when these two different types of "quantum ice cream" slow down. The scientists wanted to see how the fluid reacts, how it gets messy (turbulent), and how it handles the "friction" caused by the bowl's surface.

The Two Types of Ice Cream

The researchers compared two specific recipes:

1. The Standard Bosonic Ice Cream: This is the "normal" quantum fluid. The particles like to interact with each other in a standard way.
2. The Axionic Ice Cream: This is a special, more complex recipe. It includes "higher-order" interactions (think of it as adding a secret spice that makes the particles behave differently). This type is often used by scientists to model Dark Matter (the invisible stuff holding galaxies together) or the inside of Neutron Stars (super-dense dead stars).

The Experiment: The "Spin-Down"

In the real world, neutron stars spin down over millions of years. In this computer simulation, the scientists made them slow down very quickly to see what happens.

They also added a "crust" to the bottom of the bowl. Think of this crust like a Velcro floor. As the ice cream spins, tiny whirlpools (called vortices) get stuck to the Velcro.

When the spinning slows down, the fluid wants to keep moving, but the Velcro holds the whirlpools back. Eventually, the stress gets too high, and the whirlpolds rip free all at once. This is called "vortex depinning."

What Happened? The Great Unraveling

When the whirlpools ripped free, they didn't just disappear. They started crashing into each other, creating a chaotic mess. This is Quantum Turbulence.

Here is what the scientists found when comparing the two ice creams:

1. The Shape of the Bowl

• Standard Ice Cream: When it slowed down, it stayed relatively fluffy and spread out.
• Axionic Ice Cream: Because of its special "spice," it stayed much tighter and more compact. It was like a dense, heavy ball of dough compared to the fluffy standard ice cream.

2. The "Velcro" Effect (Vortex Retention)

This is the most important discovery.

• In the Standard case, when the whirlpools ripped free, they flew out and broke apart easily. This created a beautiful, classic pattern of chaos (called a Kolmogorov cascade), which is how turbulence usually works in nature (like water swirling down a drain).
• In the Axionic case, because the fluid was so tight and dense, the whirlpools stuck around longer. They didn't break apart as easily. They got "retained" in the center.

The Analogy: Imagine a crowd of people running out of a stadium.

• In the Standard case, the crowd runs out, spreads into the parking lot, and disperses quickly.
• In the Axionic case, the crowd is so tightly packed that they get stuck in the exit gate. They stay bunched up, refusing to spread out.

3. The Energy Crash

Because the Axionic whirlpools stayed stuck together, the "classic" turbulence pattern (the smooth energy cascade) broke down. The energy didn't flow smoothly from big swirls to small swirls. Instead, the Axionic system held onto its energy in a chaotic, stuck state for longer.

Why Does This Matter?

This isn't just about ice cream. It helps us understand the universe:

• Neutron Stars: When a neutron star spins down, it sometimes has a "glitch" (a sudden jump in speed). This happens because the superfluid inside rips free from the crust. This paper suggests that if the star's interior is made of "Axionic" matter, these glitches might look different because the whirlpools get stuck more easily.
• Dark Matter: If the universe is made of this "Axionic" dark matter, it might behave differently when galaxies collide or spin down than we previously thought.

The Bottom Line

The study shows that how the tiny particles interact (the "micro" physics) changes how the whole system behaves (the "macro" physics).

Even though both fluids started the same way, the "Axionic" fluid held onto its chaos longer. It was like a stubborn child refusing to let go of a toy, while the "Standard" fluid let go and moved on. This "stubbornness" changes how energy is dissipated in the universe, from the smallest quantum fluids to the largest stars.

Technical Summary

1. Problem Statement and Motivation

The paper investigates the turbulent spin-down dynamics of self-gravitating Bose-Einstein condensates (BECs), specifically comparing purely bosonic systems (standard two-body interactions) with axionic systems (featuring higher-order, three-body interactions).

• Context: Quantum turbulence (QT) in self-gravitating fluids is relevant to astrophysical phenomena such as pulsar glitches in neutron stars (where a superfluid core interacts with a solid crust) and fuzzy dark matter halos (composed of ultra-light axion-like particles).
• The Gap: While the transition from laminar flow to turbulence is well-studied in atomic BECs, the influence of higher-order nonlinearities (axionic terms) on vortex dynamics, energy cascades, and dissipation mechanisms in self-gravitating systems remains an open question.
• Core Question: Do bosonic and axionic condensates follow the same route to turbulence during rapid spin-down, or do their distinct microscopic nonlinearities lead to fundamentally different macroscopic dissipation pathways?

2. Methodology

The authors employ direct numerical simulations of the coupled Gross-Pitaevskii-Poisson (GPP) system in three dimensions.

• Governing Equations:
  • Gross-Pitaevskii Equation (GPE): Describes the condensate wavefunction $\psi$. It includes a standard cubic interaction term ($g_{3D}|\psi|^2$) and, for axionic systems, a quintic interaction term ($g_{3D2}|\psi|^4$).
  • Poisson Equation: Solves for the self-gravitational potential $\Phi$ generated by the condensate density $|\psi|^2$.
  • External Potentials: The system includes a centrifugal potential (due to rotation) and a static crust potential ($V_{crust}$) to model vortex pinning sites, analogous to the neutron star crust.
• Simulation Protocol:
  • Initial State: A rapidly rotating equilibrium state ($\Omega_0 = 8.0$) with vortices pinned to the crust potential.
  • Spin-Down: A linear ramp-down of the rotation frequency $\Omega(t)$ to a final frequency $\Omega_f = 1.0$ over a time $t_s \approx 2.3$. This induces a Magnus flow that eventually exceeds the critical velocity for vortex depinning.
  • Numerical Scheme: The equations are solved using a split-step Crank-Nicolson method combined with a spectral collocation approach (FFT-based) for the Poisson solver. Simulations are performed on $256^3$ and $512^3$ grids.
• Analysis Tools:
  • Decomposition of kinetic energy into incompressible (vortical), compressible (sound waves), and quantum pressure components.
  • Spectral analysis of energy spectra $E(k)$ to identify scaling laws (Kolmogorov vs. Vinen).
  • Tracking of vortex population decay and radial width evolution.

3. Key Contributions and Results

A. Static Properties and Vortex Nucleation

• Density Profiles: Axionic condensates (with attractive two-body and repulsive quintic terms) exhibit more uniform and compact density profiles compared to bosonic condensates for similar interaction strengths.
• Critical Frequency: The critical rotation frequency ($\Omega_c$) required for vortex nucleation is lower in axionic systems. The compact structure of axionic condensates facilitates easier vortex entry compared to the more extended bosonic counterparts.

B. Turbulent Transition and Energy Cascades

Upon rapid spin-down, both systems undergo vortex depinning, leading to a turbulent state characterized by two distinct regimes:

1. Kolmogorov Regime ($k^{-5/3}$): Immediately following spin-down, both systems exhibit a short-lived inertial range with Kolmogorov scaling, indicative of a direct energy cascade from large to small scales.
2. Vinen Regime ($k^{-1}$): As the vortex tangle decays, the spectrum transitions to Vinen turbulence scaling, characteristic of a sparse network of uncorrelated vortices.
3. Compressible Spectrum: The compressible energy spectrum shows thermalization ($k$ scaling) and a high-wavenumber decay ($k^{-7/2}$), driven by the expulsion of density fluctuations to the periphery.

C. The Divergence: Vortex Retention and Interaction Strength

The most significant finding is the divergence in behavior as interaction strength (and thus condensate size) increases:

• Bosonic Systems: As $g_{3D}$ increases, the condensate size grows, allowing for more vortices. The system maintains the Kolmogorov ($k^{-5/3}$) cascade relatively well, though with fluctuations.
• Axionic Systems: As $g_{3D2}$ increases, the condensate size saturates (stops growing) due to the balance of attractive and repulsive forces. Crucially, axionic systems exhibit enhanced vortex retention.
  • The compact, uniform structure and higher Magnus-flow thresholds prevent vortices from depinning and escaping efficiently.
  • This retention suppresses the vortex breakdown process necessary for the Kolmogorov cascade.
  • Result: Axionic spectra show a pronounced deviation from Kolmogorov scaling and a suppression of the $k^{-5/3}$ regime at high interaction strengths, whereas bosonic systems retain the cascade.

D. Energy Injection Mechanisms

• Quantum Pressure Dominance: The growth of incompressible energy is primarily driven by quantum pressure during vortex detachment, rather than by compressible flows.
• Exchange Channels: In hard-walled traps, compressible-to-incompressible energy exchange dominates in large systems. However, in self-gravitating systems, the lack of reflective boundaries causes compressible waves to be expelled quickly. Consequently, the quantum-pressure-driven injection remains the dominant mechanism for incompressible energy growth in both systems, though the ratio of this dominance ($\eta$) is lower in axionic systems due to enhanced compressible excitation generation.

4. Significance and Implications

• Astrophysical Relevance: The study provides a microscopic basis for understanding pulsar glitches. It suggests that the nature of the superfluid component (bosonic vs. axionic/higher-order) critically determines the efficiency of vortex unpinning and the subsequent energy dissipation.
• Dark Matter Models: The results offer insights into the turbulent dynamics of fuzzy dark matter halos. The finding that axionic interactions promote vortex retention and suppress classical turbulence cascades could influence models of structure formation and relaxation in dark matter halos.
• Universal vs. System-Specific Turbulence: The work demonstrates that while universal features (like the initial Kolmogorov cascade) exist, the microscopic interaction potentials (specifically higher-order terms) can fundamentally alter the macroscopic turbulent dissipation pathways, leading to system-specific deviations from classical scaling laws.

Conclusion

The paper establishes that while both self-gravitating bosonic and axionic condensates undergo similar initial turbulent transitions upon spin-down, their long-term turbulent dynamics diverge significantly. Axionic condensates, due to their compact structure and higher-order nonlinearities, exhibit enhanced vortex retention, which suppresses the classical Kolmogorov energy cascade. This highlights the critical role of microscopic interaction terms in governing the macroscopic dissipation and stability of quantum fluids in astrophysical environments.