{"id":"feac42eb-9537-4a65-8bf6-d38eb17c27b3","arxiv_id":"2411.18395","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Numerical simulations of dipolar Bose gas quenches show Kibble-Zurek power-law scaling of freeze-out time and correlation length, yielding critical exponents nu = 0.57(1) and z = 1.05(2).","lead":"Using computer simulations, this paper studies what happens when a gas of ultra-cold magnetic atoms is quickly tuned across a phase transition between a superfluid and a supersolid. It finds that the speed of the tuning controls how long the transition is delayed and how far the solid order extends, matching the predictions of the Kibble-Zurek mechanism.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Internal inconsistency: reported ν=0.57(1) does not follow from ζ=0.352(3) and νKZ=0.335(3) via the paper's own KZM relation; the correct value is ν≈0.52.","rationale":"The reader's weakest assumption concerned the separability/mean-field modeling, which is an external validity issue. The most load-bearing concern is internal: the paper's quoted critical exponents are not derivable from its own measured scaling laws using the standard KZM relations, indicating a likely arithmetic error or a misattribution of the correlation-length exponent. This is decidable by pure algebra and does not require new simulations. The corrected ν≈0.52 would actually sit closer to the mean-field value ν=1/2, so the qualitative thesis may survive, but the paper as written contains a central quantitative inconsistency that must be corrected before acceptance. Hence the verdict should remain CONDITIONAL, with the specific condition of fixing and re-verifying the exponent extraction.","tokens_in":24992,"tokens_out":11915,"duration_ms":106392,"concrete_test":"Recompute ν and z from the reported ζ=0.352(3) and νKZ=0.335(3) by inverting Eq. (1): z = ζ/νKZ and ν = νKZ/(1−ζ). Compare the resulting ν to Eq. (13); if it deviates by more than the quoted error, the paper's critical exponent extraction contains a numerical or attribution error. Additionally, verify which (ν,z) pair is produced using the defect-based νKZ=0.37(2) from Sec. III C.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central quantitative output is the extracted critical exponents ν=0.57(1) and z=1.05(2) in Sec. III D. These are derived from the measured KZM exponents τ̂∝τ_Q^ζ with ζ=0.352(3) and X∝τ_Q^νKZ with νKZ=0.335(3). Using the paper's own formulas in Eq. (1), ζ = zν/(1+zν) and νKZ = ν/(1+zν), the two measured exponents imply z = ζ/νKZ = 1.05(2) and ν = νKZ/(1−ζ) = 0.335/(1−0.352) = 0.517(7). This differs from the reported ν=0.57(1) by about 7σ, far beyond the stated error bars. Notably, the reported ν=0.57 would instead require νKZ≈0.37, which is the defect-density exponent from Sec. III C, not the correlation-length exponent from Sec. III B. The text explicitly attributes the extraction to subsections III A and III B, so this is an internal mathematical inconsistency in the headline result, not merely a modeling assumption. The nominal z is consistent with the stated inputs, but ν is not. Because ν and z are the paper's principal claims, this inconsistency directly undermines the central result as reported.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a truncated-Wigner extended Gross-Pitaevskii study of linear quenches of the s-wave scattering length across the continuous superfluid-to-supersolid transition of an elongated 164Dy gas. The authors define a freeze-out time from the delay in the drop of the superfluid fraction, extract a frozen correlation length from a fit to the density-density correlation function, and count crystal-phase defects. They find power-law scalings τ̂ ∝ τ_Q^{0.352(3)}, X ∝ τ_Q^{0.335(3)}, and a defect-count exponent of 0.37(2), and from these they quote critical exponents ν = 0.57(1) and z = 1.05(2), comparing them with Bogoliubov mean-field predictions ν = 1/2, z = 1. The paper also includes checks with quantum-only noise, larger system sizes, and full 3D simulations.","tokens_in":25284,"tokens_out":7158,"duration_ms":64649,"significance":"The topic is timely and the study is in a regime of current experimental interest. The main strengths are the extended dynamical range of τ_Q, the ensemble averaging over hundreds of realizations, and the explicit checks in Appendices A–C that the exponents are robust to the noise model, system size, and, within larger error bars, the dimensionality of the simulation. If the reported scalings survive a corrected analysis, this would be one of the first quantitative Kibble-Zurek exponent extractions for a dipolar supersolid transition and would provide a useful benchmark for experiments. However, the headline exponent extraction contains an internal inconsistency that must be resolved before the quantitative claims can be accepted.","major_comments":[{"comment":"The quoted exponents ν = 0.57(1) and z = 1.05(2) are not consistent with the inputs stated in the text. With ζ_KZ = 0.352(3) and ν_KZ = 0.335(3), Eq. (1) gives z = ζ_KZ/ν_KZ = 1.05(2) and ν = ν_KZ/(1 − ζ_KZ) = 0.335/(1 − 0.352) = 0.517(5). The difference from the reported 0.57 is about ten standard errors, far outside the quoted uncertainty. Because Sec. III D explicitly states that the extraction uses the results of Secs. III A and III B, this is an arithmetic inconsistency in the central result rather than an alternative fitting choice. The abstract, Sec. III D, and Sec. IV should be revised to report the corrected exponents and the correspondingly adjusted comparison with mean-field predictions.","section":"§III D and Eq. (1)"},{"comment":"The paper presents three mutually inconsistent estimates of the KZM correlation-length exponent for the same quenches: ν_KZ = 0.335(3) from g(2) in Sec. III B, 0.275(3) from the half-width of C(x) in Sec. III C, and 0.37(2) from the defect count with Δϕ > π/2 in Sec. III C. Equation (11) predicts that the defect-density exponent equals ν_KZ, and the three values are all presented as probes of the same frozen correlation scale; the spread between C(x) and g(2) is roughly 14σ on the quoted errors. The manuscript acknowledges these discrepancies only qualitatively. The authors should either provide a quantitative account of the systematic differences or present the headline exponents with error bars that reflect this systematic spread.","section":"§III C and Eq. (11)"}],"minor_comments":[{"comment":"The sentence 'Sections III D-III C presents our observations' appears to be a typo; it should read 'Sections III A-III D present' or similar, with the sections in the correct order.","section":"Sec. I"},{"comment":"The caption contains the incomplete phrase 'A giving an approximate scaling'; a verb is missing and the sentence should be rewritten.","section":"Fig. 12 caption"},{"comment":"The text refers to a rescaling of 'the g(1) correlator', but the quantity plotted and discussed in Sec. III B and in Fig. 15 is g(2); please correct the label to avoid confusion.","section":"Appendix E"},{"comment":"The open-symbol convention for quenches that extend past the supersolid phase boundary is used consistently, but the figures would be easier to read if the caption or the main text stated explicitly that open markers correspond to the modified fast-quench protocol described in Sec. III A.","section":"Figs. 3–6"}],"recommendation":"major_revision","confidential_remarks":"The arithmetic inconsistency in §III D is straightforward to fix and is the main obstacle to acceptance. If the corrected analysis yields ν ≈ 0.52(1) with z ≈ 1.05(2), the qualitative conclusion—near-mean-field exponents with clean power-law scaling—would survive, but the text, abstract, and comparison with Bogoliubov theory would need to be updated accordingly. The additional spread among the three ν_KZ estimates should also be addressed, either by identifying the systematics or by enlarging the quoted uncertainties."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a valuable numerical study and the first systematic KZM scaling analysis for the dipolar superfluid-supersolid transition. The raw results are strong: clean power laws for the freeze-out time and correlation length over more than two decades in quench time, sensible robustness checks against quantum noise, system size, and 3D simulations, and an honest discussion of defects and their Poisson statistics. The universal scaling collapse in Appendix E is a nice extra. The comparison to Bogoliubov predictions is appropriate and the conclusion that the transition appears continuous and not (1+1)D XY is reasonable.\n\nBut there is a real problem with the central quantitative claim. The paper quotes ζ = 0.352(3) from Sec. III A and νKZ = 0.335(3) from Sec. III B, and then reports ν = 0.57(1) and z = 1.05(2). Using the paper's own Eq. (1), those two measurements give z = ζ/νKZ = 1.05(2) and ν = νKZ/(1−ζ) = 0.335/(1−0.352) ≈ 0.517(11). That is several error bars away from the reported ν = 0.57(1). The z value is consistent, but ν is not. The reported ν = 0.57 matches what you get from ζ = 0.352 and νKZ = 0.37, the defect-density exponent from Sec. III C, not the correlation-length exponent from Sec. III B. Since the text explicitly says Secs. III A and III B are used, this is a genuine internal inconsistency in the headline result, not a modeling subtlety.\n\nThis matters, but it is also fixable. The corrected ν ≈ 0.52 is actually closer to the Bogoliubov value 1/2, so the overall story—KZM scaling with near-mean-field exponents—survives and arguably becomes cleaner. Still, as posted, the numbers do not follow from the data by the paper's own formulas, and readers would be misled.\n\nThe other soft spots are minor by comparison. The threshold f_s = 0.98 and the Gaussian envelope fit for g(2) are empirical, but the authors test alternatives and the scaling is consistent across them. The separability ansatz in Eq. (3) is standard for dipolar systems; the 3D check has large error bars but does not contradict the 1D results.\n\nBottom line: the paper deserves a serious referee, but the editors should ask the authors to recheck their exponent extraction. This is exactly the kind of error that a careful referee would catch. Recommend: send to peer review, with the expectation of a revision to correct ν.","headline":"Solid numerical KZM study of the dipolar superfluid-supersolid transition, but the headline critical exponents are internally inconsistent with the paper's own scaling relations.","tokens_in":25905,"tokens_out":4210,"would_cite":false,"duration_ms":35798,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulations show Kibble-Zurek scaling governs the superfluid-to-supersolid transition in an elongated dipolar gas, with critical exponents matching mean-field theory.","keywords":["Kibble-Zurek mechanism","supersolid","dipolar quantum gas","critical exponents","roton instability","quench dynamics","extended Gross-Pitaevskii equation","correlation length scaling"],"falsifier":"A direct test would repeat the same linear quenches in full three-dimensional simulations or in a uniform ring-trap experiment with $^{164}$Dy and measure the freeze-out delay and correlation length at the critical crossing over the same quench-rate range. If the power laws depart from $\\hat{\\tau}\\propto\\tau_Q^{0.35}$ and $X\\propto\\tau_Q^{0.34}$ beyond the quoted fit errors—for instance by showing exponential Berezinskii-Kosterlitz-Thouless scaling at slow rates—the central claim would be falsified. A complementary check is to measure how the roton gap closes as a function of $a_s$ in the full model and confirm independently that $z\\nu=1/2$ and $z=1$.","tokens_in":24689,"feed_emoji":"⚛️","tokens_out":11507,"duration_ms":92978,"temperature":0.7,"pith_summary":"This paper asks whether the Kibble-Zurek mechanism governs how a dipolar gas forms a supersolid when the interaction strength is swept linearly across the superfluid-to-supersolid transition. Simulating an elongated gas of dipolar atoms with a quasi-one-dimensional extended Gross-Pitaevskii equation, the authors find that the formation delay (freeze-out time) and the growing crystalline correlation length scale as powers of the quench rate over several orders of magnitude. The measured exponents give critical exponents $\\nu = 0.57(1)$ and $z = 1.05(2)$, close to the mean-field values $\\nu = 1/2$, $z = 1$ that follow from the linear closing of the roton gap. If the paper is right, a single universal scaling description—not an exponential Berezinskii-Kosterlitz-Thouless one—captures the onset of supersolid order, and the transition belongs to a near-mean-field universality class that remains to be determined.","feed_headline":"Supersolid transition obeys Kibble-Zurek scaling in simulations","feed_subtitle":"Quench-rate power laws give critical exponents ν=0.57(1), z=1.05(2), close to mean-field predictions.","key_machinery":"The load-bearing object is the Kibble-Zurek scaling ansatz for a linear quench: as the control parameter approaches its critical value, the relaxation time $\\tau$ and correlation length $\\xi$ diverge as power laws, giving a freeze-out time $\\hat{\\tau}\\propto \\tau_Q^{z\\nu/(1+z\\nu)}$ and a frozen length $\\hat{\\xi}\\propto \\tau_Q^{\\nu/(1+z\\nu)}$. In this system the underlying critical softening is the roton minimum of the Bogoliubov excitation spectrum, which closes as $\\epsilon_{\\mathrm{rot}}\\propto\\sqrt{a_s-a_c}$; the dynamics are generated with the extended Gross-Pitaevskii equation reduced to one dimension through a variational transverse profile. The paper identifies the freeze-out time with the delay before the superfluid fraction drops below 0.98 and the correlation length with the Gaussian envelope width of the density-density correlation function at that time.","core_discovery":"The paper simulates linear ramps of the s-wave scattering length from the uniform superfluid into the supersolid phase in an elongated tube of dipolar atoms, using the reduced extended Gross-Pitaevskii equation with stochastic initial noise and many independent realizations. It claims that the delay in the onset of density modulation after crossing the critical point—the Kibble-Zurek freeze-out time—and the width of the density-density correlation function at that moment both obey clean power laws in the quench time, $\\hat{\\tau} \\propto \\tau_Q^{0.352(3)}$ and $X \\propto \\tau_Q^{0.335(3)}$, over roughly three decades of quench rates. From these exponents the paper extracts $\\nu = 0.57(1)$ and $z = 1.05(2)$, compatible with the Bogoliubov mean-field values $\\nu = 1/2$, $z = 1$, and it finds that the number of crystal-phase defects at freeze-out also follows a power law when the phase jumps defining a defect are sufficiently large. The authors take this as evidence for a continuous transition whose universality class is near mean-field and not that of the (1+1)-dimensional XY model.","pith_inferences":["If the near-mean-field exponents survive a fully self-consistent many-body treatment, the transition may behave as if it were above its upper critical dimension; an exact many-body calculation of $\\nu$ and $z$ for the same Hamiltonian would settle that without relying on the quasi-1D mean-field reduction.","The two correlation-length measures in the paper—the density-envelope exponent 0.335 and the crystal-phase width exponent 0.275—may indicate separate density and phase coherence lengths; tracking both after the freeze-out time would clarify whether one or two order parameters are needed.","Because only the onset of supersolid formation is required, the same protocol could be applied to shorter-lived supersolids and ring traps, and a reverse quench that melts the supersolid should show symmetric freeze-out if the mechanism is Kibble-Zurek rather than coarsening."],"forward_implications":["The freeze-out time and correlation length of supersolid formation are predictable for any linear quench rate in this regime, up to the fast-quench breakdown where domains approach the lattice spacing.","The transition is consistent with a continuous second-order transition with near-mean-field exponents, which would rule out the (1+1)-dimensional XY/BKT universality class for these parameters.","Defects of the supersolid crystal appear as local jumps of the crystal phase; for jumps larger than $\\pi/2$ their count follows the predicted KZM power law and their statistics are Poissonian, so defect counting is a viable experimental probe.","Experimental verification needs only the onset of density modulation, not long-lived supersolid coherence, which relaxes the lifetime requirement for ring-trap experiments.","Very fast quenches break the scale separation between domain size and crystal periodicity, producing a saturation that is itself a universal breakdown of Kibble-Zurek scaling."],"supporting_citations":[{"why":"It supplies the Kibble-Zurek predictions for freeze-out time and domain size versus quench rate that the paper tests.","marker":"[24–27]"},{"why":"It provides the reduced quasi-1D model and Bogoliubov dispersion; the roton gap closes as a square root and linearizes at the zone edge, giving zν=1/2 and z=1.","marker":"[73]"},{"why":"It establishes the phase diagram and identifies the continuous superfluid-to-supersolid transition line at the linear density used for the quenches.","marker":"[15]"},{"why":"It supplies the critical-point spectrum, the Higgs mode interpretation, and the compressibility discontinuity supporting a second-order transition.","marker":"[75]"},{"why":"It provides the infinite-tube phase diagram and the full 3D critical point used in the appendix comparison to three-dimensional simulations.","marker":"[16]"},{"why":"It supplies the stochastic classical-field method used to generate the 500 independent noise realizations per quench.","marker":"[99]"},{"why":"It explains the saturation of defect scaling at fast quench rates as a universal breakdown of Kibble-Zurek scaling.","marker":"[88]"},{"why":"It supplies the Poisson defect-count statistics used to define which phase jumps count as genuine KZM defects.","marker":"[89]"}],"fun_headline_variants":["Kibble-Zurek scaling confirmed in supersolid","Dipolar gas quenches show Kibble-Zurek power laws","Supersolid freeze-out obeys Kibble-Zurek exponents","Universal scaling in supersolid transition simulations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the reduced quasi-one-dimensional extended Gross-Pitaevskii equation, with its fixed variational transverse profile and approximate quantum-fluctuation term, faithfully represents the real three-dimensional many-body dynamics of the dipolar gas across the transition.","fun_headline_variants_meta":{"raw":{"variants":["Kibble-Zurek scaling confirmed in supersolid","Dipolar gas quenches show Kibble-Zurek power laws","Supersolid freeze-out obeys Kibble-Zurek exponents","Universal scaling in supersolid transition simulations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000203,"raw_usage":{"total_tokens":1435,"prompt_tokens":1047,"completion_tokens":388,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":320}},"tokens_in":663,"tokens_out":388,"duration_ms":4378,"temperature":1.0,"reasoning_tokens":320,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:13:26.032760+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would repeat the same linear quenches in full three-dimensional simulations or in a uniform ring-trap experiment with $^{164}$Dy and measure the freeze-out delay and correlation length at the critical crossing over the same quench-rate range. If the power laws depart from $\\hat{\\tau}\\propto\\tau_Q^{0.35}$ and $X\\propto\\tau_Q^{0.34}$ beyond the quoted fit errors—for instance by showing exponential Berezinskii-Kosterlitz-Thouless scaling at slow rates—the central claim would be falsified. A complementary check is to measure how the roton gap closes as a function of $a_s$ in the full model and confirm independently that $z\\nu=1/2$ and $z=1$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the reduced quasi-1D model and Bogoliubov dispersion; the roton gap closes as a square root and linearizes at the zone edge, giving zν=1/2 and z=1."},{"cited_title":"Sinatra, C","cited_arxiv_id":null,"evidence_quote":"It supplies the stochastic classical-field method used to generate the 500 independent noise realizations per quench."},{"cited_title":"Bland, Q","cited_arxiv_id":null,"evidence_quote":"It explains the saturation of defect scaling at fast quench rates as a universal breakdown of Kibble-Zurek scaling."},{"cited_title":"Zeng, C.-Y","cited_arxiv_id":null,"evidence_quote":"It supplies the Poisson defect-count statistics used to define which phase jumps count as genuine KZM defects."}],"review_version":1}