{"id":"15118d96-4523-45e5-8efe-02f65a558cf8","arxiv_id":"2506.21670","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Simulations show a 2D Bose-Einstein condensate formed by a thermal quench fills with quantum vortices whose energy spectrum obeys Kolmogorov k^{-5/3} scaling and collapses across quench rates under Kibble-Zurek rescaling.","lead":"A numerical study shows that forming a Bose-Einstein condensate by quenching temperature creates a turbulent tangle of quantum vortices, not just a smooth superfluid. The turbulence follows both Kibble-Zurek scaling with quench speed and Kolmogorov's energy spectrum, so a single process connects defect formation to fluid turbulence.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main-text spectral fit (k^-1.65) is incompatible with the supplement's own ESS exponent, which implies k^-1.82 for the same inertial range.","rationale":"The Kibble–Zurek part of the paper is reasonably supported: vortex number and equilibration-time scalings are close to the mean-field prediction, and the collapse in Fig. 3 is a nontrivial empirical demonstration. The most serious weakness is internal to the turbulence diagnostic. The main text concludes k^-5/3 from direct fits over a narrow inertial range, but the supplement's ESS analysis yields zeta(2)=0.820±0.022, which via the standard relation E(k) ~ k^-(1+zeta(2)) implies k^-1.82 for the same data. The error bars make this a ~4σ discrepancy, and the K62 intermittency model invoked in the supplement makes the spectrum steeper, not shallower. The BKT/critical-exponent concern raised in the earlier review is real but secondary: the simulations effectively determine the KZ exponents, and the data support mean-field values; however, the theoretical justification for mean-field exponents in a 2D system still needs explicit discussion. The appropriate response is conditional acceptance: the authors must reconcile the spectral and structure-function analyses, or revise the Kolmogorov claim and the associated conclusions about spontaneous quantum turbulence.","tokens_in":19982,"tokens_out":14343,"duration_ms":165665,"concrete_test":"From the same ensemble-averaged fields used for Fig. 2(b) and Fig. S3, compute S_2(r) and E_i(k) over exactly the same inertial range (d_v ≤ r ≤ l_v, equivalently 2π/l_v ≤ k ≤ 2π/d_v) and form the compensated spectra E_i(k) k^5/3 and E_i(k) k^1.82. If E_i(k) k^1.82 is flat while E_i(k) k^5/3 decreases/increases systematically, the Kolmogorov claim fails. Additionally, fit S_3(r) directly: if its exponent differs from 1, the reported zeta(p) values are relative exponents and must be rescaled before comparing to the spectrum; the comparison must then be repeated with the absolute zeta(2).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that the incompressible energy spectrum shows Kolmogorov k^-5/3 scaling is contradicted by the paper's own extended-self-similarity analysis. For the tau_Q=350 dataset, Fig. 2(b) gives a direct spectral fit E_i(k) ~ k^-1.654±0.031, while Fig. S3(b) reports the second-order longitudinal structure function exponent zeta(2)=0.820±0.022. The supplement itself states that zeta(2)=2/3 corresponds to E_i(k) ~ k^-5/3, i.e. E(k) ~ k^-(1+zeta(2)). Using the measured zeta(2), the same inertial range should give E_i(k) ~ k^-1.820±0.022. The discrepancy is ~0.17, about 4.4 combined standard deviations. The K62 intermittency correction invoked in the supplement makes the spectrum steeper (alpha ≈ 1.78–1.82), not shallower, so it cannot reconcile the direct fit with the ESS result. Thus the two diagnostics used to establish quantum turbulence are internally inconsistent for the same data; the claim of agreement with 5/3 is not supported by the evidence as presented.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports stochastic projected Gross-Pitaevskii simulations of a homogeneous 2D Bose gas driven across condensation by a linear chemical-potential quench. It finds that the vortex number at equilibration scales as tau_Q^-0.534±0.047, the equilibration time as tau_Q^0.541±0.069, and the incompressible kinetic-energy spectrum exhibits a power-law region with fitted exponent alpha≈1.65–1.71 that the authors associate with Kolmogorov k^-5/3 scaling. It further proposes a scaling ansatz leading to collapse of spectra for different quench times under the rescalings E_i tau_Q^{3/4} vs. k tau_Q^{1/4}, and it reports a universal collapse in Fig. 3. The authors interpret these results as spontaneous quantum turbulence generated by the Kibble-Zurek mechanism.","tokens_in":20278,"tokens_out":11091,"duration_ms":124986,"significance":"The paper addresses a question of genuine current interest: whether a finite-time quench can produce quantum turbulence without external stirring. It has clear strengths: the KZ predictions are tested directly against the simulations rather than used as fitting assumptions; the vortex-number and equilibration-time power laws are checked over more than an order of magnitude in tau_Q; the spectral collapse in Fig. 3 is a nontrivial prediction; and the data set is openly available. The single-vortex derivation leading to E_i ∝ tau_Q^{-1} is simple and consistent with the numerical fits in Fig. S4. However, the central claim that the spectrum is Kolmogorov k^-5/3 is weakened by an internal inconsistency with the supplement's own extended-self-similarity analysis, as detailed below; the title and abstract claim of 'spontaneous quantum turbulence' should be evaluated after that inconsistency is resolved.","major_comments":[{"comment":"The direct spectral fit and the extended-self-similarity (ESS) fit are mutually inconsistent for the same tau_Q=350 dataset. The main text reports E_i(k) ∝ k^{-1.654±0.031} in the inertial range, while the supplement reports zeta(2)=0.820±0.022 for S_2 ∝ [S_3]^{zeta(2)} over the corresponding range. Using the paper's own relation (S_3 ∝ r in the inertial range, and E_i(k) ∝ k^{-(1+zeta(2))}), the ESS value implies E_i(k) ∝ k^{-1.820±0.022}; the difference is about 4.4 combined standard deviations. The K62 intermittency correction would make the spectrum steeper, not shallower, so it cannot reconcile the two diagnostics. Please reconcile the fit ranges and definitions, or revise the claim that the spectrum agrees with 5/3.","section":"Fig. 2(b) and Fig. S3(b)"},{"comment":"The KZ exponents and the collapse variables in Eq. (4) are derived using mean-field critical exponents nu=1/2, z=2, but the manuscript does not justify why these apply to the homogeneous 2D Bose gas, whose equilibrium transition is of Berezinskii-Kosterlitz-Thouless type with exponential, not power-law, scaling. The numerical exponents for n_v and t_eq are consistent with nu=1/2, z=2 within error, but they are not by themselves a proof of the mean-field choice over the simulated range. Please provide a justification for the effective mean-field behavior (for example, the role of the c-field cutoff or finite-size effects), or present the collapse using fitted effective exponents and show the sensitivity of the collapse to their uncertainties.","section":"Main text, 'Kibble-Zurek dynamics of the BEC transition', and Eq. (4)"},{"comment":"The derivation of the amplitude A(tau_Q) in Eq. (4) uses the total incompressible energy E_i ∝ tau_Q^{-(1+2nu)/(1+z nu)} together with the single-length-scale ansatz E_i(k; tau_Q)=A(tau_Q) F(k xi_hat). However, the spectrum also contains a k^{-3} vortex-core region whose amplitude is independent of tau_Q, and the authors themselves state that no collapse holds there. The total-energy integration therefore includes a tau_Q-independent contribution. The argument should either restrict the ansatz to the inertial window and show that the core contribution is negligible for the total-energy scaling, or the derivation should be stated as approximate for the inertial range only.","section":"Main text, 'Kibble-Zurek universality of SQT'"}],"minor_comments":[{"comment":"The fit intervals used to obtain the exponents alpha and beta in Fig. 2 are not stated explicitly; they are only indicated by vertical lines. Please state the k ranges used for each fit.","section":"Fig. 2 caption and surrounding text"},{"comment":"The reference to the supplementary material appears as 'url will be inserted by publisher'; this placeholder should be completed before publication.","section":"Reference [49]"},{"comment":"The equilibration-time fit is a two-term expression, t_eq = (3.779±0.693) tau_Q^{0.541±0.069} + (0.078±0.059) tau_Q, but the main text quotes only the power-law part. The linear term is non-negligible at the largest simulated tau_Q (about 78 at tau_Q=1000), so its potential effect on the quoted KZ exponent should be discussed.","section":"Eq. (S6) in the Supplementary Material"}],"recommendation":"major_revision","confidential_remarks":"The KZ scaling part of the manuscript is largely sound, and the spectral-collapse test is a genuine strength; the main roadblock is the internal inconsistency between the direct spectral exponent and the ESS exponent, which affects the central claim of Kolmogorov 5/3 scaling. I do not see this as an irreparable error, but it requires careful re-analysis before the paper can be accepted. A number of KZ references (refs 28, 36, 37, 68) are from the same group; this is not disqualifying in a focused subfield, but the authors should ensure that the prior-work discussion is balanced."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely new: a finite-time thermal quench of a 2D Bose gas produces a newborn condensate whose incompressible energy spectrum is close to k^{-5/3}, and spectra for different quench times collapse under the Kibble-Zurek rescaling k->k tau_Q^{1/4}, E_i->E_i tau_Q^{3/4}. No stirring, no external potential sweep—just the quench rate controls the vortex density and the inertial range. That's a clean protocol and a nice bridge between two universality concepts. The numerics look careful: 1000 realizations per point, vortex detection with Newton refinement, and the data are openly available. The fits for n_v and t_eq match KZ scaling with the expected exponents, and the energy scaling E_i ~ tau_Q^{-1} is verified at two temperatures.\n\nNow the soft spots, in proportion. The biggest is the mean-field exponents. The paper assumes nu=1/2, z=2 for a homogeneous 2D Bose gas, whose equilibrium transition is BKT. KZ scaling for BKT is not power-law in the standard way; the authors don't justify why mean-field exponents should apply to their quench rates. This matters directly: the collapse variable in Eq. (4) is built from nu/(1+z nu), so if the effective exponents are different, the collapse is a fit to a chosen form, not a test of universality.\n\nSecond, there's an internal inconsistency in the tau_Q=350 data that a referee will catch. Fig. 2(b) gives a direct fit E_i ~ k^{-1.654±0.031}. The supplement's ESS analysis reports zeta(2)=0.820±0.022, and the supplement itself states that E_i(k) ~ k^{-(1+zeta(2))}. That gives k^{-1.820±0.022}, which is ~4.4 sigma away from the direct fit. The K62 intermittency correction makes the spectrum steeper, not shallower, so it doesn't reconcile them. The ESS fit may cover a broader range than the IR used for the spectrum, but the paper doesn't say that or explain why the exponents differ. As written, the two diagnostics conflict, and the claim of 5/3 scaling is weaker than it appears.\n\nMinor: the inertial range is narrow (less than a decade), which is typical for BECs but limits the force of the spectral fit. The derivation of the E_i scaling uses a single-vortex integration with a logarithmic correction; that's fine.\n\nOverall, this deserves a serious referee and likely publication after revision, but the authors need to address the BKT/mean-field issue and the spectral inconsistency head-on. I'd recommend sending it to review with those questions explicit.","headline":"A genuinely new protocol for quench-generated quantum turbulence with a beauty KZ collapse across quench times, but the 5/3 claim is weakened by an unresolved internal inconsistency between the spectral fit and the ESS exponent, and by unexamined mean-field exponents in a 2D BKT system.","tokens_in":20794,"tokens_out":3775,"would_cite":false,"duration_ms":38991,"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":"A finite-time thermal quench alone can generate quantum turbulence in a newborn Bose-Einstein condensate, with vortex energy spectra obeying both Kolmogorov and Kibble-Zurek scaling.","keywords":["quantum turbulence","Kibble-Zurek mechanism","Bose-Einstein condensation","quantum vortices","Kolmogorov scaling","stochastic projected Gross-Pitaevskii equation","two-dimensional superfluids","universal scaling"],"falsifier":"Repeat the SPGPE quench simulations (or perform the equivalent experiment in a uniform quasi-2D BEC) over a wider range of $\\tau_Q$, extending below 10 and above 1000, and test two things: whether the inertial-range exponent stays at $5/3$ and whether the spectra still collapse under $k \\tau_Q^{1/4}$ versus $E_i \\tau_Q^{3/4}$. If the collapse breaks down at slow quenches or the fitted exponents drift with $\\tau_Q$, the combination of Kolmogorov and Kibble-Zurek scaling for SQT is not universal. A complementary test replaces the power-law correlation length with the BKT exponential scaling and asks whether the collapse variable becomes logarithmic in $\\tau_Q$.","tokens_in":2376,"feed_emoji":"🌀","tokens_out":3163,"duration_ms":142856,"temperature":0.7,"pith_summary":"The paper argues that there is a simpler route to quantum turbulence than stirring a superfluid: cool a two-dimensional Bose gas through its condensation transition at a finite rate. The quench forces the order parameter to break symmetry at the Kibble-Zurek freeze-out time, leaving a newborn condensate seeded with a random tangle of vortices and antivortices whose density is set entirely by the quench duration. Using stochastic Gross-Pitaevskii simulations, the authors find that the incompressible kinetic energy of this vortex tangle shows the Kolmogorov $k^{-5/3}$ spectrum, the hallmark of turbulence, and that spectra from different quench times collapse onto one universal curve under the Kibble-Zurek rescaling. If true, the quench rate serves as a single control knob that determines both the vortex density and the turbulent state, and two distinct universality principles, defect scaling at phase transitions and energy-spectrum scaling in turbulence, meet in one system.","feed_headline":"Quench a Bose gas and quantum turbulence appears on its own","feed_subtitle":"No stirring involved: quench speed alone sets the vortex tangle, whose energy spectrum obeys Kolmogorov 5/3 scaling.","key_machinery":"The argument is carried by three objects working together: the stochastic projected Gross-Pitaevskii equation (SPGPE), a classical-field equation of motion for the low-energy coherent modes of a Bose gas in contact with a thermal reservoir, which supplies the ensemble of quench trajectories; the canonical Helmholtz-like decomposition of the density-weighted velocity field $\\mathbf{u} = \\sqrt{\\rho}\\,\\mathbf{v}$ into divergence-free (incompressible) and curl-free (compressible) parts, which isolates the vortex contribution to the kinetic energy; and the Kibble-Zurek correlation length $\\hat{\\xi} \\propto \\tau_Q^{\\nu/(1+z\\nu)}$, which sets both the vortex density $n_v \\propto \\hat{\\xi}^{-2}$ and the characteristic momentum scale at equilibration. The load-bearing identity is the scaling ansatz $E_i(k;\\tau_Q) = A(\\tau_Q) F(k\\hat{\\xi})$ with $A(\\tau_Q) \\propto \\tau_Q^{-(1+\\nu)/(1+z\\nu)}$, which converts into the concrete collapse rule of Eq. (4), $k \\to k \\tau_Q^{\\nu/(1+z\\nu)}$ and $E_i \\to E_i \\tau_Q^{(1+\\nu)/(1+z\\nu)}$, evaluated at $\\nu=1/2$, $z=2$ so that the exponents become $1/4$ and $3/4$.","core_discovery":"On the paper's own terms, the central discovery is that spontaneous quantum turbulence (SQT) is an intrinsic outcome of finite-time Bose-Einstein condensation: a linear quench of the chemical potential through the critical point produces a condensate populated with vortices of winding number $\\pm 1$, whose number scales as $n_v \\propto \\tau_Q^{-1/2}$ and whose density-weighted superfluid velocity splits into compressible and incompressible parts with the incompressible kinetic energy spectrum obeying $E_i(k) \\propto k^{-5/3}$ over an inertial range bounded below by the vortex-core scale and above by the mean intervortex distance. The spectral shape is universal across quench rates: under the rescaling $k \\to k \\tau_Q^{1/4}$ and $E_i \\to E_i \\tau_Q^{3/4}$ (the Kibble-Zurek exponents obtained from the mean-field values $\\nu=1/2$, $z=2$), all spectra collapse onto a single curve, and the total incompressible energy follows $E_i \\propto \\tau_Q^{-1}$. The same collapse works for the compressible spectrum at low temperature, where phonons emitted by vortex-antivortex annihilation dominate, and the velocity structure functions show extended self-similarity with intermittency corrections described by the refined K62 model.","pith_inferences":["Editorial inference: the load-bearing choice of mean-field exponents $\\nu=1/2$, $z=2$ deserves scrutiny because the homogeneous 2D Bose gas has a Berezinskii-Kosterlitz-Thouless equilibrium transition with exponential rather than power-law critical scaling; re-running the collapse test with the BKT scaling form over slower and faster quenches would reveal whether the reported exponents and rescali","Editorial inference: the same mechanism should generalize to other symmetry-breaking transitions that create vortex-like defects, such as spinor condensates or higher-dimensional superfluids, with the rescaling exponents set by the relevant $\\nu$ and $z$; the paper notes a 4D study in passing but does not develop this analogy.","Editorial inference: because the inertial range spans barely a decade, the extended self-similarity analysis of velocity structure functions is what actually pins the $5/3$ law; a sharper experimental test would measure $S_p(r)$ directly in a quenched gas rather than relying on the energy spectrum alone."],"forward_implications":["No external stirring is required to create a turbulent superfluid: the quench duration $\\tau_Q$ alone sets the vortex density and, through it, the energy-injection scale and the width of the Kolmogorov inertial range.","The universal collapse of Eq. (4) gives experimenters a concrete diagnostic: measure $E_i(k)$ at equilibration for several quench times and check whether $E_i \\tau_Q^{3/4}$ plotted against $k \\tau_Q^{1/4}$ falls on one curve.","The turbulent state is transient by construction: vortex-antivortex annihilation, atom losses, and coarsening suppress the Kolmogorov scaling after equilibration, so SQT is best probed in a window around the equilibration time.","Slower quenches broaden the inertial range (larger intervortex separation at fixed core size), while quenches faster than a threshold ($\\tau_Q \\lesssim 10$) leave the Kibble-Zurek scaling regime and lose the collapse.","The total incompressible kinetic energy at equilibration obeys $E_i \\propto \\tau_Q^{-1}$ for the mean-field exponents, a scaling that follows from the same argument and is verified in the data."],"supporting_citations":[{"why":"Supplies the critical scaling $\\xi \\propto |\\varepsilon|^{-\\nu}$, $\\tau \\propto |\\varepsilon|^{-z\\nu}$ and the Kibble-Zurek freeze-out time that the paper uses to build the rescaling of Eq. (4).","marker":"[28]"},{"why":"Defines the stochastic projected Gross-Pitaevskii equation, the simulation method that produces the quenched condensates and vortex ensembles.","marker":"[43–48]"},{"why":"The experimental observation of spontaneous Kibble-Zurek vortices in Bose-Einstein condensation that anchors the vortex-formation premise.","marker":"[31]"},{"why":"The 2D BEC experiment reporting Kolmogorov scaling of the incompressible spectrum, the empirical fingerprint the paper reproduces in a quench-generated state.","marker":"[18]"},{"why":"Supplies the density-weighted velocity field and the compressible/incompressible decomposition of kinetic energy used to define $E_i(k)$.","marker":"[51–54]"},{"why":"Establishes the $k^{-3}$ vortex-core spectrum and the intervortex-distance and core-size bounds that frame the inertial range.","marker":"[55]"},{"why":"Provides the criterion (change of condensate growth from exponential to linear) used to define the equilibration time at which spectra are measured.","marker":"[34]"},{"why":"Supplies the refined K62 intermittency model used to fit the structure-function exponents measured via extended self-similarity.","marker":"[66]"},{"why":"Provides the universal-breakdown argument for Kibble-Zurek scaling in fast quenches that the paper invokes to explain deviations for $\\tau_Q \\leq 10$.","marker":"[36]"}],"fun_headline_variants":["Quench a Bose gas, and vortices self-assemble into turbulence","Kibble-Zurek quench spontaneously seeds quantum vortex turbulence","No stirring: a single quench drives quantum turbulence formation","Universal Kolmogorov scaling from quench-born vortex tangles"],"cache_read_input_tokens":22912,"weakest_assumption_plain":"The predicted exponents and the spectral collapse assume mean-field critical exponents $\\nu = 1/2$ and $z = 2$ for the 2D Bose-gas transition, even though the homogeneous 2D gas has a Berezinskii-Kosterlitz-Thouless transition with exponential rather than power-law scaling; if the effective exponents differ over the simulated quench times, the Kibble-Zurek exponents and the collapse variables change.","fun_headline_variants_meta":{"raw":{"variants":["Quench a Bose gas, and vortices self-assemble into turbulence","Kibble-Zurek quench spontaneously seeds quantum vortex turbulence","No stirring: a single quench drives quantum turbulence formation","Universal Kolmogorov scaling from quench-born vortex tangles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000305,"raw_usage":{"total_tokens":1749,"prompt_tokens":940,"completion_tokens":809,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":745}},"tokens_in":556,"tokens_out":809,"duration_ms":8850,"temperature":1.0,"reasoning_tokens":745,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:22:26.968069+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the SPGPE quench simulations (or perform the equivalent experiment in a uniform quasi-2D BEC) over a wider range of $\\tau_Q$, extending below 10 and above 1000, and test two things: whether the inertial-range exponent stays at $5/3$ and whether the spectra still collapse under $k \\tau_Q^{1/4}$ versus $E_i \\tau_Q^{3/4}$. If the collapse breaks down at slow quenches or the fitted exponents drift with $\\tau_Q$, the combination of Kolmogorov and Kibble-Zurek scaling for SQT is not universal. A complementary test replaces the power-law correlation length with the BKT exponential scaling and asks whether the collapse variable becomes logarithmic in $\\tau_Q$.","supporting_citations":[{"cited_title":"Universality of phase transition dynamics: Topological defects from symmetry breaking,","cited_arxiv_id":null,"evidence_quote":"Supplies the critical scaling $\\xi \\propto |\\varepsilon|^{-\\nu}$, $\\tau \\propto |\\varepsilon|^{-z\\nu}$ and the Kibble-Zurek freeze-out time that the paper uses to build the rescaling of Eq. (4)."},{"cited_title":"Spontaneous vortices in the formation of Bose-Einstein condensates,","cited_arxiv_id":null,"evidence_quote":"The experimental observation of spontaneous Kibble-Zurek vortices in Bose-Einstein condensation that anchors the vortex-formation premise."},{"cited_title":"Kolmogorov scaling in turbulent 2d bose-einstein condensates,","cited_arxiv_id":null,"evidence_quote":"The 2D BEC experiment reporting Kolmogorov scaling of the incompressible spectrum, the empirical fingerprint the paper reproduces in a quench-generated state."},{"cited_title":"Energy spectra of vortex distributions in two-dimensional quantum turbu- lence,","cited_arxiv_id":null,"evidence_quote":"Establishes the $k^{-3}$ vortex-core spectrum and the intervortex-distance and core-size bounds that frame the inertial range."},{"cited_title":"Defect formation beyond kibble-zurek mechanism and hologra- phy,","cited_arxiv_id":null,"evidence_quote":"Provides the criterion (change of condensate growth from exponential to linear) used to define the equilibration time at which spectra are measured."},{"cited_title":"A refinement of previous hypotheses concerning the local structure of turbulence in a viscous incompressible fluid at high reynolds number,","cited_arxiv_id":null,"evidence_quote":"Supplies the refined K62 intermittency model used to fit the structure-function exponents measured via extended self-similarity."},{"cited_title":"Universal breakdown of kibble-zurek scaling in fast quenches across a phase transition,","cited_arxiv_id":null,"evidence_quote":"Provides the universal-breakdown argument for Kibble-Zurek scaling in fast quenches that the paper invokes to explain deviations for $\\tau_Q \\leq 10$."}],"review_version":1}