Stable Quantum Vortices in Lee-Huang-Yang Dipolar Superfluids

This study numerically demonstrates that Lee-Huang-Yang corrections stabilize rotating dipolar Bose-Einstein condensates into a robust pure LHY superfluid state, revealing distinct critical frequencies and asymmetric frequency ranges for the formation of even versus odd numbers of vortices due to the system's nonlinearity and trapping geometry.

The Gist

The Big Picture: A Quantum Dance Floor

Imagine a ballroom filled with millions of dancers. In a normal crowd, everyone bumps into each other, pushes, and pulls in a chaotic mess. But in a Bose-Einstein Condensate (BEC), these dancers are so cold and calm that they stop acting like individuals and start moving as a single, giant "super-dancer." This is a superfluid: a liquid with zero friction that can flow forever without slowing down.

Now, imagine spinning this ballroom. If you spin a bucket of water, it stays flat until it spins fast enough, then a whirlpool (a vortex) forms in the middle. In a superfluid, these whirlpools are special. They are quantized, meaning they can't be just any size; they are like distinct, indivisible coins. You can have one coin, two coins, or ten coins, but never "half a coin."

This paper explores what happens when you spin this super-dancer ballroom, but with a twist: the dancers have two types of personalities (interactions) that usually fight each other, and a hidden "quantum safety net" that keeps them from falling apart.

---

The Cast of Characters

To understand the experiment, we need to know the three forces at play:

1. The Contact Push (Mean-Field): Imagine the dancers are wearing puffy jackets. If they get too close, they push each other away. This is the standard repulsion in these gases.
2. The Magnetic Pull (Dipolar Interaction): Now, imagine every dancer is also a tiny magnet. Depending on how they are oriented, they might pull toward each other or push away. The scientists can use a "remote control" (magnetic fields) to change the angle of these magnets, making them pull or push at will.
3. The Quantum Safety Net (Lee-Huang-Yang or LHY): This is the star of the show. In the past, if the dancers pushed too hard or pulled too hard, the whole group would collapse or fly apart. But quantum mechanics has a secret rule: even when things look like they should collapse, tiny "quantum fluctuations" create a repulsive force that acts like a safety net. This is the LHY correction.

The Experiment: Canceling the Chaos

The scientists wanted to create a very specific, rare state of matter: a Pure LHY Superfluid.

Think of it like balancing a seesaw.

• On one side, you have the "Contact Push" (repulsion).
• On the other side, you have the "Magnetic Pull" (attraction).

Usually, one side wins, and the system behaves normally. But the scientists used their "remote control" to perfectly balance the seesaw. They adjusted the magnets so the pull exactly canceled out the push.

The Result: The usual forces vanished! The only thing left holding the dancers together was the Quantum Safety Net (LHY). This created a fluid that exists only because of quantum fluctuations. It's like a building held up entirely by invisible, vibrating energy rather than bricks or steel.

The Discovery: The "Even Number" Rule

Once they created this fragile, pure quantum fluid, they started spinning the ballroom (increasing the rotation speed, $\Omega$) to see how many whirlpools (vortices) would form.

Here is where it gets weird and wonderful:

• In normal fluids: As you spin faster, you get 1 vortex, then 2, then 3, then 4. It's a steady, predictable climb.
• In this Pure LHY fluid: The dancers are picky.
  • Getting 1 or 3 vortices is incredibly hard. It's like trying to balance a pencil on its tip. You have to spin at a perfectly specific speed. If you speed up or slow down by a tiny fraction, the single vortex disappears or splits into two.
  • Getting 2 or 4 vortices is easy. The fluid loves even numbers. There is a wide range of speeds where 2 or 4 vortices are stable and happy.

The Analogy: Imagine trying to stack blocks.

• With normal forces, you can stack 1, 2, 3, or 4 blocks easily.
• With this Pure LHY fluid, stacking 1 block is like balancing it on a razor blade (very unstable). Stacking 2 blocks is like putting them in a sturdy box (very stable). The fluid strongly prefers to form pairs or groups of four.

