{"id":"68ea3bec-d2d2-438c-806e-8e6c7973eb31","arxiv_id":"2607.06130","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":7,"one_line_summary":"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 fraction up to 0.8.","lead":"This paper uses quantum Monte Carlo simulations to show that ultracold NaCs molecules under double microwave shielding form self-bound droplets and droplet arrays when the microwave ellipticity angle is tuned away from zero. The results matter because self-bound superfluid droplet arrays in polar molecules would be a new platform for studying supersolidity and strongly correlated quantum matter.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The supersolid claim rests on an upper bound (fs ≤ 0.8) that is fully consistent with fs = 0; no direct superfluid fraction measurement is provided.","rationale":"The reader correctly identified the superfluid fraction estimate as the weakest link, but the concern is more fundamental than 'difficult to evaluate accurately' or 'no raw data shown.' An upper bound of 0.8 from Leggett's relation is not evidence of superfluidity—it is consistent with fs = 0. The reader states the bound was 'computed from the winding number or density response,' but Leggett's bound is derived from the density profile, not from winding numbers. The deeper issue is that the paper needs a direct measurement (winding number or phase-twist response) rather than an upper bound. The self-bound droplet claim is solid and well-supported by PIGS energy calculations. The CONDITIONAL verdict is appropriate: the droplet formation result stands, but the supersolid claim requires direct evidence of superfluidity before it can be accepted. I note that the ξ→-ξ symmetry breaking for |ξ| > 6° raises a secondary concern about whether PIGS has converged to the true ground state or is stuck in a trap-history-dependent metastable state, but this affects the detailed density profiles rather than the core self-binding result. The missing simulation parameters (bead count, time step τ) limit reproducibility but are not load-bearing for the central claims given the established nature of the PIGS method.","tokens_in":8111,"tokens_out":4126,"duration_ms":153047,"concrete_test":"Compute the superfluid fraction directly at ξ = 9° using the winding number estimator (mean-square winding number) or, if open-boundary PIGS precludes this, compute the response of the ground-state energy to an imposed phase twist (twisted boundary conditions). If the direct estimate yields fs significantly above zero (e.g., > 0.1 with controlled statistical error), the supersolid claim is supported. If fs is consistent with zero, the claim should be retracted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline novel claim is a supersolid state of self-bound droplets. The evidence for superfluidity is twofold: (1) non-vanishing inter-droplet density (Fig. 5), which is necessary but not sufficient for superfluidity—quantum zero-point motion in a PIGS ground-state calculation will produce non-zero inter-droplet density regardless of phase coherence; and (2) an upper bound fs ≤ 0.8 ± 0.1 from Leggett's relation, described as 'difficult to evaluate accurately' with no methodology or raw data shown. The critical issue is that an upper bound of 0.8 on a quantity ranging from 0 to 1 does not constitute evidence of superfluidity. It merely says the density modulation is not so extreme as to rule it out. The actual superfluid fraction could be zero. The standard approach in path integral methods is to compute the superfluid fraction directly via the winding number estimator (or equivalently, the response to a phase twist / Galilean boost), which yields a direct estimate rather than an upper bound. The paper does not report this. The contrast-based estimate (~0.51) is acknowledged by the authors themselves as '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 unsupported. The self-bound droplet formation claim, by contrast, is well-supported by the negative PIGS energy calculations and is consistent with the cited experiment for positive ξ.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","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.","tokens_in":9001,"tokens_out":1174,"duration_ms":242590,"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":[{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null}],"minor_comments":[{"comment":"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.","section":null},{"comment":"Fig. 2a: the heatmaps would benefit from explicit color scale labels and consistent color ranges across panels for ξ = −6°, 0°, +6°.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null}],"recommendation":"major_revision","confidential_remarks":"The self-bound droplet formation result is solid and publishable on its own. The supersolid claim, however, is the headline novelty and is not supported by the evidence presented. If the authors can either (a) compute the superfluid fraction directly via the winding number estimator in PIGS, or (b) substantially soften the supersolid claim to a conjecture supported by the density bridge and upper bound, the paper could be acceptable. The current framing oversells what the calculations show."},"author_rebuttal":{"model":"glm-5.2","summary":"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.","responses":[{"response":"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_made":"yes","referee_comment":"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."},{"response":"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_made":"yes","referee_comment":"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."},{"response":"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_made":"yes","referee_comment":"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."}],"tokens_in":7976,"tokens_out":857,"duration_ms":104516,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The main result that holds up: PIGS Monte Carlo shows self-bound droplets of NaCs polar molecules under double microwave shielding for ellipticity angles above about 3 degrees. The energy calculations are clean — negative energy per particle, clear transition from gas to bound state — and the interaction potential is taken directly from the experimental parameters in Ref. [17], not fitted to the results. The snapshots and density profiles showing one- and two-droplet arrays are convincing. The agreement with experiment for positive ξ on droplet size and formation threshold is real. This is a legitimate first-principles calculation on a system where mean-field theory breaks down, and the application of PIGS to the fully anisotropic double-shielded molecular potential is genuinely new. Credit is due there. The group has a track record with PIGS on dipolar systems and the methodology is sound for what it claims to do. The ξ → −ξ symmetry analysis is also correct and the honest reporting of disagreement with experiment for negative ξ is appropriate. What does not hold up: the supersolid claim. The paper offers two pieces of evidence for superfluidity. First, non-vanishing inter-droplet density in Fig. 5. But PIGS is a ground-state method — quantum zero-point motion will always produce some inter-droplet density regardless of phase coherence. This is necessary but not sufficient. Second, an upper bound fs ≤ 0.8 ± 0.1 from Leggett's relation, which the authors themselves call 'difficult to evaluate accurately.' An upper bound of 0.8 on a quantity ranging from 0 to 1 is fully consistent with fs = 0. It does not constitute evidence of superfluidity. The standard tool here is the winding number estimator, which gives a direct superfluid fraction in path integral methods. The paper does not report it. The contrast-based estimate of ~0.51 is acknowledged by the authors as measuring particle overlap, not superfluidity. So the supersolid claim is unsupported as written. This is the central soft spot and it is load-bearing for the headline. Minor issues: N=1500 only, no thermodynamic limit analysis, no bead count or time step stated, no code or data shipped. These limit reproducibility but are secondary. The self-bound droplet result is solid and worth publishing. The supersolid framing needs to be either backed by a direct winding number calculation or substantially softened. I'd send this to a serious referee who can evaluate whether the winding number estimator is feasible in this geometry and push the authors to either include it or reframe. The paper deserves peer review — the droplet physics is real and the method is appropriate. But the supersolid claim as it stands should not survive review without additional evidence.","headline":"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.","tokens_in":9128,"tokens_out":632,"would_cite":false,"duration_ms":74194,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Microwave-shielded polar molecules form self-bound droplet arrays","keywords":[],"falsifier":"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.","tokens_in":8389,"feed_emoji":"💧","tokens_out":1004,"duration_ms":137146,"temperature":0.7,"pith_summary":"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.","feed_headline":"Microwave-shielded molecules form self-bound droplet arrays","feed_subtitle":"Monte Carlo simulations of NaCs polar molecules show self-binding above 3 degrees of ellipticity, with hints of a supersolid state without a","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Double microwave shielding yields trap-free molecular droplet arrays","Ultracold NaCs molecules form trap-free superfluid droplet arrays","Anisotropic microwave shielding drives self-bound droplet formation","Self-bound droplet arrays form in microwave-shielded NaCs molecules","Tunable microwave shielding creates trap-free superfluid droplet arrays"],"cache_read_input_tokens":0,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Double microwave shielding yields trap-free molecular droplet arrays","Ultracold NaCs molecules form trap-free superfluid droplet arrays","Anisotropic microwave shielding drives self-bound droplet formation","Self-bound droplet arrays form in microwave-shielded NaCs molecules","Tunable microwave shielding creates trap-free superfluid droplet arrays"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":1666,"prompt_tokens":511,"completion_tokens":1155,"prompt_tokens_details":null},"tokens_in":511,"tokens_out":1155,"duration_ms":64187,"temperature":1.0,"reasoning_tokens":1124,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T15:40:04.647029+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"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.","supporting_citations":[],"review_version":1}