Long-Range, Dipolar, and Molecular Quantum Gases

Long-range quantum gases are controlled by more than an interaction tail. Dipole orientation, short-distance regularization, transverse confinement, molecular internal states, shielding, loss, and the preparation protocol determine the effective Hamiltonian. A roton minimum, density modulation, self-bound cloud, or condensate fraction is then an observable-specific statement rather than an automatic phase identification.

Required background. Ultracold platforms, scales, and traps supplies density, confinement, LDA, loss, and observation scales. Beyond-Bogoliubov dilute-gas expansion supplies the fluctuation expansion and its gas-parameter limit.

Helpful background. Self-bound quantum droplets develops fluctuation stabilization and finite-droplet diagnostics.

Dipolar interaction and momentum-space stability

For identical polarized dipoles aligned along unit vector $\hat{\mathbf e}$,

$$V_{\mathrm{dd}}(\mathbf r) =\frac{C_{\mathrm{dd}}}{4\pi} \frac{1-3(\hat{\mathbf e}\!\cdot\!\hat{\mathbf r})^2}{r^3}.$$

The singular expression is a long-distance potential; its distributional contact part and short-range boundary condition must be matched consistently. Define

$$g=\frac{4\pi\hbar^2a_s}{m}, \qquad g_{\mathrm{dd}}=\frac{C_{\mathrm{dd}}}{3}, \qquad \epsilon_{\mathrm{dd}}=\frac{g_{\mathrm{dd}}}{g} =\frac{a_{\mathrm{dd}}}{a_s}, \qquad a_{\mathrm{dd}}=\frac{mC_{\mathrm{dd}}}{12\pi\hbar^2}.$$

In a homogeneous three-dimensional convention with the contact piece assigned as above,

$$\widetilde V(\mathbf k)= g+g_{\mathrm{dd}} \left[3(\hat{\mathbf e}\!\cdot\!\hat{\mathbf k})^2-1\right].$$

The Bogoliubov spectrum is

$$E_{\mathbf k}^2 =\varepsilon_k \left[\varepsilon_k+2n\widetilde V(\mathbf k)\right], \qquad \varepsilon_k=\frac{\hbar^2k^2}{2m}.$$

At long wavelength the most attractive direction has $\widetilde V_{\min}=g(1-\epsilon_{\mathrm{dd}})$, so the uniform mean-field gas requires $\epsilon_{\mathrm{dd}}<1$. This criterion is not a universal collapse boundary: confinement changes the momentum dependence, fluctuations shift stability, losses alter the sampled density, and a finite trap can support metastable states. Lahaye et al. 2009, §§2–5 review these Fourier, geometry, and stability conventions.

In quasi-two-dimensional confinement, integrating over a transverse oscillator wave function gives

$$\widetilde V_{\mathrm{2D}}(k_\parallel,\varphi) =g_{\mathrm{2D}} +g_{\mathrm{dd,2D}}\, \mathcal F(k_\parallel a_\perp,\varphi;\hat{\mathbf e}),$$

where the nonconstant form factor $\mathcal F$ can soften a finite wave number and create a roton minimum. The roton condition is therefore found from the complete confined dispersion, not from $\epsilon_{\mathrm{dd}}=1$ alone.

The validity diagram separates this Hamiltonian construction from the phase claim. Inspect the loss branch: shielding that lengthens lifetime changes the collision problem and must be included in the same effective interaction used to interpret the state.

Validity map for dipolar and molecular gases. A calibrated anisotropic interaction can support a roton or instability calculation, but phase identification additionally requires equilibrium, coherence, finite-size, loss, and observation controls. Original schematic, not to scale; platform evidence is bounded through 10 August 2026.

Fluctuation stabilization and its boundary

For a dilute dipolar Bose gas, a local-density Lee–Huang–Yang correction is often written

$$\delta\mu_{\mathrm{LHY}} =\frac{32}{3\sqrt\pi}\, gn\sqrt{na_s^3}\, \mathcal Q_5(\epsilon_{\mathrm{dd}}),$$

