arxiv: 2604.24446 · v1 · submitted 2026-04-27 · ✦ hep-ph

Recognition: unknown

Energy spectrum of magnetic fields from electroweak symmetry breaking

K\'aroly Seller , G\"unter Sigl

Authors on Pith no claims yet

Pith reviewed 2026-05-08 02:44 UTC · model grok-4.3

classification ✦ hep-ph

keywords electroweak symmetry breakingprimordial magnetic fieldsHiggs field configurationsanalytic spectrumstatistical isotropycausalityearly universephase transitions

0 comments

The pith

Random configurations of an inhomogeneous Higgs field produce an analytic spectrum for magnetic fields generated during electroweak symmetry breaking.

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

The paper establishes that the energy spectrum of magnetic fields from the electroweak phase transition can be derived analytically by treating the initial Higgs field as random and inhomogeneous. This exploits statistical isotropy and causality to define correlation functions that match numerical results without adjustable parameters. A new continuous-field simulation method is introduced to capture small-scale structures that lattice methods miss. If correct, the approach replaces expensive simulations with direct analytic expressions that can be compared immediately to data.

Core claim

By modeling the Higgs field during electroweak symmetry breaking as random inhomogeneous configurations and imposing statistical isotropy together with causality, the authors obtain an essentially analytic expression for the magnetic-field energy spectrum. This spectrum is obtained without lattice discretization and is shown to agree with a new continuous-field simulation framework that resolves finer scales than previous lattice runs. General correlation functions derived from isotropy and causality are fitted to the simulation outputs and compared directly to the analytic prediction.

What carries the argument

Analytic correlation functions obtained by applying statistical isotropy and causality to random inhomogeneous Higgs field configurations.

If this is right

The magnetic spectrum becomes available for immediate comparison to simulation data without post-processing adjustments.
Continuous-field simulations resolve small-scale magnetic structures inaccessible to prior lattice calculations.
General correlation functions derived from isotropy and causality can be fitted to numerical data and tested against the analytic result.
The framework removes the need for repeated costly lattice runs when exploring parameter variations in the early Universe.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The analytic method could be applied to other cosmological phase transitions where random scalar fields are present.
The continuous simulation approach may allow tracking of magnetic field evolution over longer times or larger volumes than lattice grids permit.
If the spectrum holds, it provides a cleaner starting point for evolving primordial fields to later epochs such as recombination.

Load-bearing premise

That random initial Higgs configurations plus statistical isotropy and causality are sufficient to produce a parameter-free analytic spectrum that matches simulations without further tuning.

What would settle it

A high-resolution continuous simulation with the same random initial Higgs conditions that produces a magnetic power spectrum deviating from the analytic formula at wavenumbers accessible to both methods.

Figures

Figures reproduced from arXiv: 2604.24446 by G\"unter Sigl, K\'aroly Seller.

**Figure 1.** Figure 1: FIG. 1. Energy spectrum calculated in the 4PA case on an view at source ↗

**Figure 2.** Figure 2: FIG. 2. An example slice of the magnitude of the magnetic view at source ↗

**Figure 3.** Figure 3: FIG. 3. Spatial correlation function view at source ↗

**Figure 4.** Figure 4: FIG. 4. Comparison between the simulated correlation func view at source ↗

read the original abstract

We study the magnetic fields produced in the early Universe during the electroweak symmetry breaking by considering random configurations of an inhomogeneous Higgs field. By exploiting the inherent randomness of the initial configurations the spectrum of the produced magnetic field is essentially analytic, which bypasses the need for costly lattice simulations. On the numerical side, we devise a simulation framework which results in continuous fields capable of resolving the small-scale structure of the fields that was inaccessible for the lattice-based calculation. Finally, by revisiting the effects of statistical isotropy and causality on the spectrum, we define general correlation functions that are then fitted to the simulation data and compared to the analytic results.

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.

Desk Editor's Note private letter to a colleague

Random Higgs configurations yield an analytic magnetic spectrum at electroweak breaking, backed by a continuous-field simulation, but the fit to correlation functions adds a parameter that undercuts the fully independent claim.

read the letter

