We use the Ground-State Path Integral Monte Carlo method to study a Bose-Einstein condensate of strongly interacting NaCs polar molecules under the action of a fully anisotropic double microwave shielding potential characterized by a linear and an elliptical polarization field. In particular, we analyze the ground state of the system and its structure as a function of the ellipticity angle
ξ. While for the circularly polarized case (
ξ=0) a gas phase is realized, one or more self-bound droplets are observed for small
∣ξ∣'s above a threshold value near
3∘. With increasing
ξ, the observed droplets rapidly become tightly bound and are estimated to form a superfluid array. Our results compare favorably to the experimental observations in [Zhang et al., Nature \textbf{651}, 601 (2026)] for positive
ξ, while moderate differences show up for
ξ<0 where our simulations conform to the expected symmetries of the intermolecular potential.
The strong, tunable dipolar interactions of ultracold molecules make them a powerful platform for quantum many-body physics, but reaching degeneracy by evaporative cooling requires suppressing inelastic collisions. Microwave collisional shielding has enabled the preparation of degenerate Fermi gases and Bose-Einstein condensates of polar molecules, whereas static electric field shielding has been limited to fermionic species, which are less prone to inelastic loss. We demonstrate static-field Förster shielding of bosonic
23Na
87Rb molecules, suppressing two-body loss by up to four orders of magnitude. Within a bound-state-free electric-field window, three-body loss is also strongly suppressed, enabling efficient evaporation. Evaporating with an efficiency of 2.09(9), we increase the phase-space density of the gas by two orders of magnitude, reaching degeneracy with 7200(1000) molecules. We observe self-bound droplets at the end of evaporation over a wide range of field strengths, emerging from either degenerate or non-degenerate parent gases. We perform Path Integral Monte Carlo simulations, which suggest that the observed droplets are filamentary in nature, and find good agreement with the experimentally observed droplet formation temperatures. Our results establish Förster shielding as a single-field route to prepare degenerate gases and self-bound droplets of bosonic polar molecules.
We investigate the physics of ultracold dipolar molecules using path-integral quantum Monte Carlo simulations, and construct the complete phase diagram extending from weak to strong interactions and from small to mesoscopic particle numbers. Our calculations predict the formation of self-bound quantum droplets at interaction strengths lower than previous estimates for molecular condensates. Strikingly, we observe that, for stronger interactions, the oblate quantum droplet transitions to a two-dimensional sheet or superfluid membrane with a thickness of a single molecule. As interactions are increased further the system continually loses superfluidity while correlations develop, and is eventually found to undergo a transition to a crystalline monolayer that remains self-bound without external confinement. The spontaneous formation of such two-dimensional phases from a three-dimensional quantum gas is traced back to the peculiar anisotropic form of the dipole-dipole interaction generated by microwave dressing of rotational molecular states. For sufficiently large particle numbers, crystallization takes place for comparably low interaction strengths that do not promote two-body bound states and should thus be observable in ongoing experiments without limitations from three-body recombination.