{"id":"d2e2c45d-0c65-49ee-852d-f9cf074a101c","arxiv_id":"2412.13986","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A double sine-Gordon effective model reproduces subdiffusive coarsening of a spin-1 Bose gas, with multi-well field spread distinguishing subdiffusive from diffusive scaling.","lead":"This paper shows that a simple model, the double sine-Gordon equation, captures how magnetic order coarsens in a spin-1 Bose gas after a quench, with slow subdiffusive growth when the field spreads across many potential wells. The result offers a microscopic route to classifying far-from-equilibrium coarsening in many-body systems.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The DSG reduction is not proven in the multi-well regime where subdiffusive scaling occurs; the neglected density-phase couplings and the ad hoc regularization at potential maxima are the load-bearing assumptions.","rationale":"The reader's weakest assumption is exactly the validity of the DSG mapping, and my stress-test sharpens that concern: the expansion in Eq. (B6) drops terms that may not be small when φ_s crosses the potential maxima, and the regularization in App. B.3 is explicitly non-rigorous. This is the single most load-bearing point because the paper's central claim is that the spin-1 Bose gas belongs to the DSG universality class; if the effective theory fails in the multi-well regime where subdiffusive scaling is observed, the entire identification collapses. The DSG numerics are internally consistent and convincingly show a qualitative difference between few-well and many-well dynamics, but the quantitative link to the spin-1 parameters is weakened by the hand-tuned couplings, and the experimental evidence comprises only a PDF of φ_s, not the scaling exponents. All of this supports the reader's conditional verdict: the central qualitative claim is plausible but not yet firmly established. The proposed test would either validate the neglected-term expansion in the relevant regime or demonstrate that a more careful effective theory is needed. Since no new information here moves the verdict away from conditional acceptance, I recommend no change to the reader's assessment.","tokens_in":23907,"tokens_out":4509,"duration_ms":42970,"concrete_test":"Using the truncated-Wigner trajectories of the full spin-1 gas that produced Fig. 4, compute the spacetime integrals of the neglected terms |δρ ∇δρ ∇δρ| and |δρ ∇φ_s ∇φ_s| from Eq. (B6) relative to the retained quadratic terms L2 (Eq. B10), restricting the evaluation to the space-time regions where φ_s is within one spin healing length of a potential maximum (2ν+1)π. If the ratio of the neglected terms to L2 is not small (≪1) in those regions, the Gaussian integration underlying Eq. (10) is uncontrolled precisely where the multi-well subdiffusive dynamics is claimed to reside.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation of the DSG effective theory (Eq. 10, App. B) relies on two controlled approximations: a Gaussian integration over density fluctuations after dropping terms of O(δρ∇δρ∇δρ) and O(δρ∇φ∇φ) (Eq. B6), and an expansion of the resulting nonlocal Lagrangian around φ_s = 2πZ, with the divergences at φ_s = (2ν+1)π regularized by replacing ∇^2 with -k_ξs^2 (App. B.3). The central claim concerns the subdiffusive scaling regime in which the field spreads over many minima; there φ_s necessarily passes through the potential maxima where the k→0 Green's function diverges and where the expansion around 2πZ is invalid. The paper itself states that the restoration of DSG around the maxima is 'not intended to constitute a rigorous derivation' (App. B.3). If the neglected terms in Eq. (B6) are not small on the actual spin-1 trajectories in the scaling regime, or if the k~k_ξs regularization does not correctly capture the dynamics of the kinks that carry the field across maxima, then the observed agreement between DSG and spin-1 scaling could be coincidental, and the claimed identification of the spin-1 gas with the DSG universality class is not established. The hand-tuned DSG couplings (λ = 10λ_s for the subdiffusive run, versus λ_spin-1 ≈ 5.8λ_s,spin-1 derived from Eq. 10, App. C2b) further mean that the numerical demonstration of β = 0.28 was not performed at the microscopically derived parameters, leaving the parameter independence of the exponent within DSG unresolved.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the universal subdiffusive coarsening dynamics observed in the one-dimensional easy-plane spin-1 Bose gas after a quench belongs to the double sine-Gordon (DSG) universality class. It derives a low-energy effective DSG Lagrangian for the spinor phase by integrating out density fluctuations (Sec. III.A, App. B), and numerically shows that the DSG model exhibits self-similar subdiffusive scaling with β≈0.28 in 1+1D when the field spreads over many potential wells, and diffusion-type scaling β≈0.5 when only two wells are occupied (Sec. II). These results are compared with truncated-Wigner simulations and experimental PDFs of the spinor phase (Sec. III.B), and with a 2+1D DSG simulation that gives β≈0.51 (Sec. III.C).","tokens_in":24353,"tokens_out":24894,"duration_ms":202577,"significance":"If correct, the work would establish a concrete microscopic route from the spin-1 Bose gas parameters to a DSG effective description, providing a step toward a classification of coarsening dynamics and offering spinor gases as a platform for studying sine-Gordon physics. The paper is commendable for its careful scaling analysis (residuals, reference-time dependence) and for presenting both numerical and experimental evidence. However, the central identification rests on an approximate derivation that the authors themselves state is not rigorous (App. B.3), on DSG simulations whose couplings are tuned rather than taken from the microscopic mapping (App. C.2.b), and on a fit that imposes α=dβ. These issues currently weaken the claim that the spin-1 gas belongs to the DSG universality class.","major_comments":[{"comment":"The subdiffusive DSG simulations are performed with λ=10λ_s, whereas the microscopic mapping (Eq. (10)) yields λ_{spin-1}≈5.8λ_s. The text states that the couplings 'were chosen such as to achieve reliable self-similar scaling' (App. C.2.b). This leaves open whether the DSG model at the microscopically derived couplings exhibits β≈0.28; consequently, the numerical agreement between the DSG and spin-1 exponents could reflect the parameter choice rather than universality. The authors should demonstrate that β is independent of λ/λ_s over a range including the microscopic value, or repeat the simulation at λ≈5.8λ_s.","section":"App. C.2.b / Eq. (10)"},{"comment":"The mapping to the DSG model is not controlled in the regime where the subdiffusive scaling occurs. The derivation expands around φ_s=2πZ and neglects terms of O(δρ∇δρ∇δρ, δρ∇φ∇φ) (Eq. (B6)); the k→0 Green's function diverges at φ_s=(2ν+1)π, and the regularization ∇²→-k_ξs² is stated to be 'not intended to constitute a rigorous derivation' (App. B.3). The multi-well field configurations that are central to the claimed subdiffusive class necessarily cross these maxima. Without a quantitative check of the neglected terms on the actual spin-1 trajectories, the identification of the spin-1 gas with the DSG class rests on an uncontrolled approximation.","section":"App. B.3 / Sec. III.A"},{"comment":"The scaling collapse in Fig. 1 imposes α=dβ, justified by 'conservation of the momentum integral over S.' However, the DSG equation (1) has no visible Noether symmetry that would conserve ∫S(k,t) dk; φ is a noncompact field with a potential that breaks shift symmetry. The authors should either derive this conservation law for the DSG model or fit α and β independently (as done in the 2D case in Fig. 5) to confirm that β≈0.28 is not an artifact of the imposed relation.","section":"Sec. II.A / Fig. 1"},{"comment":"The sign convention relating the effective theory to the DSG simulations should be clarified. For c_1<0, the sin²φ_s term in Eq. (10) has a negative coefficient (and the same sign appears in Eq. (B22)), whereas the DSG Lagrangian (C2) is simulated with λ_s>0. These appear to correspond to opposite signs of the sin(2φ) nonlinearity in the equations of motion. Please state explicitly how the coefficients of Eq. (10) map to λ and λ_s in Eq. (C2), and confirm that the simulated DSG model is the one obtained from the microscopic derivation.","section":"Eq. (10) / Eq. (C2)"}],"minor_comments":[{"comment":"The experimental comparison uses a Boltzmann approximation to extract V_eff(φ_s) from a far-from-equilibrium distribution; the agreement is qualitative, and both the mean-field shift and calibration offset are adjusted. The authors may wish to state this more cautiously in the main text.","section":"Sec. III.B / Fig. 4"},{"comment":"The paragraph distinguishing the observed scaling from phenomenological diffusion equations is somewhat verbose; a concise statement of the relation to Refs. [71,80] would help the reader.","section":"Sec. II.C"},{"comment":"The sentence that one may neglect terms of order ˙φ_j sin²(φ_s/2) and (∇φ_j)² sin²(φ_s/2) is a strong truncation; since these terms are quadratic in derivatives and sin², it would be useful to state why they are negligible compared to the retained terms in the scaling regime.","section":"App. B.3"},{"comment":"The definition of k_ξs in the caption, k_ξs=(2Mρ̃|c_1|)^{1/2}≈4Q, mixes dimensional and numerical statements; the units and the numerical value should be defined more clearly.","section":"Fig. 1 caption"},{"comment":"Reference [80] is a master's thesis; if a published version exists or becomes available, it would be preferable to cite that instead.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper reports a substantial experimental and numerical effort, and the scaling analysis is careful. However, the central identification of the spin-1 gas with the DSG universality class rests on a derivation that the authors admit is not rigorous and on DSG simulations with hand-tuned couplings. The apparent sign mismatch between the derived sin²φ_s term and the simulated λ_s is a serious technical point that should be carefully checked. The α=dβ assumption in the 1D fit also needs justification. These issues are addressable with additional simulations and analysis, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The genuinely new result is the derivation of a noncompact double sine-Gordon theory for the spinor phase of the easy-plane spin-1 gas, together with the claim that subdiffusive coarsening is tied to field configurations spread over many potential wells, while two-well configurations give β ≈ 1/2. The DSG numerics are careful: the scaling analysis includes reference-time stability and residual checks, and the 2D DSG results line up with the earlier spin-1 results from Schmied et al. The experimental PDF of φs matches the DSG potential shape once you allow the stated calibration offset. Credit where due: the paper is also honest about its own limits, explicitly saying the restoration of DSG around the potential maxima is not a rigorous derivation.\n\nThe soft spots are real but not fatal. First, the effective-theory derivation is controlled only near φs = 2πZ. The subdiffusive regime is exactly where the field passes through the maxima, and there the divergence is regularized by replacing ∇² with −k_ξs². That is an assumption about kink-scale physics, not a derivation. If that regularization does not capture the actual dynamics of the kinks that carry the field across maxima, the agreement with the full spin-1 simulations could be coincidental. Second, the DSG couplings in the subdiffusive run were tuned by hand: λ = 10λ_s, whereas the spin-1 mapping gives λ ≈ 5.8λ_s. So the β = 0.28 demonstration was not performed at the microscopically derived couplings, and parameter independence within DSG is not shown. Third, the abstract overstates the experimental support: the experiment provides the PDF of the spinor phase, not the scaling exponents.\n\nDoes the central claim hold up? Within the DSG model, the two-regime picture looks solid: multi-well occupation gives subdiffusive, two-well gives diffusive. What is not yet established is that the spin-1 gas actually belongs to the same universality class. The mapping is plausible and the numerical agreement is suggestive, but the maxima regularization and the coupling choice are the load-bearing assumptions that need strengthening.\n\nThis paper is for the quantum-gas and nonequilibrium-scaling community. It deserves a serious referee. I would send it to peer review with a request for two things: a coupling scan showing the subdiffusive exponent is robust in the relevant range, and a direct comparison of the effective theory to full spin-1 dynamics at the derived couplings. That would move the claim from plausible to convincing.","headline":"A useful and honest paper that maps the easy-plane spin-1 gas to a double sine-Gordon model and ties subdiffusive scaling to multi-well occupation, but the load-bearing reduction is approximate and the couplings are hand-tuned, so the universality-class claim is plausible rather than proven.","tokens_in":24852,"tokens_out":2944,"would_cite":true,"duration_ms":25022,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.75.Mn","05.70.Ln"],"model":"deepseek-v4-flash","headline":"The easy-plane spin-1 Bose gas coarsens according to the double sine-Gordon universality class, with the spinor phase as the single relevant field.","keywords":["double sine-Gordon model","spin-1 Bose gas","coarsening dynamics","universality class","nonthermal fixed point","spinor phase","subdiffusive scaling","self-similar scaling"],"falsifier":"Measure the spinor-phase structure factor in a long one-dimensional $^{87}$Rb condensate quenched from the polar into the easy-plane phase and fit the self-similar collapse: if the extracted $\\beta$ equals $0.28(3)$ while the phase PDF covers several $2\\pi$ wells, the DSG claim is supported, whereas a value near $0.5$ in the same multi-well regime would falsify it.","tokens_in":23700,"feed_emoji":"🌀","tokens_out":9952,"duration_ms":82312,"temperature":0.7,"pith_summary":"This paper aims to establish that the coarsening dynamics of a spin-1 Bose gas after a quench into the easy-plane phase belongs to the universality class of the double sine-Gordon model: a single real scalar field, the spinor phase, carries the universal scaling. By integrating out small density fluctuations, the authors derive a low-energy effective action for $\\varphi_s$ that has both $\\cos\\varphi_s$ and $\\sin^2\\varphi_s$ terms, and they show numerically that this model reproduces the subdiffusive coarsening of the full spinor gas, with $\\beta=0.28(3)$ in one dimension. The same model yields diffusion-type scaling, $\\beta\\simeq 0.5$, when the field occupies only two minima of its periodic potential, so the spread over many minima is identified as the precondition for the slow exponent. If the reduction is correct, it explains why the spinor phase alone can describe the scaling, and it connects a microscopic many-body Hamiltonian to a universality class of coarsening.","feed_headline":"Spin-1 gas coarsening falls in the double sine-Gordon class","feed_subtitle":"A spinor-phase field theory predicts subdiffusive scaling beta=0.28 in 1D, confirmed by ultracold-atom data.","key_machinery":"The central object is the double sine-Gordon model of a real scalar field, $\\ddot{\\varphi}=c_s^2\\Delta\\varphi-\\lambda\\sin\\varphi+\\lambda_s\\sin(2\\varphi)$, whose periodic potential combines a $\\cos\\varphi$ term and a $\\sin^2\\varphi$ term. The argument is carried by the Gaussian integration of density fluctuations in the spin-1 Lagrangian, which eliminates the $\\delta\\rho$ and $\\delta\\epsilon$ fields and leaves an effective action for the two phase angles; after expanding about $\\varphi_s=2\\pi\\mathbb{Z}$, the spinor phase obeys the DSG equation while the Larmor phase becomes a free massless field. The $\\sin^2\\varphi_s$ term is load-bearing: truncating it to a pure sine-Gordon model leaves the one-dimensional power spectra static. The other indispensable ingredient is the regularization of the effective potential near $\\varphi_s=(2\\nu+1)\\pi$ by momenta of order the spin healing momentum $k_{\\xi_s}$, which allows localized kinks to interpolate between adjacent minima and thus lets the field spread over many wells.","core_discovery":"The paper's central claim is that the low-energy physics of the easy-plane spin-1 Bose gas reduces to a double sine-Gordon theory for the spinor phase $\\varphi_s$, with effective Lagrangian $\\mathcal{L}_{\\mathrm{eff}}=-\\frac{1}{32c_1}\\dot{\\varphi}_s^2-\\frac{n(\\tilde{\\rho}-2n)}{4M\\tilde{\\rho}}(\\nabla\\varphi_s)^2-\\left[2c_1n(\\tilde{\\rho}-2n)-\\frac{q^2}{16c_1}\\right]\\cos\\varphi_s+\\frac{q^2}{32c_1}\\sin^2\\varphi_s$. In one spatial dimension this theory shows self-similar subdiffusive coarsening with $\\alpha=\\beta=0.28(3)$ when the unwrapped field wanders over many minima of the periodic potential, and diffusion-type scaling with $\\beta=0.52(4)$ when only two minima are populated; the contrasting exponents are therefore linked to how many wells the field configuration visits, not to domain size alone. The paper further claims that the noncompact DSG field inherits the topological content of the compact spinor phase by unwrapping, so kink and vortex information is encoded without explicit topological degrees of freedom. Numerical simulations of the full spin-1 model and experimental probability distributions of $\\varphi_s$ from a quasi-one-dimensional condensate are presented as evidence that the spinor phase indeed localizes at multiples of $2\\pi$ and spreads across several wells, as the DSG picture requires.","pith_inferences":["If the mapping is generic, one could expect a similar spinor-phase DSG reduction for other easy-plane ferromagnetic spin-$F$ condensates, with $\\beta\\simeq0.28$ in one dimension whenever the field visits many wells.","A clean experimental extension would be to tune the quench depth or quadratic Zeeman shift to control how many minima the spinor phase explores and watch $\\beta$ cross from about $0.5$ to $0.28$ within the same system.","The paper's focus on spread over many minima suggests a general criterion for coarsening universality classes: the relevant quantity may be the number of potential wells visited by the order-parameter field, not merely the rate of domain growth.","One direct test of the unwrapping picture would be to compare kink statistics in the full spin-1 gas with the kink statistics of the noncompact DSG field; if they match, the topological encoding claim is confirmed."],"forward_implications":["The spinor phase alone, without the Larmor or total phase, is enough to reproduce the subdiffusive coarsening of the full spin-1 gas in one dimension.","The number of occupied minima of the periodic potential sets the universality class: many wells give $\\beta=0.28(3)$, two wells give $\\beta\\simeq0.5$.","In two dimensions the DSG model reproduces diffusion-type scaling with $\\beta=0.51(8)$ and $\\alpha=0.98(20)$, consistent with earlier spin-1 results, so the same effective theory covers both dimensionalities.","The $\\sin^2\\varphi_s$ term is essential for the scaling: a pure sine-Gordon truncation leaves the power spectra static in one dimension.","The measured PDF of $\\varphi_s$ matching the DSG potential makes the spinor condensate a platform for studying sine-Gordon dynamics, including soliton collisions and breathers."],"supporting_citations":[{"why":"Provides the Luttinger-liquid-type effective theory and scaling formalism for phase correlations that the DSG reduction extends to topological excitations.","marker":"[9]"},{"why":"Reports the subdiffusive coarsening in the spin-1 gas structure factor that the DSG model is claimed to reproduce.","marker":"[60]"},{"why":"Found diffusion-type scaling in a two-dimensional spin-1 gas, which the two-dimensional DSG simulations corroborate with $\\beta\\approx0.5$.","marker":"[62]"},{"why":"Gives the analytic scaling analysis of sine-Gordon-type models predicting both the subdiffusive and diffusion-type exponents used here.","marker":"[71]"},{"why":"Observed space-time vortex defects and the same subdiffusive scaling in the full spin-1 model, providing the numerical target for the effective theory.","marker":"[72]"},{"why":"Supports the two-exponent scaling picture for sine-Gordon-type models through an alternative approximation to the effective action.","marker":"[80]"},{"why":"Supplies the mean-field phase structure (polar and easy-plane) of the spin-1 gas on which the mapping is built.","marker":"[81]"}],"fun_headline_variants":["Double sine-Gordon class unifies spin-1 gas coarsening dynamics","Spin-1 BEC coarsening: subdiffusive when field spans many wells","Coarsening in spinor gas: exponent beta=0.28 from double sine-Gordon","Two coarsening regimes in spin-1 gas traced to sine-Gordon wells","Multi-well phase configurations govern spin-1 gas scaling laws"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The mapping assumes density fluctuations are weak enough to be integrated out to second order and that the spinor phase can be expanded around $\\varphi_s=2\\pi\\mathbb{Z}$, with the singular behaviour at potential maxima cured by momenta near the spin healing length; the paper states the last step is not intended to constitute a rigorous derivation.","fun_headline_variants_meta":{"raw":{"variants":["Double sine-Gordon class unifies spin-1 gas coarsening dynamics","Spin-1 BEC coarsening: subdiffusive when field spans many wells","Coarsening in spinor gas: exponent beta=0.28 from double sine-Gordon","Two coarsening regimes in spin-1 gas traced to sine-Gordon wells","Multi-well phase configurations govern spin-1 gas scaling laws"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000251,"raw_usage":{"total_tokens":1620,"prompt_tokens":1070,"completion_tokens":550,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":686,"completion_tokens_details":{"reasoning_tokens":446}},"tokens_in":686,"tokens_out":550,"duration_ms":5735,"temperature":1.0,"reasoning_tokens":446,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:35:09.317986+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the spinor-phase structure factor in a long one-dimensional $^{87}$Rb condensate quenched from the polar into the easy-plane phase and fit the self-similar collapse: if the extracted $\\beta$ equals $0.28(3)$ while the phase PDF covers several $2\\pi$ wells, the DSG claim is supported, whereas a value near $0.5$ in the same multi-well regime would falsify it.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the subdiffusive coarsening in the spin-1 gas structure factor that the DSG model is claimed to reproduce."},{"cited_title":"Spitz, J","cited_arxiv_id":null,"evidence_quote":"Gives the analytic scaling analysis of sine-Gordon-type models predicting both the subdiffusive and diffusion-type exponents used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Observed space-time vortex defects and the same subdiffusive scaling in the full spin-1 model, providing the numerical target for the effective theory."},{"cited_title":"Siovitz, S","cited_arxiv_id":null,"evidence_quote":"Supports the two-exponent scaling picture for sine-Gordon-type models through an alternative approximation to the effective action."},{"cited_title":"Gliott, A","cited_arxiv_id":null,"evidence_quote":"Supplies the mean-field phase structure (polar and easy-plane) of the spin-1 gas on which the mapping is built."}],"review_version":1}