🌀 nonlinear sciences

Thermal effects on density-modulated phases in a dipolar Bose--Einstein condensate

This study demonstrates that finite-temperature corrections, modeled via a temperature-dependent extended Gross-Pitaevskii equation, significantly reshape the structural boundaries and density contrast evolution of density-modulated phases in trapped dipolar Bose-Einstein condensates, shifting phase transitions and smoothing out the sharp zero-temperature contrasts.

Original authors: Changjian Yu, Jinbin Li, Kui-Tian Xi

Published 2026-06-26

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Changjian Yu, Jinbin Li, Kui-Tian Xi

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a pot of very cold, special liquid made of atoms that act like tiny magnets. In physics, we call this a dipolar Bose-Einstein condensate (BEC). Usually, when scientists study these liquids, they look at them at absolute zero temperature (the coldest possible state). But in this paper, the researchers asked: "What happens if we warm this liquid up just a tiny bit?"

Here is a simple breakdown of what they found, using everyday analogies:

1. The Tug-of-War in the Pot

Inside this atomic liquid, there are three forces fighting for control:

The "Bump" Force: Atoms don't like to be in the same spot, so they push each other away (contact repulsion).
The "Magnet" Force: Because the atoms are magnetic, they attract or repel each other depending on their orientation (dipole-dipole interaction).
The "Jitter" Force: Even at very low temperatures, atoms wiggle and fluctuate (quantum and thermal fluctuations).

When these forces balance perfectly, the liquid doesn't just sit there as a smooth blob. Instead, it organizes itself into cool patterns, like honeycombs, maze-like stripes, or pumpkin shapes.

2. The "Warmth" Changes the Map

The researchers used a computer model to simulate what happens when you add a little bit of heat (finite temperature) to this system.

Think of the different patterns (honeycombs vs. stripes vs. smooth blobs) as different "neighborhoods" on a map.

At Absolute Zero: The borders between these neighborhoods are sharp and clear. If you change the conditions just a tiny bit, you instantly jump from one pattern to another.
With a Little Heat: The "heat" acts like a foggy lens. It doesn't just erase the patterns; it shifts the borders.
- Some patterns (like the connected honeycomb mazes) can now survive in conditions where they previously would have collapsed.
- The researchers found that to get these cool patterns, you now need to tweak the "stickiness" of the atoms (scattering length) slightly differently than you would at zero temperature.

The Analogy: Imagine a sandcastle. At zero temperature, the castle is rigid; a small wave either leaves it alone or destroys it instantly. At a slightly warmer temperature, the sand is a bit more fluid. The castle might shift its shape or require a slightly different pile of sand to hold its form, but it doesn't necessarily fall apart immediately. The "rules" for building the castle have changed.

3. The "Smooth Transition" vs. The "Snap"

One of the most interesting findings is how the liquid changes from a smooth blob to a patterned maze.

At Zero Temperature: It's like a light switch. You flip it, and the pattern appears instantly. The contrast between the high and low points of the liquid density changes sharply.
With Heat: It's like a dimmer switch. As you change the conditions, the pattern fades in or out gradually. The transition is smooth, not sudden. The researchers call this a "crossover-like evolution."

4. How "Flowy" is the Liquid? (Superfluidity)

Some of these patterned liquids are "supersolids," meaning they are solid-like (having a pattern) but can also flow without friction (superfluid).

The researchers used a mathematical tool called Leggett's bound to estimate how easily the liquid could flow.
The Result: If the liquid is a smooth blob, it flows easily (like water in a pipe). If it forms a maze or a honeycomb, the "roads" for the atoms to flow through get narrow and winding.
The Analogy: Imagine a highway. A smooth BEC is a wide, open highway where traffic moves fast. A patterned state is a city with narrow, winding streets. The more complex the pattern (like isolated droplets), the harder it is for the "traffic" (superfluid flow) to get through. The researchers calculated that for some patterns, the ability to flow drops by nearly 100 times compared to the smooth state.

Summary

This paper doesn't invent a new machine or cure a disease. Instead, it refines our understanding of how heat affects the shape of these exotic atomic liquids.

The main takeaway is that heat isn't just a destructive force that ruins these delicate patterns. Instead, a little bit of heat acts as a new control knob. It can shift the boundaries of where these patterns appear and make the transitions between them smoother, offering scientists a new way to tune and study these quantum materials.

Technical Summary: Thermal Effects on Density-Modulated Phases in a Dipolar Bose–Einstein Condensate

