Phase diagram of rotating Bose-Einstein condensates trapped in power-law and hard-wall potentials

This paper investigates the rotational phase diagram of quasi-two-dimensional Bose-Einstein condensates in power-law and hard-wall traps, revealing that while weak interactions cause discontinuous transitions between multiply-quantized vortex states and stronger interactions lead to continuous transitions to mixed states, the two confinement types exhibit distinct qualitative behaviors regarding central density stability and scaling properties.

The Gist

Imagine a giant, invisible dance floor made of super-cooled atoms. This is a Bose-Einstein Condensate (BEC), a state of matter where thousands of atoms act as a single, giant "super-atom" moving in perfect unison.

Now, imagine you start spinning this dance floor. What happens? Do the atoms just spin with the floor, or do they form little whirlpools (vortices) to handle the spin?

This paper is a theoretical study by G. M. Kavoulakis that maps out exactly how these atomic dancers behave when spun, depending on the shape of the container holding them. The author compares two very different types of containers:

1. The "Soft Bowl" (Power-Law Trap): Imagine a bowl that gets steeper the further out you go, but it's smooth and curved everywhere, like a giant, soft funnel.
2. The "Hard Box" (Hard-Wall Trap): Imagine a container with perfectly vertical, impenetrable walls, like a circular room with a hard floor and a hard ceiling.

Here is the breakdown of the findings, translated into everyday language:

1. The Two Ways the Dance Can Go Wrong

When you spin the container, the atoms want to carry the "spin" (angular momentum). They do this by forming vortices. The paper finds two distinct ways the atoms react as you spin faster or make them interact more strongly:

• The "Snap" (Discontinuous Transition): If the atoms barely talk to each other (weak interaction), they are stubborn. As you speed up the spin, they suddenly jump from having one giant whirlpool to having two, then three, then four. It's like a light switch flipping on and off. The system jumps instantly from one state to the next.
• The "Blend" (Continuous Transition): If the atoms are chatty (strong interaction), they get tired of holding one giant whirlpool. Instead of snapping to a new number, they slowly break that giant whirlpool apart into a "necklace" of many tiny, single whirlpools. This happens gradually, like melting ice cream.

2. The Big Surprise: The Shape of the Container Matters

The most exciting part of this paper is that the Soft Bowl and the Hard Box behave in completely opposite ways when the atoms start breaking apart.

In the Soft Bowl (Power-Law Trap):

• The "Donut" Effect: As the spin gets faster, the atoms are pushed outward by the centrifugal force. Because the bowl is soft and curved, the atoms flee the center completely.
• The Result: The center of the dance floor becomes empty. The atoms form a ring, like a donut or a tire. If you look at the center, there is no density there at all. The giant whirlpool splits, but the pieces stay away from the middle.

In the Hard Box (Hard-Wall Trap):

• The "Crowded Center" Effect: This is the counter-intuitive part. Even though the atoms are being spun, the hard walls force them to behave differently. When the giant whirlpool breaks apart, it doesn't just push everything to the edge.
• The Result: The atoms actually fill the center. One of the tiny whirlpools stays right in the middle of the dance floor, while the others orbit around it. The center remains crowded with atoms.

The Analogy:
Imagine spinning a bucket of water.

• In the Soft Bowl, the water climbs the sides so high that the bottom is completely dry.
• In the Hard Box, the water is forced against the walls, but the physics of the "hard wall" forces a little bit of water to stay right in the middle, creating a stable spot in the center.

3. Why Does This Happen?

The author explains this using a battle between two forces:

1. The "Single Particle" Force: This force likes to keep the atoms in a specific, organized pattern (like a single giant whirlpool).
2. The "Interaction" Force: This force (the atoms bumping into each other) wants the atoms to spread out evenly to avoid crowding.

• In the Soft Bowl: The shape of the bowl changes as you spin faster (it becomes a "Mexican Hat" shape). The lowest energy spot moves away from the center. So, the atoms naturally leave the center empty to find that new low-energy spot.
• In the Hard Box: The walls are rigid. The "lowest energy spot" is always at the edge, but the math of the hard wall forces the atoms to keep one piece of the puzzle in the center to satisfy the rules of the container.

4. Why Should We Care?

This isn't just math for math's sake. Scientists can actually build these traps in real labs using lasers and magnets.

• Testing the Theory: By spinning a cloud of atoms and taking a picture, scientists can see if the center is empty (Soft Bowl) or full (Hard Box).
• Measuring the "Stiffness": The way the atoms split tells scientists exactly how "stiff" or "soft" their trap is. It's like diagnosing a patient by how they react to a specific medicine.

Summary

This paper is a map for how super-cold atoms behave when spun. It tells us that the shape of the container is the boss.

• If you use a soft, curved trap, the atoms will run away from the center, leaving a hole in the middle.
• If you use a hard, box-like trap, the atoms will stay in the center, even while spinning.

It's a beautiful example of how the rules of the universe (quantum mechanics) play out differently depending on the "room" you put them in.

Technical Summary

1. Problem Statement

The paper investigates the rotational phase diagram of a quasi-two-dimensional, weakly-interacting Bose-Einstein condensate (BEC). While the rotational properties of BECs in harmonic traps are well-understood, this study focuses on two distinct confinement geometries that allow for rotation frequencies ($\Omega$) exceeding the trap frequency:

1. Power-law traps: Anharmonic potentials where the confining potential rises faster than quadratically ($V \propto \rho^{2p}$ with $p > 2$).
2. Hard-wall traps: Potentials that are zero within a radius $R_0$ and infinite outside, mimicking a rotating bucket.

