{"id":"f6e0fcf7-e591-4fe9-b8a7-e871ec0a8c38","arxiv_id":"1909.00681","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Stronger quenches in a spinor Bose-Einstein condensate produce finer spin-domain patterns and suppress long-wavelength excitations, as confirmed by experiment and predicted by Bogoliubov analysis.","lead":"In a cloud of ultracold sodium atoms, physicists suddenly changed the magnetic energy in the system and watched the atoms' spin patterns reorganize. The experiment maps out how the resulting spin structures change from large irregular patches to fine speckles as the quench gets stronger, providing a direct test of a 2007 and 2010 theory.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Calibration of q̃ is the load-bearing hinge: a wrong scattering length or density factor shifts the q̃~1 and q̃>2 boundaries and can erase the strong-quench regime.","rationale":"I could not find a fatal flaw in the experiment itself; the seeded-sample protocol is a direct probe of k=0 stability and the imaging data support a qualitative crossover. The load-bearing uncertainty is exactly the calibration of q̃, which the reader also flagged. Since the threshold q̃>2 sits at the edge of the accessible range and depends on a recently updated scattering length and an effective-density choice, the quantitative part of the claim is conditional on those values. A sensitivity check would settle it, and the conditional verdict is appropriate.","tokens_in":10614,"tokens_out":20418,"duration_ms":205573,"concrete_test":"Perform a two-parameter sensitivity reanalysis of Figs. 4 and 5: recompute q̃ for all shots using (i) as=0.823 a0 instead of 1.88 a0 and (ii) n̄=4/5 n0 instead of 2/3 n0. Then check whether the seeded-sample suppression onset and the t1 minimum remain at the recalibrated q̃≈2 and q̃≈1. If changing as moves the highest data point below 2, or if changing n̄ shifts the suppression boundary from 2 to 2.4, the strong-quench claim loses quantitative support; if the boundaries track the recalibrated q̃ as predicted by the Bogoliubov dispersion, the calibration concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Sec. IV.C) is that long-wavelength instability is strongly suppressed for q̃>2 and that this agrees with the Bogoliubov crossover at q̃~1 and q̃~2. All thresholds are expressed in q̃=qf/(c2 n̄), so the entire comparison rests on two calibration inputs: as=1.88 a0 in footnote 2 (replacing the prior 0.823 a0 used in Refs. [5,7]) and n̄=2/3 n0 in Sec. III. A factor-2.28 error in as would compress the quoted range 0.12<q̃<4.4 to 0.05<q̃<1.93, putting every 'strong-quench' data point below the q̃=2 boundary and removing the predicted k=0 stabilization. The density factor is a second independent hazard: for a quasi-2D Thomas-Fermi column, the relevant effective density is not automatically the linear axial average 2/3 n0; a density-squared-weighted average (∫n²dz/∫n dz = 4/5 n0 for a TF profile) would shift all q̃ values by 20% and move the suppression threshold to q̃≈2.4. The paper says n̄=2/3 was verified by spin-sound measurements, but Fig. 5 shows t1,t2 without error bars, so the quantitative comparison with Γ^{-1}(q̃) does not currently bound this calibration uncertainty. This is a correctness risk, not an internal inconsistency: the qualitative imaging data support a crossover, but the specific boundaries q̃~1 and q̃>2 are only as solid as this calibration.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports an experimental study of the early-time dynamics of a quasi-two-dimensional spin-1 antiferromagnetic 23Na Bose-Einstein condensate after a sudden quench from the easy-plane polar (EPP) phase to the easy-axis polar (EAP) phase. The authors introduce the dimensionless quench strength q̃ = qf/(c2 n̄) and observe that the spatial scale of the initially generated mz = 0 spin domains decreases with increasing q̃, with an array-like pattern appearing for q̃ around unity and a speckled pattern for q̃ > 2. They also prepare a \"seeded\" sample containing a residual mz = 0 population to selectively enhance zero-momentum spin excitations, and show that the resulting early enhancement—a collective spin oscillation and a reduced onset time t1—weakens with increasing q̃ and nearly disappears for q̃ > 2. The observations are interpreted as consistent with the Bogoliubov dynamic instability of the initial EPP state, in which the k = 0 transverse mode is unstable only for q̃ < 2 and the most unstable mode acquires finite momentum for q̃ > 1.","tokens_in":10954,"tokens_out":9635,"duration_ms":83247,"significance":"If the calibration of q̃ is sound, the experiment provides a clear qualitative confirmation of the Bogoliubov picture across three quench regimes, and the seeded-sample experiment is a direct, falsifiable test of the predicted stabilization of the k = 0 mode for q̃ > 2. A notable strength is that no free parameters are fitted to force agreement: the theory curve comes from the standard Bogoliubov dispersion, and the experimental variation is controlled by the same q̃ parameter. The result would be a useful benchmark for spinor-BEC quench dynamics and supports the idea of using quenches as a magnon thermometer. However, the central quantitative claims are threshold statements in q̃, and the comparison is only as solid as the calibration of q̃ through the scattering length and the effective density; this is where the manuscript currently needs the most work.","major_comments":[{"comment":"The quantitative quench-strength axis is q̃ = qf/(c2 n̄), and every threshold claim (the crossover at q̃ ≈ 1 and the stabilization of k = 0 for q̃ > 2) is stated in these units. The two inputs to q̃ are not accompanied by uncertainty or sensitivity analysis: as = 1.88 a0 in footnote 2 is a factor 2.28 larger than the value used in the authors' previous work, and the effective-density factor n̄ = 2/3 n0 rests on the unpublished Ref. [23]. A shift in either input moves the q̃ = 1 and q̃ = 2 boundaries and could eliminate the strong-quench regime entirely. In addition, the quoted range 0.12 < q̃ < 4.4 is internally inconsistent with c2 n0 = h × 30.7 Hz and n̄ = 2/3 n0, since qf/h = 82 Hz gives q̃ = 4.0. The authors should provide error bars for t1 and t2 in Fig. 5 and propagate calibration uncertainties through the comparison with the theoretical growth rate.","section":"Sec. III and Fig. 5"},{"comment":"The caption of Fig. 5 states that the solid green line is Γ^(-1)(q̃), but Γ(q̃) is defined in Sec. II as the maximum growth rate over momentum, which is constant for 1 < q̃ < 2. The seeded-sample comparison in Sec. IV.C is instead made against Im[Et(k = 0)] ∝ sqrt(q̃(2 − q̃)), which decreases to zero at q̃ = 2. These two quantities have different q̃-dependence between q̃ = 1 and q̃ = 2, so the legend must specify which curve is actually plotted. As printed, the reader cannot verify the claimed agreement of the seeded t1 data with the k = 0 growth rate.","section":"Sec. IV.C and Fig. 5"},{"comment":"For the seeded sample, η(t) exhibits a large-amplitude collective oscillation during the early stage, and t1 is defined as the first time at which η = 0.2. If this crossing occurs while η is still following the oscillatory pulse rather than the exponential instability growth, then t1 is not a direct measure of the k = 0 growth rate and the comparison with Im[Et(k = 0)] is not justified. The authors should show, for example by fitting the exponential envelope or by identifying that the crossing occurs in the post-pulse growth stage, that the seeded t1 values are set by the instability growth rate and not by the seed-induced oscillation.","section":"Sec. IV.C"}],"minor_comments":[{"comment":"The abstract says the quench strength is varied 'up to 4', while Sec. III states the range as 0.12 < q̃ < 4.4; these numbers should be reconciled after the calibration issue is resolved.","section":"Abstract and Sec. III"},{"comment":"The middle image in Fig. 3 is said to show 'the same data shown in Fig. 1(c)', but the corresponding data appear in Fig. 2(c); the cross-reference should be corrected.","section":"Fig. 3 caption"},{"comment":"The phrase 'free from a proportional factor' for the Γ^(-1) curve in Fig. 5 is unclear; it should be replaced by an explicit statement of the vertical normalization used when overlaying t1 on the inverse growth rate.","section":"Sec. IV.B"},{"comment":"The claim that n̄ = 2/3 n0 is 'experimentally verified' via the spin-sound measurement relies on Ref. [23], which is an unpublished preprint; the verification should be described in enough detail for the reader to judge the density-averaging procedure, or an independent calibration should be provided.","section":"Sec. III"}],"recommendation":"major_revision","confidential_remarks":"The central idea and the seeded-sample test are well suited to the journal, but the quantitative comparison hinges on two calibration inputs that are not currently supported by error analysis in the manuscript. The switch from as = 0.823 a0 to as = 1.88 a0 and the use of n̄ = 2/3 n0 from an unpublished preprint should be verified independently, and the Fig. 5 curve-definition issue should be fixed, before the threshold claims q̃ ≈ 1 and q̃ > 2 can be accepted as quantitative."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this one. The paper is a careful experimental study of the EPP-to-EAP quench in a spin-1 sodium BEC, and it does something new: it sweeps the dimensionless quench strength q̃ from 0.12 to 4.4 and shows the early-time dynamics crossing over from long-wavelength instability (irregular domains) to short-wavelength, speckled structures, with the long-wavelength instability strongly suppressed for q̃>2. The seeded-sample experiment—preparing a condensate with a deliberately enhanced m_F=0 component—is a genuinely clever way to test the stability of the k=0 mode. The amplitude of the collective spin oscillation decreases with q̃ and vanishes above q̃~2, which is a direct, qualitative confirmation of the Bogoliubov prediction.\n\nThe theory is not new (Lamacraft 2007, Matuszewski 2010), and the paper doesn't pretend otherwise. No free parameters are fitted. The imaging and the η(t) curves are consistent with the instability picture.\n\nThe soft spot is the calibration. The dimensionless quench strength is q̃=q_f/(c2 n̄), and the paper changed the sodium scattering length from 0.823 a0 (used in their earlier work) to 1.88 a0, citing Knoop et al. That's a factor 2.28. If the old value were right, every data point would sit below q̃=2 and the claimed suppression of the k=0 mode would vanish. They also use an effective density n̄=2/3 n0, the linear axial average, whereas a density-weighted average for a TF profile would give 4/5 n0—a 20% shift in all q̃ values. They say the density factor was verified by spin-sound measurements, but the reference is an unpublished arXiv preprint, and Fig. 5 shows no error bars on t1 and t2, so the comparison against Γ^{-1}(q̃) doesn't bound the calibration uncertainty. They are upfront about the imaging resolution limit in footnote 1, which is why the intermediate-regime comparison stays qualitative.\n\nNone of this is an internal contradiction. The qualitative crossover is real, and the central argument—that the early-time dynamics is governed by the momentum-dependent instability of the initial state—holds up. But the specific thresholds q̃~1 and q̃>2 are only as solid as the calibration. A referee should push for a sensitivity analysis or at least a more transparent discussion of the uncertainty in as and n̄.\n\nFor the spinor-BEC community, this is a useful, citable experimental result. I'd send it to peer review; it deserves careful refereeing, and with reasonable revisions on the calibration, it's publishable. Not a home run, but a solid advance.","headline":"A solid experimental test of the weak-to-strong quench crossover in a spinor BEC; the seeded-sample probe is genuinely clever, but the q̃ calibration (scattering length and effective density) is the hinge and deserves sharper treatment.","tokens_in":11476,"tokens_out":7507,"would_cite":true,"duration_ms":101090,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.75.Mn"],"model":"deepseek-v4-flash","headline":"Quenching a spin-1 antiferromagnetic Bose-Einstein condensate from the easy-plane to the easy-axis polar phase produces a crossover in early spin dynamics that matches the momentum-dependent Bogoliubov instability of the initial state…","keywords":["spinor Bose-Einstein condensate","quantum quench","dynamic instability","Bogoliubov theory","antiferromagnetic spin-1","easy-plane polar state","spin domain formation","quench strength"],"falsifier":"An experiment that measured the momentum-resolved growth of density modulations in the $m_z=0$ component after the quench (for example, by Fourier analysis of absorption images) and found significant low-wavenumber growth for $\\tilde{q}>2$, or a seeded sample whose early dynamics still accelerates at those quenches, would falsify the claim that the zero-momentum mode has stabilized.","tokens_in":10441,"feed_emoji":"⚛️","tokens_out":14997,"duration_ms":115053,"temperature":0.7,"pith_summary":"This paper reports experiments on a sodium-23 spin-1 antiferromagnetic Bose-Einstein condensate quenched from the easy-plane polar (EPP) phase to the easy-axis polar (EAP) phase by suddenly reversing the sign of the quadratic Zeeman energy. It argues that the early-time post-quench dynamics is set by the momentum-dependent Bogoliubov instability of the initial EPP state, and that the character of this instability changes with the dimensionless quench strength $\\tilde{q}=q_f/(c_2 \\bar{n})$: for $\\tilde{q}<1$ the most unstable mode is at zero momentum, for $1<\\tilde{q}<2$ a finite wavenumber dominates, and for $\\tilde{q}>2$ the zero-momentum mode becomes stable, suppressing long-wavelength spin excitations. Using a sample seeded with a small condensate population in the $m_z=0$ component, the authors show that the zero-momentum seed loses its influence for $\\tilde{q}>2$, confirming the predicted disappearance of the long-wavelength instability. If correct, this establishes a quantitative connection between quench dynamics and linear stability analysis in a spinor condensate, and it maps out distinct weak, intermediate, and strong quench regimes that govern defect and pattern formation in these systems.","feed_headline":"Zero-momentum spin instability disappears at quench strength 2","feed_subtitle":"Experiments confirm the Bogoliubov prediction that long-wavelength magnons stop growing in strong quenches.","key_machinery":"The load-bearing object is the transverse magnon branch of the easy-plane polar state, with energy $E_t(k)=\\sqrt{(\\epsilon_k - q)(\\epsilon_k - q + 2c_2 n)}$, where $\\epsilon_k=\\hbar^2 k^2/2m$ is the single-particle energy and $c_2 n$ the spin interaction energy. After the quench to $q>0$ this energy becomes imaginary for a range of wavenumbers, making the corresponding modes dynamically unstable; the maximum imaginary part gives the growth rate $\\Gamma(\\tilde{q})$. As $\\tilde{q}$ increases, the unstable momentum region changes topology: for $\\tilde{q}<1$ it is a disk centered at $k=0$, for $1<\\tilde{q}<2$ the most unstable mode moves to $k_m=k_s\\sqrt{\\tilde{q}-1}$, and for $\\tilde{q}>2$ the $k=0$ mode becomes stable so the unstable region is an annulus with $k_m>\\xi_s^{-1}$. This momentum-space topology change is the mechanism the paper identifies as the origin of the crossover in the observed spin dynamics.","core_discovery":"The central claim is that the crossover from weak to strong quench in the EPP-to-EAP transition is governed by the momentum dependence of the transverse magnon dispersion of the initial easy-plane polar state. For $\\tilde{q}<1$, the most unstable mode is the zero-momentum magnon, producing long-wavelength spin excitations and irregular spin domains; for $1<\\tilde{q}<2$, the most unstable wavenumber becomes $k_m=k_s\\sqrt{\\tilde{q}-1}$, and the early spin texture appears granulated or array-like; for $\\tilde{q}>2$, the $k=0$ mode is no longer unstable, the unstable region in momentum space becomes an annulus, and the pattern turns speckled. The decisive experiment uses a seeded sample with an enhanced $m_z=0$ condensate fraction, i.e., enhanced zero-momentum transverse magnons: the seed accelerates and makes coherent the early growth of $\\eta$ (the $m_z=0$ fraction) for low $\\tilde{q}$, but its effect monotonically fades with increasing $\\tilde{q}$ and nearly vanishes for $\\tilde{q}>2$, in agreement with the prediction that the zero-momentum mode has stabilized. The authors conclude that the observed crossover is consistent with the Bogoliubov description of the dynamic instability of the initial spinor condensate.","pith_inferences":["The seeded-sample technique could be developed into a momentum-resolved quench spectroscopy: by preparing seeds at different wavenumbers, one could map out the full imaginary branch of the magnon dispersion, directly imaging the unstable momentum region and its topology change at $\\tilde{q}=1$ and $\\tilde{q}=2$.","A quantitative comparison with numerical simulations that include nonlinear interactions beyond the linear-growth stage would test whether the Bogoliubov picture survives into the late-time dynamics, especially in the strong-quench regime where the most unstable mode has single-particle character.","If the suppression of long-wavelength modes at high $\\tilde{q}$ holds generally, it implies that quench dynamics in this system cannot be described by a universal scaling with a single length scale (e.g., the spin healing length); instead the dynamics selects a characteristic scale set by $k_m^{-1}$, which may have implications for the interpretation of early-time scaling laws in other symmetry-br"],"forward_implications":["For $\\tilde{q}>2$, the composite-defect nucleation picture (domain walls bounded by half-quantum vortices) that holds in the weak-quench regime no longer applies; the small domain size makes hosting $m_z=\\pm1$ cores energetically prohibitive, so defect formation in strong quenches must be described by a different mechanism.","The reverse quench (EAP to EPP) should exhibit the same crossover behavior as a function of the initial excitation energy, because the two phases share the same magnon dispersion structure; this is directly predicted by the paper.","The seeded-sample protocol provides a selective amplifier for zero-momentum magnons, offering a practical tool to probe the stability of the $k=0$ mode and to calibrate spin-interaction parameters in spinor condensates.","The three identified quench regimes (weak, intermediate, strong) give a classification that future experiments on quench-induced pattern formation in spinor gases should respect when comparing to theory."],"supporting_citations":[{"why":"Supplies the Bogoliubov spectrum of the EPP state, including the two magnon branches used in the instability analysis.","marker":"[10]"},{"why":"Predicted that a finite-wavenumber mode becomes most unstable for quenches with $\\tilde{q}>1$, motivating the search for the crossover.","marker":"[15]"},{"why":"Predicted rotonlike instability and quasi-periodic pattern formation from finite-wavenumber spin excitations, used to interpret the array-like patterns.","marker":"[16]"},{"why":"Provides the static structure factors and excitation branches for spin-1 condensates that support the dispersion and instability calculation.","marker":"[17]"},{"why":"Established the experimental platform and the observation of spin turbulence in the EAP-to-EPP quench, which this work extends to the reverse transition.","marker":"[5]"},{"why":"Reported the wall-vortex composite defects in the weak-quench EPP-to-EAP transition, which this paper extends to higher quench strengths.","marker":"[7]"},{"why":"Provides the sodium scattering length $a_s=1.88a_0$ used to compute $c_2$ and hence the quench strength $\\tilde{q}$.","marker":"[21]"},{"why":"Confirms the effective density $\\bar{n}=2n_0/3$ via spin-sound speed measurement, anchoring the $\\tilde{q}$ calibration.","marker":"[23]"}],"fun_headline_variants":["Spinor condensate quench: instability vanishes at q=2","Weak to strong quench crossover in spinor BEC","Momentum shift kills long-wavelength spin growth","Zero-momentum magnon growth stops at quench strength 2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim rests on the calibration of the quench strength $\\tilde{q}$ using a specific scattering length and an averaged density; if that calibration is off, the crossover positions shift and the quantitative match to Bogoliubov theory breaks down.","fun_headline_variants_meta":{"raw":{"variants":["Spinor condensate quench: instability vanishes at q=2","Weak to strong quench crossover in spinor BEC","Momentum shift kills long-wavelength spin growth","Zero-momentum magnon growth stops at quench strength 2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000258,"raw_usage":{"total_tokens":1598,"prompt_tokens":976,"completion_tokens":622,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":592,"completion_tokens_details":{"reasoning_tokens":551}},"tokens_in":592,"tokens_out":622,"duration_ms":6104,"temperature":1.0,"reasoning_tokens":551,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:39:26.401704+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An experiment that measured the momentum-resolved growth of density modulations in the $m_z=0$ component after the quench (for example, by Fourier analysis of absorption images) and found significant low-wavenumber growth for $\\tilde{q}>2$, or a seeded sample whose early dynamics still accelerates at those quenches, would falsify the claim that the zero-momentum mode has stabilized.","supporting_citations":[{"cited_title":"Kawaguchi and M","cited_arxiv_id":null,"evidence_quote":"Supplies the Bogoliubov spectrum of the EPP state, including the two magnon branches used in the instability analysis."},{"cited_title":"Lamacraft, Quantum Quenches in a Spinor Conden- sate, Phys","cited_arxiv_id":null,"evidence_quote":"Predicted that a finite-wavenumber mode becomes most unstable for quenches with $\\tilde{q}>1$, motivating the search for the crossover."},{"cited_title":"Matuszewski, Rotonlike Instability and Pattern For- mation in Spinor Bose-Einstein Condensates, Phys","cited_arxiv_id":null,"evidence_quote":"Predicted rotonlike instability and quasi-periodic pattern formation from finite-wavenumber spin excitations, used to interpret the array-like patterns."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the static structure factors and excitation branches for spin-1 condensates that support the dispersion and instability calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Established the experimental platform and the observation of spin turbulence in the EAP-to-EPP quench, which this work extends to the reverse transition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reported the wall-vortex composite defects in the weak-quench EPP-to-EAP transition, which this paper extends to higher quench strengths."},{"cited_title":"Knoop, T","cited_arxiv_id":null,"evidence_quote":"Provides the sodium scattering length $a_s=1.88a_0$ used to compute $c_2$ and hence the quench strength $\\tilde{q}$."},{"cited_title":"Observation of Two Sound Modes in a Binary Superfluid Gas","cited_arxiv_id":"1907.10289","evidence_quote":"Confirms the effective density $\\bar{n}=2n_0/3$ via spin-sound speed measurement, anchoring the $\\tilde{q}$ calibration."}],"review_version":1}