Dipole Condensates: Quantum Fluids & Lattice Systems

• Dipole condensates are quantum-condensed states where dipolar interactions dictate the order parameter and collective dynamics.
• They span both atomic/molecular systems with anisotropic long-range forces and lattice models featuring particle–hole dipolar condensation.
• Field-theoretic and mean-field frameworks reveal self-bound droplets, supersolids, and multipolar phase transitions in dipolar systems.

Searching arXiv for recent and relevant papers on dipole condensates to ground the article in current literature. arxiv_search query: "dipole condensates Bose-Einstein condensate dipolar droplets supersolids tilted Bose-Hubbard chains" Dipole condensates are quantum-condensed states in which dipolar degrees of freedom are central to the order parameter, collective modes, or dominant interactions. In the literature, the term is used in at least two technically distinct senses. In dipolar Bose–Einstein condensates, it refers to atomic or molecular condensates governed by long-range, anisotropic magnetic or electric dipole–dipole interactions, often supplemented by beyond-mean-field stabilization terms or conservative three-body interactions (Schmidt et al., 2021). In strongly constrained lattice systems, it denotes condensation of neutral particle–hole dipoles, characterized by two-body off-diagonal order such as $\langle b_i^\dagger b_{i+\mathbf s}\rangle\neq 0$ even when $\langle b_i\rangle=0$ (Xu et al., 2023). Both usages emphasize that dipolar objects, rather than conventional scalar condensate order alone, control the macroscopic physics.

1. Terminological scope and definitions

A spinor dipolar condensate is a spinor Bose–Einstein condensate in which the long-range, anisotropic magnetic dipole–dipole interaction plays a role on the same footing as the $s$-wave collisions. For spin-1 gases, the order parameter is $\psi_\alpha(\mathbf r,t)$ with $\alpha=+1,0,-1$, and the magnetic dipole–dipole interaction enters through

$$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$

with local spin density $$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$ (Yasunaga et al., 2010).

A different definition appears in fracton-related and tilted-lattice settings. There a dipole condensate is defined by off-diagonal long-range order in a two-body correlator,

$$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$

so that a bound particle–hole pair is the condensing object (Xu et al., 2023). In the one-dimensional dipole-conserving Bose–Hubbard chain, the natural condensed operator is $d_n\equiv b_n b_{n+1}$ or equivalently the nearest-neighbor dipole correlator $\langle b_n^\dagger b_{n+1}\rangle$, while single-boson correlators remain exponentially decaying (Lake et al., 2022).

The coexistence of these usages is not merely terminological. In the first case, the condensate is made of particles carrying permanent or induced dipole moments; in the second, the condensed object is itself a dipole. This suggests that “dipole condensate” is best treated as a family resemblance term spanning dipolar quantum fluids and condensates of dipolar quasiparticles rather than a single microscopic phase.

2. Field-theoretic and mean-field descriptions of dipolar condensates

For dipolar Bose gases, the standard starting point is an extended Gross–Pitaevskii functional. In molecular condensates with strong dipole–dipole interactions,

$\langle b_i\rangle=0$0

with

$\langle b_i\rangle=0$1

The Lee–Huang–Yang term is treated in local-density approximation as

$\langle b_i\rangle=0$2

(Schmidt et al., 2021).

Closely related formulations appear in strongly dipolar atomic condensates. In the generalized nonlocal nonlinear Schrödinger equation of Wächtler and Santos, the beyond-mean-field term is written through

$\langle b_i\rangle=0$3

with

$\langle b_i\rangle=0$4

and the dipolar mean field

$\langle b_i\rangle=0$5

(Wächtler et al., 2016).

Flattened and quasi-two-dimensional geometries admit reduced descriptions with effective interaction parameters. For harmonic confinement along one direction, Baillie and Blakie introduce

$\langle b_i\rangle=0$6

which control the axial condensate profile, roton structure, and mean-field stability (Baillie et al., 2014).

Not all stabilization mechanisms are fluctuation-driven. A separate line of work augments the dipolar Gross–Pitaevskii theory by a conservative three-body interaction $\langle b_i\rangle=0$7, giving

$\langle b_i\rangle=0$8

(Blakie, 2016). A Gaussian-state theory goes further by replacing the usual coherent-state ansatz with a multimode squeezed coherent state and adding an explicit three-body Hamiltonian

$\langle b_i\rangle=0$9

thereby treating squeezing and coherent occupation self-consistently (Wang et al., 2020).

3. Self-bound droplets, supersolids, and interaction-driven liquid behavior

