{"id":"8d29e31b-5697-4b36-b271-922c179ed234","arxiv_id":"2607.02743","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"EDQNM simulations of nonhelical MHD find that at Re or Rm ≳ 10^6 the kinematic growth rate, saturated magnetic-to-kinetic energy ratio (~0.55), and integral-scale ratio (~3) become independent of Pm, with both spectra approaching k^{-3/2} via Alfvénisation.","lead":"A turbulence closure model shows that small-scale dynamos reach universal saturated states at extreme Reynolds numbers, independent of magnetic Prandtl number. This gives concrete numbers for magnetic energy and coherence scales in galaxies and clusters that direct simulations cannot yet reach.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Asymptotic ratios rest on a fixed Cs=0.26 and isotropic EDQNM damping that is never varied; the claimed universality may be an artifact of that closure choice.","rationale":"The Reader correctly isolates the weakest link: the fixed Cs = 0.26 and the isotropy assumption that underwrite the nonlinear damping. The analytic equivalence to Kazantsev is clean and the public code makes the numerics reproducible, so an outright REJECT is unwarranted. At the same time the headline numbers (0.55 and 3) have never been shown to be stable under variation of the only free parameter that controls the nonlinear transfers. A short Cs-sensitivity suite would either confirm the asymptotes or reveal that they move, thereby deciding whether the claim can be upgraded or must remain conditional. That is precisely the load-bearing concern; everything else (intermittency, reconnection, DNS comparison) is secondary until this internal robustness check is performed.","tokens_in":22083,"tokens_out":593,"duration_ms":5981,"concrete_test":"Re-run the three highest-Re/Rm cases that define the asymptotic plateaus (Pm = 10^{-3}, 1, 10^3) with Cs varied by ±30 % (Cs = 0.18 and 0.34) while keeping all other parameters fixed; if either χ or k_M/k_V shifts by more than ~10 % the claimed universal numbers are closure-dependent rather than robust.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that, once Re ≳ 10^6 (Pm>1) or Rm ≳ 10^6 (Pm<1), nonhelical EDQNM produces Pm-independent asymptotic values χ ≃ 0.55 and k_M/k_V ≃ 3 together with a shared -3/2 inertial range. Those numbers are obtained exclusively with the damping rate of Eq. (20) that contains the single free coefficient Cs = 0.26 (taken from Grappin et al. 1982) and with the Markovian triad time of Eq. (19). Section VI acknowledges that EDQNM assumes strict spectral isotropy and cannot capture intermittency or reconnection, yet no sensitivity study of Cs, of the Alfvénic contribution inside µ_k, or of the Markovian approximation is reported. Because the nonlinear transfer terms T^M_MM and T^V_VM that set the saturated state are proportional to θ_kpq, any systematic change in the damping immediately rescales the effective nonlinear diffusion and therefore the location of the spectral peak and the energy ratio. The analytic recovery of Kazantsev (θ=const) is solid, but the quantitative asymptotes that constitute the paper’s strongest claim are therefore still conditional on an untested closure parameter.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper shows that the nonhelical EDQNM closure for incompressible MHD recovers the classical Kazantsev/Kulsrud–Anderson kinematic small-scale dynamo when the triad relaxation time is taken constant (θ_kpq = τ), both analytically (Appendix A) and numerically. It then uses the full scale-dependent EDQNM equations to explore SSD growth and saturation over a wide (Re, Rm, Pm) range inaccessible to DNS. The central claims are that, for Re ≳ 10^6 (Pm > 1) or Rm ≳ 10^6 (Pm < 1), the kinematic growth rate becomes approximately Pm-independent, the saturated magnetic-to-kinetic energy ratio converges to χ ≃ 0.55, the integral-wavenumber ratio k_M/k_V asymptotes to ≃ 3, and both spectra develop a shared −3/2 inertial range attributed to Alfvénisation.","tokens_in":22439,"tokens_out":1219,"duration_ms":17300,"significance":"If the reported asymptotic regimes hold, they supply concrete, Pm-independent targets (χ ≃ 0.55, k_M/k_V ≃ 3, shared −3/2 spectra) for modeling SSD-dominated systems such as elliptical galaxies and the ICM, and for interpreting future high-Re global simulations. Strengths include a careful analytic recovery of the Kazantsev equation (Appendix A), term-by-term transfer analysis of saturation, public release of code and data, and systematic coverage of both Pm ≫ 1 and Pm ≪ 1 at extreme Reynolds numbers. The work usefully unifies earlier kinematic models with a nonlinear spectral closure and identifies an effective Pm = 1 saturated state driven by Alfvénisation.","major_comments":[{"comment":"The quantitative asymptotes χ ≃ 0.55 and k_M/k_V ≃ 3 (Figs. 16–17, §V.C) and the Pm-independent growth-rate plateau (Fig. 13, §V.B) are obtained exclusively with the damping rate of Eq. (20) and the fixed coefficient C_s = 0.26 taken from Grappin et al. (1982). Because the nonlinear transfers T^M_MM and T^V_VM that set saturation are proportional to θ_kpq, a change in C_s (or in the Alfvénic contribution inside µ_k) rescales the effective nonlinear diffusion and can shift the spectral peak and energy ratio. No sensitivity scan of C_s, of the Alfvénic term, or of the Markovian approximation is reported. A modest variation of C_s (or an explicit statement that the numbers are closure-dependent) is needed before the claimed universality can be treated as robust.","section":null},{"comment":"The thresholds Re_asym ∼ 10^6 and Rm_asym ∼ 10^6 are stated as sharp asymptotic onsets (§V.B–C and abstract), yet the manuscript provides neither resolution-variation checks at those extremes beyond the Fres = 16 vs 32 comparison in Appendix A.3 nor uncertainty estimates on the fitted plateaus. Given that the logarithmic grid and the choice of k_max = 4 max(k_η, k_ν) control the dissipative range, the precise numerical values of the thresholds and of χ, k_M/k_V should be accompanied by a short resolution/parameter-uncertainty discussion so that readers can judge how firmly the asymptotes are established.","section":null},{"comment":"Section VI correctly notes that EDQNM assumes spectral isotropy and cannot capture intermittency or reconnection, yet the abstract and conclusions present the asymptotic ratios as guidance for astrophysical systems without quantifying how strongly those assumptions may bias the numbers. A clearer separation between (i) the robust qualitative existence of high-Re/Rm plateaus and (ii) the precise numerical values that are closure- and isotropy-dependent would strengthen the paper’s claims and avoid over-interpretation by the community.","section":null}],"minor_comments":[{"comment":"Fig. 12: the scaling R_cr_m ∝ Re^{−5/4} in the viscous regime is stated but not overlaid on the plot; adding a guide line would help the reader.","section":null},{"comment":"Notation: both M(k) and M_k are used for the magnetic spectrum; a single convention would improve readability.","section":null},{"comment":"Eq. (26) defines k_M/k_V with the inverse-moment convention of Monin & Yaglom; a brief reminder that this is not the peak-wavenumber ratio would avoid confusion with DNS literature that often quotes k_peak.","section":null},{"comment":"Typographical: “Alfv´enisation” / “Alfv`enisation” appear with inconsistent accents; standardize to Alfvénisation throughout.","section":null},{"comment":"Table I lists only two Kazantsev runs; a short note on how many full-EDQNM runs underlie Figs. 13–17 would aid reproducibility.","section":null}],"recommendation":"major_revision","confidential_remarks":"The analytic Kazantsev recovery and the public code are genuine strengths. The main risk is that the headline numbers (0.55, 3, 10^6) will be quoted as universal without the closure caveats. Requiring a short C_s sensitivity test (or explicit framing as EDQNM-specific) is proportionate and fixable; I would not reject on that basis. Scope is appropriate for a theory/astro-plasma journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The paper’s real payload is two things: a clean analytic and numerical demonstration that EDQNM with constant triad time recovers the classic Kazantsev kinematic equation, and a set of previously inaccessible asymptotic plateaus (Pm-independent growth rate, χ≃0.55, kM/kV≃3, shared –3/2 spectra) once Re or Rm exceeds ~10^6. Those numbers are new because DNS cannot reach them, and the authors ship the code, so the numerics are checkable.\n\nWhat they do well is the bookkeeping. Appendix A walks through the transfer integrals carefully; the kinematic growth-rate match and the self-similar Kazantsev spectrum are solid. The nonlinear term-by-term analysis (T^M_MM as hyperdiffusion at large k and enhanced diffusion at small k, Alfvénisation flattening the kinetic spectrum) is transparent and matches earlier analytic pictures. They also map Rcr_m(Pm) and show the expected viscous-to-inertial transition. For people who need concrete saturation ratios for non-rotating plasmas or for sub-grid models, this is usable guidance.\n\nThe soft spot is exactly the one the stress-test flags, and it is real but not fatal. All the quantitative asymptotes sit on the fixed Cs=0.26 and the isotropic Markovian damping of Eqs. 19–20; no sensitivity scan is reported. Because the nonlinear transfers that set the saturated state are proportional to θ_kpq, a different eddy-damping choice can move χ and kM/kV. The authors themselves note the isotropy and intermittency limitations in §VI, so they are not overselling. The thresholds ~10^6 are also resolution-dependent in the usual logarithmic-grid sense. None of this undoes the kinematic equivalence or the qualitative picture of Alfvénisation-driven saturation; it just means the precise numbers should be treated as EDQNM numbers, not universal constants.\n\nThis is for dynamo theorists and people building global MHD models of clusters or ellipticals who need high-Re guidance. The math and citation pattern look careful; the free parameters are few and inherited. I would send it to referees. Worth reading and citing for the asymptotes and the public code, with the usual closure caveat attached.","headline":"Clean Kazantsev recovery plus new high-Re/Rm asymptotic plateaus inside EDQNM; the numbers are useful but still conditional on the fixed Cs=0.26 damping that was never varied.","tokens_in":23047,"tokens_out":558,"would_cite":true,"duration_ms":6753,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"At extreme Reynolds numbers the small-scale dynamo saturates with a magnetic-to-kinetic energy ratio near 0.55 and an integral-scale ratio near 3, independent of magnetic Prandtl number.","keywords":["small-scale dynamo","EDQNM closure","magnetic Prandtl number","nonlinear saturation","Alfvénisation","MHD turbulence","Reynolds number asymptotics"],"falsifier":"A high-resolution direct numerical simulation of nonhelical MHD turbulence at Re or Rm of order 10^6 that measures a saturated magnetic-to-kinetic energy ratio substantially different from 0.55, or an integral-scale ratio substantially different from 3, would falsify the claimed asymptotes.","tokens_in":23014,"feed_emoji":"🌌","tokens_out":971,"duration_ms":7224,"temperature":0.7,"pith_summary":"Small-scale dynamos are the main candidate for generating magnetic fields in turbulent plasmas that lack strong rotation, such as galaxy clusters and elliptical galaxies. Their early, kinematic growth is well mapped, but the later nonlinear saturation is hard to reach with direct simulations at realistic Reynolds numbers. This paper shows that a classical spectral closure for magnetohydrodynamic turbulence recovers the classic kinematic theory and then carries it into the nonlinear regime across a wide range of magnetic Prandtl numbers. In highly turbulent regimes the growth rate, the saturated energy ratio and the ratio of magnetic to kinetic integral scales all become independent of Prandtl number, while both spectra settle into the same -3/2 inertial range through Alfvénic coupling. The result supplies concrete asymptotic targets for future simulations and for modelling magnetised astrophysical plasmas that live far beyond present computational reach.","feed_headline":"Small-scale dynamos saturate at 55% energy, scale ratio 3","feed_subtitle":"Extreme turbulence erases Prandtl-number dependence, giving universal targets for galaxy-cluster fields","key_machinery":"The eddy-damped quasi-normal Markovian (EDQNM) closure for incompressible MHD turbulence: its triad relaxation time reduces exactly to the classic Kazantsev kinematic dynamo equation when magnetic back-reaction is neglected and the relaxation time is held constant, and the same equations with scale-dependent damping then evolve the nonlinear saturation.","core_discovery":"When fluid Reynolds number exceeds roughly 10^6 (for Pm > 1) or magnetic Reynolds number exceeds roughly 10^6 (for Pm < 1), nonhelical EDQNM simulations of the small-scale dynamo yield a Prandtl-number-independent kinematic growth rate, a saturated magnetic-to-kinetic energy ratio that converges to approximately 0.55, and a magnetic-to-kinetic integral-wavenumber ratio that asymptotes to approximately 3; both spectra share a -3/2 inertial range produced by Alfvénisation.","pith_inferences":["The same asymptotic ratios may serve as fixed points for any spectral model that includes Alfvénic damping, not only EDQNM.","If the -3/2 range is observed in Faraday-rotation or synchrotron spectra of clusters, it would be indirect evidence that the plasma has already entered the Prandtl-independent regime.","The critical magnetic Reynolds number’s mild rise at low Pm suggests that planetary and stellar interiors may still host efficient small-scale dynamos once Rm exceeds a few tens.","Extensions that restore weak helicity could test whether the same energy and scale ratios survive when a large-scale dynamo is also present."],"forward_implications":["Astrophysical models of galaxy-cluster or elliptical-galaxy magnetisation can adopt a universal saturation efficiency of order 55 percent once the plasma is sufficiently turbulent.","The saturated magnetic integral scale is predicted to lie only a factor of three below the kinetic integral scale, independent of microscopic Prandtl number.","Saturated kinetic and magnetic spectra should both display a -3/2 inertial range once Alfvénisation has equalised the effective viscous and resistive cut-offs.","Global simulations that cannot reach extreme Reynolds numbers can still calibrate their sub-grid dynamo models against these asymptotic ratios."],"fun_headline_variants":["SSD saturates at 0.55 energy ratio and scale ratio 3 across Pm","High-Re SSD energy ratio converges to 0.55 independent of Pm","EDQNM shows Pm-independent SSD saturation energy 0.55 scale 3","Nonlinear SSDs asymptote to energy ratio 0.55 via Alfvénisation","Extreme Re yields universal SSD targets: 55% energy scale ratio 3"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The fixed eddy-damping coefficient and the assumption that the turbulence remains spectrally isotropic continue to hold all the way through nonlinear saturation; if either fails, the reported asymptotic ratios can shift.","fun_headline_variants_meta":{"raw":{"variants":["SSD saturates at 0.55 energy ratio and scale ratio 3 across Pm","High-Re SSD energy ratio converges to 0.55 independent of Pm","EDQNM shows Pm-independent SSD saturation energy 0.55 scale 3","Nonlinear SSDs asymptote to energy ratio 0.55 via Alfvénisation","Extreme Re yields universal SSD targets: 55% energy scale ratio 3"]},"model":"grok-4.5","effort":"low","cost_usd":0.005896,"raw_usage":{"total_tokens":1666,"prompt_tokens":929,"num_sources_used":0,"completion_tokens":91,"cost_in_usd_ticks":58960000,"prompt_tokens_details":{"text_tokens":929,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":646,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":929,"tokens_out":91,"duration_ms":7222,"temperature":1.0,"reasoning_tokens":646,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-12T07:19:49.369808+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A high-resolution direct numerical simulation of nonhelical MHD turbulence at Re or Rm of order 10^6 that measures a saturated magnetic-to-kinetic energy ratio substantially different from 0.55, or an integral-scale ratio substantially different from 3, would falsify the claimed asymptotes.","supporting_citations":[],"review_version":1}