where $\mathcal Q_5$ is an angular average of the Bogoliubov branches in the chosen stable continuation. Lima and Pelster 2011, main text derive the dipolar fluctuation correction in the homogeneous dilute regime. The control parameter is $\sqrt{na_s^3}\ll1$, together with a separation between the healing, confinement, range, and droplet-surface scales. Near a mean-field instability, simply retaining the real part of a formally complex local expression is a modeling prescription, not a derivation. Nonlocal fluctuation terms, finite-range physics, three-body correlations, and density-dependent losses can compete.

A self-bound droplet requires negative bulk pressure balanced by repulsive corrections and positive surface energy; the fluctuation-stabilized mixture mechanism of Petrov 2015, main text provides a clean example of this balance outside the dipolar setting. A density-modulated supersolid requires both translational order and global phase coherence. A single in-situ modulation can also arise from trap modes, interference, fragmentation, or imaging transfer; a long-lived compact cloud can be dynamically arrested rather than an equilibrium droplet.

Molecules, internal states, and present evidence

Polar molecules add rotational and hyperfine states. Microwave or electric dressing can generate a tunable effective dipole while suppressing short-range loss, but the effective interaction depends on detunings, Rabi frequencies, polarization impurities, avoided crossings, and collision energy. A scalar $C_{\mathrm{dd}}$ is insufficient when several dressed channels are nearby.

Shi et al. 2026, main text and Methods reported Bose–Einstein condensation of ground-state sodium–rubidium molecules under dual microwave shielding, together with tunable dipolar behavior and a gas-to-droplet observation. The experiment establishes a platform and bounded observables; it does not by itself certify every equilibrium phase assignment in other species, dimensionalities, or dressing regimes. Condensate number, lifetime, expansion dynamics, internal-state purity, and interaction calibration remain part of the claim.

The source and platform assessment here is current through 10 August 2026. Later lifetime records, molecular phase evidence, revised collision models, and corrections belong in the Quantum Matter and Emergence Research synthesis. The canonical cold-atom and synthetic-matter claim test matrix places interaction, scale hierarchy, preparation, resolution, loss, discrepancy, and evidence ceiling together. A reproducible verification workflow propagates these quantities through a selected stability or observable calculation.

Exercise

Angular stability of a homogeneous dipolar gas. For $\widetilde V(\theta)=g[1+\epsilon_{\mathrm{dd}}(3\cos^2\theta-1)]$ with $g>0$, find its minimum over $\theta$ and the long-wavelength mean-field stability condition. What changes if the polarization is tilted in a quasi-two-dimensional trap?

Solution

Since $3\cos^2\theta-1$ ranges from $-1$ at $\theta=\pi/2$ to $2$ at $\theta=0$, for $\epsilon_{\mathrm{dd}}>0$ the minimum is

$$\widetilde V_{\min}=g(1-\epsilon_{\mathrm{dd}}).$$

The sound-speed squared in that direction is proportional to $n\widetilde V_{\min}/m$, so the homogeneous long-wavelength mean-field condition is $\epsilon_{\mathrm{dd}}<1$. With transverse confinement and tilt, the effective interaction contains a momentum-dependent form factor and in-plane angular dependence. The first unstable mode can occur at finite $k_\parallel$ and depends on $a_\perp$, tilt, contact coupling, and density; the three-dimensional $\epsilon_{\mathrm{dd}}<1$ criterion is no longer sufficient.

References

• Lahaye, T., Menotti, C., Santos, L., Lewenstein, M., and Pfau, T. (2009). “The physics of dipolar bosonic quantum gases.” Reports on Progress in Physics 72, 126401. doi:10.1088/0034-4885/72/12/126401.
• Lima, A. R. P., and Pelster, A. (2011). “Quantum fluctuations in dipolar Bose gases.” Physical Review A 84, 041604(R). doi:10.1103/PhysRevA.84.041604.
• Petrov, D. S. (2015). “Quantum mechanical stabilization of a collapsing Bose–Bose mixture.” Physical Review Letters 115, 155302. doi:10.1103/PhysRevLett.115.155302.
• Shi, Z., Huang, Z., Deng, F., Jin, W.-J., Yi, S., Shi, T., and Wang, D. (2026). “Bose–Einstein condensate of ultracold sodium–rubidium molecules with tunable dipolar interactions.” Nature Physics, published 9 July 2026. doi:10.1038/s41567-026-03362-9.