Phase identity, inter-droplet coherence, and stabilization mechanism of self-bound droplets in strongly dipolar microwave-dressed molecular gases

Self-bound droplets and droplet arrays were observed in a Bose-Einstein condensate of microwave-dressed sodium-cesium molecules (Zhang et al., Nature 651, 601 (2026)), realizing the first strongly dipolar quantum phases of ultracold molecules, with densities up to two orders of magnitude above the parent condensate and a dipolar length about four times the interparticle spacing. The experiment left the physical character of these phases undetermined. This problem asks for a validated characterization of the droplet phases: whether the strongly interacting single droplet is a superfluid quantum liquid or a crystal-like solid (tested by vorticity, collective excitations, and shear response); whether droplet arrays are phase-coherent and hence supersolid, or incoherent metastable clusters; and what the microscopic stabilization mechanism and the nature of the transitions between the weakly dipolar BEC, single-droplet, and droplet-array phases are — including whether the arrays are equilibrium phases or non-equilibrium states fixed by the preparation protocol — together with the auxiliary question of whether many-body correlations bias the loss-based density inference on which the strong-interaction claims rest. As of August 2026 no experiment has probed droplet phase identity or inter-droplet coherence, and microscopic theory remains quantitatively unsettled, with ab initio methods disagreeing at the factor-of-two level on the droplet-formation boundary.

Background

Ultracold gases of polar molecules carry permanent electric dipole moments, and the resulting dipole-dipole interaction — long-ranged (falling off as $1/r^3$) and anisotropic (attractive head-to-tail, repulsive side-by-side) — has been predicted for two decades to stabilize quantum phases with no short-range-interaction counterpart: self-bound quantum droplets, supersolids, and dipolar crystals. Two obstacles long blocked this program: rapid collisional loss destroyed molecular samples before degeneracy, and achievable dipoles were weak. Both were removed by microwave dressing: coupling the rotational ground state of a molecule to an excited rotational state with near-resonant microwave fields induces a controllable effective dipole-dipole interaction while simultaneously shielding molecules from inelastic collisions ('collisional shielding'). In double microwave shielding, the molecule is simultaneously dressed by a $\pi$ field, linearly polarized along the $z$ axis, and a $\sigma$ field whose elliptical polarization rotates in the perpendicular $xy$ plane; this dressing shields molecules from inelastic collisions, while tuning the ellipticity $\xi$ of the $\sigma$ field away from circular polarization induces strong in-plane dipolar interactions whose anisotropy is set by $\xi$. The residual van der Waals interaction acts as a steep repulsive hardcore potential at short range (radius $\sim 2{,}000\,a_0$ for NaCs, where $a_0$ is the Bohr radius), capping the density while suppressing loss. This enabled, successively, the first Bose-Einstein condensate (BEC) of dipolar molecules (Bigagli et al., Nature 631, 289 (2024)) and the observation of self-bound droplets (Zhang et al., Nature 651, 601 (2026)).

A self-bound quantum droplet is a liquid-like many-body state held together by its own interactions, without external confinement. In ultracold atoms, droplets were first observed in bosonic mixtures and in gases of strongly magnetic atoms (dysprosium, erbium); there, mean-field attraction is arrested by the Lee-Huang-Yang (LHY) correction — the leading beyond-mean-field contribution of quantum fluctuations to the energy of a dilute Bose gas — and the resulting droplets are dilute and weakly interacting. When such droplets form a regular array whose members retain quantum phase coherence, the state is a supersolid: it breaks continuous translational symmetry through density modulation while supporting dissipationless flow, and supersolids were realized in dipolar magnetic atoms in 2019. Whether analogous, or qualitatively new, phases appear when the dipoles are much stronger is the current frontier.

