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.