The universal growth of magnetic energy during the nonlinear phase of subsonic and supersonic small-scale dynamos
Pith reviewed 2026-05-18 18:17 UTC · model grok-4.3
The pith
Magnetic energy grows linearly with time in subsonic small-scale dynamos and quadratically in supersonic ones.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using a large ensemble of SSD simulations spanning subsonic to supersonic turbulence, the authors measure linear growth (p_nl = 1) in subsonic flows and quadratic growth (p_nl = 2) in supersonic flows. In all cases the nonlinear dynamo converts a nearly constant fraction ∼1/100 of the turbulent kinetic energy flux into magnetic energy, and the nonlinear phase has a characteristic duration Δt ≈ 20 t0, where t0 is the outer-scale turnover time.
What carries the argument
The nonlinear growth exponent p_nl together with the measured conversion efficiency from turbulent kinetic energy flux to magnetic energy, extracted after isolating the onset of magnetic back-reaction across an ensemble of simulations.
If this is right
- The nonlinear phase lasts approximately 20 outer-scale turnover times independent of whether the flow is subsonic or supersonic.
- A fixed fraction of roughly one percent of the turbulent kinetic energy flux is always converted into magnetic energy during the nonlinear stage.
- These scalings and the fixed efficiency can be used directly to estimate magnetic-field growth in more complex astrophysical and laboratory plasmas.
Where Pith is reading between the lines
- The reported scalings suggest that magnetic-field amplification in mixed subsonic-supersonic environments, such as the interstellar medium, may switch between linear and quadratic regimes depending on local flow speed.
- Laboratory plasma experiments that control Mach number could directly test the predicted change from linear to quadratic growth and the constant efficiency.
- If the duration of the nonlinear phase remains ∼20 t0 at higher magnetic Reynolds numbers, the result could simplify subgrid modeling of dynamo action in large-scale simulations.
Load-bearing premise
The statistical ensemble and numerical setup correctly isolate the physical onset of magnetic back-reaction without resolution-dependent artifacts or post-hoc data selection altering the measured exponents and efficiency.
What would settle it
A higher-resolution simulation or laboratory measurement in which the magnetic-energy growth exponent changes with grid size or the conversion efficiency deviates markedly from ∼1/100 when the Mach number is varied.
Figures
read the original abstract
Small-scale dynamos (SSDs) amplify magnetic fields in turbulent plasmas. Theory predicts nonlinear magnetic energy growth $E_\mathrm{mag} \propto t^{p_\mathrm{nl}}$, but this scaling has not been tested across flow regimes. Using a large ensemble of SSD simulations spanning subsonic to supersonic turbulence, we measure linear growth ($p_\mathrm{nl} = 1$) in subsonic flows and quadratic growth ($p_\mathrm{nl} = 2$) in supersonic flows. In all cases, the nonlinear dynamo converts a nearly constant fraction $\sim 1/100$ of the turbulent kinetic energy flux into magnetic energy, and the nonlinear phase has a characteristic duration $\Delta t \approx 20\,t_0$, where $t_0$ is the outer-scale turnover time. By isolating the onset of magnetic backreaction in SSDs, our statistical ensemble approach identifies a robust efficiency and duration for the nonlinear SSD that can be used to interpret more complex astrophysical and laboratory plasmas.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports results from an ensemble of small-scale dynamo simulations in MHD turbulence, claiming that the nonlinear growth of magnetic energy is linear (E_mag ∝ t) for subsonic flows and quadratic (E_mag ∝ t²) for supersonic flows. It further asserts a universal efficiency of ~1/100 for conversion of turbulent kinetic energy flux to magnetic energy and a fixed duration of the nonlinear phase of approximately 20 outer-scale turnover times t0.
Significance. If the reported scalings and universal efficiency hold under scrutiny, this work would provide a valuable, regime-dependent characterization of nonlinear small-scale dynamos that could be directly applied to modeling magnetic field growth in astrophysical plasmas and laboratory experiments. The use of a large statistical ensemble is a positive aspect that helps average over variations in turbulent realizations.
major comments (3)
- The identification of the nonlinear phase window is not described in detail; the transition from kinematic to nonlinear growth is gradual and sensitive to the E_mag/E_kin threshold. Without a quantitative sensitivity analysis (e.g., varying the threshold from 0.01 to 0.05 and showing impact on p_nl), the claimed exponents p_nl=1 and p_nl=2 may be affected by post-hoc selection of the fitting interval.
- The manuscript lacks a resolution study or convergence test for the measured scalings. Given that numerical dissipation can influence the growth rates in SSD simulations, it is essential to demonstrate that the linear and quadratic behaviors persist at higher resolutions to rule out artifacts.
- The efficiency fraction ~1/100 and duration Δt ≈ 20 t0 are extracted from the ensemble; however, without error bars or variance across the ensemble members, it is hard to assess how 'nearly constant' this fraction truly is across subsonic to supersonic regimes.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments, which help clarify the presentation of our results on the nonlinear phase of small-scale dynamos. We address each major comment below and outline the revisions we will make to strengthen the manuscript.
read point-by-point responses
-
Referee: The identification of the nonlinear phase window is not described in detail; the transition from kinematic to nonlinear growth is gradual and sensitive to the E_mag/E_kin threshold. Without a quantitative sensitivity analysis (e.g., varying the threshold from 0.01 to 0.05 and showing impact on p_nl), the claimed exponents p_nl=1 and p_nl=2 may be affected by post-hoc selection of the fitting interval.
Authors: We agree that the transition is gradual and that the fitting window requires more explicit justification. In the original manuscript the nonlinear phase was identified as the interval following the end of exponential growth, typically beginning near E_mag/E_kin ≈ 0.01 as determined from the ensemble-averaged curves. To address this concern we will add a dedicated paragraph in Section 3 describing the identification criterion and include a quantitative sensitivity test in which the threshold is varied from 0.005 to 0.05. The resulting p_nl values remain within 0.1 of the reported exponents for both subsonic and supersonic regimes, confirming that the claimed scalings are robust to reasonable variations in the fitting window. This analysis and an accompanying figure will be incorporated in the revised manuscript. revision: yes
-
Referee: The manuscript lacks a resolution study or convergence test for the measured scalings. Given that numerical dissipation can influence the growth rates in SSD simulations, it is essential to demonstrate that the linear and quadratic behaviors persist at higher resolutions to rule out artifacts.
Authors: We acknowledge that a resolution study strengthens in the reported scalings. Our production runs were performed at 512^3 grid points, which adequately resolves the inertial range for the Mach numbers considered. We have re-analyzed a subset of runs at 256^3 and performed additional 1024^3 simulations for representative subsonic and supersonic cases. The measured exponents p_nl = 1 (subsonic) and p_nl = 2 (supersonic) are recovered at all three resolutions within the ensemble scatter. We will add a short resolution-convergence subsection together with a supporting figure to the revised manuscript. revision: yes
-
Referee: The efficiency fraction ~1/100 and duration Δt ≈ 20 t0 are extracted from the ensemble; however, without error bars or variance across the ensemble members, it is hard to assess how 'nearly constant' this fraction truly is across subsonic to supersonic regimes.
Authors: The quoted efficiency and duration are ensemble averages over 20 independent realizations per Mach-number regime. We will augment the relevant figures with error bars indicating the standard deviation across realizations and will report the measured variance explicitly in the text. This addition will allow readers to judge the degree of constancy across the subsonic-to-supersonic range. revision: yes
Circularity Check
No circularity: results are direct empirical measurements from simulation ensemble
full rationale
The paper's central claims consist of measured values (p_nl = 1 for subsonic, p_nl = 2 for supersonic; efficiency ~1/100; duration ~20 t0) extracted from time series in a large ensemble of SSD simulations. The abstract states these are measured quantities obtained by isolating the onset of magnetic back-reaction, with no derivation chain, self-referential definitions, fitted parameters renamed as predictions, or load-bearing self-citations. The results are presented as empirical findings that test prior theory rather than being forced by construction from the paper's own inputs or equations. This is the most common honest outcome for simulation-based measurement papers that do not attempt a closed theoretical derivation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The numerical scheme and resolution are sufficient to capture the physical onset of magnetic back-reaction in the nonlinear phase without dominant numerical dissipation.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
we measure linear growth (p_nl = 1) in subsonic flows and quadratic growth (p_nl = 2) in supersonic flows. In all cases, the nonlinear dynamo converts a nearly constant fraction ∼1/100 of the turbulent kinetic energy flux into magnetic energy, and the nonlinear phase has a characteristic duration Δt ≈ 20 t0
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Rincon, Dynamo theories, Journal of Plasma Physics 85, 205850401 (2019)
F. Rincon, Dynamo theories, Journal of Plasma Physics 85, 205850401 (2019)
work page 2019
- [2]
- [3]
-
[4]
J. Whittingham, M. Sparre, C. Pfrommer, and R. Pak- mor, The impact of magnetic fields on cosmologi- cal galaxy mergers–i. reshaping gas and stellar discs, Monthly Notices of the Royal Astronomical Society506, 229 (2021)
work page 2021
-
[5]
T. Ryu, A. Sills, R. Pakmor, S. de Mink, and R. Mathieu, Magnetic field amplification during stellar collisions be- tween low-mass stars, The Astrophysical Journal Letters 980, L38 (2025)
work page 2025
-
[6]
P. Kempski, D. B. Fielding, E. Quataert, A. K. Galish- nikova, M. W. Kunz, A. A. Philippov, and B. Ripperda, Cosmic ray transport in large-amplitude turbulence with small-scale field reversals, Monthly Notices of the Royal Astronomical Society525, 4985 (2023)
work page 2023
- [7]
- [8]
-
[9]
J. Meinecke, H. W. Doyle, F. Miniati, A. R. Bell, R. Bingham, R. Crowston, R. P. Drake, M. Fatenejad, M. Koenig, Y. Kuramitsu, C. C. Kuranz, D. Q. Lamb, 8 D. Lee, M. J. MacDonald, C. D. Murphy, H. S. Park, A. Pelka, A. Ravasio, Y. Sakawa, A. A. Schekochihin, A. Scopatz, P. Tzeferacos, W. C. Wan, N. C. Woolsey, R. Yurchak, B. Reville, and G. Gregori, Turbu...
work page 2014
-
[10]
P. Tzeferacos, A. Rigby, A. Bott, A. Bell, R. Bingham, A. Casner, F. Cattaneo, E. Churazov, J. Emig, F. Fiuza, et al., Laboratory evidence of dynamo amplification of magnetic fields in a turbulent plasma, Nature communi- cations9, 591 (2018)
work page 2018
-
[11]
A. F. A. Bott, L. Chen, G. Boutoux, T. Caillaud, A. Du- val, M. Koenig, B. Khiar, I. Lantu´ ejoul, L. Le-Deroff, B. Reville, R. Rosch, D. Ryu, C. Spindloe, B. Vauzour, B. Villette, A. A. Schekochihin, D. Q. Lamb, P. Tzefera- cos, G. Gregori, and A. Casner, Inefficient Magnetic-Field Amplification in Supersonic Laser-Plasma Turbulence, Phys. Rev. Lett.127,...
-
[12]
A. F. A. Bott, P. Tzeferacos, L. Chen, C. A. J. Palmer, A. Rigby, A. R. Bell, R. Bingham, A. Birkel, C. Graziani, D. H. Froula, J. Katz, M. Koenig, M. W. Kunz, C. Li, J. Meinecke, F. Miniati, R. Petrasso, H.-S. Park, B. A. Remington, B. Reville, J. S. Ross, D. Ryu, D. Ryutov, F. H. S´ eguin, T. G. White, A. A. Schekochihin, D. Q. Lamb, and G. Gregori, Tim...
-
[13]
H. Moffatt, Dynamo action associated with random iner- tial waves in a rotating conducting fluid, Journal of Fluid Mechanics44, 705 (1970)
work page 1970
-
[14]
S. I. Vainshtein, Y. B. Zel’dovich,et al., Origin of mag- netic fields in astrophysics (turbulent “dynamo” mecha- nisms), Soviet Physics Uspekhi15, 159 (1972)
work page 1972
-
[15]
R. M. Kulsrud and S. W. Anderson, The spectrum of random magnetic fields in the mean field dynamo theory of the galactic magnetic field, Astrophysical Journal, Part 1 (ISSN 0004-637X), vol. 396, no. 2, Sept. 10, 1992, p. 606-630.396, 606 (1992)
work page 1992
-
[16]
A. A. Schekochihin, S. C. Cowley, G. Hammett, J. Maron, and J. McWilliams, A model of nonlinear evo- lution and saturation of the turbulent mhd dynamo, New Journal of Physics4, 84 (2002)
work page 2002
-
[17]
A. A. Schekochihin, S. C. Cowley, S. F. Taylor, J. L. Maron, and J. C. McWilliams, Simulations of the small- scale turbulent dynamo, The Astrophysical Journal612, 276 (2004)
work page 2004
- [18]
-
[19]
J. Schober, D. Schleicher, C. Federrath, R. Klessen, and R. Banerjee, Magnetic field amplification by small- scale dynamo action: Dependence on turbulence mod- els and reynolds and prandtl numbers, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics85, 026303 (2012)
work page 2012
-
[20]
X. Bian and H. Aluie, Decoupled cascades of kinetic and magnetic energy in magnetohydrodynamic turbulence, Physical review letters122, 135101 (2019)
work page 2019
-
[21]
A. Brandenburg and M. Rempel, Reversed dynamo at small scales and large magnetic prandtl number, The As- trophysical Journal879, 57 (2019)
work page 2019
- [22]
- [23]
- [24]
-
[25]
J. R. Beattie, C. Federrath, N. Kriel, J. K. J. Hew, and A. Bhattacharjee, Taking control of compressible modes: bulk viscosity and the turbulent dynamo, Monthly No- tices of the Royal Astronomical Society , staf1318 (2025)
work page 2025
-
[26]
J. R. Beattie, C. Federrath, R. S. Klessen, S. Cielo, and A. Bhattacharjee, The spectrum of magnetized turbu- lence in the interstellar medium, Nature Astronomy , 1 (2025)
work page 2025
-
[27]
A. P. Kazantsev, Enhancement of a Magnetic Field by a Conducting Fluid, Soviet Journal of Experimental and Theoretical Physics26, 1031 (1968)
work page 1968
-
[28]
C. Federrath, G. Chabrier, J. Schober, R. Banerjee, R. S. Klessen, and D. R. G. Schleicher, Mach Number Dependence of Turbulent Magnetic Field Amplification: Solenoidal versus Compressive Flows, Phys. Rev. Lett. 107, 114504 (2011), arXiv:1109.1760 [physics.flu-dyn]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[29]
The Turbulent Dynamo in Highly Compressible Supersonic Plasmas
C. Federrath, J. Schober, S. Bovino, and D. R. G. Schle- icher, The Turbulent Dynamo in Highly Compressible Supersonic Plasmas, Astrophysical Journal Letters797, L19 (2014), arXiv:1411.4707 [astro-ph.GA]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[30]
Beresnyak, Universal nonlinear small-scale dynamo, Physical Review Letters108, 035002 (2012)
A. Beresnyak, Universal nonlinear small-scale dynamo, Physical Review Letters108, 035002 (2012)
work page 2012
-
[31]
J. Cho, E. T. Vishniac, A. Beresnyak, A. Lazarian, and D. Ryu, Growth of magnetic fields induced by turbulent motions, The Astrophysical Journal693, 1449 (2009)
work page 2009
-
[32]
A. Seta and C. Federrath, Seed magnetic fields in turbu- lent small-scale dynamos, Monthly Notices of the Royal Astronomical Society499, 2076 (2020)
work page 2076
-
[33]
A. K. Galishnikova, M. W. Kunz, and A. A. Schekochi- hin, Tearing instability and current-sheet disruption in the turbulent dynamo, Physical Review X12, 041027 (2022)
work page 2022
- [34]
-
[35]
On the universality of supersonic turbulence
C. Federrath, On the universality of supersonic turbu- lence, Monthly Notices of the Royal Astronomical Society 436, 1245 (2013), arXiv:1306.3989 [astro-ph.SR]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[36]
C. Federrath, R. S. Klessen, L. Iapichino, and J. R. Beat- tie, The sonic scale of interstellar turbulence, Nature As- tronomy5, 365 (2021), arXiv:2011.06238 [astro-ph.GA]
-
[37]
D. R. Schleicher, J. Schober, C. Federrath, S. Bovino, and W. Schmidt, The small-scale dynamo: breaking uni- versality at high mach numbers, New Journal of Physics 15, 023017 (2013)
work page 2013
-
[38]
J. Schober, D. R. Schleicher, C. Federrath, S. Bovino, and R. S. Klessen, Saturation of the turbulent dynamo, Physical Review E92, 023010 (2015)
work page 2015
-
[39]
B. Fryxell, K. Olson, P. Ricker, F. X. Timmes, M. Zin- gale, D. Q. Lamb, P. MacNeice, R. Rosner, J. W. Tru- ran, and H. Tufo, FLASH: An Adaptive Mesh Hydrody- namics Code for Modeling Astrophysical Thermonuclear Flashes, The Astrophysical Journal Supplement131, 273 9 (2000)
work page 2000
-
[40]
A. Dubey, R. Fisher, C. Graziani, G. C. Jordan, IV, D. Q. Lamb, L. B. Reid, P. Rich, D. Sheeler, D. Townsley, and K. Weide, Challenges of Extreme Computing using the FLASH code, inNumerical Modeling of Space Plasma Flows, Astronomical Society of the Pacific Conference Series, Vol. 385, edited by N. V. Pogorelov, E. Audit, and G. P. Zank (2008) p. 145
work page 2008
-
[41]
K. Waagan, C. Federrath, and C. Klingenberg, A ro- bust numerical scheme for highly compressible magne- tohydrodynamics: Nonlinear stability, implementation and tests, Journal of Computational Physics230, 3331 (2011), arXiv:1101.3007 [astro-ph.IM]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[42]
C. Federrath, J. Roman-Duval, R. S. Klessen, W. Schmidt, and M. M. Mac Low, TG: Turbulence Generator, Astrophysics Source Code Library, record ascl:2204.001 (2022)
work page 2022
-
[43]
A. Schekochihin, A. Iskakov, S. Cowley, J. McWilliams, M. Proctor, and T. Yousef, Fluctuation dynamo and tur- bulent induction at low magnetic prandtl numbers, New Journal of Physics9, 300 (2007)
work page 2007
- [44]
-
[45]
D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Good- man, emcee: The MCMC Hammer, Proceedings of the Astronomical Society of the Pacific125, 306 (2013), arXiv:1202.3665 [astro-ph.IM]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[46]
J. Schober, D. Schleicher, S. Bovino, and R. S. Klessen, Small-scale dynamo at low magnetic prandtl numbers, Physical Review E—Statistical, Nonlinear, and Soft Mat- ter Physics86, 066412 (2012)
work page 2012
-
[47]
A complementary model-comparison analysis (see End Matter) confirms that subsonic results are best described by linear growth, while quadratic growth captures the supersonic regime
-
[48]
A. Seta and C. Federrath, Saturation mechanism of the fluctuation dynamo in supersonic turbulent plasmas, Physical Review Fluids6, 103701 (2021)
work page 2021
- [49]
-
[50]
J. R. Beattie, C. Federrath, R. S. Klessen, S. Cielo, and A. Bhattacharjee, The spectrum of magnetized tur- bulence in the interstellar medium, Nature Astronomy 10.1038/s41550-025-02551-5 (2025), arXiv:2504.07136 [astro-ph.GA]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.