A central result of the dipolar-condensate literature is that mean-field collapse need not be the endpoint of strong dipolar attraction. On the mean-field level a purely dipolar condensate with $s$0 would collapse, but the LHY term provides a repulsive $s$1 pressure which can stabilize a finite-size self-bound droplet (Schmidt et al., 2021). In molecular Bose–Einstein condensates, the critical molecule number obeys a universal curve

$s$2

which diverges as $s$3 and falls rapidly for larger $s$4. The droplet central density saturates as

$s$5

and numerically follows $s$6 with $s$7–$s$8 over a wide range, while the droplet sizes satisfy $s$9 in the large-$\psi_\alpha(\mathbf r,t)$0, saturated regime (Schmidt et al., 2021).

In trapped pancake geometries the same framework yields a phase diagram in $\psi_\alpha(\mathbf r,t)$1 with three main regimes: dilute BEC for $\psi_\alpha(\mathbf r,t)$2, roton-softened gauge-broken states with “blood-cell” shapes once a roton gap closes, and symmetry-broken supersolids and isolated droplet arrays at still larger $\psi_\alpha(\mathbf r,t)$3. The empirical boundaries are

$\psi_\alpha(\mathbf r,t)$4

with constants set by trap geometry. In pancake geometry the roton minimum appears for $\psi_\alpha(\mathbf r,t)$5 when $\psi_\alpha(\mathbf r,t)$6, and density modulation may be quantified by

$\psi_\alpha(\mathbf r,t)$7

The Leggett bound gives

$\psi_\alpha(\mathbf r,t)$8

and numerical solutions yield $\psi_\alpha(\mathbf r,t)$9–$\alpha=+1,0,-1$0 deep in the supersolid regime (Schmidt et al., 2021).

For dysprosium parameters, arrested collapse produces filament-like droplets. Wächtler and Santos identify a minimal droplet atom number $\alpha=+1,0,-1$1 at $\alpha=+1,0,-1$2, and a larger threshold $\alpha=+1,0,-1$3 where the internal energy passes through zero. Numerical solutions give axial half-length $\alpha=+1,0,-1$4m, transverse radius $\alpha=+1,0,-1$5m, aspect ratio $\alpha=+1,0,-1$6, and peak density $\alpha=+1,0,-1$7–$\alpha=+1,0,-1$8 (Wächtler et al., 2016).

The stabilization mechanism is not unique across models. In the trapped theory with conservative three-body repulsion, two distinct fluid states appear: a low-density phase and a high-density droplet phase, separated by a first-order transition terminating at a critical point (Blakie, 2016). In Gaussian-state theory for dysprosium, self-bound states are instead split into an expanding gas, a self-bound gas, and a self-bound liquid, with the self-bound gas–self-bound liquid boundary marked by an abrupt jump in the coherent fraction $\alpha=+1,0,-1$9 (Wang et al., 2020).

Interaction anisotropy can also change droplet dimensionality. In a doubly dipolar condensate with both magnetic and electric dipole moments, varying the angle $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$0 between the two dipoles produces a crossover from quasi-1D cigar droplets elongated along the effective polarization axis to quasi-2D pancake droplets flat in the plane perpendicular to the soft direction. For dysprosium at typical fields this occurs near $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$1 (Mishra et al., 2019). At still stronger dipolarity, a harmonically confined condensate with a central Gaussian hill can host a hollow cylindrical quasi-one-dimensional metastable droplet with ring topology, with a ring branch appearing as the second excited stationary state for $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$2, $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$3, and $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$4 between $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$5 and $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$6 (Adhikari, 14 Jul 2025).

4. Geometry, polarization, and interaction engineering

