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.

In this paper we investigate large non-gaussianity in axion N-flation models, taking account while dynamically a large number of axions begin away from the hilltop region(come down from the hill) and so serve only to be the source of the Hubble rate. Therefore the single field stays closest to the hilltop sources the non-Gaussianity. In this case most of axions can be replaced by a single effective field with a quadratic potential. So our potential will contain two fields. The full cosine is responsible for the axion closest to hilltop and quadratic term which is a source for Hubble rate [4]. We obtain power spectrum, spectral index and non-gaussianity parameter, then we impose conditions from WMAP for power spectrum and spectral index and see how large on non-gaussianity parameter it is possible to achieve with such conditions. Finally we swap quadratic term to {\lambda}{\phi}^4 and see whether this makes it harder or easier to achieve large non-gaussianity.We find large non-gaussianity is achievable by imposing data from WMAP conditional on axion decay constant f has reasonable value in connection with Planck mass and by requiring number of e-folds must be bounded between 40-60.When we swap to {\lambda}{\phi}^4 we find that it is harder to achieve non-Gaussianity, because we are imposed to investigate only {\lambda}{\phi}^4 domination in consistency with WMAP data for spectral index. Although, in this case we find that large non-gaussianity is still achievable. Finally we verify imposing the condition for spectral index to be nearly one and find acceptable and detectable value for non-gaussianity typically of order of 10 and 100 depending on value of decay constant f.Swapping to {\lambda}{\phi}^4 in this case does not give us any significant relation. In this paper we consider models which can generate large non-gaussianity. We restrict ourselves to axion N-flation models.