Approximate vortex lattices of atomic Fermi superfluid on a spherical surface

This paper characterizes approximate vortex structures of atomic Fermi superfluids on a spherical surface under an effective monopole field by employing two Ginzburg-Landau-based constructions—geometric scaffolding and free-energy minimization—which yield vortex configurations that converge to the planar Abrikosov lattice limit as the number of vortices increases.

The Gist

Imagine you are a master architect trying to build a perfect city of skyscrapers. On a flat piece of land (like a standard city), you can lay out a grid of streets and buildings in a perfect, repeating pattern forever. This is what happens in flat superfluids: when you spin them or apply a magnetic field, they form a perfect, repeating grid of tiny whirlpools called vortices. This is known as an Abrikosov lattice.

But now, imagine you have to build that same perfect city on the surface of a giant, smooth beach ball.

The Problem: The "Beach Ball" Limit

The paper tackles a fundamental geometric problem: You cannot tile a sphere perfectly with more than 20 points.

Think of it like trying to cover a soccer ball with hexagons (like a honeycomb). No matter how hard you try, you eventually hit a spot where the shape doesn't fit. You are forced to introduce "defects"—places where you have pentagons instead of hexagons. In math, this is why you can't have a perfect, infinite crystal lattice on a sphere. If you have more than 20 vortices, the lattice must be imperfect.

So, the question the authors ask is: If we can't build a perfect city on a sphere, what does the "best possible" imperfect city look like?

The Setup: The Atomic Bubble

The researchers are studying atomic Fermi superfluids. Imagine a cloud of ultra-cold atoms (so cold they act like a single quantum fluid) trapped in a thin, hollow shell, like a soap bubble.

To make these atoms swirl and form vortices, they use a "synthetic magnetic monopole."

• The Metaphor: Imagine a magnet that has only a North pole and no South pole, sitting right in the center of the bubble.
• The Effect: This creates a magnetic field that radiates outward in all directions, like the spokes of a wheel or the rays of the sun. This forces the atoms to spin and create tiny whirlpools (vortices) on the surface of the bubble.

The Two Approaches: "The Blueprint" vs. "The Sculptor"

The authors used two different methods to figure out where these whirlpools should sit to be as happy (energetically stable) as possible.

1. The Geometric Approach (The Blueprint)

Instead of calculating every single atom, they used pre-existing patterns as scaffolds (blueprints) to place the whirlpools. They tried three types of patterns:

• Random: Throwing darts at the sphere. (Predictably messy).
• Geodesic Dome: Like the structure of a geodesic dome or a soccer ball. This is based on dividing triangles. It works great for small numbers of whirlpools but gets messy and uneven as you add more.
• Fibonacci Lattice: This is the "Golden Ratio" pattern. Imagine a spiral that winds from the North Pole to the South Pole, placing points at a specific irrational angle (the Golden Ratio). This creates a pattern that looks incredibly uniform, almost like a honeycomb, even on a sphere.

The Result: They built a quantum wave function that had "zeros" (empty spots) exactly where these blueprints said to put them. They checked the physics and confirmed: Yes, these are real whirlpools with current spinning around them.

2. The Minimization Approach (The Sculptor)

This method didn't use a blueprint. Instead, it used a computer to act like a sculptor.

• The computer started with a guess and then constantly tweaked the position and strength of the whirlpools.
• It asked: "If I move this whirlpool a tiny bit, does the total energy of the system go down?"
• It kept adjusting until it found the absolute lowest energy state (the most stable arrangement).

The Big Discovery: Convergence

Here is the most exciting part of the paper.

When they compared the Fibonacci Blueprint (the geometric method) with the Computer Sculptor (the minimization method), they found something beautiful:

1. For small numbers of vortices: The "Geodesic Dome" (soccer ball style) was the most efficient.
2. For large numbers of vortices: The Fibonacci Lattice became almost identical to the perfect Computer Sculptor solution.

