Pith. sign in

REVIEW 3 major objections 7 minor 31 references

Microwave-shielded polar molecules form self-bound droplet arrays

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · glm-5.2

2026-07-08 15:40 UTC pith:ZAGCMMDO

load-bearing objection Self-bound droplets of NaCs molecules are well-supported; the supersolid claim is not — it rests on an upper bound consistent with zero superfluid fraction. the 3 major comments →

arxiv 2607.06130 v1 pith:ZAGCMMDO submitted 2026-07-07 cond-mat.quant-gas

Self-Bound Droplets of Ultracold Dipolar Molecules under Tunable Double Microwave Shielding

Roger Melero , Jordi Boronat , Ferran Mazzanti This is my paper
classification cond-mat.quant-gas
keywords dropletsdoublemicrowavemoleculesobservedpotentialself-boundshielding
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper uses ground-state Path Integral Monte Carlo (PIGS) to simulate 1500 NaCs polar molecules cooled to quantum degeneracy under a double microwave shielding potential combining linearly (pi) and elliptically (sigma) polarized fields. The ellipticity angle xi controls the anisotropy of the intermolecular interaction. At xi = 0 (circular polarization), the dipolar interaction vanishes and the system stays in a gas phase. For |xi| above a threshold near 3 degrees, attractive anisotropic channels open, the energy per particle turns negative, and the system self-organizes into one or more elongated, self-bound droplets that persist without external confinement. At larger |xi| (around 9 degrees), two droplets form with non-vanishing inter-droplet density, and a rough estimate of the superfluid fraction yields fs <= 0.8 +/- 0.1, suggesting the array may be superfluid and thus a candidate supersolid state. The simulations reproduce the xi -> -xi, phi -> -phi spatial symmetries of the interaction potential, which the experiment appears to break, and the authors attribute this discrepancy to possible metastable trapping in the experimental protocol.

Core claim

The central finding is that the fully anisotropic double-microwave-shielded NaCs interaction potential, when simulated ab initio with PIGS for N = 1500 molecules, produces self-bound droplet states for ellipticity angles |xi| > ~3 degrees and multi-droplet arrays at larger |xi|, with the energy per particle turning sharply negative. The non-vanishing density between droplets in the two-droplet regime at xi = 9 degrees, combined with a superfluid fraction upper bound near 0.8, points toward a self-bound supersolid array of polar molecules that survives without external trapping.

What carries the argument

The PIGS (Ground-State Path Integral Monte Carlo) method with an O(tau^6) short-time propagator and a trivial trial wave function Phi_T = 1, applied to the NaCs double-microwave-shielded interaction potential (dipolar plus anisotropic short-range van der Waals). The ellipticity angle xi is the single tuning parameter. Leggett's upper bound on the superfluid fraction provides the supersolidity diagnostic.

Load-bearing premise

The supersolid claim rests on a single rough estimate of the superfluid fraction using Leggett's upper bound, which the authors themselves describe as difficult to evaluate accurately, and no raw convergence data for this quantity is shown.

What would settle it

If a more precise superfluid fraction calculation (or an experimental measurement of inter-droplet phase coherence) yields a value near zero, the supersolid interpretation collapses despite the self-bound droplet array remaining real.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the supersolid interpretation holds, polar-molecule systems would achieve self-bound supersolidity without external trapping, unlike magnetic dipolar atom experiments that require confinement to prevent evaporation.
  • The xi -> -xi symmetry of the interaction potential provides a falsifiable prediction: ground-state droplet arrays at +xi and -xi should be rotated by 90 degrees but otherwise identical, which could be tested experimentally with improved trap-release protocols.
  • The threshold ellipticity angle near 3 degrees for self-binding identifies a sharp phase boundary that could be mapped experimentally by scanning xi with fine resolution.
  • The discrepancy between simulation (2 droplets) and experiment (3 droplets for positive xi) suggests the experimental system may be in a metastable state rather than the true ground state, motivating studies of relaxation dynamics.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the superfluid fraction estimate could be refined beyond Leggett's bound, the supersolid claim would either strengthen or weaken substantially; the current single rough number is the load-bearing evidence and a more direct superfluidity measurement would be decisive.
  • The N = 1500 simulation size may not capture finite-size effects relevant to experimental ensembles or the thermodynamic limit; scaling studies with larger N could reveal whether inter-droplet coherence persists or whether additional droplets appear.
  • The inter-droplet density contrast (~0.51) and the superfluid fraction (~0.8) are related but distinct quantities; a systematic study of how both vary with xi could reveal whether there is a continuous transition from isolated droplets to a coherent supersolid array.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 7 minor

