{"id":"cefd9caa-bc5e-44ff-ade7-24f22e041796","arxiv_id":"2607.26620","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Cooling a 3D Bose gas exactly to its BEC critical momentum yields a distinct non-thermal fixed point (α=3/2, β=3/4) with superdiffusive spreading, making the equilibrium threshold a dynamical critical point.","lead":"This paper argues that cooling a 3D Bose gas across the Bose–Einstein condensation threshold organizes its far-from-equilibrium dynamics into three universal regimes, including a new 'critical' fixed point exactly at threshold. If correct, the equilibrium condensation point also acts as a dynamical phase transition, giving a unified map for quench experiments in quantum gases.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Effective-coupling closure (Eq. 2) is unbenchmarked for the threshold quench; if it fails, the critical fixed point (3/2,3/4,2) is an artifact.","rationale":"The reader's weakest assumption correctly identifies the uncontrolled effective-coupling closure as the load-bearing premise. All three dynamical regimes—thermal, coarsening, and especially the new critical fixed point—are computed from Eq. (1)+(2). The thermal regime is low-risk because it is checked against known equilibrium properties (Eq. 3, Fig. 4a). The coarsening regime is partially supported by independent physics (vortex-line tension argument, p^-4 tail). But the critical fixed point has no independent support: its novel exponent β=3/4 comes only from the same closure, the same numerical data, and an ad hoc front criterion. A classical-field GPE simulation is the natural, decisive test: it makes no use of the closure and directly tests whether the threshold quench produces superdiffusive (β=3/4) rather than diffusive (β=1/2) front propagation. Such a test would distinguish a genuine non-perturbative fixed point from a truncation artifact. Because the paper is explicitly conditional on this closure and no independent benchmark is supplied, the appropriate verdict remains CONDITIONAL; my stress-test does not change that verdict.","tokens_in":16974,"tokens_out":7302,"duration_ms":81269,"concrete_test":"Run classical-field Gross–Pitaevskii simulations in 3D for the same quench protocol: a uniform dilute Bose gas with 4π^2ρξ^3=10^3, initial distribution n_p(0)∝exp(−p^4/p_c^4) with p0=pc≈7.7/ξ, and evolve to times ≳10^4/(gρ). Extract n_0(t) and the cascade front p_f(t) using Eq. (B2), varying λ from 0.3 to 0.8. If the GPE yields an asymptotic front exponent β consistent with 3/4 and clearly distinct from the subcritical β=1/2, the closure and the fixed point are supported; if the GP data instead show β≈1/2 or a λ-dependent crossover, the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—a distinct critical fixed point at p0=pc with exponents (α,β,γ)=(3/2,3/4,2)—rests entirely on the quantum-kinetic Eq. (1) with the non-perturbative effective coupling Eq. (2) evaluated using the instantaneous distribution n_k(t). This closure is imported from Ref. [7], not derived here, and is the only regularization of the Boltzmann finite-time singularity. Its accuracy is never benchmarked against ab initio Gross–Pitaevskii simulations or experiments specifically for threshold quenches. In the regime beyond t*, where the paper claims new physics, occupations are large and g_eff deviates strongly from g (Fig. 7), so uncontrolled errors in Eq. (2) could produce an apparent new attractor that is not present in the underlying Bose gas dynamics. The distinction between coarsening (β=1/2) and critical (β=3/4) is delicate: both regimes share α=3/2, so the entire claim hinges on a single exponent extracted through the ad hoc front definition (B2) with λ=0.6, without error bars or sensitivity analysis. The statement that β=3/4 is 'independently confirmed' by Appendix B is not an independent check, because that appendix uses the same front-fitting procedure. No code or data are provided, so these numerics cannot be re-run.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies far-from-equilibrium dynamics of a 3D Bose gas after a cooling quench, using a quantum kinetic equation with a non-perturbative effective coupling geff (Eq. (2)). It claims that the equilibrium BEC threshold pc acts as a dynamical critical point: quenches above pc relax to a thermal fixed point; quenches below pc first show weak turbulence and then coarsening with exponents (α,β,γ)=(3/2,1/2,4); quenches exactly to pc are controlled by a new critical non-thermal fixed point with exponents (α,β,γ)=(3/2,3/4,2), characterized by superdiffusive propagation of critical fluctuations. The evidence is numerical: scaling collapses of the momentum distribution, zero-mode growth, and the inverse-cascade front p_f(t).","tokens_in":17345,"tokens_out":3954,"duration_ms":42015,"significance":"If correct, the paper would establish a unified dynamical phase diagram for cooling quenches across the BEC transition and identify a previously unknown critical fixed point. This is a conceptually appealing and potentially influential result for quantum-gas experiments and non-equilibrium statistical mechanics. The paper also proposes a physically motivated resolution of the long-standing Boltzmann finite-time singularity. However, the central claim rests on a single numerical framework with an imported closure (Eq. (2)) and on exponent extractions that lack error quantification; the result is therefore not yet established at the level required for a definitive claim.","major_comments":[{"comment":"The central distinction between the critical fixed point (β=3/4) and the coarsening fixed point (β=1/2) rests on the front exponent β, extracted via the arbitrary threshold λ=0.6 in Eq. (B2). No error bars, sensitivity analysis to λ, or alternative front definitions are given. The phrase 'independently confirmed by an analysis of the cascade front in Appendix B' is misleading: Appendix B uses the same definition, same data, and same fitting procedure. Since α=3/2 is shared by both regimes, all the weight is on this single exponent.","section":"Critical fixed point; Eq. (B2)"},{"comment":"The entire study is a numerical solution of Eq. (1) with the non-perturbative closure geff=g/|1+gΠ|. This closure is imported from Ref. [7] and is not benchmarked for the threshold quench. For p0=pc, Fig. 7 shows geff deviates strongly from g at long times, so the new critical dynamics lives precisely in the regime where the closure is most consequential. A failure of Eq. (2) would make the critical fixed point an artifact. The authors should validate against Gross–Pitaevskii simulations (at least for a range of p0 including pc) or against existing experiments; without such a benchmark the claim remains conditional.","section":"Eq. (2); Model"},{"comment":"The selection of z=2/3 for the critical fixed point is circular: it uses the observed α/β=2, which itself comes from the fitted α and β of the same numerical data. Equation (5) is then used to infer z. This is a consistency check, not a derivation. To make the claim load-bearing, the scaling of geff should be measured directly from the numerical solution (e.g., from the kernel in Eq. (A2)) and compared with the prediction z=2/3, or an independent analytical argument should be supplied.","section":"Appendix C; exponent relation Eq. (5)"},{"comment":"The collapse in Fig. 5c is presented as 'almost perfect', but no quantitative comparison is made with the alternative coarsening collapse (α=3/2, β=1/2). Given that both fixed points share α, the reader cannot judge whether the data actually discriminate between β=3/4 and β=1/2. A quantitative goodness-of-fit or residual analysis, and preferably a plot showing that the (3/2,1/2) collapse fails, is required.","section":"Fig. 5c; 'Critical fixed point'"},{"comment":"No code or data are provided, and the numerical solver details, while extensive, are not sufficient to reproduce the results without significant effort. For a paper whose conclusions are purely numerical, the absence of public code/data is a serious deficiency. At minimum, the authors should provide the extracted exponent values with uncertainty estimates and make the code available.","section":"Appendices A–D; reproducibility"}],"minor_comments":[{"comment":"The caption writes n0(t→∞)∼t^{3/2} for p0<pc; this should be n0(t)∼t^{3/2} for long times, not t→∞ (the asymptotic state is not actually reached in finite time). Similar wording appears around Eq. (4) and in the coarsening section.","section":"Fig. 1 caption"},{"comment":"The text says the coarsening dynamics is 'no longer self-similar' because α/β=3≠γ=4. But Eq. (4) with f(x)~x^{-4} still represents a self-similar collapse (n_p=t^{3/2} f(p t^{1/2})). The terminology needs clarification: if the authors mean that the power-law exponent differs from the self-similarity ratio, this should be stated explicitly and reconciled with the use of Eq. (4) for the collapse.","section":"Fig. 4c, right panel"},{"comment":"The fit giving ν=1 is quoted without an error bar or a plot showing the numerical ζ vs. (p0-pc). The inset in Fig. 4a is mentioned, but a quantitative comparison or at least the number of points and residuals should be given.","section":"Eq. (3)"},{"comment":"The critical row lists γ=2=α/β; this is fine, but for coarsening the table lists γ=4 while α/β=3, which again highlights the 'self-similar' ambiguity noted above. The table would benefit from a column defining whether the scaling function is a pure power law.","section":"Table I"},{"comment":"Some references appear incorrectly formatted or incomplete (e.g., Ref. [19] uses arXiv:2605.23600, Ref. [41] lacks a journal volume/page, Ref. [44] has a comma in the author list). Please check all references for completeness.","section":"General"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper before reading it. First, it is a numerical study of a quantum kinetic equation with a self-consistent effective coupling (recycled from Chantesana et al.), and it claims a third universal attractor — a critical fixed point — for quenches landing exactly on the BEC threshold. Second, that claim is not yet supported by anything beyond the numerics that define it. The exponents are fitted from the same simulation data that are used to establish the collapse, and the effective coupling is never checked against an independent calculation. So treat the critical NTFP as an interesting conjecture, not a result.\n\nThe paper does have real merits. It gives a clean phase diagram for cooling quenches and shows convincingly that below threshold the system moves through a weak-turbulence stage and then into coarsening, consistent with earlier work. The thermal fixed point above threshold reproduces the expected equilibrium critical behavior, which is a good consistency check. The scaling ansatz and the exponent relation (2+z)β−α=1/2 are set out clearly, and the flow diagram in Fig. 5d is a nice way to visualize the regimes. The authors are also honest that Eq. (2) is imported from Ref. [7] and not derived here.\n\nThe soft spots are concentrated on the new branch. The entire case for the critical fixed point rests on the closure geff=g/|1+gΠ| evaluated with the instantaneous distribution n_k(t). In the long-time, high-occupation regime where the paper claims new physics, geff deviates strongly from g, so the approximation is exactly where it matters. There is no comparison to Gross–Pitaevskii simulations or experiment for the threshold quench, and the stress-test note is right: if the closure fails there, the critical point is an artifact. The distinction from coarsening is also delicate: both regimes have α=3/2, so the whole claim hangs on β=3/4 vs 1/2. That β is extracted with the front definition Eq. (B2) using λ=0.6, with no error bars and no sensitivity analysis. And the 'independent confirmation' in Appendix B is not independent; it is the same front-fitting procedure. Finally, no code or data are shipped, so the numerics cannot be re-run. None of these are fatal — they are fixable — but they are real.\n\nThis paper is for people working on NTFPs and quench dynamics in Bose gases. It gives them a map of regimes and a specific, testable conjecture. I would not cite it as evidence, but I would bring it to a reading group, because the closure problem is a good starting point for discussion. A serious editor should send this to peer review: the question is important, the work is technically careful, and a referee can demand the missing benchmarks and data/code. My own verdict is conditional, leaning skeptical on the new branch.","headline":"Plausible new branch in a known kinetic framework, but the critical fixed point hinges on an unbenchmarked closure and self-fitted exponents.","tokens_in":17884,"tokens_out":2720,"would_cite":false,"duration_ms":27520,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.30.Jp","03.75.Kk","05.70.Ln"],"model":"deepseek-v4-flash","headline":"The equilibrium condensation threshold also acts as a dynamical critical point, separating three universal relaxation regimes in a quenched Bose gas.","keywords":["Bose-Einstein condensation","non-thermal fixed point","dynamical phase transition","quantum kinetic theory","quench dynamics","coarsening","critical dynamics","superdiffusion"],"falsifier":"Compute the long-time dynamics of a uniform 3D Bose gas quenched precisely to the condensation threshold using the Gross–Pitaevskii equation (or measure it in an experiment) and extract the inverse-cascade front scale p_f(t). If p_f(t) does not scale as t^(−3/4) at long times while n_0(t) ~ t^(3/2), or if the rescaled distributions do not collapse onto a single curve with f(x) ~ x^(−2), the critical fixed point is not realized.","tokens_in":1572,"feed_emoji":"⚛️","tokens_out":5509,"duration_ms":68220,"temperature":0.7,"pith_summary":"The paper claims that the Bose-Einstein condensation temperature is not only a thermodynamic boundary but also the organizing center of nonequilibrium relaxation in a 3D Bose gas after a cooling quench. Using a non-perturbative quantum kinetic equation, the authors show that quenches above the threshold relax to a thermal fixed point, quenches below pass through weak turbulence and then a vortex-coarsening fixed point, and quenches directly to the threshold are controlled by a distinct critical fixed point where fluctuations spread superdiffusively. The central result is a new set of dynamical exponents (α, β, γ) = (3/2, 3/4, 2) for the critical quench, which differs from the coarsening exponents (3/2, 1/2, 4). If correct, this establishes a far-from-equilibrium counterpart of the condensation phase transition, with the equilibrium critical point shaping long-time universal dynamics.","feed_headline":"Quenching a Bose gas to its critical point creates a new universal regime","feed_subtitle":"At the condensation threshold, critical fluctuations spread superdiffusively, a distinct dynamical class between coarsening and thermal.","key_machinery":"The central object is the non-perturbative effective interaction g_eff(ε,p,t) = g / |1 + g Π(ε,p,t)|, a ladder-resummed coupling that replaces the bare contact interaction in the quantum Boltzmann equation. Evaluated with the instantaneous momentum distribution n_k(t), it regularizes the finite-time singularity of the perturbative Boltzmann equation and produces algebraic momentum dependencies that select the different fixed points. At the critical fixed point, g_eff ~ p^(2/3), which the authors show is the only exponent compatible with an inverse cascade and self-similar scaling, and it indicates that the dynamics is dominated by on-shell scattering processes.","core_discovery":"The central discovery is that the post-quench energy density selects the asymptotic nonequilibrium attractor, and at the critical energy the system approaches a previously unidentified non-thermal fixed point. At this critical fixed point, the inverse-cascade front grows as p_f(t) ~ t^(-3/4), faster than the diffusive t^(-1/2) of coarsening, and the momentum distribution obeys n_p(t) = t^(3/2) f(p t^(3/4)) with f(x) ~ x^(-2) at intermediate momenta. The exponents are selected uniquely by combining the general scaling relation (2+z)β − α = 1/2 with the self-similarity constraint α/β = 2 and the effective coupling exponent z = 2/3. This establishes a distinct universality class for critical qu","pith_inferences":["If the critical fixed point is genuine, critical quenches should exhibit aging and universal two-time correlation functions with a new aging exponent, analogous to classical critical quenches; two-time measurements in cold gases could test this prediction.","The dynamical phase diagram suggests a sharp observable signature: a small change in quench depth across p_c changes the long-time front exponent from 1/2 to 3/4, which can be resolved in time-of-flight imaging of the momentum distribution.","The z = 2/3 effective coupling at criticality might be detectable through momentum-resolved Bragg scattering or by measuring collision rates at low momenta, since it implies a nontrivial suppression of the interaction relative to the bare coupling.","The same non-perturbative kinetic framework could be adapted to other transitions with a threshold energy, such as the 2D BKT transition, where the predicted exponents would differ and provide a sharp test of the framework's generality."],"forward_implications":["Quenches above the condensation threshold relax to the equilibrium thermal fixed point, with the correlation length diverging as (p_0 − p_c)^(−1), reproducing conventional critical behavior.","Quenches below threshold exhibit two sequential regimes: a transient weak-turbulence state with exponents (α,β,γ) ≈ (−2.7, −1.1, 2.45), followed by a coarsening fixed point with exponents (3/2, 1/2, 4) governed by vortex-line recombination.","Quenches precisely at the threshold are controlled by a distinct critical fixed point with exponents (3/2, 3/4, 2), characterized by superdiffusive spreading of critical fluctuations and scale invariance that already emerges on microscopic timescales.","The long-time effective coupling scales as p^z with z = 2 below the threshold and z = 2/3 at the threshold, uniquely selecting the fixed points through the exponent relation.","For deep subcritical quenches the weak-turbulence regime disappears entirely, and the system evolves directly toward the coarsening fixed point, while also developing a direct cascade at large momenta."],"fun_headline_variants":["Critical quench reveals a new non-thermal fixed point","Bose gas critical quench: superdiffusive universality class","New dynamical class at Bose condensation threshold","Superdiffusive spreading at critical quench fixed point","Critical quench uncovers a distinct non-equilibrium attractor"],"cache_read_input_tokens":19072,"weakest_assumption_plain":"The entire classification rests on the self-consistent effective coupling g/(1 + gΠ) evaluated with the instantaneous momentum distribution being accurate beyond the Boltzmann blow-up time and especially at the threshold quench, since no ab initio simulation or experimental comparison is provided for that case.","fun_headline_variants_meta":{"raw":{"variants":["Critical quench reveals a new non-thermal fixed point","Bose gas critical quench: superdiffusive universality class","New dynamical class at Bose condensation threshold","Superdiffusive spreading at critical quench fixed point","Critical quench uncovers a distinct non-equilibrium attractor"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000446,"raw_usage":{"total_tokens":2092,"prompt_tokens":748,"completion_tokens":1344,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":492,"completion_tokens_details":{"reasoning_tokens":1265}},"tokens_in":492,"tokens_out":1344,"duration_ms":8990,"temperature":1.0,"reasoning_tokens":1265,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T12:10:30.233526+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the long-time dynamics of a uniform 3D Bose gas quenched precisely to the condensation threshold using the Gross–Pitaevskii equation (or measure it in an experiment) and extract the inverse-cascade front scale p_f(t). If p_f(t) does not scale as t^(−3/4) at long times while n_0(t) ~ t^(3/2), or if the rescaled distributions do not collapse onto a single curve with f(x) ~ x^(−2), the critical fixed point is not realized.","supporting_citations":[],"review_version":1}