Large-scale dynamo action precedes turbulence in shearing box simulations of the magnetorotational instability

We study the dynamo generation (exponential growth) of large-scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal (x-y) averaging, we also demonstrate the presence of large-scale fields when vertical (y-z) averaging is employed instead. By computing space- time planar averaged fields and power spectra, we find large-scale dynamo action in the early MRI growth phase - a previously unidentified feature. Non-axisymmetric linear MRI modes with low horizontal wavenumbers and vertical wavenumbers near that of expected maximal growth, amplify the large-scale fields exponentially before turbulence and high wavenumber fluctuations arise. Thus the large-scale dynamo requires only linear fluctuations but not nonlinear turbulence (as defined by mode-mode coupling). Vertical averaging also allows for monitoring the evolution of the large-scale vertical field and we find that a feedback from horizontal low wavenumber MRI modes provides a clue as to why the large-scale vertical field sustains against turbulent diffusion in the non-linear saturation regime. We compute the terms in the mean field equations to identify the individual contributions to large-scale field growth for both types of averaging. The large-scale fields obtained from vertical averaging are found to compare well with global simulations and quasi-linear analytical analysis from a previous study by Ebrahimi & Blackman. We discuss the potential implications of these new results for understanding the large-scale MRI dynamo saturation and turbulence.