As the number of vortices grew huge (approaching infinity), both methods converged to the exact same efficiency as the perfect flat-grid city. Even though the surface is curved, if you have enough tiny whirlpools, the sphere looks flat locally, and the "Golden Ratio" spiral is the best way to arrange them.

Why Does This Matter?

This isn't just about math puzzles.

• Real Experiments: Scientists are currently building "bubble traps" for atoms (even on the International Space Station!). They want to see these vortex patterns in real life.
• Detection: In these atomic bubbles, the whirlpools are hard to see because the atoms don't disappear completely inside them (unlike in water). The authors suggest ways to "quench" the system (change the temperature suddenly) to make the whirlpools visible as dark spots.
• New Physics: It shows us how nature balances geometry and physics. Even when the rules of the universe (curved space) prevent perfection, nature finds a "good enough" solution that mimics the perfect flat world when things get crowded.

The Takeaway

Imagine trying to arrange a crowd of people on a giant beach ball so everyone has the same amount of personal space.

• If there are only 20 people, you can arrange them perfectly (like the vertices of a soccer ball).
• If there are 1,000 people, you can't make a perfect grid.
• But, if you arrange them in a Golden Spiral (Fibonacci), they end up spaced out almost perfectly, just as efficiently as if they were standing on a flat floor.

This paper proves that the Fibonacci spiral is the secret code nature uses to organize quantum whirlpools on a sphere, bridging the gap between curved geometry and the perfect order of flat space.

Technical Summary

1. Problem Statement

The paper addresses the challenge of characterizing vortex lattice structures in atomic Fermi superfluids confined to a spherical surface under an effective magnetic monopole field.

• Geometric Constraint: Unlike planar geometries where Abrikosov vortex lattices form perfect periodic structures, a spherical surface forbids perfect lattices with more than 20 vertices due to the non-existence of regular polyhedra with $>20$ vertices.
• Physical Context: While vortex lattices in Bose-Einstein condensates (BECs) on spheres have been studied, the behavior of Fermi superfluids (governed by BCS theory) in this curved geometry remains less explored.
• Core Question: How do vortex lattices arrange themselves on a sphere to minimize free energy, and what are their structural properties as the number of vortices ($N_v$) increases?

2. Methodology

The authors employ the Ginzburg-Landau (GL) theory adapted for a spherical surface with a central magnetic monopole. The radial flux from the monopole acts as the perpendicular magnetic field equivalent to the planar case.

Theoretical Framework

• Linearized GL Equation: The system is described by the eigenvalue equation $H_{ang}\psi = E\psi$, where the ground states are monopole harmonics ($Y_{\mu, j, m}$).
• Order Parameter Construction: The superfluid order parameter $\psi(\theta, \phi)$ is constructed as a superposition of degenerate monopole harmonics:
  $$\psi(\theta, \phi) = \sum_{m=-N}^{N} C_m Y_{N,N,m}(\theta, \phi)$$
  where $N$ is related to the total magnetic flux quanta, and $N_v = 2N$ is the number of vortices.
• Verification: The authors verify that the zeros of $\psi$ correspond to physical vortices by calculating the gauge-invariant current density ($J_s$) to ensure circulating currents exist around the cores.
• Stability Metric: The Abrikosov parameter ($\beta_A = \langle|\psi|^4\rangle / \langle|\psi|^2\rangle^2$) is used to quantify the uniformity of the vortex distribution. Lower $\beta_A$ indicates a more energetically favorable, uniform lattice.

Two Construction Approaches

The paper utilizes two distinct methods to determine the coefficients $C_m$ and locate the vortices:

