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.