By using the first-principles approach, we derive a system of three quantum kinetic equations governing the production and evolution of charged scalar particles by electric field in expanding Universe. Analyzing the ultraviolet asymptotic behavior of the kinetic functions, we found the divergent parts of the electric current and the energy-momentum tensor of produced particles and determined the corresponding counterterms. The renormalized system of equations is used to study the generation of electromagnetic fields during and after inflation in the kinetic coupling model $\mathcal{L}_{\rm EM}=-(1/4)f^{2}(\phi)F_{\mu\nu}F^{\mu\nu}$ with the Ratra coupling function $f=\exp(\beta\phi/M_{p})$. It is found that the electric current of created particles is retarded with respect to the electric field. This leads to an oscillatory behavior of both quantities in agreement with the results obtained previously in phenomenological kinetic and hydrodynamical approaches.

We study magnetogenesis in axionlike inflation driven by a pseudoscalar field $\phi$ coupled axially to the electromagnetic (EM) field $(\beta/M_{p})\phi F_{\mu\nu}\tilde{F}^{\mu\nu}$ with dimensionless coupling constant $\beta$. A set of equations for the inflaton field, scale factor, and expectation values of quadratic functions of the EM field is derived. These equations take into account the Schwinger effect and the backreaction of generated EM fields on the Universe expansion. It is found that the backreaction becomes important when the EM energy density reaches the value $\rho_{\rm EM}\sim (\sqrt{2\epsilon}/\beta)\rho_{\rm inf}$ ($\epsilon$ is the slow-roll parameter and $\rho_{\rm inf}$ is the energy density of the inflaton) slowing down the inflaton rolling and terminating magnetogenesis. The Schwinger effect becomes relevant when the electric energy density exceeds the value $\rho_{E}\sim \alpha_{\rm EM}^{-3} (\rho_{\rm tot}^{2}/M_{p}^{4})$, where $\rho_{\rm tot}=3H^{2}M_{p}^{2}$ is the total energy density and $\alpha_{\rm EM}$ is the EM coupling constant. For large $\beta$, produced charged particles could constitute a significant part of the Universe energy density even before the preheating stage. Numerically studying magnetogenesis in the $\alpha$-attractor model of inflation, we find that it is possible to generate helical magnetic fields with the maximal strength $10^{-15}\,{\rm G}$, however, only with the correlation length of order $1\,{\rm pc}$ at present.

We study the generation of electromagnetic fields during inflation when the conformal invariance of Maxwell's action is broken by the kinetic coupling $f^{2}(\phi)F_{\mu\nu}F^{\mu\nu}$ of the electromagnetic field to the inflaton field $\phi$. We consider the case where the coupling function $f(\phi)$ decreases in time during inflation and, as a result, the electric component of the energy density dominates over the magnetic one. The system of equations which governs the joint evolution of the scale factor, inflaton field, and electric energy density is derived. The backreaction occurs when the electric energy density becomes as large as the product of the slow-roll parameter $\epsilon$ and inflaton energy density, $\rho_{E}\sim \epsilon \rho_{\rm inf}$. It affects the inflaton field evolution and leads to the scale-invariant electric power spectrum and the magnetic one which is blue with the spectral index $n_{B}=2$ for any decreasing coupling function. This gives an upper limit on the present-day value of observed magnetic fields below $10^{-22}\,{\rm G}$. It is worth emphasizing that since the effective electric charge of particles $e_{\rm eff}=e/f$ is suppressed by the coupling function, the Schwinger effect becomes important only at the late stages of inflation when the inflaton field is close to the minimum of its potential. The Schwinger effect abruptly decreases the value of the electric field, helping to finish the inflation stage and enter the stage of preheating. It effectively produces the charged particles, implementing the Schwinger reheating scenario even before the fast oscillations of the inflaton. The numerical analysis is carried out in the Starobinsky model of inflation for the powerlike $f\propto a^{\alpha}$ and Ratra-type $f=\exp(\beta\phi/M_{p})$ coupling functions.

By assuming the kinetic coupling $f^{2}(\phi)FF$ of the effective inflaton field $\phi$ with the electromagnetic field, we explore magnetogenesis during the inflation and preheating stages in the Starobinsky model. We consider the case of the exponential coupling function $f(\phi)=\exp(\alpha\phi/M_{p})$ and show that for $\alpha \sim 12-15$ it is possible to generate the large scale magnetic fields with strength $\gtrsim 10^{-15}$ Gauss at the present epoch. The spectrum of generated magnetic fields is blue with the spectral index $n=1+s$, $s>0$. We have found that for the relevant values of the coupling parameter, $\alpha=12-15$, model avoids the back-reaction problem for all relevant modes.

By making use of a simple model that captures the key features of the anomalous Maxwell equations, we study the role of inhomogeneities on the evolution of magnetic fields in a chiral plasma. We find that inhomogeneities of the chiral asymmetry by themselves do not prevent the anomaly-driven inverse cascade and, as in the homogeneous case, the magnetic helicity is transferred from shorter to longer wavelength helical modes of the magnetic field. However, we also find that the evolution appears to be sensitive to the effects of diffusion. In the case when diffusion is negligible, the inverse cascade slows down considerably compared to the homogeneous scenario. In the case of the primordial plasma, though, we find that the diffusion is substantial and efficiently suppresses chiral asymmetry inhomogeneities. As a result, the inverse cascade proceeds practically in the same way as in the chirally homogeneous model.

Particle production in the background of an external classical oscillating field is a key process describing the stage of preheating after inflation. For sufficiently strong couplings between the inflaton and matter fields, this process is known to proceed non-perturbatively. Parametric resonance plays crucial role for bosonic fields in this case, and the evolution of the occupation numbers for fermions is non-perturbative as well. In the Minkowski space, parametric resonance for bosons and non-perturbative effects for fermions would still persist even in the case of weak coupling. In particular, the energy density of created bosons would grow exponentially with time. However, the situation is quite different in the expanding universe. We give a simple demonstration how the conditions of the expanding universe, specifically, redshift of the field modes, lead to the usual perturbative expressions for particle production by an oscillating inflaton in the case of weak couplings. The results that we obtain are relevant and fully applicable to the Starobinsky model of inflation.