1. Geometric Construction (Scaffold Method):

   • The zeros of the wavefunction are prescribed to lie on specific lattice scaffolds.
   • Scaffolds used:
     • Random Lattice: Points distributed uniformly at random.
     • Geodesic-Dome Lattice: Based on icosahedral symmetry (Class I, II, III), subdividing icosahedral faces. These allow for specific "magic numbers" of vertices but introduce topological defects (disclinations).
     • Fibonacci Lattice: An analytical construction using the golden ratio to place $N_v$ points quasi-uniformly on the sphere for any integer $N_v$.
   • The coefficients $C_m$ are solved by setting $\psi(\theta_n, \phi_n) = 0$ at the scaffold points.

2. Numerical Minimization:

   • The coefficients $C_m$ are optimized numerically to minimize the Abrikosov parameter $\beta_A$ directly, without pre-defining vortex positions.
   • Algorithms used: L-BFGS-B and Trust-Region constrained optimization (SciPy).
   • Multiple random initial guesses are used to avoid local minima.

3. Key Contributions

• Extension of Abrikosov Theory: Successfully extends the theory of Abrikosov lattices from 2D planes to spherical geometries for Fermi superfluids.
• Validation of Vortex Nature: Unlike previous minimization studies that only found zeros, this work explicitly verifies that the zeros correspond to physical vortices with circulating gauge-invariant currents.
• Comparative Analysis: Provides a rigorous comparison between geometric approximations (Fibonacci/Geodesic) and fully numerical energy minimization.
• Asymptotic Behavior: Demonstrates that despite the curvature-induced defects, the system converges to the planar triangular lattice limit as $N_v \to \infty$.

4. Key Results

• Small $N_v$ ($<20$):
  • Perfect lattices are possible.
  • Geodesic-dome constructions yield the lowest $\beta_A$ because they align with the natural symmetry of the sphere (e.g., the icosahedron for $N_v=12$).
• Large $N_v$ ($>20$):
  • Geodesic-dome lattices suffer from non-uniformity around 5-fold disclinations, causing $\beta_A$ to increase rapidly for non-ideal vertex counts.
  • Fibonacci lattices provide a highly uniform distribution. Their $\beta_A$ values are comparable to the numerical minimization solutions and significantly better than random or imperfect geodesic lattices.
  • Numerical Minimization: Produces the lowest possible $\beta_A$ for a given $N_v$. As $N_v$ increases, the minimization solution visually and structurally converges to the Fibonacci lattice configuration.
• Asymptotic Limit ($N_v \to \infty$):
  • Both the Fibonacci lattice and the numerical minimization solutions extrapolate to $\beta_A \approx 1.16$.
  • This value matches the Abrikosov parameter for the triangular vortex lattice on a 2D plane, confirming that local flatness dominates at high vortex densities.
• Defect Structure:
  • Contrary to predictions from elasticity theory (Thomson problem) which suggest 5-7 dislocation pairs to recover planar energy, the Fermi superfluid vortex lattices in this study do not show clear 5-7 dislocation pairs. The system appears quasi-crystalline without an obvious network of such dislocations.

5. Significance and Implications

• Experimental Relevance: The findings are directly applicable to recent experiments with ultracold atomic bubbles (e.g., on the International Space Station) and phase-separated shell structures on Earth.
• Detection Strategies: The paper discusses that while vortices in Fermi superfluids are harder to detect via density dips than in BECs (due to unpaired fermions), techniques like interaction quenches to the BEC regime could make them visible.
• Theoretical Insight: The work highlights the interplay between geometry, topology, and many-body interactions. It suggests that while curvature forces topological defects (12 five-fold disclinations are required by Euler's theorem), the nonlinear interactions in the GL free energy lead to a different structural organization than purely repulsive charge systems (Thomson problem).
• Practical Utility: The Fibonacci lattice is identified as a computationally efficient and analytically tractable approximation for spherical vortex lattices, offering a reliable alternative to heavy numerical minimization for large systems.

In conclusion, the paper establishes that while perfect spherical lattices are impossible, approximate vortex lattices formed by Fermi superfluids naturally evolve toward the planar triangular lattice structure as the system size increases, with the Fibonacci lattice serving as an excellent analytical proxy for these complex curved-space configurations.