Electrostriction in a Bose-Einstein Condensate of Dipolar Molecules

Abstract: The recent creation of a Bose-Einstein condensate (BEC) of dipolar molecules has opened a new frontier for many-body quantum systems in which dipolar interactions can drive novel self-organization phenomena. Here, we observe electrostriction in a molecular BEC, an elliptical deformation driven by anisotropic dipolar interactions. We use double microwave dressing, involving $σ$- and $π$-polarized fields, to control non-axially symmetric dipolar interactions. We compare the experimental observations of electrostriction to a model based on an extended Gross-Pitaevskii equation and find excellent agreement in the regime of weak to moderate interactions. Using electrostriction, we demonstrate that the molecular BEC can be torqued by dynamically changing the orientation of the elliptical $σ$ microwave field. This provides a route to setting molecular quantum gases into rotation, opening opportunities to probe vorticity, superfluidity, and supersolidity in strongly dipolar matter.

• The paper observes interaction-driven electrostriction in a microwave-dressed NaCs molecular BEC, with anisotropic expansion increasing as the dipolar length reaches approximately 6900 a₀.
• The experiments agree with extended Gross–Pitaevskii equation simulations to about 10% in cloud size across weak-to-moderate interactions, showing that the Lee–Huang–Yang correction prevents the collapse predicted by standard mean-field theory.
• The paper demonstrates electrostirring by rotating the microwave polarization ellipse, producing a 47(2)° phase lag and centrifugal aspect-ratio dynamics that suggest—but do not yet prove—vortex formation.

Overview

Kwak, Stevenson, and colleagues report the observation of electrostriction—the elliptical deformation of a condensate along the attractive axis of anisotropic dipole-dipole interactions—in a Bose-Einstein condensate (BEC) of microwave-dressed NaCs molecules (2608.19180). The work extends the magnetostriction phenomenology established in magnetic-atom gases to molecular condensates and, more importantly, provides a direct benchmark of the extended Gross-Pitaevskii equation (eGPE) in a system where the contact-interaction approximation is far less obviously justified than for magnetic atoms. The authors also demonstrate that electrostriction can be exploited as an actuator: rotating the polarization ellipse of the microwave dressing field exerts a torque on the condensate and spins it up.

Interaction control via double microwave dressing

The molecules are dressed simultaneously by a circularly polarized $\sigma$ field and a linearly polarized $\pi$ field, which both suppress inelastic two-body loss through a repulsive shield and allow independent tuning of the multipole components of the dipolar interaction. Setting the $|m_l|=0$ and $|m_l|=1$ dipolar lengths $a_{\mathrm{d}0}$ and $a_{\mathrm{d}1}$ near zero leaves a purely non-axially symmetric interaction controlled by the ellipticity $\xi$ of the $\sigma$ field:

$$V_{\rm dd}(r) = \frac{\sqrt{3}\hbar^2}{mr^3}\, a_{\mathrm{d}2} \sin^2\theta\cos 2\varphi,$$

attractive along one in-plane axis and repulsive along the orthogonal axis. At short range ($r \lesssim 2000\,a_0$), the dressed potential is approximated by a contact term with s-wave scattering length $\pi$0, computed from coupled-channels calculations. Two features distinguish this platform from magnetic atoms: the scattering lengths (~$\pi$1–$\pi$2) and hard-core radius constitute a significant fraction of the interparticle spacing (~$\pi$3), and $\pi$4 itself depends on $\pi$5, decreasing monotonically from $\pi$6 at zero ellipticity toward negative values beyond $\pi$7. The eGPE comparison is therefore a nontrivial test rather than a foregone conclusion.

Observation of electrostriction and agreement with the eGPE

The experiments begin from BECs of roughly $\pi$8 molecules (condensate fraction above 70%) in an oblate trap round in the $\pi$9 plane, with peak densities around $|m_l|=0$0. Because the trap is isotropic in-plane, any observed anisotropy is attributable to the interactions. After ramping $|m_l|=0$1 to values between $|m_l|=0$2 and $|m_l|=0$3 (equivalently, $|m_l|=0$4 from ~$|m_l|=0$5 to ~$|m_l|=0$6), time-of-flight absorption imaging reveals strongly elliptical expansion: for $|m_l|=0$7 the aspect ratio grows continuously during expansion rather than showing the inversion characteristic of a contact-interacting condensate released from an anisotropic trap.