In the NaCs experiment, ramping the microwave-induced interactions across ellipticity $\xi$ produced three regimes: a weakly dipolar BEC for $|\xi|<3.5^{\circ}$; a single self-bound droplet for $\xi<-3.5^{\circ}$; and droplet arrays for $\xi>3.5^{\circ}$ (robust one-dimensional chains at slow interaction ramps, fluctuating two-dimensional structures at fast ramps spanning four orders of magnitude in ramp rate). Droplet mean densities reach $\sim 10^{14}\,\mathrm{cm}^{-3}$, about one hundred times the BEC density. At $\xi=-7.1(5)^{\circ}$ the mean interparticle spacing, $3{,}600(900)\,a_0$, approaches the hardcore radius $\sim 2{,}000\,a_0$, while the dipolar length (the characteristic length scale of the dipole-dipole interaction strength, here tunable with the dressing ellipticity), $13{,}000(900)\,a_0$, is about four spacings — the strongly interacting regime, quantified by $\langle n\rangle|a_{\mathrm{d2}}|^3>1$ (density times the cube of the dipolar length exceeding unity, meaning interaction energy dominates kinetic energy in the degenerate gas). The droplets are smaller than the imaging resolution, so densities were inferred indirectly from the decay of near-universal two-body loss — two-body loss at a rate coefficient close to the universal limit fixed by short-range collision physics — after switching off the dressing fields within $10\,\mu\mathrm{s}$, assuming many-body correlations do not modify the loss; lifetimes exceed $100$ ms. The theory that worked for atomic droplets does not simply transfer: the extended Gross-Pitaevskii equation (eGPE — mean-field theory supplemented by a local-density LHY correction) develops imaginary, internally inconsistent fluctuation contributions at the large relative dipolarities reachable with molecules, and ab initio path-integral Monte Carlo methods (PIMC, and its zero-temperature variant path-integral ground-state, PIGS — sampling the many-body density matrix by mapping quantum particles onto classical ring polymers) currently disagree at about the factor-of-two level on when molecular droplets form. The experiment's conclusion therefore explicitly poses the questions formalized below.

Problem Statement

Determine the physical character of the self-bound droplet phases of strongly dipolar, microwave-dressed molecular Bose gases — as observed in the NaCs experiment — and the mechanism by which they form and are stabilized. The problem has three linked components serving this single objective, plus one auxiliary component:

(i) Phase identity of the single droplet. Is the strongly interacting single droplet (mean density up to $\sim 10^{14}\,\mathrm{cm}^{-3}$, mean interparticle spacing $3{,}600(900)\,a_0$ approaching the hardcore radius $\sim 2{,}000\,a_0$, dipolar length about four spacings) a superfluid quantum liquid or a crystal-like solid? Discriminating diagnostics, as proposed by the source experiment, are vorticity under rotation, the collective excitation spectrum, and the shear response.

(ii) Coherence of droplet arrays. Do the droplets in an array share a common quantum phase — i.e., is the array a supersolid — or are they incoherent clusters, and can a phase-coherent molecular supersolid be reached in this platform? This requires phase-sensitive probes of coherence across droplets.

(iii) Stabilization mechanism and transitions. What is the microscopic stabilization mechanism of strongly dipolar droplets of microwave-dressed molecules (hardcore repulsion versus quantum-fluctuation contributions versus interaction-anisotropy effects), and what is the nature of the changes between the weakly dipolar BEC, single-droplet, and droplet-array regimes — equilibrium phase transitions (for example driven by roton-like softening, i.e., the instability of a condensate excitation mode at finite momentum) or non-equilibrium states fixed by the preparation protocol?

(iv) Auxiliary. Does the near-universal two-body-loss model used to infer droplet densities remain valid when many-body correlations at the highest measured densities are accounted for, or is the inferred density scale biased?

Components (i)–(iii) are jointly necessary and serve one objective — a validated account of what the observed phases are and why they are stable: component (iii) explains why the states of (i) and (ii) exist and whether the arrays of (ii) are equilibrium supersolids or metastable quenched states, while auxiliary component (iv) secures the density evidence on which the strong-interaction claims of (i) rest. The problem is complete when every sub-question of components (i)–(iii) — the single droplet's phase character, liquid or crystalline (component (i)); the arrays' coherence and whether a phase-coherent molecular supersolid is reachable in this platform (component (ii)); and the microscopic stabilization mechanism together with whether the regime changes are equilibrium or non-equilibrium (component (iii)) — is established by measurement and/or by microscopic theory quantitatively consistent with the experimentally measured critical ellipticities (approximately $\pm 3.5^{\circ}$), densities, and lifetimes; auxiliary component (iv) supplies supporting evidence for the density consistency that component (i) demands and is not a separate completion requirement.

