Pith. sign in

REVIEW 3 major objections 1 minor 36 references

Isochrones in primordial magnetic field evolution

T0 review · 3 major / 1 minor · reviewed 2026-06-27 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read For any initial primordial magnetic field, a proper time adjustment makes early and late isochrones parallel.

desk verdict A fitted proper-time offset from Alfvén time makes early and late isochrones parallel in 2D decaying MHD simulations for varied initial peaks, but the anastrophy conservation basis lacks reported verification. read the letter →

arxiv 2606.10863 v1 pith:UQPSPSNY submitted 2026-06-09 astro-ph.CO

Axel Brandenburg , Mattia Cielo , Oksana Iarygina , Franco Vazza This is my paper · ORCID
classification astro-ph.CO
keywords primordialmagneticfieldsisochronesMHDturbulenceinversecascadeanastrophyturbulentdecayradiation-dominatedera
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper establishes that primordial magnetic fields from different magnetogenesis mechanisms can be placed on universal isochrones by determining a suitable proper time for each simulation. This matters because it allows different initial length scales and field strengths to follow the same evolutionary tracks despite independent nominal generation times. In two-dimensional simulations of decaying MHD turbulence during the radiation-dominated era, the isochrones traced by the largest processed eddy become parallel after the adjustment. The evolution follows from conservation of anastrophy, and a fit to the Alfvén time also estimates both the decay-time factor and the proper-time offset. An initial velocity field of realistic strength produces straighter tracks.

What carries the argument

The isochrone: the track of the size of the largest processed eddy, which scales with the Alfvén speed and moves toward larger scales with increasing time.

What would settle it

A simulation in which anastrophy is not conserved would produce non-parallel isochrones even after the proper-time adjustment.

Watch

Extended reading notes

Core claim

The central claim is that for any initial field, a proper time can be determined such that the isochrones at early times are parallel to those at late times. The magnetic field parameters lie on universal isochrones even for early times. This is demonstrated in two-dimensional numerical simulations of decaying MHD turbulence by varying the initial position of the peak of the magnetic energy spectrum, where the evolution is governed by the conservation of anastrophy. A fit to the Alfvén time yields an accurate estimate of the factor by which the decay time is longer than the Alfvén time, while the offset provides an estimate of the proper time added to the nominal time since the beginning of

Load-bearing premise

The evolution in the simulations is governed by the conservation of anastrophy.

Editorial extensions

If this is right

  • Varying the initial peak position of the magnetic energy spectrum still yields parallel isochrones after proper-time adjustment.
  • The fit to the Alfvén time accurately estimates the factor by which decay time exceeds Alfvén time.
  • The offset in that fit estimates the proper time to add to nominal simulation time.
  • An initial velocity field of realistic strength produces a straighter evolutionary track.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • This unification could map magnetic field evolution from any generation mechanism onto one diagram independent of nominal start time.
  • Extending the result to three-dimensional simulations would test whether universality survives without strict anastrophy conservation.
  • The same proper-time construction might apply to other inverse-cascade systems in astrophysical fluids.
  • Intergalactic magnetic field observations could be compared directly against these adjusted universal tracks.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 1 minor

Summary. The paper claims that in the radiation-dominated era, primordial magnetic fields undergoing turbulent decay exhibit universal isochrones once a proper time offset (extracted from an Alfvén-time fit) is applied, making early-time isochrones parallel to late-time ones regardless of the initial peak position in the magnetic energy spectrum. This is demonstrated via 2D MHD simulations of decaying turbulence, attributed to anastrophy conservation, with an additional note that realistic initial velocity fields produce straighter tracks.

Significance. If substantiated, the result would allow different magnetogenesis scenarios to be mapped onto a single set of universal isochrones, simplifying theoretical predictions for the evolution of primordial fields and their observational consequences. The approach of varying initial spectral peaks in simulations is a positive step toward generality, though the current support rests on unverified assumptions about conserved quantities.