The central quantitative result is the agreement between measured cloud shapes and eGPE simulations across the full explored range, from nearly non-dipolar conditions ($|m_l|=0$8) up to $|m_l|=0$9. The aspect ratio agrees essentially within experimental uncertainty; absolute cloud sizes agree to about 10%, with expansion slightly faster in experiment—a residual discrepancy attributed plausibly to higher-order fluctuations or thermal-cloud coupling. A standard GPE without the LHY term predicts mean-field collapse above $|m_l|=1$0 and fails to converge beyond $|m_l|=1$1, while the experiment shows stable, increasingly deformed clouds. This directly establishes that quantum-fluctuation stabilization, captured by the LHY correction, is essential even in this weak-to-moderate interaction regime.

Two caveats attach to this conclusion. First, the dominant theoretical uncertainty is the $|m_l|=1$2 calibration of residual $|m_l|=1$3, translating into roughly $|m_l|=1$4 uncertainty in the predicted aspect ratio. Second, the simulations show density-modulated structures—two or three droplets connected by bridges—for $|m_l|=1$5, which the experiment cannot resolve owing to finite imaging resolution; the reported experimental aspect ratio of 2.0(2) versus a simulated value reaching 22.5 at $|m_l|=1$6 reflects this resolution limit rather than disagreement, but it means the droplet regime remains unverified by direct imaging here.

Treatment of the LHY term

For the non-axially symmetric interaction, the LHY coefficient involves the angular integral $|m_l|=1$7, which acquires an imaginary part when the Bogoliubov spectrum develops instabilities. Rather than discarding the imaginary contribution conventionally, the authors introduce a direction-dependent momentum cutoff that excludes unstable modes, yielding a real-valued modified $|m_l|=1$8. Within the explored range the two prescriptions differ only marginally—the imaginary part reaches at most ~7% of the real part at $|m_l|=1$9—but it grows rapidly beyond, so the validity of any LHY-based treatment at larger ellipticity is explicitly left open.

Electrostirring and rotation

Rotating the orientation of the $a_{\mathrm{d}0}$0-field polarization ellipse at $a_{\mathrm{d}0}$1 Hz while holding $a_{\mathrm{d}0}$2 torques the electrostricted condensate. After ~100 ms of spinning, the cloud follows the rotating ellipse with a phase lag of $a_{\mathrm{d}0}$3, and its aspect ratio shows an initial centrifugal overshoot to 2.3 followed by relaxation to ~1.5—qualitatively consistent with vortex nucleation in rotating superfluids. Direct evidence of vortices is not obtainable in the current setup: healing-length-scale cores (~$a_{\mathrm{d}0}$4m in situ) expand only to ~$a_{\mathrm{d}0}$5m during time of flight, below the imaging resolution. The rotational dynamics therefore suggest, but do not demonstrate, superfluid behavior in the molecular condensate.

Limitations and open questions

The paper's claims are carefully bounded. The eGPE validation covers only weak-to-moderate interactions ($a_{\mathrm{d}0}$6); beyond this, droplet formation sets in and the growing imaginary component of the conventional LHY integral undermines confidence in mean-field-plus-LHY descriptions. For $a_{\mathrm{d}0}$7, $a_{\mathrm{d}0}$8 becomes negative, yet stable molecular gases have been observed in that regime—an outright failure mode of the contact-interaction eGPE that no theory presented here addresses. The intermediate window $a_{\mathrm{d}0}$9 is identified as the natural testbed for discriminating between the eGPE, variational approaches, and quantum Monte Carlo methods. Experimentally, whether electrostirring nucleates quantized vortices, and whether rotating molecular droplets or supersolids behave as predicted, remain unobserved pending improved imaging resolution.

Conclusion

This work demonstrates quantitative control and understanding of anisotropic dipolar deformation in a molecular BEC and validates the eGPE—including an LHY term derived for non-axially symmetric interactions—as the correct effective description in the weak-to-moderate coupling regime, despite short-range physics that departs substantially from the magnetic-atom paradigm. It further introduces electrostirring as a precision tool for rotational excitation of molecular quantum gases. The principal open problems are theoretical: extending accurate descriptions into the droplet-forming and negative-scattering-length regimes where the current framework demonstrably breaks down.