Completing the gravitational-wave prediction for inflationary gauge-field sourcing: post-inflationary magnetic fields, reheating transitions, tensor non-Gaussianity, and renormalization

Teuscher, Durrer, Martineau and Barrau (JCAP 03 (2026) 043) computed the gravitational-wave background sourced during slow-roll inflation by an Abelian gauge field coupled to the inflaton through kinetic and axial couplings, obtaining a nearly scale-invariant, strongly polarized, non-Gaussian tensor spectrum that can exceed the standard vacuum background. The computation is deliberately conservative: the source is assumed to shut off at the end of slow-roll (instantaneous reheating), so post-inflationary tensor sourcing by the surviving helical magnetic field — including its inverse-cascade evolution — is neglected; era transitions are instantaneous; vacuum contributions are excised by a time-dependent ultraviolet cutoff rather than an explicit renormalization; the analysis is restricted to $0\le\gamma_1<4$ with analytic super-horizon asymptotics that degrade near $\gamma_1\simeq0$ and $\gamma_1\simeq1$; and the tensor bispectrum is not derived (only an expected scaling $\sim\mathcal P_T^{3/2}$). The open problem is to complete this prediction: quantify the post-inflationary magnetic-field contribution to the present-day spectrum, incorporate gradual transitions and a detailed reheating phase, derive the quantitative bispectral signature, and place the result on formal footing via explicit renormalization, controlled asymptotics in the special $\gamma_1$ regimes, and an infrared-cutoff treatment for $\gamma_1\ge4$.

Background

During cosmological inflation — an early phase of accelerated expansion — the metric itself carries transverse-traceless perturbations $h_{ij}$, the gravitational waves (GWs). Their observable content today is summarized by the stochastic GW energy density per logarithmic wavenumber interval, $\Omega_{\rm GW}(k)$, together with the spectral tilt, the degree of circular polarization, and non-Gaussian statistics such as the bispectrum (the connected three-point function, which vanishes for a Gaussian field). In the minimal scenario these tensors are vacuum fluctuations of the gravitational field; in non-standard scenarios additional fields act as sources through their anisotropic stress, the traceless spatial part of the energy-momentum tensor.

A leading motivation for such sources is magnetogenesis: magnetic fields are observed in galaxies, clusters, and filaments, with indirect evidence for fields even in cosmic voids, suggesting a primordial origin. In four space-time dimensions the Maxwell action is conformally invariant, so cosmic expansion alone cannot excite the electromagnetic vacuum; generating magnetic fields during inflation requires non-minimal couplings. Under mild assumptions the most general such couplings are a kinetic coupling $(1+i_1(\phi))F_{\mu\nu}F^{\mu\nu}/4$ (the Ratra-type coupling, which rescales the effective electric charge and is subject to a strong-coupling problem: $1+i_1$ must never become small, as emphasized by Demozzi, Mukhanov and Rubinstein) and an axial coupling $i_2(\phi)F_{\mu\nu}\tilde F^{\mu\nu}/4$, where $\phi$ is the inflaton and $\tilde F$ the dual field tensor. The axial coupling violates parity and acts with opposite signs on the two circular polarizations, exponentially amplifying one of them (growth $\sim e^{\pi\gamma_2/2}$) and producing a helical, parity-odd field. Helical fields are of special interest because an inverse cascade in the later conducting plasma can transfer their magnetic power to large scales.

Teuscher, Durrer, Martineau and Barrau analyzed this system for couplings that are power laws in conformal time $\tau$ (time in comoving units), parametrized by the kinetic exponent $\gamma_1$ and the axial strength $\gamma_2$ over $0\le\gamma_1<4$ (the upper bound avoids infrared divergences of the electric-field spectrum). The gauge-field anisotropic stress sources tensor modes; solving the tensor equation during slow-roll inflation (the regime in which the inflaton rolls slowly enough that the Hubble rate $H$ is nearly constant), with sub-horizon vacuum modes excluded from the source by a time-dependent ultraviolet cutoff $\Lambda(\tau)=k_h(\tau)=\gamma_3\mathcal H$ (an "electromagnetic horizon", $\mathcal H$ being the comoving Hubble rate), yields a nearly scale-invariant spectrum $\mathcal P_T(k)=\left(\frac{H}{M_{\rm Pl}}\right)^4\left(\frac{k}{k_*}\right)^{-4\epsilon}\mathcal F(\gamma_1,\gamma_2),$ where $M_{\rm Pl}$ is the Planck mass, $\epsilon$ the first slow-roll parameter, and $\mathcal F$ grows like $e^{2\pi\gamma_2}$. The tilt $n_T=-4\epsilon$ differs from the vacuum value $-2\epsilon$; the background is strongly circularly polarized and non-Gaussian, since the source is quadratic in a nearly Gaussian gauge field. The authors also provide an analytic matching formalism that propagates the sourced tensors through an arbitrary number of post-inflationary eras of constant equation-of-state parameter $w$ (pressure-to-energy-density ratio), assuming instantaneous transitions and instantaneous reheating, and they chart the parameter region where this secondary background exceeds the standard inflationary one while backreaction (the feedback of the gauge-field energy on the inflationary dynamics) stays negligible.

