pith. machine review for the scientific record. sign in

arxiv: 2604.00454 · v2 · submitted 2026-04-01 · 🌌 astro-ph.CO · gr-qc

Recognition: 2 theorem links

· Lean Theorem

Relic Magnetic Fields from Non-Adiabatic Photon Freeze-Out at Recombination

Hyeong-Chan Kim
Authors on Pith no claims yet

Pith reviewed 2026-05-13 22:30 UTC · model grok-4.3

classification 🌌 astro-ph.CO gr-qc
keywords primordial electromagnetic relicrecombinationThomson relaxationnon-adiabatic squeezingmagnetic field spectrumphoton freeze-outcosmic magnetic fields
0
0 comments X

The pith

Finite Thomson relaxation at recombination produces a frozen electromagnetic relic peaked at 10-20 Mpc scales today.

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

The paper examines how the cosmic photon gas, treated as an open system coupled to the electron plasma, departs from thermal equilibrium because the Thomson relaxation rate remains finite. This departure triggers non-adiabatic mode squeezing that then freezes out once the relaxation rate drops too low to maintain tracking. The resulting non-equilibrium electromagnetic relic projects onto the magnetic sector with a spectrum whose peak is set by the combination k³S_k rather than the squeezing parameter alone. In representative cases the peak today corresponds to scales of 10-20 Mpc, although the field amplitude is extremely small. A reader would care because the mechanism generates a primordial electromagnetic relic using only the standard recombination process without additional fields or particles.

Core claim

Treating the photon sector as an open system coupled to the electron plasma, a finite Thomson relaxation rate generates a departure from instantaneous thermal equilibrium that produces non-adiabatic mode squeezing. As the relaxation rate rapidly decreases across recombination the system loses the ability to amplify the deviation further, so the squeezing freezes out at a small but finite value. This dynamics is recast by a canonical transformation into a forced oscillator with a smooth effective potential. Projecting the relic onto the magnetic sector yields a spectrum controlled by the weighted combination k³S_k whose characteristic peak today lies at scales of order 10-20 Mpc while the com

What carries the argument

A canonical transformation that recasts the reduced evolution equation into a forced oscillator with a smooth effective potential, thereby capturing the origin of the squeezing and selecting the relic scale through the combination k³S_k.

Load-bearing premise

The Thomson relaxation rate decreases rapidly enough across recombination to produce a narrow transition layer between adiabatic tracking and post-relaxation freeze-out.

What would settle it

A measurement of the intergalactic magnetic power spectrum at comoving scales around 10-20 Mpc that either detects or places a firm upper limit below the predicted tiny amplitude and specific k³S_k shape.

read the original abstract

We propose a new mechanism for generating a primordial electromagnetic relic during the recombination--decoupling transition, based on the rate-dependent thermodynamics of the cosmic photon gas. Treating the photon sector as an open system coupled to the electron plasma, we show that a finite Thomson relaxation rate generates a departure from instantaneous thermal equilibrium, leading to non-adiabatic mode squeezing. As this relaxation rate rapidly decreases across recombination, the system quickly loses the ability to further amplify the deviation, and the squeezing freezes out at a small but finite value. This dynamics is naturally described as a narrow transition layer between an adiabatic tracking regime and a post-relaxation freeze-out regime. By a canonical transformation, the reduced evolution equation is recast into a forced oscillator with a smooth effective potential, clarifying the origin of the squeezing and the selection of the relic scale. Projecting the resulting non-equilibrium electromagnetic relic onto the magnetic sector, we derive the corresponding spectrum and show that its characteristic peak is controlled not by the squeezing parameter alone but by the weighted combination \(k^3\mathscr S_k\). In representative realizations, the peak corresponds today to scales of order \(10\)--\(20\) Mpc, while the present-day field amplitude remains extremely small. The mechanism is therefore better viewed as a source of a frozen non-equilibrium electromagnetic relic than as a complete explanation of the observed cosmic magnetic fields.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes a mechanism for generating a primordial electromagnetic relic during the recombination-decoupling transition by treating the photon sector as an open system coupled to the electron plasma. A finite Thomson relaxation rate produces a departure from instantaneous thermal equilibrium and non-adiabatic mode squeezing; as the rate decreases rapidly, the squeezing freezes out at a small finite value. The dynamics is recast via canonical transformation into a forced oscillator with a smooth effective potential. Projecting the resulting relic onto the magnetic sector yields a spectrum whose peak is controlled by the weighted combination k³S_k and corresponds today to scales of order 10–20 Mpc, although the present-day amplitude remains extremely small. The mechanism is presented as a source of frozen non-equilibrium relics rather than a complete explanation of observed cosmic magnetic fields.