major comments (3)
  1. [abstract] Abstract: the assertion that 'the evolution is governed by the conservation of anastrophy' is invoked to explain the parallelism after proper-time adjustment, yet no supporting diagnostics (e.g., time series of anastrophy for different initial peaks or with added velocity fields) are referenced. Without this verification the physical basis for universality independent of initial peak location remains untested.
  2. [abstract] Abstract (paragraph on 2D simulations): the proper-time offset is obtained from a fit to the Alfvén time performed on the same simulation data subsequently claimed to lie on universal isochrones. This procedure risks circularity; an independent determination of the offset (theoretical or from a separate diagnostic) is needed to establish that the parallelism is not an artifact of the fitting choice.
  3. [abstract] Abstract: the claim that 'the magnetic field parameters lie on universal isochrones even for early times' is load-bearing for the central result, but the manuscript provides no quantitative measure (e.g., scatter or slope comparison before/after offset) of how well the early-time tracks align across the varied initial peaks once the offset is applied.
minor comments (1)
  1. [abstract] The abstract mentions varying the initial position of the peak but does not specify the range of peak wavenumbers or the number of runs performed; adding this detail would clarify the scope of the universality claim.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We address each major comment below. The requested diagnostics and quantitative measures will be added to strengthen the manuscript.

read point-by-point responses

  1. Referee: [abstract] Abstract: the assertion that 'the evolution is governed by the conservation of anastrophy' is invoked to explain the parallelism after proper-time adjustment, yet no supporting diagnostics (e.g., time series of anastrophy for different initial peaks or with added velocity fields) are referenced. Without this verification the physical basis for universality independent of initial peak location remains untested.

    Authors: We agree that explicit verification is needed. In the revised manuscript we will add time-series plots of anastrophy (and its conservation) for the full set of runs with varied initial spectral peaks, including the cases with realistic initial velocity fields. revision: yes

  2. Referee: [abstract] Abstract (paragraph on 2D simulations): the proper-time offset is obtained from a fit to the Alfvén time performed on the same simulation data subsequently claimed to lie on universal isochrones. This procedure risks circularity; an independent determination of the offset (theoretical or from a separate diagnostic) is needed to establish that the parallelism is not an artifact of the fitting choice.

    Authors: The offset is extracted from the late-time asymptotic regime where the Alfvén-time scaling is already established and is then applied to test early-time behavior; this is not strictly circular. Nevertheless, to strengthen the result we will also provide an independent theoretical estimate of the offset derived from the initial anastrophy and Alfvén time and compare it with the fitted value. revision: partial

  3. Referee: [abstract] Abstract: the claim that 'the magnetic field parameters lie on universal isochrones even for early times' is load-bearing for the central result, but the manuscript provides no quantitative measure (e.g., scatter or slope comparison before/after offset) of how well the early-time tracks align across the varied initial peaks once the offset is applied.

    Authors: We will add quantitative diagnostics in the revised manuscript: root-mean-square scatter of the isochrone positions across initial-peak runs, and direct slope comparisons, both before and after the time offset is applied. revision: yes

Circularity Check

1 steps flagged · score 6.0 of 10

Proper-time offset fitted from simulation data then used to claim universal isochrones

  1. fitted input called prediction [abstract]
    "A fit to the Alfvén time yields an accurate estimate of the factor by which the decay time is longer than the Alfvén time, while the offset in the fit provides an additional estimate of the proper time that needs to be added to the nominal time since the beginning of each simulation. [...] The magnetic field parameters lie on universal isochrones even for early times."

    The offset that defines the proper time is extracted by fitting the Alfvén time directly to the simulation outputs for different initial peaks. Adding this fitted offset is then said to produce parallel (universal) isochrones on the same data, so the universality result is statistically forced by the choice of fit rather than predicted from first principles or external constraints.

full rationale

The paper determines the proper time offset via a fit performed on the 2D MHD simulation data (varying initial spectral peak), then adds that offset to assert that the same data lie on universal isochrones at early times. This adjustment is therefore a direct output of the fitting procedure applied to the data being analyzed, reducing the claimed universality to a consequence of the fit rather than an independent derivation. No self-citations, ansatzes, or other load-bearing reductions are exhibited in the provided text.