The prediction is deliberately conservative, and the paper itself enumerates what is missing: after inflation the electric field is rapidly damped by plasma currents and small-scale magnetic power diffuses away, but the late-time anisotropic stress of the surviving magnetic field — including the inverse cascade that moves its helical power to large scales — is not included and could modify the final amplitude; the era matching assumes instantaneous transitions and neglects model-dependent GW production during reheating; the vacuum subtraction is implemented by a physical cutoff rather than an explicit renormalization; the super-horizon analytic asymptotics lose accuracy near $\gamma_1\simeq0$ and $\gamma_1\simeq1$ where logarithmic terms matter; the restriction $\gamma_1<4$ excludes the $\gamma_1\simeq6$ case that would give scale-invariant magnetic fields, whose treatment requires an infrared cutoff set by the beginning of inflation; and the bispectrum of the sourced tensor field is not computed, only an expected order $\mathcal P_T^{3/2}$. Completing this prediction is the open problem.

Problem Statement

The problem is to complete, for the inflationary gauge-field model of Teuscher, Durrer, Martineau and Barrau — an Abelian $U(1)$ field coupled to the inflaton through a kinetic coupling $1+i_1$ and an axial coupling $i_2$, both power laws in conformal time parametrized by $(\gamma_1,\gamma_2)$ with $0\le\gamma_1<4$ and $\gamma_2\in\mathbb{R}$, during slow-roll inflation — the prediction of the present-day stochastic gravitational-wave background, which the source paper computes only for the inflationary era under conservative assumptions. A complete answer requires all four of the following results, which jointly serve this single objective:

(a) Post-inflationary sourcing: compute the gravitational waves sourced after inflation by the magnetic field that survives the transition (the electric field being damped by plasma currents, small-scale magnetic power by diffusion), including the effect of the helical inverse cascade in the conducting primordial plasma, and determine whether and how it modifies the predicted present-day $\Omega_{\rm GW}(k)$ and its polarization.

(b) Era transitions and reheating: extend the propagation beyond instantaneous transitions between constant-equation-of-state eras, either by incorporating gradual transitions or by modeling a detailed reheating phase (including any model-dependent GW production during reheating), and quantify the resulting changes in the mapping to present-day observables.

(c) Non-Gaussianity: derive the quantitative non-Gaussian statistics of the sourced tensor field — in particular its bispectrum, with amplitude, momentum-space shape, and polarization dependence — beyond the order-of-magnitude expectation $\sim\mathcal P_T^{3/2}$, at the level needed for observational predictions.

(d) Formal robustness: replace the time-dependent ultraviolet cutoff with an explicit renormalization of vacuum contributions to the anisotropic stress (which must vanish in the free-Maxwell limit so that the vacuum sources nothing); provide controlled super-horizon asymptotics near $\gamma_1\simeq0$ and $\gamma_1\simeq1$, where logarithmic and higher-order series terms dominate, at minimum by numerical confirmation; and treat the infrared regime $\gamma_1\ge4$ using a cutoff determined by the beginning of inflation, characterizing its imprint on the GW spectrum and the associated sensitivity to the total duration of inflation.

A complete answer should state which conclusions of the conservative treatment — the near scale invariance, the tilt $n_T=-4\epsilon$, the strongly polarized and non-Gaussian character, and the parameter window in $(\gamma_1,\gamma_2)$ and the inflation scale where the sourced background is observable while backreaction stays negligible — survive once these effects are included. The source paper additionally asks whether the results generalize to non-Abelian gauge fields; that generalization is recorded here as a related but independently separable question, not as part of the completion boundary above.

References

