We use the $f^{2}FF$ model to study the generation of primordial magnetic fields (PMF) in the context of large field inflation (LFI), described by the potential, $V \sim M \phi^{p}$. We compute the magnetic and electric spectra for all possible values of the model parameters under de Sitter and power law expansion. We show that scale invariant PMF are not obtained in LFI to first order in the slow roll approximation, if we impose the constraint $V(\phi=0)\sim 0$. Alternatively, if these constraints are relaxed, the scale invariant PMF can be generated. The associated electric field energy can fall below the energy density of inflation, $\rho_{\rm{Inf}}$ for the ranges of comoving wavenumbers, $ k > 8 \times 10^{-7} \rm{Mpc^{-1}}$ and $ k > 4 \times 10^{-6} \rm{Mpc^{-1}}$ in de Sitter and power law (PL) expansion. Further, it can drop below $\rho_{\rm{Inf}}$ on the ranges, e-foldings $N > 51$, $p<1.66$, $p >2.03$, $l_0 > 3 \times 10^5 {M_{\rm{Pl}}}^{-1} (H_i < 3.3 \times 10^{-6} M_{\rm{Pl}})$, and $M > 2.8 \times 10^{-3} M_{\rm{Pl}}$. All of the above ranges fit with the observational constraints.

We study the simple gauge invariant model ${f^2}FF$ as a way to generate primordial magnetic fields (PMF) in Natural Inflation (NI). We compute both magnetic and electric spectra generated by the ${f^2}FF$ model in NI for different values of model parameters and find that both de Sitter and power law expansion lead to the same results at sufficiently large number of e-foldings. We also find that the necessary scale invariance property of the PMF cannot be obtained in NI in first order of slow roll limits under the constraint of inflationary potential, $V\left( 0 \right) \simeq 0$. Furthermore, if this constraint is relaxed to achieve scale invariance, then the model suffers from the backreaction problem for the co-moving wave number, $k \lesssim 8.0\times 10^{-7} \rm{Mpc^{-1}}$ and Hubble parameter, $H_i \gtrsim 1.25\times 10^{-3} \rm{M_{\rm{Pl}}}$. The former can be considered as a lower bound of $k$ and the later as an upper bound of $H_i$ for a model which is free from the backreaction problem. Further, we show that there is a narrow range of the height of the potential $\Lambda $ around ${\Lambda _{\min }} \approx 0.00874{M_{{\rm{Pl}}}}$ and of $k$ around ${k_{\min }} \sim 0.0173{\rm{Mp}}{{\rm{c}}^{ - 1}}$, at which the energy of the electric field can fall below the energy of the magnetic field. The range of $k$ lies within some observable scales. However, the relatively short range of $k$ presents a challenge to the viability of this model.