{"id":"741c40ff-7ffe-4290-b418-9ed43388e98b","arxiv_id":"2506.23820","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Elliptically polarized microwaves induce anisotropic dipolar interactions in microwave-shielded NaCs molecules, and path-integral Monte Carlo simulations show a supersolid phase appears at experimentally accessible parameters.","lead":"A theoretical proposal shows that adding an elliptically polarized microwave to an ultracold gas of sodium-cesium molecules creates an anisotropic interaction that can stabilize a supersolid, where molecules order in a crystal while remaining superfluid. Quantum Monte Carlo simulations map the phase diagram and find the supersolid in a parameter range close to current experiments, suggesting a new molecular platform for supersolidity.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The phase-diagram calculation is internally consistent, but the abstract's 'accessible to current experiments' claim is unsecured: no inelastic two-body or three-body loss rate is computed for the proposed dual-microwave parameters, so the lifetime extension is only an assumption.","rationale":"I read the manuscript as a theoretical proposal with a well-posed many-body calculation: the effective potential is derived from a Floquet treatment, checked against adiabatic potentials, and cross-validated with multichannel scattering, and the PIMC with worm algorithm is an unbiased method for the stated Hamiltonian. The supersolid phase in the xi-omega_y plane is therefore a concrete theoretical prediction. The weakest step is the jump from this loss-free pairwise Hamiltonian to the abstract's claim of experimental accessibility. The paper itself states that lower density extends the lifetime, but it never quantifies two-body inelastic loss or three-body recombination for the proposed dual-microwave parameters. SM3 introduces beta_inel and the elastic-to-inelastic ratio gamma but does not give values at the operating point, so the authors' own machinery could settle the question. I also noted the textual condition attached to Eq. (1), 'for |xi| > 15 degrees', while all simulations use xi between about 3 and 5 degrees. Read literally this is a contradiction; however, the sin(2 xi) dependence and the validations at xi = 3 and 6 degrees strongly suggest the intended condition is an upper bound, likely '|xi| less than about 15 degrees', and this should be corrected editorially rather than treated as a scientific flaw. The arbitrary 1% thresholds and missing particle number are secondary because they affect phase-boundary placement, not the qualitative existence of a coherent density-modulated phase. Overall, the reader's conditional verdict is the right weighting: the theoretical supersolid is plausible, but the experimental accessibility claim needs a quantitative loss check.","tokens_in":16691,"tokens_out":10806,"duration_ms":130730,"concrete_test":"Run the Floquet multichannel scattering calculation of SM3 at the manuscript parameters (Omega_sigma=2pi x 7.9 MHz, delta_sigma=-2pi x 8 MHz, Omega_pi=2pi x 6.5 MHz, delta_pi=-2pi x 10 MHz) for xi = 3.9, 4.0, and 4.2 degrees, and report beta_el, beta_inel, and gamma = beta_el / beta_inel. Combine these with the three-body loss coefficient K3 measured for NaCs at n = 1.5 x 10^12 cm^-3 to estimate the lifetime of the supersolid state. If the estimated lifetime is comparable to or shorter than the experimental preparation and imaging time, the 'accessible to current experiments' claim should be weakened or made conditional on additional loss suppression.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim has two parts: (i) the PIMC phase diagram shows a supersolid for a Hamiltonian with the pairwise Veff of Eq. (1), and (ii) this regime is accessible to current experiments. Part (i) is reasonably supported: Eq. (1) is cross-validated against the adiabatic potential and multichannel scattering in SM2 and SM3, and the PIMC is exact for the stated Hamiltonian. The load-bearing weak point is part (ii). Equation (2) contains no loss terms, and SM3 defines the inelastic rate beta_inel but reports no value for the proposed parameters (Omega_sigma=2pi x 7.9 MHz, delta_sigma=-2pi x 8 MHz, Omega_pi=2pi x 6.5 MHz, delta_pi=-2pi x 10 MHz, xi around 4 degrees). The 'significantly extended lifetime' claim is inferred only from the low peak density of about 1.5 x 10^12 cm^-3. However, the proposed dual-microwave configuration with C3,0 approximately 0 is not the demonstrated single-microwave NaCs shielding setup, so the two-body inelastic scattering rate and the three-body recombination rate at this density and temperature are unknown; a larger beta_inel or K3 could make the supersolid lifetime too short to observe. The existence of the thermodynamic supersolid does not depend on this, but the abstract's central promise of experimental accessibility does.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":16919,"tokens_out":10597,"duration_ms":108621,"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":[{"comment":"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.","section":"Eq. (2), SM3"},{"comment":"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.","section":"Many-body quantum phases"}],"minor_comments":[{"comment":"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 ξ').","section":"Effective inter-molecular interaction"},{"comment":"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.","section":"Many-body quantum phases"},{"comment":"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.","section":"Many-body quantum phases"},{"comment":"There are several typographical errors, including 'anistropy' in the Fig. 4 caption and 'ebahnced' in SM3; a careful proofread is recommended.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid theoretical proposal. The main blocker is the unquantified loss rate for the dual-microwave scheme; I would not reject on this basis because the authors have the tools (SM3) to compute β_inel, but they must either supply a number or weaken the abstract's accessibility claim. The missing particle number N is a simple but essential fix. I also noticed the 'for |ξ| > 15°' statement, which looks like a typo; the editor may want to have the authors verify the domain of validity of Eq. (1)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe short version: this is a credible theory proposal with a genuinely new mechanism. Replacing the sigma-polarized microwave with an elliptically polarized one introduces a C3,2 anisotropy in the effective dipolar interaction, and with a pi-field tuned to make C3,0 ~ 0 the PIMC worm simulations for NaCs find a supersolid in a sizeable region of the xi-omega_y plane. The effective potential in Eq. (1) is derived from Floquet theory and cross-checked against the adiabatic potential and multichannel scattering; the agreement is good away from the shape resonance, and the supersolid sits away from it. The phase diagram is an honest PIMC calculation, not a fit. That part holds up.\n\nThe weak point is the experimental claim. The abstract says the supersolid is 'accessible to current experiments' and the paper argues the lower density gives a 'significantly extended lifetime.' But this rests on two unsecured things: the pi-microwave Rabi frequency Omega_pi = 6.5 MHz is assumed, not demonstrated, and the inelastic two-body rate beta_inel for this dual-frequency configuration is never computed. They define it in SM3 but don't report a number for these parameters. So the lifetime extension is an inference from the low peak density, not a computed result. A referee should push on this. If two-body loss is much worse in the non-cylindrical setup, the supersolid may be unobservable in practice. To be clear, this doesn't affect the equilibrium many-body physics; it's an overstatement in the abstract.\n\nMinor points: the phase-boundary criteria (S > 1%, g_pp < 1%) are threshold choices, and the boundaries have no error bars. The coherence curve drops monotonically, so the phase assignments look reasonable, but I'd want particle numbers and extrapolation details in the supplement. The citation pattern is fine; the self-citations are to the group's prior dual-shielding work, which is relevant.\n\nWho's this for? Anyone working on dipolar molecules or supersolid proposals. The anisotropy mechanism is a useful addition, and the phase diagram is a benchmark for future work. It deserves a serious referee. My recommendation: send it to review, but require either a loss-rate estimate or a toned-down claim of experimental accessibility.","headline":"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.","tokens_in":17521,"tokens_out":4930,"would_cite":true,"duration_ms":48245,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["supersolid","polar molecules","microwave shielding","dipolar interaction","path-integral Monte Carlo","NaCs","quantum phases","elliptically polarized microwave"],"falsifier":"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.","tokens_in":16414,"feed_emoji":"❄","tokens_out":3900,"duration_ms":46103,"temperature":0.7,"pith_summary":"The paper proposes that dressing a gas of microwave-shielded polar molecules with an elliptically polarized microwave breaks the cylindrical symmetry of the inter-molecular dipolar interaction, turning it repulsive along one in-plane axis and attractive along the other. Using unbiased path-integral Monte Carlo simulations of trapped NaCs molecules, it maps a finite-temperature phase diagram in the ellipticity–trap-frequency plane and finds a supersolid phase with density modulation along the weak trap axis and measurable inter-droplet coherence. The supersolid appears near an ellipticity of about 4 degrees and for y-direction trap frequencies between about 60 and 120 Hz, parameters within reach of present experiments. If correct, this gives a molecular platform for studying supersolidity at peak densities near 1.5×$10^{12}$ $cm^{-3}$, much lower than atomic dipolar supersolids, which the authors argue should extend the lifetime of the molecular gas.","feed_headline":"Microwave ellipticity creates a molecular supersolid at low density","feed_subtitle":"Simulations find a supersolid window near 4 degrees ellipticity and 60-120 Hz trap frequency, within reach of today's NaCs experiments.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the experimental NaCs condensate setup and microwave parameters used to fix the concrete system in the simulations.","marker":"[5]"},{"why":"Provides the earlier result that microwave-shielded polar molecule condensates have reduced condensate fraction and antibunched correlations, motivating the need for a low-density supersolid scheme.","marker":"[6]"},{"why":"Previous path-integral Monte Carlo study of dual-microwave shielded molecules that this work extends by adding ellipticity.","marker":"[8]"},{"why":"Shows that strong single-microwave shielding produces a self-bound monolayer crystal with vanishing superfluid fraction, the contrast case the present scheme aims to avoid.","marker":"[23]"},{"why":"Lays the foundation of continuous-time worldline Monte Carlo for exact finite-temperature quantum many-body statistics.","marker":"[24]"},{"why":"Introduces the worm algorithm for continuous-space path-integral Monte Carlo, the numerical method used to map the phase diagram.","marker":"[25]"},{"why":"Contains the derivation of the effective potential and its validation against adiabatic and multichannel scattering results, the basis for the many-body Hamiltonian.","marker":"[32]"},{"why":"Provides the three-body recombination measurements of microwave-shielded NaCs that define the loss process the low-density supersolid aims to suppress.","marker":"[33]"},{"why":"Gives the density of an atomic dipolar supersolid, about 5×10^14 cm^-3, used as the benchmark showing the molecular supersolid operates at much lower density.","marker":"[18]"}],"fun_headline_variants":["Elliptical microwave unlocks molecular supersolid","Supersolid predicted in microwave-shielded polar gas","Polar molecules get supersolid via elliptical microwaves","Low-density supersolid from elliptical microwave control","Microwave ellipticity steers polar molecules into supersolid"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Elliptical microwave unlocks molecular supersolid","Supersolid predicted in microwave-shielded polar gas","Polar molecules get supersolid via elliptical microwaves","Low-density supersolid from elliptical microwave control","Microwave ellipticity steers polar molecules into supersolid"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000272,"raw_usage":{"total_tokens":1559,"prompt_tokens":800,"completion_tokens":759,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":416,"completion_tokens_details":{"reasoning_tokens":685}},"tokens_in":416,"tokens_out":759,"duration_ms":8678,"temperature":1.0,"reasoning_tokens":685,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:30:34.217111+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Bigagli, W","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental NaCs condensate setup and microwave parameters used to fix the concrete system in the simulations."},{"cited_title":"Bose-Einstein condensates of microwave-shielded polar molecules","cited_arxiv_id":"2406.06412","evidence_quote":"Provides the earlier result that microwave-shielded polar molecule condensates have reduced condensate fraction and antibunched correlations, motivating the need for a low-density supersolid scheme."},{"cited_title":"Quantum Phases for Finite-Temperature Gases of Bosonic Polar Molecules Shielded by Dual Microwaves","cited_arxiv_id":"2503.02644","evidence_quote":"Previous path-integral Monte Carlo study of dual-microwave shielded molecules that this work extends by adding ellipticity."},{"cited_title":"Self-bound monolayer crystals of ultracold polar molecules","cited_arxiv_id":"2504.02682","evidence_quote":"Shows that strong single-microwave shielding produces a self-bound monolayer crystal with vanishing superfluid fraction, the contrast case the present scheme aims to avoid."},{"cited_title":"Prokof’Ev, B","cited_arxiv_id":null,"evidence_quote":"Lays the foundation of continuous-time worldline Monte Carlo for exact finite-temperature quantum many-body statistics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contains the derivation of the effective potential and its validation against adiabatic and multichannel scattering results, the basis for the many-body Hamiltonian."},{"cited_title":"Stevenson, S","cited_arxiv_id":null,"evidence_quote":"Provides the three-body recombination measurements of microwave-shielded NaCs that define the loss process the low-density supersolid aims to suppress."}],"review_version":1}