1. Martin Teuscher, Ruth Durrer, Killian Martineau, Aurélien Barrau, "Gravitational waves sourced by gauge fields during inflation", JCAP 03 (2026) 043, DOI 10.1088/1475-7516/2026/03/043, arXiv:2510.00869, https://arxiv.org/abs/2510.00869
2. Martin Teuscher, Ruth Durrer, Killian Martineau, Aurélien Barrau, "A universal scaling law for gravitational waves induced during inflation", Phys. Rev. D 113 (2026) 103545, DOI 10.1103/sn8s-d1sg, arXiv:2512.14670, https://arxiv.org/abs/2512.14670
3. Chiara Caprini, Lorenzo Sorbo, "Adding helicity to inflationary magnetogenesis", JCAP 10 (2014) 056, DOI 10.1088/1475-7516/2014/10/056, arXiv:1407.2809, https://arxiv.org/abs/1407.2809
4. Chiara Caprini, Ruth Durrer, Elisa Fenu, "Can the observed large scale magnetic fields be seeded by helical primordial fields?", JCAP 11 (2009) 001, DOI 10.1088/1475-7516/2009/11/001, arXiv:0906.4976, https://arxiv.org/abs/0906.4976
5. Tina Kahniashvili, Axel Brandenburg, Ruth Durrer, Alexander G. Tevzadze, Winston Yin, "Scale-invariant helical magnetic field evolution and the duration of inflation", JCAP 12 (2017) 002, DOI 10.1088/1475-7516/2017/12/002, arXiv:1610.03139, https://arxiv.org/abs/1610.03139
6. Subhasis Maiti, Debaprasad Maity, Rohan Srikanth, "Probing reheating phase via nonhelical magnetogenesis and secondary gravitational waves", Phys. Rev. D 112 (2025) 063552, DOI 10.1103/4n86-9nsc, arXiv:2505.13623, https://arxiv.org/abs/2505.13623
7. Bill Atkins, Debika Chowdhury, Alisha Marriott-Best, Gianmassimo Tasinato, "Inflationary magnetogenesis beyond slow-roll and its induced gravitational waves", Phys. Rev. D 112 (2025) 063534, DOI 10.1103/btjz-pqgv, arXiv:2507.01772, https://arxiv.org/abs/2507.01772
8. Arko Bhaumik, Theodoros Papanikolaou, Anish Ghoshal, "Vector induced gravitational waves sourced by primordial magnetic fields", JCAP 08 (2025) 054, DOI 10.1088/1475-7516/2025/08/054, arXiv:2504.10477, https://arxiv.org/abs/2504.10477
9. Arko Bhaumik, Theodoros Papanikolaou, Anish Ghoshal, "Magnetically assisted primordial scalar perturbations: Scalar-Induced Gravitational Waves", arXiv:2607.17875, https://arxiv.org/abs/2607.17875
10. Alberto Roper Pol, Chiara Caprini, Andrii Neronov, Dmitri Semikoz, "Gravitational wave signal from primordial magnetic fields in the Pulsar Timing Array frequency band", Phys. Rev. D 105 (2022) 123502, DOI 10.1103/PhysRevD.105.123502, arXiv:2201.05630, https://arxiv.org/abs/2201.05630
11. Maresuke Shiraishi, Toyokazu Sekiguchi, "First observational constraints on tensor non-Gaussianity sourced by primordial magnetic fields from cosmic microwave background", Phys. Rev. D 90 (2014) 103002, DOI 10.1103/PhysRevD.90.103002, arXiv:1304.7277, https://arxiv.org/abs/1304.7277
12. Konstantinos Dimopoulos, Anish Ghoshal, Theodoros Papanikolaou, "Is inflationary magnetogenesis sensitive to the postinflationary history?", Phys. Rev. D 111 (2025) 083550, DOI 10.1103/PhysRevD.111.083550, arXiv:2412.10367, https://arxiv.org/abs/2412.10367
13. Neil Barnaby, Ryo Namba, Marco Peloso, "Phenomenology of a Pseudo-Scalar Inflaton: Naturally Large Nongaussianity", JCAP 04 (2011) 009, DOI 10.1088/1475-7516/2011/04/009, arXiv:1102.4333, https://arxiv.org/abs/1102.4333
14. Peter Adshead, Mark Wyman, "Natural Inflation on a Steep Potential with Classical Non-Abelian Gauge Fields", Phys. Rev. Lett. 108 (2012) 261302, DOI 10.1103/PhysRevLett.108.261302, arXiv:1202.2366, https://arxiv.org/abs/1202.2366
15. A. Maleknejad, M. M. Sheikh-Jabbari, J. Soda, "Gauge Fields and Inflation", Phys. Rept. 528 (2013) 161-261, DOI 10.1016/j.physrep.2013.03.003, arXiv:1212.2921, https://arxiv.org/abs/1212.2921
16. Ema Dimastrogiovanni, Matteo Fasiello, Martino Michelotti, Lucas Pinol, "Primordial gravitational waves in non-minimally coupled chromo-natural inflation", JCAP 02 (2024) 039, DOI 10.1088/1475-7516/2024/02/039, arXiv:2303.10718, https://arxiv.org/abs/2303.10718
17. O. O. Sobol, A. V. Lysenko, E. V. Gorbar, S. I. Vilchinskii, "Gradient expansion formalism for magnetogenesis in the kinetic coupling model", Phys. Rev. D 102 (2020) 123512, DOI 10.1103/PhysRevD.102.123512, arXiv:2010.13587, https://arxiv.org/abs/2010.13587
18. Chiara Caprini, Daniel G. Figueroa, "Cosmological Backgrounds of Gravitational Waves", Class. Quant. Grav. 35 (2018) 163001, DOI 10.1088/1361-6382/aac608, arXiv:1801.04268, https://arxiv.org/abs/1801.04268
19. V. Demozzi, V. Mukhanov, H. Rubinstein, "Magnetic fields from inflation?", JCAP 08 (2009) 025, DOI 10.1088/1475-7516/2009/08/025, arXiv:0907.1030, https://arxiv.org/abs/0907.1030