Summary. The manuscript uses ground-state Path Integral Monte Carlo (PIGS) to study N=1500 NaCs polar molecules under double microwave shielding with tunable ellipticity angle ξ. The interaction potential (Eqs. 1–3) is taken from the experimental parameters of Ref. [17] without fitting. The authors find a gas-to-droplet transition at |ξ|≈3°, report energy-per-particle calculations showing self-bound states, and observe one- or two-droplet arrays at larger |ξ|. They compare their results with the experiment of Zhang et al. [17], finding reasonable agreement for ξ>0 but discrepancies for ξ<0. The paper concludes with an estimate of the superfluid fraction via Leggett's upper bound, suggesting a possible supersolid state.

Significance. The self-bound droplet formation claim is well-supported by the PIGS energy calculations (Fig. 3) and is a non-trivial ab-initio result for polar molecules with a parameter-free interaction potential. The use of PIGS with a high-order propagator (Ref. [28]) and Φ_T=1 is methodologically sound. The comparison with the experimental data of Ref. [17] adds value. However, the supersolid claim, which is the most novel headline claim, rests on an upper bound that is consistent with zero superfluid fraction and lacks supporting raw data or methodology. The paper does not provide machine-checked proofs, reproducible code, or falsifiable predictions beyond the existing experimental comparison.