Problem Statement
Dipolar Bose–Einstein condensates (dBECs) exhibit complex spatially modulated structures, such as droplets, honeycomb networks, and labyrinthine patterns, arising from the competition between contact repulsion, long-range anisotropic dipole–dipole interactions, and beyond-mean-field quantum fluctuations. While zero-temperature structural maps have been established, the role of finite temperature remains less understood. Thermal fluctuations are typically associated with decoherence and depletion, yet recent work suggests they may also modify the stability of modulated states. This paper addresses the specific question of how finite-temperature corrections modify stationary density patterns and their structural boundaries within a trapped dipolar condensate, moving beyond the assumption that temperature solely acts as a source of instability.

Methodology
The authors investigate a single-component dipolar BEC of $^{162}$ Dy atoms confined in a cylindrically symmetric harmonic trap. The study employs a temperature-dependent extended Gross–Pitaevskii equation (TeGPE) within the local-density approximation (LDA). The framework incorporates two key fluctuation corrections:

Quantum Fluctuations: Modeled via the Lee–Huang–Yang (LHY) correction, which provides a repulsive contribution stabilizing the gas against mean-field collapse.
Thermal Fluctuations: Included via a temperature-dependent term derived from the thermal Bogoliubov contribution. This term requires an infrared cutoff ( $k_c \simeq \pi/2\xi$ , where $\xi$ is the local healing length) to handle cases where the dipolar interaction strength $\epsilon_{dd} > 1$ .

The stationary states are obtained via imaginary-time propagation on a Cartesian grid using Fourier pseudospectral discretization. The system is normalized to a fixed condensate atom number $N$ . The authors analyze the resulting density morphologies across a parameter space of $s$ -wave scattering length ( $a_s$ ) and atom number ( $N$ ) at various temperatures ( $T = 0, 50, 100, 150, 200$ nK).

Key Contributions and Results

Reshaping of Structural Boundaries: The finite-temperature correction does not uniformly shift all phase boundaries in the same direction. Instead, it reshapes the structural map in a state-dependent manner. Specifically, the thermal term shifts certain boundaries toward larger scattering lengths. For instance, at $N = 10^4$ and $T = 100$ nK, the boundary between density-modulated and BEC-like profiles shifts from $a_s \approx 87.5 a_0$ (at $T=0$ ) to higher values, with the modulated region extending to $\approx 88.87 a_0$ . Conversely, the critical atom number required for the transition between connected density-modulated and labyrinthine states increases with temperature.
Density Contrast Evolution: The paper analyzes the density contrast $C = (\rho_{max} - \rho_{min})/(\rho_{max} + \rho_{min})$ $C = (ρ_{ma x} - ρ_{min}) / (ρ_{ma x} + ρ_{min})$ as a function of $a_s$ $a_{s}$ and $N$ $N$ .
- At $T=0$ , the contrast exhibits a sharp, discontinuous-like drop as the system transitions from a modulated state to a uniform BEC.
- At finite temperatures (e.g., $T=200$ nK), this transition becomes significantly smoother, occurring over a broader range of scattering lengths. The authors characterize this finite-temperature behavior as a "smooth crossover-like evolution" within the finite trapped system, rather than a sharp thermodynamic phase transition.
Superfluid Response Estimation: Using Leggett's upper bound, the authors estimate the geometric limit on the superfluid fraction for representative density profiles.
- A nearly uniform BEC-like state yields a high bound ( $f_S \leq 0.942$ ).
- A connected density-modulated state shows a reduced bound ( $f_S \leq 0.346$ ) due to transport constraints imposed by low-density channels.
- An isolated droplet-array state exhibits a drastically suppressed bound ( $f_S \leq 0.0146$ ), indicating that the density profile strongly inhibits flow along the chosen direction.
- The authors emphasize that these values serve as geometric upper limits on superfluid response, not standalone proofs of supersolidity.

Significance and Claims
The paper concludes that finite-temperature corrections play a constructive role in the pattern formation of dipolar quantum gases, acting as an additional control parameter rather than merely a source of decoherence. The thermal term modifies the local equation of state, altering the balance between contact, dipolar, and fluctuation energies, which in turn changes the parameter regions where connected density-modulated states are stable.

The authors maintain a modest stance regarding their findings:

The structural map presented is a classification of stationary density morphologies within the TeGPE-LDA model, not a complete thermodynamic phase diagram.
The "smooth crossover" description applies to the finite trapped system studied; further thermodynamic evidence would be required to confirm the nature of the transition in the thermodynamic limit.
The results are most reliable for dilute, weakly depleted gases where the LDA treatment of fluctuations is controlled. The authors acknowledge that near structural boundaries, in dilute tails, or at higher depletion, coupling to the thermal cloud and higher-order correlations may become significant.

The study suggests that experimental verification could be pursued through in-situ density imaging, observing the temperature dependence of density contrast, and analyzing static structure factors or time-of-flight images.