Significance. If the central derivation holds, the work supplies a first-principles, parameter-free route to a relic electromagnetic spectrum whose characteristic scale is selected by the freeze-out dynamics of the Thomson rate rather than by external parameters or initial conditions. The canonical transformation to a forced oscillator clarifies the origin of the squeezing and provides an explicit prediction for the peak location via k³S_k. However, the extremely small present-day amplitude limits the mechanism’s relevance to explaining observed intergalactic fields, so its primary value lies in identifying a previously overlooked non-equilibrium relic from standard recombination physics.

major comments (2)
  1. [Abstract] The assumption that the Thomson relaxation rate decreases rapidly enough across recombination to produce a narrow transition layer between adiabatic tracking and post-relaxation freeze-out (Abstract) is load-bearing for the canonical transformation to a forced oscillator with a smooth effective potential. If the decrease is only gradual on the timescale of the relevant modes, the effective potential is no longer smooth, the squeezing parameter S_k acquires additional time dependence, and the claimed control of the peak location by k³S_k is no longer guaranteed.
  2. [Abstract] The projection onto the magnetic sector that yields the spectrum controlled by k³S_k (Abstract) is stated without the intermediate steps that relate the forced-oscillator solution to the electromagnetic field components. Explicit equations showing how the weighting k³ arises and how the present-day amplitude is obtained are required to verify that the 10–20 Mpc peak is robust rather than an artifact of the representative realizations.
minor comments (2)
  1. [Abstract] The abstract refers to “representative realizations” without specifying the range of parameters varied or the criterion used to select them; a brief description of the parameter space explored would clarify the robustness of the 10–20 Mpc scale.
  2. [Abstract] The statement that the present-day field amplitude “remains extremely small” is given without a numerical estimate or comparison to current observational upper bounds on primordial magnetic fields.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for identifying two points where the presentation in the abstract and main text requires strengthening. Both comments concern the justification of the rapid-transition approximation and the explicit steps linking the forced-oscillator solution to the magnetic spectrum. We have revised the manuscript to supply the missing timescale comparison and the intermediate projection equations; the central physical picture remains unchanged.

read point-by-point responses

  1. Referee: [Abstract] The assumption that the Thomson relaxation rate decreases rapidly enough across recombination to produce a narrow transition layer between adiabatic tracking and post-relaxation freeze-out (Abstract) is load-bearing for the canonical transformation to a forced oscillator with a smooth effective potential. If the decrease is only gradual on the timescale of the relevant modes, the effective potential is no longer smooth, the squeezing parameter S_k acquires additional time dependence, and the claimed control of the peak location by k³S_k is no longer guaranteed.

    Authors: We agree that the rapidity of the Thomson-rate drop is essential for the smoothness of the effective potential. In the original manuscript we stated that the rate falls by several orders of magnitude over a narrow interval, but we did not quantify the comparison with the mode oscillation period. The revised version adds a dedicated paragraph in Section II together with a new figure that plots the Thomson rate against the Hubble rate and against the oscillation frequency k/a for the wave-numbers that dominate the k³S_k peak. The comparison shows that, for the relevant modes (today’s scales 10–20 Mpc), the transition layer is narrower than one oscillation period, validating the sudden-freeze-out approximation and the smoothness of the potential. For much longer wavelengths the approximation weakens, but those modes contribute negligibly to the peak; we have added a brief caveat to this effect. revision: yes

  2. Referee: [Abstract] The projection onto the magnetic sector that yields the spectrum controlled by k³S_k (Abstract) is stated without the intermediate steps that relate the forced-oscillator solution to the electromagnetic field components. Explicit equations showing how the weighting k³ arises and how the present-day amplitude is obtained are required to verify that the 10–20 Mpc peak is robust rather than an artifact of the representative realizations.

    Authors: We accept that the projection step was presented too concisely. The revised manuscript inserts a new Appendix B that starts from the canonical transformation, writes the explicit solution for the quadrature operators, projects onto the transverse magnetic field via the standard electromagnetic decomposition, and derives the factor k³ from the mode normalization (the magnetic energy density scales as k² times the squared amplitude). The present-day amplitude is obtained by red-shifting the frozen squeezing parameter through the subsequent matter-dominated expansion; the resulting expression for B(k) is given in closed form. With these steps the location of the k³S_k maximum is now traceable to the freeze-out dynamics rather than to any ad-hoc choice of representative realizations. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation starts from the open-system treatment of the photon gas with finite Thomson relaxation rate, identifies non-adiabatic squeezing during the rapid drop across recombination, and applies a canonical transformation to obtain the forced-oscillator form with smooth effective potential. The spectrum peak is then stated to be set by the combination k³S_k where S_k is the freeze-out squeezing value obtained from that dynamics. No equation or step in the provided text reduces S_k or the peak location to a fitted input, a self-citation, or a redefinition of the output; the chain remains independent of the final relic spectrum it produces. The narrow-transition assumption is an explicit physical premise rather than a hidden tautology.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard cosmological background evolution plus the modeling choice of treating the photon sector as an open system with a finite, rapidly falling Thomson relaxation rate. No new particles or forces are introduced.

