We study the onset of the reheating epoch at the end of axion-driven inflation where the axion is coupled to an Abelian, $U(1)$, gauge field via a Chern-Simons interaction term. We focus primarily on $m^2\phi^2$ inflation and explore the possibility that preheating can occur for a range of coupling values consistent with recent observations and bounds on the overproduction of primordial black holes. We find that for a wide range of parameters preheating is efficient. In certain cases the inflaton is seen to transfer all its energy to the gauge fields within a few oscillations. In most cases, we find that the gauge fields on sub-horizon scales end preheating in an unpolarized state due to the existence of strong rescattering between the inflaton and gauge-field modes. We also present a preliminary study of an axion monodromy model coupled to $U(1)$ gauge fields, seeing a similarly efficient preheating behavior as well as indications that the coupling strength has an effect on the creation of oscillons.

In this work we compute the production of magnetic fields in models of axion inflation coupled to the hypercharge sector of the Standard Model through a Chern-Simons interaction term. We make the simplest choice of a quadratic inflationary potential and use lattice simulations to calculate the magnetic field strength, helicity and correlation length at the end of inflation. For small values of the axion-gauge field coupling strength the results agree with no-backreaction calculations and estimates found in the literature. For larger couplings the helicity of the magnetic field differs from the no-backreaction estimate and depends strongly on the comoving wavenumber. We estimate the post-inflationary evolution of the magnetic field based on known results for the evolution of helical and non-helical magnetic fields. The magnetic fields produced by axion inflation with large couplings to $U(1)_Y$ can reach $B_{\rm eff} \gtrsim 10^{-16}\, {\rm G}$, exhibiting a field strength $B_{\rm phys} \approx 10^{-13}\, {\rm G}$ and a correlation length $\lambda_{\rm phys}\approx10\, {\rm pc}$. This result is insensitive to the exact value of the coupling, as long as the coupling is large enough to allow for instantaneous preheating. Depending on the assumptions for the physical processes that determine blazar properties, these fields can be found consistent with blazar observations based on the value of $B_{\rm eff}$. Finally, the intensity of the magnetic field for large coupling can be enough to satisfy the requirements for a recently proposed baryogenesis mechanism, which utilizes the chiral anomaly of the Standard Model.