The main advance is deriving an essentially analytic energy spectrum for the magnetic fields by assuming random initial configurations of an inhomogeneous Higgs field. They combine this with a new simulation framework that produces continuous fields instead of discrete lattices, which lets them resolve smaller scales than earlier work. Revisiting statistical isotropy and causality to set up general correlation functions is also a useful step. These pieces together offer a quicker route than pure lattice runs for this part of magnetogenesis calculations. The continuous simulation method is a clear technical improvement and deserves credit for accessing structure that standard approaches miss. The analytic spectrum itself is the part that could reduce reliance on heavy numerics if it holds up. The soft spot sits in the comparison. They fit the general correlation functions to the simulation data and then check against the analytic prediction. Even with one parameter, this step means the agreement is not a pure test of the randomness assumption alone. The abstract presents the analytic result as coming directly from the initial conditions and symmetries, yet the fit introduces data-driven adjustment that needs to be shown as non-essential. If the shapes only line up after calibration, the bypass of simulations is less complete than claimed. This paper is for cosmologists and particle physicists working on primordial magnetic fields and electroweak phase transitions. Readers who need faster analytic estimates or improved small-scale numerics will find concrete value in the spectrum formula and the simulation technique. The central idea shows clear thinking and honest engagement with the literature, so it deserves a serious referee rather than a desk rejection. I would recommend sending it for peer review, with the main request being a clearer demonstration that the analytic result stands without the fit doing the heavy lifting.

Referee Report

1 major / 0 minor

Summary. The manuscript claims that random configurations of an inhomogeneous Higgs field during electroweak symmetry breaking, combined with statistical isotropy and causality, permit an essentially analytic derivation of the magnetic field energy spectrum. This approach is said to bypass costly lattice simulations. The authors introduce a continuous-field simulation framework to resolve small-scale structures inaccessible to lattice methods and define general correlation functions (from isotropy and causality) that are fitted to the simulation data for comparison with the analytic predictions.

Significance. If the analytic spectrum can be shown to follow directly from the random Higgs configurations without dependence on fitted parameters, the result would be significant for early-Universe cosmology. It would enable efficient, potentially parameter-free predictions of primordial magnetic fields from the electroweak phase transition, reducing reliance on full numerical simulations and aiding constraints on magnetogenesis models. The continuous-field simulation method is a constructive technical contribution that addresses resolution limits of lattice approaches.

major comments (1)

[Abstract] Abstract: The claim that the spectrum is 'essentially analytic' by exploiting the inherent randomness of initial configurations (thereby bypassing lattice simulations) is in tension with the subsequent step of fitting general correlation functions to simulation data before comparing them to the analytic results. This fit introduces free parameters, so the manuscript must demonstrate (in the derivation section) that the final analytic spectrum is independent of any such calibration and follows solely from the randomness, isotropy, and causality assumptions. Without this clarification the central claim that the approach is essentially analytic and simulation-independent is not yet established.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive major comment. We address it point by point below and have revised the manuscript to improve clarity on the analytic derivation.

read point-by-point responses

Referee: [Abstract] Abstract: The claim that the spectrum is 'essentially analytic' by exploiting the inherent randomness of initial configurations (thereby bypassing lattice simulations) is in tension with the subsequent step of fitting general correlation functions to simulation data before comparing them to the analytic results. This fit introduces free parameters, so the manuscript must demonstrate (in the derivation section) that the final analytic spectrum is independent of any such calibration and follows solely from the randomness, isotropy, and causality assumptions. Without this clarification the central claim that the approach is essentially analytic and simulation-independent is not yet established.

Authors: We thank the referee for identifying this point of potential confusion. The analytic derivation in Section 3 proceeds directly from the random inhomogeneous Higgs configurations together with the requirements of statistical isotropy and causality. It yields an explicit functional form for the magnetic-field energy spectrum expressed solely in terms of the Higgs two-point correlation functions; no free parameters or calibration enter at any stage of this derivation. The general correlation functions themselves are constrained by isotropy and causality alone. The fitting step described in Section 4 is performed only after the analytic expression has been obtained: it determines the concrete shape of those correlation functions from the continuous-field simulations and then compares the resulting spectrum against the analytic prediction for validation. The spectrum expression therefore remains independent of the fit; the simulations supply the input functions that could in principle be obtained by other means. To eliminate the tension in the abstract, we have revised the wording to state that the spectrum is derived analytically from the stated assumptions, with the continuous-field simulations used exclusively for validation and to extract the correlation functions. We have also added an explicit paragraph in the derivation section that walks through the steps showing the absence of any calibration parameters. revision: yes

Circularity Check

0 steps flagged

Analytic spectrum from random Higgs configurations is independent; fitting is only for post-derivation validation against simulations.

full rationale