axioms (2)
  • domain assumption Standard FLRW expansion and recombination physics with known Thomson scattering rate
    The paper invokes the usual recombination transition without deriving it.
  • domain assumption Photon gas can be treated as an open quantum system whose relaxation rate controls adiabaticity
    This is the key modeling step that enables the non-adiabatic squeezing description.

pith-pipeline@v0.9.0 · 5542 in / 1465 out tokens · 73200 ms · 2026-05-13T22:30:51.154979+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

41 extracted references · 41 canonical work pages · 1 internal anchor

  1. [1]

    Ionization fractionx e(R) The recombination history is characterized by a rapid crossover from an almost fully ionized state to a partially neutral plasma, followed by a slow relaxation toward a residual ionization fraction. We model this behavior by a smooth interpolation, xe(R)≃x tail e (R) + xpre e −x tail e (R∗) 1 2 1−tanh R−R ∗ δ ,(B1) whereR ∗ denot...

    work page

  2. [2]

    Bath-branch coefficientg(R) In the main text, the reduced Gaussian state is parametrized as ρ(t)∝exp " − ˆI(t) T0 # ,(B6) whereT 0 is a fixed reference temperature. We chooseT 0 such that, at the initial reference epoch, the thermal Gaussian state satisfies ˆI= ˆH.(B7) Equivalently, if the initial state is thermal with respect to the mode Hamiltonian, ρin...

    work page

  3. [3]

    Low-frequency regime When the effective potentialQ(R) is small and the mode varies slowly, the dynamics is dominated by the source term. In this regime one obtains the quasi-static approximation yk ≃ S(R) Q(R) , u k ≃ F(R) Q(R) .(D2) Thus the deviation from equilibrium is a forcing-induced lag rather than an oscillatory response. Depending on parameters,Q...

    work page

  4. [4]

    High-frequency regime WhenQ(R)>0 and varies slowly, the solution admits a WKB form yk ∼Q −1/4 exp ±i Z Rp Q dR ! .(D3) The physical deviation is then uk ∼R 1/4 exp " − Z R λ dR # Q−1/4 exp ±i Z Rp Q dR ! .(D4) In this regime, the dynamics is dominated by oscillatory freeze-out, while the source contribution is parametrically suppressed, u(p) k ∼ F(R) Q(R)...

    work page

  5. [5]

    Durrer and A

    R. Durrer and A. Neronov,Cosmological magnetic fields: their generation, evolution and observation, Astron. Astrophys. Rev.21, 62 (2013)

    work page 2013

  6. [6]

    Subramanian,The origin, evolution and signatures of primordial magnetic fields, Rep

    K. Subramanian,The origin, evolution and signatures of primordial magnetic fields, Rep. Prog. Phys.79, 076901 (2016)

    work page 2016

  7. [7]

    L. M. Widrow,Origin of galactic and extragalactic magnetic fields, Rev. Mod. Phys.74, 775 (2002)

    work page 2002

  8. [8]

    Neronov and I

    A. Neronov and I. Vovk,Evidence for strong extragalactic magnetic fields from Fermi observations of TeV blazars, Science 328, 73 (2010)

    work page 2010

  9. [9]

    Tavecchioet al.,The intergalactic magnetic field constrained by Fermi/LAT observations of the TeV blazar 1ES 0229+200, Mon

    F. Tavecchioet al.,The intergalactic magnetic field constrained by Fermi/LAT observations of the TeV blazar 1ES 0229+200, Mon. Not. R. Astron. Soc.406, L70 (2010)

    work page 2010

  10. [10]

    Ackermannet al.(Fermi-LAT),Search for signatures of the cosmic-ray sea, blazar jets, and magnetized cosmic voids in the isotropic gamma-ray background, Astrophys

    M. Ackermannet al.(Fermi-LAT),Search for signatures of the cosmic-ray sea, blazar jets, and magnetized cosmic voids in the isotropic gamma-ray background, Astrophys. J.857, 49 (2018). 16

    work page 2018

  11. [11]

    Beck,Cosmic magnetic fields: Observations and prospects, Space Sci

    R. Beck,Cosmic magnetic fields: Observations and prospects, Space Sci. Rev.166, 215 (2012)

    work page 2012

  12. [12]

    Vernstromet al.,Evidence for large-scale magnetic fields in the filaments of the cosmic web, Mon

    T. Vernstromet al.,Evidence for large-scale magnetic fields in the filaments of the cosmic web, Mon. Not. R. Astron. Soc. 505, 4178 (2021)

    work page 2021

  13. [13]

    M. S. Turner and L. M. Widrow,Inflation-produced, large-scale magnetic fields, Phys. Rev. D37, 2743 (1988)

    work page 1988

  14. [14]

    Vachaspati,Magnetic fields from cosmological phase transitions, Phys

    T. Vachaspati,Magnetic fields from cosmological phase transitions, Phys. Lett. B265, 258 (1991)

    work page 1991

  15. [15]

    G. Sigl, A. V. Olinto, and K. Jedamzik,Primordial magnetic fields from cosmological first order phase transitions, Phys. Rev. D55, 4582 (1997)

    work page 1997

  16. [16]

    Biermann, ¨Uber den Ursprung der Magnetfelder auf Sternen und im interstellaren Raum, Z

    L. Biermann, ¨Uber den Ursprung der Magnetfelder auf Sternen und im interstellaren Raum, Z. Naturforsch. A5, 65 (1950)

    work page 1950

  17. [17]

    Matarreseet al.,Large-scale magnetic fields from density perturbations, Phys

    S. Matarreseet al.,Large-scale magnetic fields from density perturbations, Phys. Rev. D71, 043502 (2005)

    work page 2005

  18. [18]

    Fenuet al.,The seed magnetic field generated during recombination, J

    E. Fenuet al.,The seed magnetic field generated during recombination, J. Cosmol. Astropart. Phys.2011, 017 (2011)

    work page 2011

  19. [19]

    Ichikiet al.,Generation of magnetic fields at recombination epoch, Science311, 827 (2006)

    K. Ichikiet al.,Generation of magnetic fields at recombination epoch, Science311, 827 (2006)

    work page 2006

  20. [20]

    Sagaet al.,Cosmological magnetic field from second-order perturbations, Phys

    S. Sagaet al.,Cosmological magnetic field from second-order perturbations, Phys. Rev. D91, 123510 (2015)

    work page 2015

  21. [21]

    P. J. E. Peebles,Recombination of the Primeval Plasma, Astrophys. J.153, 1 (1968)

    work page 1968

  22. [22]

    Seager, D

    S. Seager, D. Sasselov, and D. Scott,A new calculation of the recombination epoch, Astrophys. J.523, L1 (1999)

    work page 1999

  23. [23]

    Seager, D

    S. Seager, D. Sasselov, and D. Scott,How exactly did the Universe become neutral?, Astrophys. J. Suppl.128, 407 (2000)

    work page 2000

  24. [24]

    Hu and N

    W. Hu and N. Sugiyama,Small scale perturbations in a general MDM cosmology, Astrophys. J.471, 542 (1996)

    work page 1996

  25. [25]

    Aghanimet al.),Planck 2018 results

    Planck Collaboration (N. Aghanimet al.),Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys.641, A6 (2020)

    work page 2018

  26. [26]

    Akramiet al.),Planck 2018 results

    Planck Collaboration (Y. Akramiet al.),Planck 2018 results. X. Constraints on inflation, Astron. Astrophys.641, A10 (2020)

    work page 2018

  27. [27]

    D. J. Fixsenet al.(COBE/FIRAS),The cosmic microwave background spectrum from the full COBE/FIRAS data set, Astrophys. J.473, 576 (1996)

    work page 1996

  28. [28]

    R. A. Sunyaev and Ya. B. Zel’dovich,Small-scale fluctuations of relic radiation, Astrophys. Space Sci.7, 3 (1970)

    work page 1970

  29. [29]

    Chluba and R

    J. Chluba and R. A. Sunyaev,The evolution of CMB spectral distortions in the early Universe, Mon. Not. R. Astron. Soc. 419, 1294 (2012)

    work page 2012

  30. [30]

    Gorini, A

    V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, J. Math. Phys.17, 821 (1976)

    work page 1976

  31. [31]

    Lindblad, Commun

    G. Lindblad, Commun. Math. Phys.48, 119 (1976)

    work page 1976

  32. [32]

    Breuer and F

    H.-P. Breuer and F. Petruccione,The Theory of Open Quantum Systems(Oxford University Press, Oxford, 2002)

    work page 2002

  33. [33]

    A. O. Caldeira and A. J. Leggett,Quantum tunnelling in a dissipative system, Ann. Phys.149, 374 (1983)

    work page 1983

  34. [34]

    Parker,Particle creation in expanding universes, Phys

    L. Parker,Particle creation in expanding universes, Phys. Rev. Lett.21, 562 (1968)

    work page 1968

  35. [35]

    L. P. Grishchuk and Y. V. Sidorov,Squeezed quantum states of relic gravitons, Phys. Rev. D42, 3413 (1990)

    work page 1990

  36. [36]

    Polarski and A

    D. Polarski and A. A. Starobinsky,Semiclassicality and decoherence of cosmological perturbations, Class. Quantum Grav. 13, 377 (1996)

    work page 1996

  37. [37]

    H. C. Kim and Y. Lee,Rate-Dependent Internal Energy from Detailed-Balance Relaxation,[arXiv:2603.00189 [cond- mat.stat-mech]]

    work page internal anchor Pith review Pith/arXiv arXiv

  38. [38]

    Caprini, R

    C. Caprini, R. Durrer, and T. Kahniashvili,Constraints on the power spectrum of primordial magnetic fields, Phys. Rev. D69, 063006 (2004)

    work page 2004

  39. [39]

    Kogutet al.,The Primordial Inflation Explorer (PIXIE): A nulling polarimeter for cosmic microwave background observations, J

    A. Kogutet al.,The Primordial Inflation Explorer (PIXIE): A nulling polarimeter for cosmic microwave background observations, J. Cosmol. Astropart. Phys.2011, 025 (2011)

    work page 2011

  40. [40]

    Allyset al.(LiteBIRD Collaboration),Probing cosmic inflation with the LiteBIRD cosmic microwave background po- larization survey, Prog

    E. Allyset al.(LiteBIRD Collaboration),Probing cosmic inflation with the LiteBIRD cosmic microwave background po- larization survey, Prog. Theor. Exp. Phys.2023, 042F01 (2023)

    work page 2023

  41. [41]

    Hyeong-Chan Kim and Youngone Lee,Nonadiabaticity of quantum harmonic oscillators,Phys. Lett. A,430, 127974 (2022). https://doi.org/10.1016/j.physleta.2022.127974

    work page doi:10.1016/j.physleta.2022.127974 2022