References

1. Siwei Zhang, Weijun Yuan, Niccolò Bigagli, Haneul Kwak, Tijs Karman, Ian Stevenson, Sebastian Will, 'Observation of self-bound droplets of ultracold dipolar molecules', Nature 651, 601–606 (2026). DOI: 10.1038/s41586-026-10245-9; arXiv:2507.15208, https://arxiv.org/abs/2507.15208
2. Niccolò Bigagli, Weijun Yuan, Siwei Zhang, Boris Bulatovic, Tijs Karman, Ian Stevenson, Sebastian Will, 'Observation of Bose–Einstein condensation of dipolar molecules', Nature 631, 289–293 (2024). DOI: 10.1038/s41586-024-07492-z, https://doi.org/10.1038/s41586-024-07492-z
3. Matteo Ciardi, Kasper Rønning Pedersen, Tim Langen, Thomas Pohl, 'Self-Bound Superfluid Membranes and Monolayer Crystals of Ultracold Polar Molecules', Physical Review Letters 135, 153401 (2025). DOI: 10.1103/v7gw-xy36; arXiv:2504.02682, https://arxiv.org/abs/2504.02682
4. Matteo Ciardi, Andreas Schindewolf, Tim Langen, Thomas Pohl, 'Equilibrium and non-equilibrium phases of microwave-dressed polar molecules beyond rotational symmetries', arXiv:2606.30589 (2026), https://arxiv.org/abs/2606.30589
5. Tim Langen, Jordi Boronat, Juan Sánchez-Baena, Raúl Bombín, Tijs Karman, Ferran Mazzanti, 'Dipolar Droplets of Strongly Interacting Molecules', Physical Review Letters 134, 053001 (2025). DOI: 10.1103/PhysRevLett.134.053001, https://doi.org/10.1103/PhysRevLett.134.053001
6. D. Baillie, 'Symmetry and self-bound droplets in dipolar molecular gases', Physical Review Research 8, 023219 (2026). DOI: 10.1103/l9wk-v25v; arXiv:2510.04634, https://arxiv.org/abs/2510.04634
7. Xinyi Tang, Tianmiao Zhang, Zibin Zhao, Guilong Li, Zhaopin Chen, Bin Liu, Boris A. Malomed, Yongyao Li, 'Formation and dynamics of self-bound droplets in dipolar molecular condensate', arXiv:2606.22503 (2026), https://arxiv.org/abs/2606.22503
8. Tiziano Arnone Cardinale, Thomas Bland, Stephanie M. Reimann, 'Exploring supersolids of single-microwave shielded molecules via exact and mean-field theories', Communications Physics 9, 191 (2026). DOI: 10.1038/s42005-026-02706-4; arXiv:2509.18051, https://arxiv.org/abs/2509.18051
9. Andreas Schindewolf, Jens Hertkorn, Ian Stevenson, Matteo Ciardi, Phillip Groß, Dajun Wang, Tijs Karman, Goulven Quéméner, Sebastian Will, Thomas Pohl, Tim Langen, 'Colloquium: Strongly dipolar molecular Bose–Einstein condensates: From few- to many-body physics', Reviews of Modern Physics 98, 031002 (2026). DOI: 10.1103/r55l-f93m; arXiv:2512.14511, https://arxiv.org/abs/2512.14511
10. Haneul Kwak, Ian Stevenson, Weijun Yuan, Siwei Zhang, Asaf Toprakci, Lin Su, Tijs Karman, Sebastian Will, 'Electrostriction in a Bose–Einstein Condensate of Dipolar Molecules', arXiv:2608.19180 (2026), https://arxiv.org/abs/2608.19180