The paper's central derivation exploits the randomness of initial inhomogeneous Higgs field configurations together with statistical isotropy and causality to obtain an essentially analytic magnetic field spectrum, explicitly bypassing lattice simulations for the spectrum itself. Simulations are introduced separately to resolve small-scale structure and to fit general correlation functions (defined from isotropy/causality) for comparison against the already-derived analytic results. Because the analytic form is obtained directly from the randomness assumption without incorporating fitted parameters or simulation outputs, the load-bearing step does not reduce to its own inputs by construction. No self-citations, uniqueness theorems, or ansatzes from prior author work are invoked in a load-bearing way. The fit serves only as an external check and does not alter or presuppose the analytic prediction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Based on the abstract only, the central claim rests on the assumption that random initial configurations permit an analytic treatment and that statistical isotropy plus causality define usable correlation functions; no explicit free parameters or new entities are named.

free parameters (1)

fit parameters for general correlation functions
The abstract states that general correlation functions are fitted to simulation data before comparison with analytic results.

axioms (1)

domain assumption Statistical isotropy and causality constrain the form of the magnetic field spectrum
The abstract explicitly says the authors revisit these effects to define the correlation functions used in the comparison.

pith-pipeline@v0.9.0 · 5398 in / 1244 out tokens · 57092 ms · 2026-05-08T02:44:48.336548+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

30 extracted references · 17 canonical work pages

[1]

admissible

Thus at leading order the power spectrum is isotropic. The power spec- trum is: ⟨|b(k)|2⟩= N 2 δx2(kL)2 h 1 +O (kL)2 i .(25) Then from (15) the energy spectrum for small wave- numbers becomes 8π2 EM(k) k4 =δx 1 +O (kL)2 .(26) This finding is consistent with the causality argument of ref. [23], and completes our proof of the initial mag- netic energy spect...
[2]

Brandenburg and K

A. Brandenburg and K. Subramanian, Astrophysical magnetic fields and nonlinear dynamo theory, Phys. Rept.417, 1 (2005), arXiv:astro-ph/0405052

work page arXiv 2005
[3]

Schober, D

J. Schober, D. Schleicher, C. Federrath, S. Glover, R. S. Klessen, and R. Banerjee, The small-scale dynamo and non-ideal magnetohydrodynamics in primordial star for- mation, The Astrophysical Journal754, 99 (2012)

2012
[4]

Schober, D.R.G

J. Schober, D. R. G. Schleicher, and R. S. Klessen, Mag- netic Field Amplification in Young Galaxies, Astron. Astrophys.560, A87 (2013), arXiv:1310.0853 [astro- 13 ph.GA]

work page arXiv 2013
[5]

Evidence for strong extragalactic magnetic fields from Fermi observations of TeV blazars

A. Neronov and I. Vovk, Evidence for strong extra- galactic magnetic fields from Fermi observations of TeV blazars, Science328, 73 (2010), arXiv:1006.3504 [astro- ph.HE]

work page Pith review arXiv 2010
[6]

Unger and G

M. Unger and G. R. Farrar, The Coherent Magnetic Field of the Milky Way, Astrophys. J.970, 95 (2024), arXiv:2311.12120 [astro-ph.GA]

work page arXiv 2024
[7]

Vachaspati, Magnetic fields from cosmological phase transitions, Phys

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

1991
[8]

The Standard Model cross-over on the lattice

M. D’Onofrio and K. Rummukainen, Standard model cross-over on the lattice, Phys. Rev. D93, 025003 (2016), arXiv:1508.07161 [hep-ph]

work page Pith review arXiv 2016
[9]

Vachaspati, Progress on cosmological magnetic fields, Rept

T. Vachaspati, Progress on cosmological magnetic fields, Rept. Prog. Phys.84, 074901 (2021), arXiv:2010.10525 [astro-ph.CO]

work page arXiv 2021
[10]

Banerjee and K

R. Banerjee and K. Jedamzik, The Evolution of cosmic magnetic fields: From the very early universe, to recom- bination, to the present, Phys. Rev. D70, 123003 (2004), arXiv:astro-ph/0410032

work page arXiv 2004
[11]

Jedamzik and G

K. Jedamzik and G. Sigl, The Evolution of the Large- Scale Tail of Primordial Magnetic Fields, Phys. Rev. D 83, 103005 (2011), arXiv:1012.4794 [astro-ph.CO]

work page arXiv 2011
[12]

Vachaspati and A

T. Vachaspati and A. Brandenburg, Spectra of magnetic fields from electroweak symmetry breaking, Phys. Rev. D111, 043541 (2025), arXiv:2412.00641 [astro-ph.CO]

