{"id":"2776af7f-dcd4-4433-bbe8-df8ae3aab5d5","arxiv_id":"2504.14578","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Numerical simulations predict triangular meron and square skyrmion topological supersolids in a spin-orbit coupled dipolar BEC, with PT-driven current oscillations whose frequency tracks the superfluid fraction.","lead":"This paper uses computer simulations of an ultracold magnetic-atom gas to predict two topological supersolid phases, in which atoms arrange in a crystal while still flowing without friction. It also shows that oscillations driven by balanced gain and loss could reveal how supersolid the gas is, offering a new way to measure a property that has been hard to probe.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper asserts, but never derives, that the driven-current oscillation frequency ω_d maps uniquely to the equilibrium Leggett superfluid fraction; Fig.","rationale":"The reader's weakest-assumption analysis correctly identifies the load-bearing gap: the paper's novelty rests on a quantitative bridge between an equilibrium superfluid fraction and a nonequilibrium observable, but that bridge is only demonstrated as a numerical correlation along one parameter path. I agree with the CONDITIONAL verdict: the equilibrium phase diagram and the driven-dissipative simulations are plausible and based on standard mean-field methods, but the central measurement claim is not established. I considered whether a more fundamental issue exists, such as the finite-trap validity of spontaneous translational symmetry breaking or the use of a non-Hermitian PT drive in a dissipative condensate, but those are secondary: even granting the equilibrium phases and the driven dynamics, the paper still lacks a demonstrated one-to-one relation between ω_d and f_NCRI. The concrete test of varying f_NCRI along an independent parameter path would directly settle whether the correlation in Fig. 3(d) reflects a causal mapping or merely a coincidental co-variation with κ. No independent support such as analytic derivation, machine-checked proof, or released code is present to compensate for this gap. Therefore the reader's CONDITIONAL verdict is appropriate and should not be changed.","tokens_in":11576,"tokens_out":6136,"duration_ms":62272,"concrete_test":"Using the same GP solver, compute equilibrium f_NCRI and the driven-current ω_d for at least two independent paths in (κ, εdd) that pass through the same f_NCRI values, for example the path in Fig. 3(d) (εdd=0.02, κ varied) and a second path with κ fixed while εdd is varied so that f_NCRI takes matching values. Overlay all points on a single f_NCRI-versus-ω_d plot; if points from different paths do not collapse onto one curve, then ω_d is not a function of f_NCRI alone and the mapping in Fig. 3(d) is non-unique. As a secondary control, repeat one supersolid simulation with γ0 varied by a factor of 10; if ω_d shifts, the frequency is drive-dependent rather than a superfluid-fraction diagnostic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central measurable claim is that the PT-driven current oscillation frequency ω_d is 'directly linked' to the equilibrium f_NCRI of Eq. (2), enabling quantitative supersolidity probing. The paper provides no derivation of this link. Equation (3) introduces a localized gain/loss pair (σ=0.1, x=±0.8, γ0=0.02), and Eq. (4) defines the current, but no equation connects the response pole to the angle-averaged Leggett integral. Figure 3(d) plots f_NCRI against ω_d for a single parameter path (εdd=0.02, κ varied), so the apparent relationship could equally be a correlation between ω_d and κ, the lattice constant, or the crystalline order parameter C, all of which vary together. The text itself concedes that damping 'depends not only on the superfluid fraction ... but also on the lattice symmetry,' indicating that the nonequilibrium response encodes more than f_NCRI. Without a control in which f_NCRI is varied while lattice stiffness, trap, and gain/loss geometry are held fixed—or an analytic linear-response derivation—the claim that ω_d quantifies supersolidity is not established. The absence of code or numerical details further prevents an independent check of Fig. 3(d).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies a two-component spin-orbit-coupled dipolar Bose-Einstein condensate in a quasi-2D harmonic trap. Using Gross-Pitaevskii (GPE) energy minimization, the authors identify a phase diagram with a single-skyrmion superfluid at weak spin-orbit coupling, a triangular meron supersolid at weak dipolar interaction, and a square skyrmion supersolid at strong dipolar interaction. They report a first-order superfluid-to-meron-supersolid transition and a second-order superfluid-to-skyrmion-supersolid transition, characterized by crystalline order and the nonclassical rotational inertia fraction f_NCRI. They then introduce parity-time-symmetric gain/loss terms and observe damped current oscillations, claiming that the oscillation frequency ω_d is directly linked to f_NCRI and can be used to quantitatively probe supersolidity. The manuscript also discusses experimental feasibility with magnetic atoms such as 52Cr, 164Dy, and 168Er.","tokens_in":11861,"tokens_out":4255,"duration_ms":40407,"significance":"If the central mapping between the driven current oscillation frequency and the equilibrium superfluid fraction were established, this work would offer a new nonequilibrium probe of supersolidity in dipolar quantum gases, complementing recent experiments on superfluid fraction measurement. The phase diagram with two distinct topological supersolids, including a square skyrmion lattice that is uncommon in single-component dipolar systems, is a valuable contribution. The paper makes specific, falsifiable predictions for transition order and for the behavior of current oscillations. However, the central quantitative claim about ω_d and f_NCRI is currently only a numerical correlation, and the numerical evidence lacks convergence tests and code; these weaknesses limit the present significance.","major_comments":[{"comment":"The claim that ω_d is 'directly linked' to f_NCRI is not derived. Equations (3) and (4) introduce the gain/loss profile and the current, but no equation or argument connects the response frequency to the equilibrium Leggett formula in Eq. (2). Figure 3(d) plots f_NCRI versus ω_d for a single parameter path (ε_dd=0.02 with κ varied), so ω_d, κ, the lattice constant, and the crystalline order C all vary together. The text itself concedes that damping 'depends not only on the superfluid fraction ... but also on the lattice symmetry,' which indicates that the nonequilibrium response encodes more than f_NCRI. Without a control in which f_NCRI is varied while lattice stiffness, trap, and gain/loss geometry are held fixed, or an analytic linear-response derivation, the central abstract and conclusion claim that 'mapping the oscillation frequency to the superfluid fraction' quantifies supersolidity is not established.","section":"Nonequilibrium probing via driven-oscillation (Fig. 3)"},{"comment":"The order of the two phase transitions is inferred from discontinuous or continuous changes in C and f_NCRI, but no metastability, energy crossing, or finite-size scaling is shown to support the first-order assignment, and the claimed 'vanishing cubic term' explaining the second-order transition is only asserted with a citation [58]. A definitive assignment of transition order for a GPE calculation would require either a demonstration of hysteresis in the variational parameters or an analysis of energy barriers; the present evidence is suggestive rather than conclusive.","section":"Topological supersolid phase transitions (Fig. 1)"},{"comment":"No numerical details are provided: the grid size, time step, imaginary-time propagation scheme, convergence criteria, particle number N (or the way gN is set), and the procedure for extracting ω_d are all absent. This prevents an independent check of the phase boundaries and of Fig. 3(d). The omission of the Lee-Huang-Yang correction is justified only by a footnote stating that it 'has been verified numerically' [reference 47], but no verification data are shown; this is particularly relevant because the phase diagram extends to ε_dd=0.4, which is not obviously the 'weak dipolar interaction' regime mentioned in the footnote.","section":"Numerical methods and reproducibility"}],"minor_comments":[{"comment":"The definitions of the averages in Eq. (2), ⟨1/ρ_s(r)⟩ and ⟨ρ(r)⟩, are not specified; moreover, the formula f_NCRI = 1/(⟨1/ρ_s(r)⟩⟨ρ(r)⟩) appears dimensionally inconsistent as written, since the left-hand side is dimensionless but the right-hand side involves products of densities with possibly different powers. Please clarify the integration domains and the normalization of ρ_s(r).","section":"Eq. (2)"},{"comment":"The text describes J_x(t) in the singly-skyrmion phase as 'undamped', but the vertical scale in Fig. 3(a) is much smaller than in (b) and (c) and the trace shows a slight decay; please quantify the damping rate or adjust the wording.","section":"Fig. 3(a) and text"},{"comment":"Reference [47] is a footnote embedded in the reference list rather than a standard citation; if the LHY correction is asserted to be negligible based on numerical checks, the check should be shown in the main text or a supplement, not only asserted in a footnote.","section":"Reference [47]"},{"comment":"The phase diagram in Fig. 1(a) does not indicate how the critical κ_c values were determined from the numerical data; a short description of the criterion (e.g., where C first exceeds a threshold) would improve reproducibility.","section":"Fig. 1(a)"}],"recommendation":"major_revision","confidential_remarks":"The core difficulty is the unsupported mapping between ω_d and f_NCRI; this is a central claim that should be either derived or backed by a decisive numerical control. The absence of numerical details and convergence tests is concerning for a paper whose main evidence is computational. If the authors can provide the missing derivation or control, and add the required numerical documentation, the work could be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new content here is the square skyrmion supersolid phase and the proposed non-equilibrium probe of supersolidity via PT-driven current oscillations. The phase diagram is plausible: the distinction between the triangular meron lattice and the square skyrmion lattice is clearly drawn, and the second-order character of the skyrmion-supersolid transition is supported by a vanishing cubic term. The non-equilibrium setup is concrete and experimentally grounded.\n\nThe soft spot is the central mapping. The abstract and conclusion state that ω_d is 'directly linked' to f_NCRI, but no derivation is given. Figure 3(d) is a single parameter path at fixed ε_dd with κ varied, so ω_d could just as well be tracking lattice constant or crystalline order. The authors concede damping depends on lattice symmetry as well as f_NCRI, meaning the response is not a clean readout of the superfluid fraction. Without a control that varies f_NCRI while holding lattice stiffness, trap, and gain/loss geometry fixed—or a linear-response argument—the measurement claim is not established. This is not a circularity problem: f_NCRI is computed from an equilibrium density integral and ω_d from real-time dynamics, so the correlation is not manufactured, but it also is not a derivation.\n\nProportionately minor issues: no code, convergence tests, or error estimates, so the numerics can't be checked independently. The novelty statement could be sharper—the model is shared with Ref. [41] and the triangular meron phase is in Ref. [48]; the square skyrmion supersolid is the distinct new ground state and deserves explicit billing. The LHY omission is noted with a claim of numerical verification in Ref. [47], but without details. Citation pattern is fine, and the experimental parameters are realistic.\n\nI think this deserves a serious referee. The right referee will push for a derivation or a much better control calculation, code release, and more careful wording about what the oscillation frequency actually measures. I would not cite the non-equilibrium claim in its current form, but the square skyrmion phase is likely to be useful in related work.","headline":"Plausible new square-skyrmion supersolid phase; the claimed ω_d–f_NCRI mapping is a numerical correlation, not a derivation, and needs controls before the probe can be trusted.","tokens_in":12394,"tokens_out":3345,"would_cite":false,"duration_ms":29103,"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":"Spin-orbit coupled dipolar condensates host topological supersolids whose superfluid fraction can be read from driven current oscillations.","keywords":["topological supersolid","spin-orbit coupling","dipolar Bose-Einstein condensate","meron lattice","skyrmion lattice","nonclassical rotational inertia","PT-symmetric dissipation","nonequilibrium dynamics"],"falsifier":"Run the same numerical calculation at fixed superfluid fraction f_NCRI while changing the gain-loss positions, widths, or amplitude, or while changing the trap anisotropy; if the oscillation frequency ω_d shifts by any amount while f_NCRI stays constant, the claimed one-to-one mapping between ω_d and the superfluid fraction is not unique and the proposed probe would require a separate calibration for each geometry.","tokens_in":11379,"feed_emoji":"🌀","tokens_out":3476,"duration_ms":32264,"temperature":0.7,"pith_summary":"This paper argues that a two-component spin-orbit coupled dipolar Bose-Einstein condensate can host topological supersolids, phases where crystalline order, superfluidity, and chiral spin textures coexist. It identifies two distinct transitions: a first-order transition from a single-skyrmion superfluid to a triangular meron supersolid, and a second-order transition from that superfluid to a square skyrmion supersolid. The paper further claims that under parity-time symmetric gain and loss, the condensate develops an oscillating current whose frequency tracks the equilibrium superfluid fraction, so that a nonequilibrium measurement can quantify supersolidity. This matters because directly measuring the superfluid fraction in supersolids has been an outstanding experimental challenge, and because the predicted skyrmion and meron textures could be useful for spintronic or topological devices.","feed_headline":"Oscillating currents could measure supersolidity in dipolar condensates","feed_subtitle":"Simulations show skyrmion and meron supersolids whose driven current frequency tracks the superfluid fraction, offering a direct probe.","key_machinery":"The central machinery is the mean-field energy functional of a Rashba spin-orbit coupled dipolar condensate, together with the angle-averaged Leggett formula for the nonclassical rotational inertia fraction, f_NCRI = 1/(⟨1/ρ_s(r)⟩⟨ρ(r)⟩), and a PT-symmetric driven-dissipative term that creates localized gain and loss. The spin-orbit coupling twists the spin texture into skyrmions or merons, the dipolar interaction selects the lattice symmetry, and the driven current J(t) responds to the resulting superfluid fraction through oscillations whose frequency ω_d is extracted numerically. The claimed connection between ω_d and f_NCRI is carried by the numerical results in Fig. 3(d), which show a monotone relation between the two quantities across the supersolid phases.","core_discovery":"The paper's central claim is that tuning the dipolar-to-contact interaction ratio and spin-orbit coupling strength in a two-component dipolar BEC produces a single-skyrmion superfluid at weak spin-orbit coupling, then either a triangular meron supersolid with alternating meron-antimeron pairs or a square skyrmion supersolid with unit topological charge as the coupling grows. The superfluid-to-meron-supersolid transition is first order, while the superfluid-to-skyrmion-supersolid transition is second order, with both transitions visible in the nonclassical rotational inertia fraction and in the nonequilibrium current response. Under PT-symmetric localized gain and loss, the current oscillation frequency ω_d shows discontinuous jumps at the phase boundaries and correlates with the superfluid fraction f_NCRI, leading the authors to propose that mapping oscillation frequency to superfluid fraction via driven nonequilibrium currents provides a new way to measure supersolidity.","pith_inferences":["If the f_NCRI-to-ω_d mapping is robust, the same nonequilibrium protocol could be adapted to other supersolid platforms, including single-component dipolar supersolids, giving a generic tool for measuring superfluid fraction from sloshing or oscillating currents.","Because the Leggett formula is angle-averaged, the mapping may depend on lattice orientation and symmetry; a systematic scan of gain-loss position and trap anisotropy could reveal whether ω_d carries extra geometric information beyond f_NCRI.","The second-order transition to the square skyrmion lattice is tied to orthogonal pairs of modulation wavevectors; an analogous construction in three dimensions or with different spin-orbit couplings might stabilize other square or rectangular topological supersolids.","A concrete extension would be to test whether the oscillation frequency remains fixed when f_NCRI is held constant but the gain-loss amplitude, width, or positions are varied; if it changes, the proposed measurement would need a calibration curve for each experimental geometry."],"forward_implications":["The oscillation frequency and its discontinuities could serve as an in-situ, nondestructive probe of both the superfluid-to-supersolid transition and the transition between meron and skyrmion supersolids.","A measured relation between ω_d and f_NCRI would give a quantitative estimate of the superfluid fraction in a dipolar supersolid without requiring direct rotational-inertia experiments.","The damping of the current oscillations depends on lattice symmetry as well as superfluid fraction, so the time trace of J(t) could distinguish triangular meron lattices from square skyrmion lattices.","The predicted phases appear accessible with ultracold magnetic atoms such as chromium, dysprosium, or erbium under Raman-induced spin-orbit coupling, making the proposal experimentally testable with current technology."],"supporting_citations":[{"why":"Defines the nonclassical rotational inertia fraction used to quantify supersolidity in Eq. (2).","marker":"[4]"},{"why":"Reports an experimental measurement of nonclassical rotational inertia in a dipolar supersolid, the challenge this paper aims to address with a nonequilibrium probe.","marker":"[19]"},{"why":"Provides a recent nonequilibrium method for probing supersolid properties whose approach the paper directly parallels.","marker":"[20]"},{"why":"Predicts a hexagonal skyrmion lattice in spin-orbit coupled condensates without dipolar interactions, the baseline for the meron and skyrmion supersolids found here.","marker":"[48]"},{"why":"Supplies the two-component dipolar condensate model and the dipole-dipole interaction potential used in the energy functional.","marker":"[40]"},{"why":"Gives the angle-averaged density formula underlying the NCRIF calculation in Eq. (2).","marker":"[57]"},{"why":"Provides the non-Hermitian driven-dissipative Hamiltonian form used to introduce PT-symmetric gain and loss.","marker":"[65]"}],"fun_headline_variants":["Current oscillations gauge supersolid fraction","Non-equilibrium currents map supersolid fraction","Driven currents reveal superfluid fraction in dipolar BEC","Meron and skyrmion supersolids probed by damped currents","Probing supersolidity with driven currents in spin-orbit BEC"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the frequency of the driven current oscillation is controlled by the equilibrium superfluid fraction, so that measuring that frequency gives a direct readout of supersolidity; if the frequency instead responds mainly to lattice stiffness, trap confinement, or the particular arrangement of gain and loss, the proposed measurement would not uniquely quantify the superfluid fraction.","fun_headline_variants_meta":{"raw":{"variants":["Current oscillations gauge supersolid fraction","Non-equilibrium currents map supersolid fraction","Driven currents reveal superfluid fraction in dipolar BEC","Meron and skyrmion supersolids probed by damped currents","Probing supersolidity with driven currents in spin-orbit BEC"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000842,"raw_usage":{"total_tokens":3666,"prompt_tokens":943,"completion_tokens":2723,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":559,"completion_tokens_details":{"reasoning_tokens":2642}},"tokens_in":559,"tokens_out":2723,"duration_ms":16394,"temperature":1.0,"reasoning_tokens":2642,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:45:54.844063+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same numerical calculation at fixed superfluid fraction f_NCRI while changing the gain-loss positions, widths, or amplitude, or while changing the trap anisotropy; if the oscillation frequency ω_d shifts by any amount while f_NCRI stays constant, the claimed one-to-one mapping between ω_d and the superfluid fraction is not unique and the proposed probe would require a separate calibration for each geometry.","supporting_citations":[{"cited_title":"Tanzi, J","cited_arxiv_id":null,"evidence_quote":"Reports an experimental measurement of nonclassical rotational inertia in a dipolar supersolid, the challenge this paper aims to address with a nonequilibrium probe."},{"cited_title":"Biagioni, N","cited_arxiv_id":null,"evidence_quote":"Provides a recent nonequilibrium method for probing supersolid properties whose approach the paper directly parallels."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Predicts a hexagonal skyrmion lattice in spin-orbit coupled condensates without dipolar interactions, the baseline for the meron and skyrmion supersolids found here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-component dipolar condensate model and the dipole-dipole interaction potential used in the energy functional."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the non-Hermitian driven-dissipative Hamiltonian form used to introduce PT-symmetric gain and loss."}],"review_version":1}