Why Does This Matter?

1. Robustness: The paper shows that this "Pure LHY" state is surprisingly strong. Even though it relies on a delicate balance of forces, once formed, it creates a very stable state of matter with a deep energy minimum. It's a "super-stable" quantum state.
2. New Physics: It proves that we can create fluids where the only thing keeping them together is the weirdness of quantum mechanics (fluctuations), not the usual atomic forces.
3. The "Magic" Angle: The study highlights how tuning the magnetic angle (the "magic angle") allows scientists to turn the usual forces off and leave only the quantum safety net active.

The Takeaway

The scientists successfully built a "ghost" fluid—a superfluid held together only by quantum whispers (fluctuations) because they canceled out the loud shouting of normal atomic forces.

When they spun this ghost fluid, they discovered it has a strange preference: it hates being alone (1 vortex) or in a trio (3 vortices) but loves being in pairs (2 vortices) or groups of four. This reveals a new, hidden rule of nature for how quantum matter organizes itself when pushed to its limits.

In short: They turned off the lights in the quantum ballroom, leaving only the "safety net" active, and found that the dancers only want to dance in even-numbered circles.

Technical Summary

1. Problem Statement

The paper addresses the stability and dynamics of quantized vortices in rotating dipolar Bose-Einstein condensates (BECs). While vortices in standard BECs (governed by mean-field contact interactions) and dipolar BECs (governed by long-range dipole-dipole interactions, DDI) are well-studied, the specific regime where quantum fluctuations dominate is less explored.

In standard mean-field (MF) theory, attractive interactions often lead to collapse. However, the Lee-Huang-Yang (LHY) correction, which accounts for quantum fluctuations beyond the MF approximation, can stabilize these systems, forming "quantum droplets." The authors investigate a specific, idealized scenario: a pure LHY superfluid. In this state, the attractive and repulsive mean-field terms (contact and dipolar interactions) are tuned to cancel each other out, leaving the LHY correction as the sole source of nonlinearity. The study aims to determine how vortices nucleate, their stability, and their critical rotational frequencies in this unique quantum state.

2. Methodology

The authors employ a theoretical and numerical approach based on the Gross-Pitaevskii (GP) equation extended to include DDI and LHY terms.

• Model System: A quasi-two-dimensional (quasi-2D) dipolar BEC confined in a harmonic-oscillator trap with a large aspect ratio ($\lambda = \omega_z/\omega_\perp = 20$), creating a "pancake" geometry.
• Governing Equation: The time-dependent GP equation includes:
  • Kinetic energy and trap potential.
  • Rotation term ($-\Omega L_z$).
  • Cubic MF term (contact interaction, $g|\psi|^2$).
  • Integral term for long-range DDI ($g_{dd}$).
  • Quintic LHY term ($g_{LHY}|\psi|^3$), representing quantum fluctuations.
• Dimensional Reduction: The 3D wavefunction is factorized into a fixed Gaussian in the $z$-direction and an effective 2D wavefunction $\psi(\rho, t)$. The DDI kernel is reduced to 2D using Fourier transforms and angular averaging.
• Pure LHY Regime Strategy: To isolate the LHY effect, the authors assume the cancellation of MF terms. By tuning the contact interaction (via Feshbach resonance) and the DDI orientation angle ($\phi$), they set the effective MF nonlinearity to zero. The system is then governed solely by the repulsive quintic LHY term.
• Numerical Simulation:
  • The 2D GP equation is solved using a split-step Crank-Nicholson method combined with Fast Fourier Transform (FFT).
  • A grid size of $512 \times 512$ is used with a time step of $\Delta t = 0.01$.
  • Simulations track the evolution of density profiles, phase patterns, energy per particle ($E$), chemical potential ($\mu$), and angular momentum ($\langle L_z \rangle$) as the rotational frequency ($\Omega$) is increased.

3. Key Contributions