Assumptions & free parameters 2 free parameters · 1 assumptions · 0 invented entities

The central claim depends on two fitted quantities extracted from the simulation data and on the assumption that anastrophy is conserved in the 2D runs; no new physical entities are introduced.

free parameters (2)
  • factor by which decay time exceeds Alfvén time
    Obtained from a fit to the Alfvén time in the simulations to quantify how much longer the actual decay takes.
  • proper-time offset
    Read from the intercept of the same fit; added to nominal time so that early and late isochrones become parallel.
assumptions (1)
  • domain assumption Conservation of anastrophy governs the decay of magnetic energy in two-dimensional MHD turbulence
    Invoked to explain why shifting the initial spectral peak still yields parallel isochrones after the offset is applied.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Isochrones in primordial magnetic field evolution." pith.science (2026). https://pith.science/paper/UQPSPSNY

@misc{pith2026260610863,
  author       = {Pith},
  title        = {Pith review of: Isochrones in primordial magnetic field evolution},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UQPSPSNY}},
  note         = {Machine review of arXiv:2606.10863}
}
read the original abstract

During the radiation-dominated era of the Universe, a primordial magnetic field undergoes a turbulent decay while its length scale increases due to an inverse cascade. At later times, the size of the largest processed eddy scales with the Alfv\'en speed and it describes an isochrone that moves toward larger scales with increasing time. Different magnetogenesis mechanisms produce different initial length scales and field strengths, independently of the nominal generation time. However, we show that for any initial field, a proper time can be determined such that the isochrones at early times are parallel to those at late times. We use two-dimensional numerical simulations of decaying MHD turbulence and vary the initial position of the peak of the magnetic energy spectrum. In this case, the evolution is governed by the conservation of anastrophy. A fit to the Alfv\'en time yields an accurate estimate of the factor by which the decay time is longer than the Alfv\'en time, while the offset in the fit provides an additional estimate of the proper time that needs to be added to the nominal time since the beginning of each simulation. We also find that the presence of an initial velocity field of realistic strength helps producing a more straight track. The magnetic field parameters lie on universal isochrones even for early times.

Figures

Figures reproduced from arXiv: 2606.10863 by the authors.