work page arXiv 2025
[13]

T. W. B. Kibble, Topology of Cosmic Domains and Strings, J. Phys. A9, 1387 (1976)

1976
[14]

Patel and T

T. Patel and T. Vachaspati, Kibble mechanism for elec- troweak magnetic monopoles and magnetic fields, JHEP 01, 059, arXiv:2108.05357 [hep-ph]

work page arXiv
[15]

Nambu, String-Like Configurations in the Weinberg- Salam Theory, Nucl

Y. Nambu, String-Like Configurations in the Weinberg- Salam Theory, Nucl. Phys. B130, 505 (1977)

1977
[16]

Patel and T

T. Patel and T. Vachaspati, Structure of elec- troweak dumbbells, Phys. Rev. D107, 093010 (2023), arXiv:2302.04886 [hep-ph]

work page arXiv 2023
[17]

Patel and T

T. Patel and T. Vachaspati, Annihilation of electroweak dumbbells, JHEP02, 164, arXiv:2311.00026 [hep-ph]

work page arXiv
[18]

E. J. Weinberg,Classical solutions in quantum field the- ory: Solitons and Instantons in High Energy Physics, Cambridge Monographs on Mathematical Physics (Cam- bridge University Press, 2012)

2012
[19]

Vachaspati and A

T. Vachaspati and A. Achucarro, Semilocal cosmic strings, Phys. Rev. D44, 3067 (1991)

1991
[20]

Hindmarsh, Existence and stability of semilocal strings, Phys

M. Hindmarsh, Existence and stability of semilocal strings, Phys. Rev. Lett.68, 1263 (1992)

1992
[21]

G. W. Gibbons, M. E. Ortiz, F. Ruiz Ruiz, and T. M. Samols, Semilocal strings and monopoles, Nucl. Phys. B 385, 127 (1992), arXiv:hep-th/9203023

work page arXiv 1992
[22]

’t Hooft, Magnetic Monopoles in Unified Gauge The- ories, Nucl

G. ’t Hooft, Magnetic Monopoles in Unified Gauge The- ories, Nucl. Phys. B79, 276 (1974)

1974
[23]

Corrigan, D

E. Corrigan, D. I. Olive, D. B. Fairlie, and J. Nuyts, Mag- netic Monopoles in SU(3) Gauge Theories, Nucl. Phys. B106, 475 (1976)

1976
[24]

Durrer and C

R. Durrer and C. Caprini, Primordial magnetic fields and causality, JCAP11, 010, arXiv:astro-ph/0305059

work page arXiv
[25]

Depending on the specific particle physics model at hand,δxmay take any value up to a maximum at the Hubble radius, which sets a causality bound

Alternatively, one can keepδxas a free parameter and think of it as the correlation length of the BEH field after the EWPT. Depending on the specific particle physics model at hand,δxmay take any value up to a maximum at the Hubble radius, which sets a causality bound
[26]

H. S. Abdul-Rahman, M. A. Gdeisat, D. R. Burton, M. J. Lalor, F. Lilley, and C. J. Moore, Fast and ro- bust three-dimensional best path phase unwrapping al- gorithm, Appl. Opt.46, 6623 (2007)

2007
[27]

Wahba,Spline Models for Observational Data(So- ciety for Industrial and Applied Mathematics, 1990) https://epubs.siam.org/doi/pdf/10.1137/1.9781611970128

G. Wahba,Spline Models for Observational Data(So- ciety for Industrial and Applied Mathematics, 1990) https://epubs.siam.org/doi/pdf/10.1137/1.9781611970128

work page doi:10.1137/1.9781611970128 1990
[28]

The remaining steps are es- sentially trivial

The main point of the derivation is to use Kummer’s transformation for the hypergeometric function which factors out an exponential from the sum and terminates the originally infinite series. The remaining steps are es- sentially trivial
[29]

Saveliev, K

A. Saveliev, K. Jedamzik, and G. Sigl, Time Evolution of the Large-Scale Tail of Non-Helical Primordial Magnetic Fields with Back-Reaction of the Turbulent Medium, Phys. Rev. D86, 103010 (2012), arXiv:1208.0444 [astro- ph.CO]

work page arXiv 2012
[30]

Saveliev, K

A. Saveliev, K. Jedamzik, and G. Sigl, Evolution of He- lical Cosmic Magnetic Fields as Predicted by Magneto- hydrodynamic Closure Theory, Phys. Rev. D87, 123001 (2013), arXiv:1304.3621 [astro-ph.CO]

work page arXiv 2013