Dipolar condensates are unusually sensitive to geometry because the kernel $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$7 couples density or spin structure directly to trap anisotropy and polarization orientation. In a quasi-two-dimensional condensate with dipoles tilted in the $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$8 plane by angle $$V_{dd}(\mathbf r-\mathbf r')=\frac{\mu_0(g_e\mu_B)^2}{4\pi}\, \frac{F_i(\mathbf r)F_j(\mathbf r')(\delta_{ij}-3e_ie_j)}{|\mathbf r-\mathbf r'|^3},$$9, the effective 2D momentum-space kernel is

$$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$0

where

$$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$1

The long-wavelength effective coupling becomes

$$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$2

At the magic angle $$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$3, one has $$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$4; for $$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$5, phonon instability sets in, while for $$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$6 just below $$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$7 an anisotropic roton minimum at finite $$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$8 appears (Mishra et al., 2016).

The post-instability dynamics depends on which branch softens. Post-phonon-instability dynamics yields stripe patterns with dislocations and eventual breakup into anisotropic bright solitons, whereas post-roton-instability dynamics yields transient, defect-free stripes followed by local collapse at stripe peaks. Slowly ramping parameters back into the stable roton region can make the stripes persist for much longer, and weak dissipation causes the stripe amplitude slowly to decay while stripe–stripe hopping is observed before the condensate returns to uniformity (Mishra et al., 2016).

Multilayer and disconnected-condensate settings provide additional control knobs. In a multilayer stack formed by a strong optical lattice, effective intra- and interlayer dipole–dipole interaction potentials can be derived analytically; the interlayer term changes the transverse aspect ratio of the ground state in central layers and reduces the excitation energy of local perturbations, affecting the development of a roton minimum (Rosenkranz et al., 2012). In physically disconnected dipolar condensates, a trapped control condensate can engineer the density of a target condensate through inter-condensate dipole–dipole interactions. With two control condensates separated by $$F_i(\mathbf r)=\Psi^\dagger(\mathbf r)\hat F_i\Psi(\mathbf r)$$9, a critical spacing $$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$0 marks where the central minimum at $$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$1 turns into a small local maximum, producing single-peak, three-peak, and four-peak target densities; a coherence measure

$$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$2

diagnoses whether adjacent peaks are phase locked or effectively isolated (Nayak et al., 2023).

The dipolar interaction itself can be induced and enhanced. In an off-resonantly illuminated neutral-atom condensate, azimuthal averaging of the retarded induced interaction gives

$$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$3

which reduces to the familiar $$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$4 form for $$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$5 and crosses over to an oscillatory $$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$6 interaction for $$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$7. The full calculation yields an enhancement by up to thirty-fold at short wavelengths (Graham et al., 2016).

An alternative engineering route exploits geometry rather than optical retardation. For cigar-shaped condensates loaded into a one-dimensional lattice, projection onto Wannier–BEC modes yields intersite dipolar matrix elements

$$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$8

so that the continuum dilute-fluctuation limit becomes a Calogero–Sutherland Hamiltonian with inverse-square interactions (Yu et al., 2011).

5. Spinor dipole condensates, magnetization dynamics, and vortical textures

In spinor condensates, dipolar interactions generate magnetization dynamics with direct analogies to condensed-matter ferromagnets. Starting from the spin-1 Gross–Pitaevskii equations, Yasunaga and Tsubota derive for the local spin density

$$\langle D_{i,\mathbf s}\rangle \equiv \langle b_i^\dagger b_{i+\mathbf s}\rangle \neq 0 \quad (\mathbf s\neq 0), \qquad \langle b_i\rangle=0,$$9

where $d_n\equiv b_n b_{n+1}$0 is a quantum-pressure or spin-diffusion term and

$d_n\equiv b_n b_{n+1}$1

Under the single-mode approximation this reduces to

$d_n\equiv b_n b_{n+1}$2

and linearization about $d_n\equiv b_n b_{n+1}$3 gives

$d_n\equiv b_n b_{n+1}$4

For a spherical condensate $d_n\equiv b_n b_{n+1}$5, so $d_n\equiv b_n b_{n+1}$6, whereas elliptic, pancake, and rod geometries shift the resonance through the shape dependence of $d_n\equiv b_n b_{n+1}$7 (Yasunaga et al., 2010).

When rotation is added, dipolar spinor condensates support an extensive defect taxonomy. Simula and collaborators compute phase diagrams in the plane $d_n\equiv b_n b_{n+1}$8 and identify regions with no defects, spontaneous spin currents in a single polar-core spin vortex $d_n\equiv b_n b_{n+1}$9, a single Mermin–Ho skyrmion $\langle b_n^\dagger b_{n+1}\rangle$0, a staggered skyrmion lattice, a staggered half-quantum vortex lattice, vortex sheets of $\langle b_n^\dagger b_{n+1}\rangle$1 defects, and a mixed lattice of $\langle b_n^\dagger b_{n+1}\rangle$2 and $\langle b_n^\dagger b_{n+1}\rangle$3 vortices, with collapse beyond $\langle b_n^\dagger b_{n+1}\rangle$4 (Simula et al., 2010).

Even without imposed rotation, semiclassical and dual descriptions predict spontaneous textures. A three-component real order parameter $\langle b_n^\dagger b_{n+1}\rangle$5 with magnetization density $\langle b_n^\dagger b_{n+1}\rangle$6 yields ground states including a flare state, a spin-vortex state, and in elongated traps a spin helix state with local azimuthal angle $\langle b_n^\dagger b_{n+1}\rangle$7 and wavevector $\langle b_n^\dagger b_{n+1}\rangle$8 (Huhtamäki et al., 2010). The same texture landscape can be reformulated through a Dirac string gas with negative string tension,

$\langle b_n^\dagger b_{n+1}\rangle$9

which favors closed strings, small curvature, and widely separated residual monopoles. In that language, three-dimensional condensates prefer a meron-like vortex texture, quasi-one-dimensional condensates prefer the axially polarized flare texture, and quasi-two-dimensional condensates exhibit either a meron texture or an in-plane polarized texture (Lian, 2015).

In dominantly dipolar scalar condensates, vortices themselves acquire dipolar signatures. With the extended Gross–Pitaevskii functional including

$\langle b_i\rangle=0$00

Bland and collaborators find that tilting the magnetic field induces magnetostriction and permits vortex nucleation by rotating the polarization direction. For $\langle b_i\rangle=0$01, $\langle b_i\rangle=0$02, and $\langle b_i\rangle=0$03, vortices appear via a dynamical instability at

$\langle b_i\rangle=0$04

for tilt angle $\langle b_i\rangle=0$05. They identify three dynamical regimes: $\langle b_i\rangle=0$06 as a stable Thomas–Fermi branch, $\langle b_i\rangle=0$07 as a “Goldilocks window” for vortex penetration and stripe-lattice formation, and $\langle b_i\rangle=0$08 as a “droplet-like stability” regime where strong head-to-tail attraction suppresses the instability (Bland et al., 2023).

6. Dipole condensation in constrained lattice systems and multipolar generalizations

In tilted Bose–Hubbard systems, the word “dipole” refers to a mobile particle–hole composite rather than to a permanent magnetic or electric moment. For a strongly tilted chain, a Schrieffer–Wolff expansion yields an effective dipole-conserving Hamiltonian

$\langle b_i\rangle=0$09

which commutes with both total boson number and total dipole moment (Lake et al., 2022). Its low-energy field theory is a quantum Lifshitz theory with Gaussian Lagrangian

$\langle b_i\rangle=0$10

Unlike the conventional Bose–Hubbard model, the phase diagram contains no compressible phases and is instead dominated by various types of exotic dipolar condensates (Lake et al., 2022).

The diagnostic of these phases is the mismatch between one-body and two-body coherence. In the dipole-condensed phase,

$\langle b_i\rangle=0$11

while

$\langle b_i\rangle=0$12

and the charge compressibility vanishes. The phase is therefore gapless but incompressible, with power-law order in neutral dipole operators rather than in single-boson operators (Lake et al., 2022).

A broader claim is advanced in the multipolar-condensate literature: dipole condensation can arise generically in bosonic systems through a self-proximity effect. In the Mott regime of the conventional Bose–Hubbard model,

$\langle b_i\rangle=0$13

so that

$\langle b_i\rangle=0$14

This “self-proximity effect” is then argued to guarantee $\langle b_i\rangle=0$15 for any nonzero $\langle b_i\rangle=0$16, even when $\langle b_i\rangle=0$17 itself vanishes (Xu et al., 2023). The same work proposes a dipolar Josephson effect with

$\langle b_i\rangle=0$18

so that supercurrents of dipoles arise in the absence of particle flows, and extends the mechanism to a hierarchy of multipolar condensates in which kinetics of the $\langle b_i\rangle=0$19th multipole induces condensation of the $\langle b_i\rangle=0$20th (Xu et al., 2023).

Synthetic gauge-field constructions sharpen the electrodynamic aspect of this picture. In a two-component ultracold-atom chain with spin-dependent tilt,

$\langle b_i\rangle=0$21

so that the dipole Hamiltonian takes the form

$\langle b_i\rangle=0$22

The resulting responses include an atomic analog of perfect Coulomb drag, in which increasing intercomponent interactions leads to equal and opposite displacements of the centers of mass of the two spin components, and a dipolar Josephson relation

$\langle b_i\rangle=0$23

after a short pulse of the rank-2 electric field imprints a dipole phase twist (Zhang et al., 17 Sep 2025).

These lattice results do not replace the dipolar-BEC meaning of dipole condensate. Rather, they delineate a second, conceptually independent branch of the subject in which the condensing object is a dipole operator and the central observables are two-body coherence, dipole transport, and multipolar Josephson effects.