Figure 1
Figure 1. Sketch of the three initial values of vA and ξM and their anticipated tracks. time t = 0, cannot lie on the same isochrone, if we expect the isochrone to lie on lines vA = ξM/t∗. Since IA ∼ v 2 Aξ 2 M, points A and C have the same value of IA, while point B has a 100 times larger value, because ξ 2 M is a hundred times larger. The initial points A, B, and C could correspond to different generation mechanisms, or eve… view at source ↗
Figure 2
Figure 2. Note that for panels (b) and (c), p(Jz/Jrms z ) de￾velops stretched exponential tails. 3.2. Evolutionary tracks In [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 2
Figure 2. Grayscale visualization of the current density along the line of sight, Jz, for Run C4u at t = 105 . The white square marks the inner part of the domain for which colorscale insets show (a) Run C4u at t = 5 × 103 , (b) Run C4u at t = 105 , and (c) Run A4 at t = 2.3 × 105 . The corresponding positions in the ξM vs. vA diagram are shown in the lower left inset for (a) and (b) on the track of Run C4u in orange, and for… view at source ↗
Figures from the paper (8 more)
Figure 3
Figure 3. Figure 3: Probability density functions p(Jz/Jrms z ) for the three insets shown in [PITH_FULL_IMAGE:figures/full_fig_p006_3.png]
Figure 4
Figure 4. Figure 4: Evolutionary tracks for Runs A1 (black), B1 (blue), C1 (red), as well as Runs A4 (green) and C4u (orange). The symbols are logarithmically spaced in time with 6 open symbols per decade and each decade is marked with a filled symbol. The dashed-dotted lines are isochron…
Figure 7
Figure 7. Figure 7: Dependence of Cξ and CE for Runs A1 (black), B1 (blue), C1 (red), as well as Runs A4 (green), C4 (yellow), and C4u (orange). Brandenburg et al. (2024) found Cξ ≈ 0.13 and CE ≈ 15 for their cases of 2D decay MHD, and Brandenburg et al. (2025) have compared earlier resul…
Figure 8
Figure 8. Figure 8: (b) it is shifted to the right. Therefore, Run B1 can￾not be made to agree with Runs A1 or C1 by a shift in the time axis. Instead, looking at [PITH_FULL_IMAGE:figures/full_fig_p007_8.png]
Figure 9
Figure 9. Figure 9: Evolution in the pq diagram for the same runs as in [PITH_FULL_IMAGE:figures/full_fig_p007_9.png]
Figure 10
Figure 10. Figure 10: Degree of anastrophy conservation quantified by the decay exponent s as a function of Lu. larger values of Lu are much closer to that line. All runs start in the lower left corner at p = q = 0, and then settle around the Alfvén line where q = 1/2 and p = 1. This is an…
Figure 12
Figure 12. Figure 12: Dependence of vA(t) vs. ξM(t) for Run B3 with u0 = 0 (black) and 0.05 (red), both with B0 = 0.1, as well as B0 = 10−2 (orange) and B0 = 10−3 again with u0 = 0.05 (blue). The vertical arrows indicate the increase of vA from the initial value to the maximum. exerted by …
Figure 13
Figure 13. Figure 13: Dependence of vA(t) vs. ξM(t) for Run C4 with α ine B = −2 (the usual case, black) and α ine B = −20 (red). The inset shows the build-up of the inertial range at times t = 0.2 (solid line), t = 100 and 300 (dotted lines), and t = 1000 (dashed line). (2025). Here, the …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 2 canonical work pages

  1. [1]

    & Jedamzik, K

    Banerjee, R. & Jedamzik, K. 2004, PhRvD, 70, 123003

    2004

  2. [2]

    Bhat, P., Zhou, M., & Loureiro, N. F. 2021, MNRAS, 501, 3074

    2021

  3. [3]

    2022, A& A, 660, A80

    Bondarenko, K., Boyarsky, A., Korochkin, A., et al. 2022, A& A, 660, A80

    2022

  4. [4]

    2020, ApJ, 901, 18

    Brandenburg, A. 2020, ApJ, 901, 18

    2020

  5. [5]

    1996, PhRvD, 54, 1 291

    Brandenburg, A., Enqvist, K., & Olesen, P. 1996, PhRvD, 54, 1 291

    1996

  6. [6]

    & Kahniashvili, T

    Brandenburg, A. & Kahniashvili, T. 2017, PhRvL, 118, 055102

    2017

  7. [7]

    2017, PhRvD, 96, 123528

    Brandenburg, A., Kahniashvili, T., Mandal, S., et al. 2017, PhRvD, 96, 123528

    2017

  8. [8]

    & Kumar, V

    Brandenburg, A. & Kumar, V. 2026, JPlPh, to be submitted

    2026

Show all 36 references
  1. [9]

    2024, A&A, 687, A18 6

    Brandenburg, A., Neronov, A., & Vazza, F. 2024, A&A, 687, A18 6

    2024

  2. [10]

    & Ntormousi, E

    Brandenburg, A. & Ntormousi, E. 2025, ApJ, 990, 223

    2025

  3. [11]

    2023, JPlPh, 8 9, 905890606

    Brandenburg, A., Sharma, R., & Vachaspati, T. 2023, JPlPh, 8 9, 905890606

    2023

  4. [12]

    2025, JPlPh, 91, E113

    Brandenburg, A., Yi, L., & W u, X. 2025, JPlPh, 91, E113

    2025

  5. [13]

    & Vazza, F

    Carretti, E. & Vazza, F. 2025, Universe, 11, 164

    2025

  6. [14]

    2001, Ph RvE, 64, 056405

    Christensson, M., Hindmarsh, M., & Brandenburg, A. 2001, Ph RvE, 64, 056405

    2001

  7. [15]

    Comisso, L., Grasso, D., & W aelbroeck, F. L. 2015, PhPl, 22, 0 42109 de Souza, R. S. & Opher, R. 2008, PhRvD, 77, 043529

    2015

  8. [16]

    2010, MNRAS, 401, 47

    Donnert, J., Dolag, K., Brunetti, G., Cassano, R., & Bonafed e, A. 2010, MNRAS, 401, 47

    2010

  9. [17]

    & Caprini, C

    Durrer, R. & Caprini, C. 2003, JCAP, 2003, 010

    2003

  10. [18]

    1998, PhRvD, 58, 12 3004

    Durrer, R., Kahniashvili, T., & Yates, A. 1998, PhRvD, 58, 12 3004

    1998

  11. [19]

    & Neronov, A

    Durrer, R. & Neronov, A. 2013, A&ARv, 21, 62

    2013

  12. [20]

    & Montgomery, D

    Fyfe, D. & Montgomery, D. 1976, JPlPh, 16, 181

    1976

  13. [21]

    2026, PhRvD, 113, 023523

    Ghosh, O., Brandenburg, A., Caprini, C., Neronov, A., & Vazz a, F. 2026, PhRvD, 113, 023523

    2026

  14. [22]

    & Rubinstein, H

    Grasso, D. & Rubinstein, H. R. 2001, Phys. Rep., 348, 163

    2001

  15. [23]

    Hosking, D. N. & Schekochihin, A. A. 2021, PhRvX, 11, 041005

    2021

  16. [24]

    Hosking, D. N. & Schekochihin, A. A. 2023, NatCo, 14, 7523

    2023

  17. [25]

    G., Brandenburg, A., & Neron ov, A

    Kahniashvili, T., Tevzadze, A. G., Brandenburg, A., & Neron ov, A. 2013, PhRvD, 87, 083007

    2013

  18. [26]

    & Brandenburg, A

    Kumar, V. & Brandenburg, A. 2026, JPlPh, submitted, arXiv:2605.18946

    2026 arXiv

  19. [27]

    2024 , arXiv e-prints, arXiv:2412.14825 Pencil Code Collaboration, Brandenburg, A., Johansen, A., et al

    Neronov, A., Vazza, F., Mtchedlidze, S., & Carretti, E. 2024 , arXiv e-prints, arXiv:2412.14825 Pencil Code Collaboration, Brandenburg, A., Johansen, A., et al. 2021, JOSS, 6, 2807

    2024

  20. [28]

    1993, Les Houches Session XL VII, 139

    Pouquet, A. 1993, Les Houches Session XL VII, 139

    1993

  21. [29]

    1976, JFM, 77, 321

    Pouquet, A., Frisch, U., & Leorat, J. 1976, JFM, 77, 321

    1976

  22. [30]

    & Banerjee, R

    Reppin, J. & Banerjee, R. 2017, PhRvE, 96, 053105

    2017

  23. [31]

    2016, Rep

    Subramanian, K. 2016, Rep. Prog. Phys., 79, 076901

    2016

  24. [32]

    2024, ApJ, 963, 135

    Tjemsland, J., Meyer, M., & Vazza, F. 2024, ApJ, 963, 135

    2024

  25. [33]

    2021, Rept

    Vachaspati, T. 2021, Rept. Prog. Phys., 84, 074901

    2021

  26. [34]

    2025, A&A, 696, A5 8

    Vazza, F., Gheller, C., Zanetti, F., et al. 2025, A&A, 696, A5 8

    2025

  27. [35]

    M., Ryu, D., Schleicher, D

    Widrow, L. M., Ryu, D., Schleicher, D. R. G., et al. 2012, SSR, 166, 37

    2012

  28. [36]

    2022, JPlPh., 88, 905 880602 Article number, page 11

    Zhou, H., Sharma, R., & Brandenburg, A. 2022, JPlPh., 88, 905 880602 Article number, page 11

    2022

Pith tools

Reviewed June 27, 2026 · model on record in the stance chip above.