• First Detailed Study of Vortices in Pure LHY Superfluids: While LHY-stabilized droplets have been observed experimentally, this paper provides the first systematic theoretical analysis of vortex nucleation and stability specifically in a regime where MF interactions are nullified.
• Identification of Critical Frequencies: The authors determine the exact critical rotational frequency ($\Omega_c$) required to create the first single vortex in a pure LHY superfluid, a value that is extremely sensitive to parameter tuning.
• Discovery of Stability Asymmetry: A novel finding is the stark difference in stability between even and odd numbers of vortices. The system exhibits large frequency ranges for stable even numbers of vortices (2, 4, 8) but extremely narrow windows for odd numbers (1, 3).
• Healing Length Analysis: The paper derives and compares the healing length (core size) of vortices in LHY fluids versus standard MF fluids, showing that LHY vortices have larger cores.

4. Key Results

A. Vortex Nucleation and Stability

• Critical Frequency ($\Omega_c$): For the pure LHY case (with LHY strength $\eta = 200$), the critical frequency for the first single vortex is found to be $\Omega_c \approx 0.6401$. This requires extremely fine-tuning; a change of $0.0001$ in $\Omega$ can prevent vortex formation or immediately generate two vortices.
• Even vs. Odd Vortex Stability:
  • Even numbers (2, 4, 8): These states are robust and exist over wide intervals of $\Omega$. For example, the 2-vortex state exists for $0.6410 \le \Omega \le 0.767$.
  • Odd numbers (1, 3): These states are highly unstable and exist only in very narrow frequency windows. The 3-vortex state, for instance, persists only between $\Omega = 0.768$ and $0.778$.
  • Interpretation: This behavior is attributed to the quartic nonlinearity of the LHY term combined with the trap's aspect ratio, which favors symmetric (even) vortex configurations.

B. Energetics and Angular Momentum

• Energy and Chemical Potential: As $\Omega$ increases, the energy per particle ($E$) and chemical potential ($\mu$) decrease, indicating that vortex formation is energetically favorable at higher rotation rates.
• Energy Cost: Generating vortices in a pure LHY fluid is more "costly" than in a standard MF contact fluid. At $\Omega=0.7$, the energy cost to create a pair of vortices in the LHY case is $\approx 22\%$ higher than the non-rotating state, compared to $\approx 14\%$ for pure MF contact interactions.
• Feynman's Relation: The angular momentum per particle ($\langle L_z \rangle$) shows a linear dependence on the number of vortices ($N_v$), consistent with Feynman's prediction for superfluids ($\langle L_z \rangle \propto N_v$).

C. Healing Length

• The healing length (core radius) in the pure LHY fluid is derived as $\xi_{LHY} \approx \hbar/\sqrt{3m\mu_{LHY}}$.
• Comparing a 4-vortex state in pure MF vs. pure LHY, the authors find that $\xi_{LHY} \approx 1.10 \times \xi_{MF}$. The vortices in the LHY superfluid have cores about 10% larger than those in standard contact-interaction BECs.

5. Significance

• Experimental Guidance: The paper provides specific targets (critical frequencies, interaction strengths, and stability windows) for experimentalists working with dipolar atoms (e.g., Dysprosium, Erbium) to realize and observe pure LHY superfluids with vortices.
• Fundamental Physics: It demonstrates that quantum fluctuations (LHY) can not only stabilize matter against collapse but also support complex topological structures (vortex lattices) with unique stability properties distinct from mean-field systems.
• Nonlinearity Effects: The results highlight how the order of nonlinearity (cubic in MF vs. quintic in LHY) fundamentally alters the phase diagram of rotating superfluids, specifically the preference for even-numbered vortex clusters.
• Robustness: The identification of the "pure LHY superfluid" as a remarkably robust state of quantum matter with deep energy minima suggests new avenues for creating stable quantum droplets and supersolids in rotating frames.

In conclusion, this work establishes that while the pure LHY regime is numerically challenging to access due to the narrow stability windows for odd vortices, it supports stable, robust vortex lattices with distinct physical properties (larger cores, higher energy costs) compared to traditional BECs.