major comments (3)
  1. The supersolid claim rests on an upper bound fs ≤ 0.8 ± 0.1 (Leggett's relation, Refs. [29–31]) that is fully consistent with fs = 0. An upper bound of 0.8 on a quantity ranging from 0 to 1 does not constitute positive evidence of superfluidity. The standard approach in path integral methods is to compute the superfluid fraction directly via the winding number estimator or the response to a phase twist. The paper does not report such a direct measurement. The authors acknowledge the Leggett bound is 'difficult to evaluate accurately' and that the contrast-based estimate (~0.51) is 'more a measure of particle overlap rather than a direct measure of a purely superfluid behavior.' Without a direct measurement or a lower bound on fs, the supersolid claim is conjectured, not established.
  2. The non-vanishing inter-droplet density shown in Fig. 5 is cited as supporting evidence for superfluidity, but quantum zero-point motion in a PIGS ground-state calculation will produce non-zero inter-droplet density regardless of phase coherence. This is necessary but not sufficient for superfluidity. The text should clarify this distinction and avoid implying that the density bridge alone is evidence of superfluid behavior.
  3. No methodology, raw data, or convergence study is provided for the Leggett bound calculation. The paper states fs ≤ (0.8 ± 0.1) for ξ = 9° but does not describe how the bound was evaluated, what input quantities entered the calculation, or how the uncertainty was determined. This makes the result impossible to assess or reproduce.
minor comments (7)
  1. The abstract states the droplets are 'estimated to form a superfluid array.' Given that the evidence is an upper bound consistent with zero, this phrasing overstates the result. Consider 'may form' or 'are candidates for' a superfluid array.
  2. Fig. 2a: the heatmaps would benefit from explicit color scale labels and consistent color ranges across panels for ξ = −6°, 0°, +6°.
  3. The sentence beginning 'Notably, the C̃6 coefficient...' contains a fragment: 'most affecting the behavior at x=y=0. at certain directions and ellipticity angles compared to the fully anisotropic description.' This is grammatically broken and difficult to parse.
  4. The number of droplets observed (at most two) is fewer than the three reported in the experiment for ξ ≳ 4.7°. The authors attribute this to uncertainty in molecule number but state this is 'fairly unlikely.' A more thorough discussion of possible causes (finite-size effects, trap geometry differences, timescale of equilibration) would help assess this explanation.
  5. The ξ → −ξ symmetry breaking observed in the simulations for |ξ| > 6° is attributed to 'a reminiscent dependence on the initial competition between the different anisotropic behaviors displayed by V(r) and the harmonic trap.' This is plausible but not demonstrated. Some analysis of whether this is a physical metastable state or a sampling artifact would strengthen the discussion.
  6. Reference [12] is listed as an arXiv preprint (arXiv:2512.14511). If the manuscript is published by the time of revision, the reference should be updated.
  7. The Conclusions state that the droplets exhibit 'superfluid behavior as estimated using Leggett's upper bound.' This conflates an upper bound with evidence of superfluid behavior; rewording is needed.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for a careful and constructive report. The referee's three major comments all concern the superfluidity/supersolidity claim, and we find them largely correct. We agree that (1) the Leggett upper bound alone does not establish superfluidity, (2) the non-vanishing inter-droplet density is necessary but not sufficient, and (3) the methodology for the Leggett bound calculation was insufficiently documented. We will revise the manuscript to soften the supersolid claim to a conjecture, clarify the distinction between density bridges and phase coherence, and add methodological details for the Leggett bound estimate. We also note that a direct winding-number measurement is in principle possible within PIGS and discuss the technical obstacles.

read point-by-point responses

  1. Referee: The supersolid claim rests on an upper bound fs ≤ 0.8 ± 0.1 (Leggett's relation) that is fully consistent with fs = 0. An upper bound of 0.8 on a quantity ranging from 0 to 1 does not constitute positive evidence of superfluidity. The standard approach is to compute fs directly via the winding number estimator or phase twist response. Without a direct measurement or a lower bound, the supersolid claim is conjectured, not established.

    Authors: The referee is correct on all counts. An upper bound of 0.8 on a quantity in [0,1] does not constitute positive evidence of superfluidity, and we should not have presented it as such. We will revise the manuscript to explicitly state that the supersolid claim is a conjecture, not an established result. The Leggett bound is reported as a consistency check, not as proof. Regarding a direct winding-number measurement: in PIGS (ground-state, T=0), the winding number estimator is not directly available in the same form as in finite-temperature PIMC, because the polymer chains are open rather than periodic. A phase-twist response or a projected estimator would require additional methodological development that is beyond the scope of this Letter. We acknowledge this limitation explicitly in the revised text. revision: yes

  2. Referee: The non-vanishing inter-droplet density in Fig. 5 is cited as supporting evidence for superfluidity, but quantum zero-point motion in a PIGS ground-state calculation will produce non-zero inter-droplet density regardless of phase coherence. This is necessary but not sufficient for superfluidity. The text should clarify this distinction.

    Authors: We agree. Non-vanishing inter-droplet density is a necessary condition for phase coherence between droplets but is not sufficient, as the referee correctly notes: zero-point motion alone can produce density bridges without superfluid connectivity. We will revise the text to state this distinction explicitly and remove any implication that the density bridge alone is evidence of superfluid behavior. The revised wording will describe the non-vanishing density as a necessary precondition for — but not evidence of — superfluidity. revision: yes

  3. Referee: No methodology, raw data, or convergence study is provided for the Leggett bound calculation. The paper states fs ≤ (0.8 ± 0.1) for ξ = 9° but does not describe how the bound was evaluated, what input quantities entered, or how the uncertainty was determined. This makes the result impossible to assess or reproduce.

    Authors: This is a fair criticism. The current manuscript does not describe the evaluation procedure for the Leggett bound. In the revised version, we will add a description of the method: the bound was estimated from the one-body density matrix and the condensate fraction extracted from the PIGS simulations, following the formulation in Refs. [29–31]. The uncertainty was estimated from the statistical noise in the density matrix elements. We will also add a note on the convergence with respect to the number of beads and projection time. We acknowledge that the estimate is rough, as already stated in the manuscript, and will frame it accordingly. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation chain is self-contained with external inputs

full rationale

The paper's derivation chain is straightforward and non-circular. (1) The interaction potential (Eqs. 1-3) uses parameters C3, C6, and A_l,m taken from the external experimental paper [17] (Zhang et al., Nature 2026) and a coupled-channels calculation [18]. These are not fitted to the droplet results claimed in this paper. (2) The PIGS Monte Carlo method [25-28] is an ab-initio ground-state quantum simulation; the energy per particle, density profiles, and droplet configurations are outputs of the simulation, not assumptions fed into it. (3) The self-bound droplet claim rests on negative energy per particle (Fig. 3), which is a direct PIGS computation result. (4) The superfluid fraction estimate fs ≤ 0.8 ± 0.1 uses Leggett's upper bound [29-31], which is an external theoretical result; while the estimate is rough and the supersolid claim is arguably overclaimed (an upper bound of 0.8 is consistent with fs=0), this is a correctness/overclaiming concern, not circularity—the result is not defined in terms of its own inputs. (5) Self-citations [16] (prior mean-field study by overlapping authors) and [26] (PIGS methodology) are contextual and methodological, not load-bearing for the central physics claims. No step in the derivation reduces to its inputs by construction.

Axiom & Free-Parameter Ledger

7 free parameters · 3 axioms · 0 invented entities

All parameters are taken from the experimental paper [17] or derived from physical constants. No new entities are postulated. The axioms are standard domain assumptions for quantum Monte Carlo, though convergence is not explicitly demonstrated.

free parameters (7)
  • C3 = 53100*sqrt(3) a0 hbar^2/m
    Dipolar interaction coefficient, taken from experimental Ref. [17], not fitted in this work.
  • C6 = (3200 a0)^4 hbar^2/m
    Short-range interaction coefficient, taken from experimental Ref. [17], not fitted in this work.
  • A_2,0 = 0.037
    Anisotropy coefficient, ratio from Ref. [17].
  • A_2,2 = 0.153
    Anisotropy coefficient, ratio from Ref. [17].
  • A_4,0 = 0.003
    Anisotropy coefficient, ratio from Ref. [17].
  • A_4,2 = 0.027
    Anisotropy coefficient, ratio from Ref. [17].
  • A_4,4 = 0.124
    Anisotropy coefficient, ratio from Ref. [17].
axioms (3)
  • domain assumption The intermolecular potential is well approximated by the sum of a dipolar term and a short-range van der Waals term (Eq. 1-3).
    This is a perturbative approximation from coupled-channels calculations (Ref. [18]); the paper notes the short-range contribution becomes comparable to the dipolar part at very short distances (~10^-2 r0).
  • domain assumption The PIGS method with the O(tau^6) propagator of Ref. [28] and Phi_T=1 converges to the exact ground state with the number of beads used.
    Standard assumption for PIGS, but the paper does not state the number of beads or demonstrate convergence explicitly.
  • domain assumption Leggett's relation provides a valid upper bound for the superfluid fraction in this anisotropic, finite-N system.
    Leggett's formula (Refs. [29-31]) is applied to estimate superfluidity; the authors note it is 'difficult to evaluate accurately' in this context.

pith-pipeline@v1.1.0-glm · 11135 in / 2668 out tokens · 410216 ms · 2026-07-08T15:40:04.647029+00:00 · methodology

0 comments
read the original abstract

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 $\xi$. While for the circularly polarized case ($\xi=0$) a gas phase is realized, one or more self-bound droplets are observed for small $|\xi|$'s above a threshold value near $3^\circ$. With increasing $\xi$, 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 $\xi$, while moderate differences show up for $\xi<0$ where our simulations conform to the expected symmetries of the intermolecular potential.

Figures

Figures reproduced from arXiv: 2607.06130 by Ferran Mazzanti, Jordi Boronat, Roger Melero.

Figure 1
Figure 1. Figure 1: FIG. 1: Schematic of the ultracold gas of NaCs molecules [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Energy per particle [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Characteristic snapshots of an ensemble of [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: One-dimensional density profiles for ellipticity [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages · 3 internal anchors

  1. [1]

    Lahaye, C

    T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. Pfau, Reports on Progress in Physics72, 126401 (2009)

    work page 2009

  2. [2]

    Schmitt, M

    M. Schmitt, M. Wenzel, F. B¨ ottcher, I. Ferrier-Barbut, and T. Pfau, Nature539, 259–262 (2016)

    work page 2016

  3. [3]

    Chomaz, S

    L. Chomaz, S. Baier, D. Petter, M. J. Mark, F. W¨ achtler, L. Santos, and F. Ferlaino, Phys. Rev. X6, 041039 (2016)

    work page 2016

  4. [4]

    Baillie and P

    D. Baillie and P. B. Blakie, Phys. Rev. A101, 043606 (2020)

    work page 2020

  5. [5]

    B¨ ottcher, M

    F. B¨ ottcher, M. Wenzel, J.-N. Schmidt, M. Guo, T. Lan- gen, I. Ferrier-Barbut, T. Pfau, R. Bomb´ ın, J. S´ anchez- Baena, J. Boronat, and F. Mazzanti, Phys. Rev. Res.1, 033088 (2019)

    work page 2019

  6. [6]

    M. Guo, F. B¨ ottcher, J. Hertkorn, J.-N. Schmidt, M. Wenzel, H. P. B¨ uchler, T. Langen, and T. Pfau, Na- ture574, 386 (2019)

    work page 2019

  7. [7]

    B¨ ottcher, J.-N

    F. B¨ ottcher, J.-N. Schmidt, M. Wenzel, J. Hertkorn, M. Guo, T. Langen, and T. Pfau, Physical Review X 9, 011051 (2019)

    work page 2019

  8. [8]

    Chomaz, D

    L. Chomaz, D. Petter, P. Ilzh¨ ofer, G. Natale, A. Traut- mann, C. Politi, G. Durastante, R. Van Bijnen, A. Patscheider, M. Sohmen,et al., Physical Review X 9, 021012 (2019)

    work page 2019

  9. [9]

    Macia, D

    A. Macia, D. Hufnagl, F. Mazzanti, J. Boronat, and R. E. Zillich, Phys. Rev. Lett.109, 235307 (2012)

    work page 2012

  10. [10]

    A.-C. Lee, D. Baillie, and P. B. Blakie, Phys. Rev. A109, 023323 (2024)

    work page 2024

  11. [11]

    J. Lin, G. Chen, M. Jin, Z. Shi, F. Deng, W. Zhang, G. Qu´ em´ ener, T. Shi, S. Yi, and D. Wang, Phys. Rev. X 13, 031032 (2023)

    work page 2023

  12. [12]

    Schindewolf, J

    A. Schindewolf, J. Hertkorn, I. Stevenson, M. Ciardi, P. Gross, D. Wang, T. Karman, G. Quemener, S. Will, T. Pohl,et al., arXiv preprint arXiv:2512.14511 (2025)

    work page arXiv 2025

  13. [13]

    Langen, J

    T. Langen, J. Boronat, J. S´ anchez-Baena, R. Bomb´ ın, T. Karman, and F. Mazzanti, Phys. Rev. Lett.134, 053001 (2025)

    work page 2025

  14. [14]

    Bigagli, W

    N. Bigagli, W. Yuan, S. Zhang, B. Bulatovic, T. Karman, I. Stevenson, and S. Will, Nature631, 289 (2024)

    work page 2024

  15. [15]

    Z. Shi, Z. Huang, F. Deng, W.-J. Jin, S. Yi, T. Shi, and D. Wang, arXiv e-prints , arXiv:2508.20518 (2025), arXiv:2508.20518 [cond-mat.quant-gas]

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  16. [16]

    S´ anchez-Baena, G

    J. S´ anchez-Baena, G. Pascual, R. Bomb´ ın, F. Mazzanti, and J. Boronat, Phys. Rev. Res.7, 033080 (2025)

    work page 2025

  17. [17]

    Zhang, W

    S. Zhang, W. Yuan, N. Bigagli, H. Kwak, T. Karman, I. Stevenson, and S. Will, Nature651, 601 (2026)

    work page 2026

  18. [18]

    Karman, N

    T. Karman, N. Bigagli, W. Yuan, S. Zhang, I. Stevenson, and S. Will, PRX Quantum6, 020358 (2025)

    work page 2025

  19. [19]

    Johnson, Journal of Computational Physics13, 445 (1973)

    B. Johnson, Journal of Computational Physics13, 445 (1973)

    work page 1973

  20. [20]

    Deng, X.-Y

    F. Deng, X.-Y. Chen, X.-Y. Luo, W. Zhang, S. Yi, and T. Shi, Phys. Rev. Lett.130, 183001 (2023)

    work page 2023

  21. [21]

    Supersolid Phases in Ultracold Gases of Microwave Shielded Polar Molecules

    W. Zhang, H. Liu, F. Deng, K. Chen, S. Yi, and T. Shi, arXiv e-prints 10.48550/arXiv.2506.23820 (2025)

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2506.23820 2025

  22. [22]

    Xu and G

    P. Xu and G. Chen, Phys. Rev. A112, 043306 (2025)

    work page 2025

  23. [23]

    F. Deng, X. Hu, W.-J. Jin, S. Yi, and T. Shi, Nature Communications16, 11219 (2025)

    work page 2025

  24. [24]

    Symmetry and Self-Bound Droplets in Dipolar Molecular Gases

    D. Baillie, arXiv preprint arXiv:2510.04634 (2025), arXiv:2510.04634 [cond-mat.quant-gas]

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  25. [25]

    Sarsa, K

    A. Sarsa, K. E. Schmidt, and W. R. Magro, The Journal of Chemical Physics113, 1366 (2000)

    work page 2000

  26. [26]

    R. Rota, J. Casulleras, F. Mazzanti, and J. Boronat, Phys. Rev. E81, 016707 (2010)

    work page 2010

  27. [27]

    Feynman,Statistical Mechanics: A Set Of Lectures (CRC press, 2018)

    R. Feynman,Statistical Mechanics: A Set Of Lectures (CRC press, 2018)

    work page 2018

  28. [28]

    S. A. Chin and C. R. Chen, The Journal of Chemical Physics117, 1409–1415 (2002)

    work page 2002

  29. [29]

    A. J. Leggett, Phys. Rev. Lett.25, 1543 (1970)

    work page 1970

  30. [30]

    Leggett, Journal of statistical physics93, 927 (1998)

    A. Leggett, Journal of statistical physics93, 927 (1998)

    work page 1998

  31. [31]

    P´ erez-Cruz, G

    D. P´ erez-Cruz, G. E. Astrakharchik, and P. Massignan, Physical Review A111, L011302 (2025)

    work page 2025