REVIEW 2 major objections 4 minor 4 cited by
Supersolid Phases in Ultracold Gases of Microwave Shielded Polar Molecules
T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Adding ellipticity to a microwave dressing field makes trapped polar molecules form a supersolid at densities an order of magnitude lower than atomic dipolar supersolids, in a parameter window current experiments can reach.
desk verdict A credible new molecular-supersolid mechanism via elliptical microwave polarization; the phase diagram holds up, but the 'accessible to current experiments' claim needs a loss-rate estimate. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the effective two-body potential Veff(r) derived to second order in the microwave coupling, Eq. (1), which contains a 1/$r^{3}$ dipolar term with an isotropic C3,0 part and an in-plane anisotropic C3,2 part, plus a 1/$r^{6}$ shielding potential parametrized by six independent C6 coefficients. The ellipticity ξ controls C3,2, which is proportional to sin 2ξ and makes the interaction repulsive along the x axis and attractive along the y axis. The load-bearing numerical machinery is path-integral Monte Carlo with the worm algorithm, which provides an unbiased finite-temperature treatment of the trapped gas without mean-field approximations; the phases are identified through column densities, the structure factor S(kx,0,0) for crystalline order, and the reduced peak-to-peak correlation g_pp for coherence.
What would settle it
In a trapped NaCs gas with dual microwave shielding, tune the ellipticity to about 4 degrees and set the y-trap frequency to about 100 Hz at roughly 4 nK, then measure the in-plane structure factor and the peak-to-peak correlation function: if no density-modulation side peaks appear in S(kx,0,0) or if the inter-peak coherence g_pp remains below 1 percent, the predicted supersolid window is not realized.
Extended reading notes
Core claim
The central claim is that replacing the σ-polarized microwave with an elliptically polarized one introduces an additional anisotropic term, C3,2 $sin^{2}$θ cos2φ, in the effective dipolar interaction of microwave-shielded polar molecules, and that the interplay of this interaction anisotropy with the trapping anisotropy produces a supersolid phase. In the ellipticity–trap-frequency plane, the authors identify three phases—expanding gas, supersolid, and self-bound droplet—and show by path-integral Monte Carlo that the supersolid has multiple density peaks along the x direction while retaining significant peak-to-peak coherence. They further show that the supersolid phase sits in a parameter window accessible to current NaCs experiments and that it develops at substantially lower density than previously considered molecular or atomic supersolid schemes.
Load-bearing premise
The many-body simulations use an effective two-body potential derived to second order in the microwave coupling and assume that three-body recombination and the finite lifetime of the dressed molecular state are negligible at the densities reached in the supersolid window.
Editorial extensions
If this is right
- A supersolid of polar molecules should appear at peak densities around 1.5×10^12 cm^-3, about two orders of magnitude lower than dysprosium droplet supersolids, which the authors argue alleviates the three-body-loss problem.
- The phase diagram on the ellipticity–trap-frequency plane predicts tunable transitions: increasing ellipticity at fixed trap drives the gas from an expanding gas through the supersolid to a self-bound droplet, while increasing the y-trap frequency at fixed ellipticity drives it from droplet to supersolid to expanding gas.
- The inter-droplet coherence, measured by g_pp, decreases monotonically with ellipticity and increases with trap frequency, giving clear experimental signatures of where the supersolid phase melts.
- The superfluid fraction is anisotropic, with the x-component dominant, and the finite y- and z-components signal phase coherence between separate density peaks, a hallmark of supersolidity.
- The scheme relies only on the permanent dipole moment and microwave shielding, so it should extend to other bialkali polar molecules beyond NaCs.
Reading between the lines
- If the two-body effective potential is accurate but three-body recombination is not fully quenched at the predicted densities, the experimental lifetimes may be shorter than the density argument suggests; a direct loss-rate measurement at the supersolid parameters would settle this.
- A sharper test of supersolidity than g_pp>1% would be a time-of-flight interference measurement of the entire droplet array, which would reveal long-range phase coherence directly.
- The same anisotropy mechanism might stabilize other modulated phases in two-dimensional traps, such as stripe or checkerboard supersolids, where the sign of C3,2 can be flipped by changing the sign of the ellipticity.
- The predicted phase boundaries could be mapped experimentally by quenching ellipticity from zero to about 4 degrees and monitoring the growth of side peaks in the structure factor.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript proposes a scheme for realizing a supersolid phase in a trapped gas of microwave-shielded NaCs polar molecules using a combination of a π-polarized and an elliptically polarized microwave. The authors derive an analytic effective two-body potential (Eq. (1)) whose in-plane anisotropy is controlled by the ellipticity ξ, and validate it against adiabatic and multichannel scattering calculations (SM2, SM3). The many-body Hamiltonian (Eq. (2)) is then studied by path-integral Monte Carlo with the worm algorithm at T = 4 nK for a harmonic trap with ωx/2π = 20 Hz and ωz/2π = 80 Hz, mapping the (ξ, ωy) plane into expanding gas (EG), supersolid (SS), and self-bound droplet (SBD) phases. The central finding is a supersolid phase with density modulation along the x axis and finite peak-to-peak coherence g_pp for ξ around 4° and ωy/2π between about 60 and 120 Hz, at peak densities of order 1.5×10^12 cm^-3. The authors argue that this density is substantially lower than in atomic dipolar supersolids, which they claim leads to a significantly extended molecular gas lifetime.
Significance. If the central claim holds, the paper provides an important new platform for supersolidity: a molecular system with a tunable, fully characterized anisotropic interaction and a numerically exact many-body treatment. The strengths are the careful analytic derivation of Veff and its cross-validation against adiabatic and multi-channel scattering results, the use of unbiased PIMC for the stated Hamiltonian, and the explicit finite-temperature phase diagram with a concrete experimental parameter set. The main weakness is that the promise of 'accessibility to current experiments' relies on unquantified inelastic loss rates and on simulation parameters (notably N) that are not stated. Provided these points are addressed, the paper would be a significant contribution to the field.
major comments (2)
- [Eq. (2), SM3] The abstract and conclusion assert that the supersolid phase is 'in the parameter regime accessible to current experiments' and that the low density 'leads to a significantly extended lifetime of molecular gases.' These claims are not supported by the calculation as presented. The many-body Hamiltonian in Eq. (2) contains only the kinetic, trap, and pairwise Veff terms; no inelastic two-body or three-body loss terms are included, and SM3 defines the inelastic rate β_inel but reports no value for the dual-microwave parameters used here (Ωσ/2π = 7.9 MHz, δσ/2π = −8 MHz, Ωπ/2π = 6.5 MHz, δπ/2π = −10 MHz, ξ ≈ 4°). Because this C3,0 ≈ 0 configuration is not the single-microwave NaCs setup demonstrated in Ref. [5], the assumption that losses are negligible at the peak density ≈1.5×10^12 cm^-3 and T = 4 nK is unverified. I request a quantitative estimate of β_inel (and, where possible, the three-body loss rate K3) for the proposed parameters, or a careful justification based on the existing shielding literature. If such an estimate is not possible, the 'accessible' and 'extended lifetime' claims should be tempered accordingly.
- [Many-body quantum phases] The main text repeatedly refers to 'a trapped gas of N NaCs molecules' but never states the value of N used in the PIMC simulations, nor does it give the PIMC discretization parameters (e.g., imaginary-time step, interaction cutoff, number of beads). The quoted peak density and the phase boundaries in Fig. 2 depend on N; without this information the results are not reproducible and the 'low-density' claim cannot be verified. Please state N explicitly and provide a brief finite-size analysis (e.g., results for two or three N values) to show that the phase regions do not shift qualitatively with N.
minor comments (4)
- [Effective inter-molecular interaction] The sentence preceding Eq. (1) states that the effective potential takes this form 'for ∣ξ∣ > 15○'. This conflicts with the rest of the paper, where all results are for ξ between 0° and about 5°, and with the SM validations for ξ = 0°, 3°, 6°. Please correct the domain of validity (likely 'ξ > 0' or 'small ξ').
- [Many-body quantum phases] The phase boundaries are defined by the cutoffs S(ksub,0,0) > 1% and g_pp ≤ 1%. Since the system is finite and trapped, these are crossover criteria rather than phase transitions; a brief discussion of how the SS region shifts when these thresholds are varied (e.g., 0.5% or 2%) would place the phase diagram on firmer footing.
- [Many-body quantum phases] The definition of ksub as 'the momentum at the sub-dominated peak along the x-direction' could be made more precise: it would be useful to state explicitly that ksub is the largest nonzero peak of S(kx,0,0) for kx ≠ 0.
- [General] There are several typographical errors, including 'anistropy' in the Fig. 4 caption and 'ebahnced' in SM3; a careful proofread is recommended.
Circularity Check
No circularity: the effective potential is derived and independently cross-validated, and the PIMC phase diagram is an unbiased output for that Hamiltonian.
full rationale
The paper's derivation chain is self-contained. The effective potential Veff(r) of Eq. (1) is obtained from Floquet perturbation theory in the Supplemental Material: the first-order coefficients C3,0 and C3,2 are given explicitly in Eqs. (SM14)-(SM15), and the second-order C6 terms are collected from Eq. (SM16). No many-body phase information is used as input. The potential is then checked against two independent two-body calculations: the adiabatic Born-Oppenheimer potential (SM2, Figs. SM2b-c) and full multichannel scattering in Floquet space (SM3). The comparison includes the s-wave scattering length versus ellipticity, with the shape-resonance position agreeing within about 0.4 degrees (Fig. SM2d). This is a cross-validation, not a fit to the target phase diagram. The many-body Hamiltonian in Eq. (2) uses this fixed pairwise potential, and the PIMC worm-algorithm simulation is unbiased for that Hamiltonian; the phase boundaries and coherence diagnostics (S(k)>1% and g_pp<=1%) are criteria applied to computed observables, not fitted parameters. The self-citations, including refs. [8], [21], [32], and [35], refer to the same group's earlier formalism, but the present manuscript actually supplies the derivation instead of merely invoking an unverified prior result. The claim that the parameters are 'accessible to current experiments' and the inferred lifetime extension are not quantitatively supported by inelastic-loss or three-body-loss rate calculations for the proposed parameters, but that is an external-validity and correctness gap, not a circularity of the derivation. No step was found in which a prediction reduces to its own input by construction.
Assumptions & free parameters
free parameters (3)
- Omega_pi (pi-polarized microwave Rabi frequency) =
2π × 6.5 MHz
- S(k_sub) phase-boundary threshold =
1%
- g_pp phase-boundary threshold =
1%
assumptions (6)
- domain assumption Born-Oppenheimer approximation: the internal molecular state adjusts instantly to the inter-molecular distance, so the interaction can be described by an adiabatic effective potential.
- domain assumption Pairwise additivity: the many-body Hamiltonian contains only two-body interactions; three-body forces and losses are neglected.
- domain assumption The gas is prepared in the highest dressed state |+> with a lifetime long enough that the equilibrium phase is relevant.
- standard math Floquet theory and second-order perturbation theory provide a converged description of the two-body interaction for the ellipticities used (|ξ| > 15 degrees, and ξ around 4 degrees in the many-body simulations).
- domain assumption The universal van der Waals interaction CvdW/r^6 is negligible compared to the microwave shielding C6/r^6 in the many-body regime.
- domain assumption The trap is perfectly harmonic with the stated frequencies; no anharmonicities or optical lattice effects are present.
Cite this review
Pith. "Pith review of Supersolid Phases in Ultracold Gases of Microwave Shielded Polar Molecules." pith.science (2026). https://pith.science/paper/N55PVGRD
@misc{pith2026250623820,
author = {Pith},
title = {Pith review of: Supersolid Phases in Ultracold Gases of Microwave Shielded Polar Molecules},
year = {2026},
howpublished = {\url{https://pith.science/paper/N55PVGRD}},
note = {Machine review of arXiv:2506.23820}
}
read the original abstract
We propose a novel scheme to realize the supersolid phase in ultracold gases of microwave-shielded polar molecules by engineering an additional anisotropy in inter-molecular dipolar interaction via an elliptically polarized microwave. It is shown through quantum Monte-Carlo calculations that the interplay of the anisotropies between the interaction and trapping potential gives rise to rich quantum phases. Particularly, it is found that the supersolid phase emerges in the parameter regime accessible to current experiments. Our study paves the way for exploring the properties of supersolid phases in ultracold gases of polar molecules.
Figures
Forward citations
Cited by 4 Pith papers
-
Ro-vibrational van der Waals interaction between ultracold polar molecules
Molecules in (v,j)=(0,1) and (1,0) feel a strong ro-vibrational van der Waals repulsion, C6=d_e^4/(9α_e), that can cut collisional loss by orders of magnitude.
-
Effective potentials for polar molecules under non-orthogonal dual microwave fields
Floquet theory yields an analytic effective potential for polar molecules in non-orthogonal dual microwave fields, showing moderate shielding reduction yet strongly suppressed inelastic losses and extra interaction tu...
-
Self-Bound Droplets of Ultracold Dipolar Molecules under Tunable Double Microwave Shielding
Ground-state path integral Monte Carlo simulations of NaCs polar molecules with double microwave shielding show self-bound droplets forming above an ellipticity threshold of ~3 degrees, with an estimated superfluid fr...
-
Symmetry and Self-Bound Droplets in Dipolar Molecular Gases
A D3 symmetry tiles the two-parameter interaction plane of microwave-dressed molecules, and this classification yields the phase diagram and scaling laws for self-bound molecular droplets.
Reference graph
Works this paper leans on
-
[5]
N. Bigagli, W. Yuan, S. Zhang, B. Bulatovic, T. Kar- man, I. Stevenson, and S. Will, Observation of bose- einstein condensation of dipolar molecules, Nature 631, 289 (2024)
work page 2024
-
[1]
Schindewolf, R
A. Schindewolf, R. Bause, X.-Y. Chen, M. Duda, T. Kar- man, I. Bloch, and X.-Y. Luo, Evaporation of microwave- shielded polar molecules to quantum degeneracy, Nature 607, 677 (2022)
2022
-
[2]
X.-Y. Chen, A. Schindewolf, S. Eppelt, R. Bause, M. Duda, S. Biswas, T. Karman, T. Hilker, I. Bloch, and X.-Y. Luo, Field-linked resonances of polar molecules, Nature 614, 59 (2022)
work page 2022
-
[3]
X.-Y. Chen, S. Biswas, S. Eppelt, A. Schindewolf, F. Deng, T. Shi, S. Yi, T. A. Hilker, I. Bloch, and X.-Y. Luo, Ultracold field-linked tetratomic molecules, Nature 626, 283 (2024)
work page 2024
-
[4]
J. Lin, G. Chen, M. Jin, Z. Shi, F. Deng, W. Zhang, G. Qu´ em´ ener, T. Shi, S. Yi, and D. Wang, Microwave shielding of bosonic narb molecules, Phys. Rev. X 13, 031032 (2023)
work page 2023
-
[6]
W.-J. Jin, F. Deng, S. Yi, and T. Shi, Bose-einstein con- densates of microwave-shielded polar molecules (2024), arXiv:2406.06412 [cond-mat.quant-gas]
work page Pith review arXiv 2024
- [7]
-
[8]
Quantum Phases for Finite-Temperature Gases of Bosonic Polar Molecules Shielded by Dual Microwaves
W. Zhang, K. Chen, S. Yi, and T. Shi, Quantum phases for finite-temperature gases of bosonic polar molecules shielded by dual microwaves (2025), arXiv:2503.02644 [cond-mat.quant-gas]
work page Pith review arXiv 2025
Show all 40 references
-
[9]
E. P. Gross, Unified theory of interacting bosons, Phys. Rev. 106, 161 (1957)
1957
-
[10]
A. F. Andreev and I. M. Lifshitz, Quantum theory of defects in crystals, Sov. Phys. JETP 29, 1107 (1969)
1969
-
[11]
A. J. Leggett, Can a solid be ”superfluid”?, Phys. Rev. Lett. 25, 1543 (1970)
1970
-
[12]
G. V. Chester, Speculations on bose-einstein condensa- tion and quantum crystals, Phys. Rev. A 2, 256 (1970)
1970
-
[13]
Boninsegni and N
M. Boninsegni and N. V. Prokof’ev, Colloquium: Super- solids: What and where are they?, Rev. Mod. Phys. 84, 759 (2012)
2012
-
[14]
D. Y. Kim and M. H. W. Chan, Absence of supersolidity in solid helium in porous vycor glass, Phys. Rev. Lett. 109, 155301 (2012)
2012
-
[15]
J.-R. Li, J. Lee, W. Huang, S. Burchesky, B. Shteynas, F. C. Top, A. O. Jamison, and W. Ketterle, A stripe phase with supersolid properties in spin–orbit-coupled bose–einstein condensates, Nature 543, 91 (2017)
2017
-
[16]
L´ eonard, A
J. L´ eonard, A. Morales, P. Zupancic, T. Don- ner, and T. Esslinger, Monitoring and manip- ulating higgs and goldstone modes in a super- solid quantum gas, Science 358, 1415 (2017), https://www.science.org/doi/pdf/10.1126/science.aan2608
2017 doi
-
[17]
M. Guo, F. Boettcher, J. Hertkorn, J.-N. Schmidt, 6 M. Wenzel, H. P. Buechler, T. Langen, and T. Pfau, The low-energy goldstone mode in a trapped dipolar super- solid, Nature 574, 386 (2019)
2019
-
[18]
Tanzi, E
L. Tanzi, E. Lucioni, F. Fam` a, J. Catani, A. Fioretti, C. Gabbanini, R. N. Bisset, L. Santos, and G. Modugno, Observation of a dipolar quantum gas with metastable supersolid properties, Phys. Rev. Lett. 122, 130405 (2019)
2019
-
[19]
B¨ ottcher, J.-N
F. B¨ ottcher, J.-N. Schmidt, M. Wenzel, J. Hertkorn, M. Guo, T. Langen, and T. Pfau, Transient supersolid properties in an array of dipolar quantum droplets, Phys. Rev. X 9, 011051 (2019)
2019
-
[20]
Chomaz, D
L. Chomaz, D. Petter, P. Ilzh¨ ofer, G. Natale, A. Traut- mann, C. Politi, G. Durastante, R. M. W. van Bijnen, A. Patscheider, M. Sohmen, M. J. Mark, and F. Ferlaino, Long-lived and transient supersolid behaviors in dipolar quantum gases, Phys. Rev. X 9, 021012 (2019)
2019
-
[21]
F. Deng, X. Hu, W.-J. Jin, S. Yi, and T. Shi, Two- and many-body physics of ultracold molecules dressed by dual microwave fields (2025), arXiv:2501.05210 [cond- mat.quant-gas]
2025 arXiv
-
[22]
Karman, N
T. Karman, N. Bigagli, W. Yuan, S. Zhang, I. Steven- son, and S. Will, Double microwave shielding (2025), arXiv:2501.08095 [cond-mat.quant-gas]
2025
-
[23]
Ciardi, K
M. Ciardi, K. R. Pedersen, T. Langen, and T. Pohl, Self- bound monolayer crystals of ultracold polar molecules (2025), arXiv:2504.02682 [cond-mat.quant-gas]
2025 arXiv
-
[24]
Prokof’Ev, B
N. Prokof’Ev, B. Svistunov, and I. Tupitsyn, Exact, complete, and universal continuous-time worldline monte carlo approach to the statistics of discrete quantum sys- tems, Journal of Experimental and Theoretical Physics 87, 310 (1998)
1998
-
[25]
Boninsegni, N
M. Boninsegni, N. Prokof’ev, and B. Svistunov, Worm algorithm for continuous-space path integral monte carlo simulations, Phys. Rev. Lett. 96, 070601 (2006)
2006
-
[26]
Boninsegni, N
M. Boninsegni, N. V. Prokof’ev, and B. V. Svistunov, Worm algorithm and diagrammatic monte carlo: A new approach to continuous-space path integral monte carlo simulations, Phys. Rev. E 74, 036701 (2006)
2006
-
[27]
K. Chen, L. Liu, Y. Deng, L. Pollet, and N. Prokof’ev, Universal conductivity in a two-dimensional superfluid- to-insulator quantum critical system, Phys. Rev. Lett. 112, 030402 (2014)
2014
-
[28]
K. Chen, Y. Huang, Y. Deng, A. B. Kuklov, N. V. Prokof’ev, and B. V. Svistunov, Deconfined criticality flow in the heisenberg model with ring-exchange interac- tions, Phys. Rev. Lett. 110, 185701 (2013)
2013
-
[29]
Boninsegni, N
M. Boninsegni, N. Prokof’ev, and B. Svistunov, Super- glass phase of 4He, Phys. Rev. Lett. 96, 105301 (2006)
2006
-
[30]
Kuklov, L
A. Kuklov, L. Pollet, N. Prokof’ev, and B. Svistunov, Transverse quantum superfluids, Annual Review of Con- densed Matter Physics https://doi.org/10.1146/annurev- conmatphys-042924-103908 (2024)
2024 doi
-
[31]
Zhang, B
C. Zhang, B. Capogrosso-Sansone, M. Boninsegni, N. V. Prokof’ev, and B. V. Svistunov, Superconducting tran- sition temperature of the bose one-component plasma, Phys. Rev. Lett. 130, 236001 (2023)
2023
-
[32]
See Supplemental Material for the time-dependent Hamiltonian for two molecules subjected to both an ellip- tically polarized and a linearly polarized microwave field, the effective interaction potential between two molecules in the highest dressed state, the justification of t...
-
[33]
Stevenson, S
I. Stevenson, S. Singh, A. Elkamshishy, N. Bigagli, W. Yuan, S. Zhang, C. H. Greene, and S. Will, Three- body recombination of ultracold microwave-shielded po- lar molecules (2024), arXiv:2407.04901 [cond-mat.quant- gas]
2024
-
[34]
Karman and J
T. Karman and J. M. Hutson, Microwave shielding of ultracold polar molecules, Phys. Rev. Lett. 121, 163401 (2018)
2018
-
[35]
Deng, X.-Y
F. Deng, X.-Y. Chen, X.-Y. Luo, W. Zhang, S. Yi, and T. Shi, Effective potential and superfluidity of microwave-shielded polar molecules, Phys. Rev. Lett. 130, 183001 (2023)
2023
-
[36]
F. Deng, X. Chen, X. Luo, W. Zhang, S. Yi, and T. Shi, Formation and dissociation of field-linked tetramers (2024), arXiv:2405.13645 [Quantum Physics]
2024 arXiv
-
[37]
Idziaszek and P
Z. Idziaszek and P. S. Julienne, Universal rate constants for reactive collisions of ultracold molecules, Phys. Rev. Lett. 104, 113202 (2010)
2010
-
[38]
Lepers, R
M. Lepers, R. Vexiau, M. Aymar, N. Bouloufa-Maafa, and O. Dulieu, Long-range interactions between polar alkali-metal diatoms in external electric fields, Phys. Rev. A 88, 032709 (2013)
2013
-
[39]
A. L. Fetter, Rotating trapped bose-einstein condensates, Rev. Mod. Phys. 81, 647 (2009)
2009
-
[40]
Floquet-Fourier
D. M. Ceperley, Path integrals in the theory of condensed helium, Rev. Mod. Phys. 67, 279 (1995). 1 Supplemental Materials This Supplemental Material is structured as follows. In Section SM 1, we present a model for polar molecules subjected to both an elliptically polarized a...
1995
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.