The central problem is to determine how the interplay between the rotation frequency ($\Omega$) and the interaction strength ($g$) governs the stability of vortex states. Specifically, the study seeks to understand the mechanisms driving transitions between multiply-quantized vortex states and "mixed" states (combinations of singly and multiply-quantized vortices), and to identify qualitative differences between power-law and hard-wall confinement.

2. Methodology

The author employs a theoretical framework based on the Gross-Pitaevskii mean-field theory adapted for rotating frames.

• Hamiltonian and Model: The system is described by a Hamiltonian including kinetic energy, the trapping potential $V(r)$, and contact interactions ($U_0$). The condensate is treated as quasi-2D, with motion along the rotation axis ($z$) frozen.
• Basis States:
  • For Power-law traps, the author uses the Lowest-Landau-Level (LLL) states of a harmonic potential as a basis, treating the anharmonic term ($\lambda$) and interactions ($g$) perturbatively.
  • For Hard-wall traps, the basis consists of Bessel function eigenstates ($J_m$) satisfying the boundary condition $\psi(R)=0$.
• Energy Minimization:
  • Discontinuous Transitions: For weak interactions, the system is assumed to occupy a single angular momentum state $m$. The transition frequency $\Omega_{m, m+1}$ is calculated by equating the energies of adjacent states in the rotating frame.
  • Continuous Transitions: As interaction strength increases, multiply-quantized states ($m_0$) become unstable. The author analyzes the stability of a state $\psi_{m_0}$ against a "mixed" triplet state: $\Psi = c_{m_0-q}\psi_{m_0-q} + c_{m_0}\psi_{m_0} + c_{m_0+q}\psi_{m_0+q}$.
  • Stability Analysis: A Hessian matrix is constructed from the second derivatives of the energy with respect to the coefficients $c_m$. The onset of instability occurs when the lowest eigenvalue of this matrix becomes negative, determining the critical coupling $g$ for a given $\Omega$.

3. Key Contributions

• Qualitative Distinction in Instability Mechanisms: The paper establishes a fundamental difference in how BECs destabilize in power-law vs. hard-wall traps.
  • Hard-wall: The leading instability always involves the state with zero angular momentum ($m=0$), resulting in a mixed state with nonzero density at the trap center.
  • Power-law: For sufficiently high rotation, the effective potential develops a "Mexican-hat" shape. The leading instability eventually excludes the $m=0$ state, causing the density at the trap center to vanish (forming a hole).
• Universal Scaling Properties:
  • Power-law: The phase diagram is invariant under a simultaneous rescaling of parameters $(1-\Omega)$, $\lambda$, and $g$.
  • Hard-wall: The phase diagram depends only on the combination $\Omega R^2$ and $g$, reflecting the intrinsic length scale $R$.
• Mapping of Instability Channels: The study systematically maps the parameter $q$ (the angular momentum difference in the mixed state) as a function of the initial vortex charge $m_0$ and trap parameters ($p, \lambda, R$).

4. Key Results

A. Weak Interaction Regime (Discontinuous Transitions)

• In both trap types, increasing $\Omega$ causes discontinuous jumps between multiply-quantized vortex states ($m \to m+1$).
• Power-law: The spacing between transition lines ($\Delta \Omega$) increases with $m$ for $p > 2$. The transition lines in the $\Omega-g$ plane have a negative slope.
• Hard-wall: The spacing between transition lines becomes constant for large $m$. The transition lines have a positive slope.

B. Strong Interaction Regime (Continuous Transitions)

• As $g$ increases, multiply-quantized vortices split into mixed states to lower interaction energy.
• Power-law Behavior:
  • For small $m_0$, the instability involves $q=m_0$ (splitting into $2m_0$ singly quantized vortices, with one at the center).
  • For large $m_0$, $q$ becomes less than $m_0$ (e.g., $q \approx 4$ for $p=4$). The central density vanishes because the effective potential minimum moves radially outward.
• Hard-wall Behavior:
  • The instability always involves $q=m_0$. The mixed state is a superposition of $(\psi_0, \psi_{m_0}, \psi_{2m_0})$.
  • Because $\psi_0$ (the $m=0$ state) is included, the density at the center remains nonzero. The system pays a kinetic energy cost to maintain a uniform density distribution favored by repulsive interactions.

C. Phase Diagram Structure

• Figure 1 & 3 (Power-law): Show phase boundaries where straight lines represent discontinuous transitions and curves represent continuous transitions. The curves shift upward as $m_0$ increases.
• Figure 6 (Hard-wall): Shows that phase diagrams for different radii $R$ are identical when the x-axis is rescaled by $R^2$.
• Density Profiles: Visualizations (Figs. 2, 4, 5, 7) confirm that power-law traps develop a central hole at high rotation, whereas hard-wall traps maintain central density even in the mixed state.

5. Significance and Experimental Relevance

• Experimental Feasibility: The predicted phenomena are within reach of current cold-atom experiments. Techniques such as Feshbach resonances allow tuning of $g$, while optical potentials can create hard-wall or anharmonic traps.
• Diagnostic Tool: The presence or absence of a central density hole serves as a direct experimental signature to distinguish between power-law and hard-wall confinement.
• Vortex Splitting: The paper predicts specific "necklace" configurations of singly quantized vortices resulting from the splitting of multiply quantized vortices, which can be observed via absorption imaging or time-of-flight expansion.
• Theoretical Insight: The work clarifies the competition between the curvature of the single-particle spectrum (favoring high-$m$ states) and repulsive interactions (favoring spatial redistribution). It demonstrates that the geometry of the trap fundamentally alters the mechanism by which a superfluid accommodates angular momentum.

In conclusion, this paper provides a comprehensive theoretical framework for understanding rotating BECs in non-harmonic traps, highlighting that the confinement geometry dictates not just the energy scales, but the fundamental topology of the vortex states and the nature of phase transitions.