the transport of mass and angular momentum in accretion disks remains one of the least understood processes in modern astrophysics .
the standard accretion disk model @xcite is based on the assumption that turbulence provides an efficient mechanism for enabling accretion but magnetic fields , thought to be crucial for driving the turbulence , do not play an explicit role .
it is currently believed that the magnetorotational instability ( mri ; @xcite ) is responsible for driving the magnetohydrodynamic ( mhd ) turbulence required for efficient angular momentum transport in astrophysical disks .    since the appreciation of the relevance of the mri to accretion disks ,
significant progress has been made in understanding the physics of the instability and in characterizing the ensuing turbulent state .
a large set of numerical simulations ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) have provided insight into the turbulent mhd flows under various physical conditions .
however , the fundamental processes that determine the strength of the turbulence in the non - linear regime are yet to be deciphered . at present
there have only been a handful of studies addressing the theoretical aspects of this problem ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
thus , understanding the mechanisms that lead to the saturation of the mri constitutes a fundamental problem in modern accretion physics .
building models that capture the essential physics of the mri and its saturation is crucial in order to build angular momentum transport models , and eventually global disk models , beyond the standard viscous accretion disk ( see , e.g. , * ? ? ?
* ; * ? ? ?
the identification , an eventual understanding , of correlations and scaling laws born out of the synergy between numerical simulations , analytical , and semi - analytical work @xcite is likely to provide important insight toward this goal .
this could also prove to be a fruitful approach toward building sub - grid models for related microphysical processes which are unfeasible to incorporate in global simulations in a direct way ( see , e.g. , * ? ? ?
in a recent paper , @xcite provided a summary of a parametric study of incompressible , mri - parasitic instabilities in a viscous , resistive plasma in dissipative regimes accessible to current numerical simulations .
they stated that the fastest growing modes are related to kelvin - helmholtz and tearing mode instabilities . in this work
, we provide the details of the solutions to the differential equations involved and solve for the fastest growing modes in order to provide support for these assertions .
we identify the existence of a critical elsasser number of order unity and show that kelvin - helmholtz and tearing modes dominate in the quasi - ideal and resistive mhd regimes , respectively . by means of a systematic study of the parameter space involved
, we reveal scaling laws that describe the behavior of the fastest growing parasites in these limits .
the joint analysis of the asymptotic behavior of the mri and the parasitic modes that feed off it provides insight into the characteristics of a viable saturation mechanism in regimes that are relevant to astrophysical ( see , e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) as well as experimental environments ( see , e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
consider a cylindrical background characterized by an angular velocity profile @xmath6 , threaded by a vertical magnetic field @xmath7 .
we work in the incompressible limit which is relevant when the magnetic fields involved are so weak that the saturation of the mri occurs at sub - equipartition fields and consider the effects of constant kinematic viscosity @xmath8 and resistivity @xmath9 .
the equations governing the dynamics of this mhd fluid in the shearing box approximation are given by @xmath10 where @xmath11 is the pressure , @xmath12 is the ( constant ) density , and the factor @xmath13 parameterizes the magnitude of the local shear at the fiducial radius @xmath14 .
the continuity equation reduces to @xmath15 and there is no need for an equation of state since the pressure can be determined from this condition .
the set of equations  ( [ eq : euler ] ) and ( [ eq : induction ] ) possesses exact solutions of the form @xmath16 against perturbations parallel to the background magnetic field . for a given wavenumber @xmath17 perpendicular to the disk mid - plane , the growth rate @xmath18 satisfies the dispersion relation @xmath19 where @xmath20 , @xmath21 , @xmath22 is the epicyclic frequency , @xmath23 is the alfvn speed , and @xmath24 is the local keplerian frequency .
the ratio of the amplitudes @xmath25 of the vectors characterizing the mri velocity and magnetic fields , @xmath26 as well as their directions @xmath27 and @xmath28 , see figures  [ fig : mri_representation_3d ] and [ fig : mri_representation_2d ] , are known functions of @xmath29 , see @xcite for a detailed discussion .
the growth rate @xmath18 has a unique maximum , @xmath30 , at @xmath31 .    unless otherwise mentioned , we work with dimensionless variables defined in terms of the characteristic length and time scales set by the background alfvn speed and the local angular frequency : @xmath32 and @xmath33 .
we subsume the scales related to viscosity and resistivity into the dimensionless quantities @xmath34 whose ratio is the magnetic prandtl number , @xmath35 . the quantity @xmath0 is known as the elsasser number , while its viscous counterpart @xmath36 is related to the reynolds number . in our dimensionless variables , magnetic field strengths are defined relative to the background field @xmath37 , while @xmath38 and @xmath39 .
the exact equations for the dynamical evolution of the velocity and magnetic fields corresponding to secondary instabilities @xmath40 and @xmath41 affecting an exact mri mode are obtained by substituting in equations  ( [ eq : euler ] ) and ( [ eq : induction ] ) the ansatz @xmath42 this procedure leads to a set of partial differential equations in space and time whose solution is beyond the scope of the present paper . in order to make the problem of the stability of the non - ideal mri modes against parasitic instabilities more tractable we adopt the same set of assumptions and approximations stated in @xcite , see also @xcite , which are similar to those adopted by @xcite for ideal magnetohydrodynamics ( mhd ) . in brief , we assume that the amplitudes of the primary modes are large enough that we can neglect the influence of the weak vertical background field , the coriolis force , and the background shear flow on the dynamics of the secondary modes . the source of free energy for the secondary instabilities should increase with the amplitudes of the primary modes .
therefore , the secondary growth rates should eventually outrun the primary mri modes .
we thus assume that we can neglect the temporal variation of the sinusoidal velocity and magnetic fields of the mri modes .      under the assumptions stated above
, the equations governing the dynamics of the secondary instabilities are @xmath43 where @xmath44 , @xmath45 stands for the pressure perturbation , and we have defined the time - independent amplitude of an unstable mri mode with wavenumber @xmath17 as @xmath46    we seek solutions to equations  ( [ eq : euler_pi ] ) and ( [ eq : induction_pi ] ) of the form @xmath47 } \ , , \\
\label{eq : b_sol } \delta { \mbox{\boldmath{$b$}}}({\mbox{\boldmath{$x$}}},t ) & = & \delta{{\mbox{\boldmath{$b$}}}}_0(z ) \ , \exp{[s t - i{\mbox{\boldmath{$k$}}}{{\mbox{\boldmath{$\cdot$}}}}{\mbox{\boldmath{$x$ } } } ] } \,,\end{aligned}\ ] ] here , @xmath48 , where @xmath49 is a horizontal wavevector , see figures  [ fig : mri_representation_3d ] and [ fig : mri_representation_2d ] , and @xmath50 is a parameter with @xmath51 .
the eigenvalue @xmath52 determines the nature of the temporal evolution of a given parasitic mode and it must be solved for together with the amplitudes @xmath53 and @xmath54 , which are @xmath55-periodic functions in @xmath56 .
the set of six first order differential equations  ( [ eq : euler_pi ] ) and ( [ eq : induction_pi ] ) can in principle be solved for the secondary velocity and magnetic fields ( @xmath45 can be eliminated using @xmath57 ) by requiring that the `` boundary conditions '' @xmath58 be satisfied for all @xmath56 .
we follow an alternative procedure that leads to a set of higher order differential equations for @xmath59 and @xmath60 .     and
@xmath61 , respectively . for a given primary mri mode with wavenumber @xmath17 , the orientation of the planes containing @xmath62 and @xmath63 , as well as the ratio of the amplitudes @xmath25 , are functions of the viscosity and resistivity .
the wavevector @xmath64 characterizes the horizontal wavelength of a given parasitic mode and @xmath65 denotes the angle between this vector and the radial direction , see also fig .
[ fig : mri_representation_2d ] . ]    substituting expressions ( [ eq : v_sol ] ) and ( [ eq : b_sol ] ) into equations  ( [ eq : euler_pi ] ) and ( [ eq : induction_pi ] ) we obtain for the @xmath56-component of the momentum and induction equations @xmath66 @xmath67 the horizontal component of the momentum equation reads @xmath68 where @xmath69 and @xmath70 . taking the scalar product of this equation with @xmath71 , and using the divergence - less nature of the secondary velocity and magnetic fields to write @xmath72 and @xmath73
, we obtain @xmath74 we can eliminate the pressure perturbations @xmath45 by taking the derivative of equation  ( [ eq : euler_h ] ) and substituting for @xmath75 from equation  ( [ eq : euler_z ] ) .
the result reads @xmath76 ( \partial_z^2 - { \mbox{\boldmath{$k$}}}_{\rm h}^2)\delta v_z \nonumber \\ & + & i ( { \mbox{\boldmath{$k$}}}_{\rm h}{{\mbox{\boldmath{$\cdot$}}}}\partial_z^2\delta{{\mbox{\boldmath{$v$}}}})\delta v_z - i ( { \mbox{\boldmath{$k$}}}_{\rm h}{{\mbox{\boldmath{$\cdot$}}}}\partial_z^2\delta{{\mbox{\boldmath{$b$}}}})\delta b_z \nonumber \\ & + & i({\mbox{\boldmath{$k$}}}_{\rm h}{{\mbox{\boldmath{$\cdot$}}}}\delta{{\mbox{\boldmath{$b$}}}})(\partial_z^2   - { \mbox{\boldmath{$k$}}}_{\rm   h}^2)\delta b_z = 0 \,.\end{aligned}\ ] ] using the explicit form of@xmath77 and @xmath78 from equations  ( [ eq : delta_v_mri ] ) and ( [ eq : delta_b_mri ] ) , we obtain a complete set of equations for the @xmath56-components of the secondary amplitudes @xmath59 and @xmath60 @xmath79 ( \partial_z^2 - { \mbox{\boldmath{$k$}}}_{\rm h}^2)\delta v_z \nonumber \\ & - & ik^2 ( { \mbox{\boldmath{$k$}}}_{\rm h}{{\mbox{\boldmath{$\cdot$}}}}\delta{{\mbox{\boldmath{$v$}}}})\delta v_z + ik^2 ( { \mbox{\boldmath{$k$}}}_{\rm h}{{\mbox{\boldmath{$\cdot$}}}}\delta{{\mbox{\boldmath{$b$}}}})\delta
b_z \nonumber \\ & + & i({\mbox{\boldmath{$k$}}}_{\rm h}{{\mbox{\boldmath{$\cdot$}}}}\delta{{\mbox{\boldmath{$b$}}}})(\partial_z^2   - { \mbox{\boldmath{$k$}}}_{\rm   h}^2)\delta b_z = 0\,,\end{aligned}\ ] ] @xmath80\delta b_z + i({\mbox{\boldmath{$k$}}}_{\rm h}{{\mbox{\boldmath{$\cdot$}}}}\delta{{\mbox{\boldmath{$b$}}}})\delta v_z = 0 \ , .
\nonumber \\\end{aligned}\ ] ]     plane .
the vectors @xmath81 and @xmath82 represent the projection of the mri velocity and magnetic fields , with associated direction angles @xmath83 and @xmath84 .
the versors @xmath85 and @xmath86 characterize , respectively , the directions parallel and perpendicular to the horizontal wavevector @xmath64 of a given parasitic mode . ]    finally , defining the differential operator @xmath87 and regrouping terms , we obtain the equations of motion for the parasites presented in @xcite : @xmath88 @xmath89 thus , the set of assumptions stated in the beginning of this section allows us to reduce the problem of analyzing the stability of the exact mri modes against secondary perturbations in terms of a set of linear , ordinary differential equations with periodic coefficients .
the set of equations  ( [ eq : linear_sys_nu_1 ] ) and ( [ eq : linear_sys_eta_1 ] ) can in principle be integrated along the @xmath56-coordinate subject to the boundary conditions ( [ eq : bc_v ] ) and ( [ eq : bc_b ] ) . note , however , that the terms proportional to the viscosity @xmath8 and resistivity @xmath9 constitute singular perturbations in the limit @xmath90 , i.e. , they multiply the highest derivatives in the differential equations .
it is thus convenient to work in fourier space and seek solutions of the form @xmath91    presuming that @xmath92 , the operator @xmath93 is hermitian and positive definite on functions @xmath94 that are sufficiently well behaved for @xmath95 : specifically , @xmath96 and its first derivative are square - integrable .
when @xmath96 represents @xmath59 or @xmath60 , these mathematical conditions correspond to the physical requirement that the energy in the parasitic mode be finite : @xmath97 \ , dz \ , .
\nonumber \\\end{aligned}\ ] ] the series ( [ eq : fourier_v ] ) and ( [ eq : fourier_b ] ) have infinite energy because they do not vanish as @xmath98 ; however , superpositions of such solutions with a small range of @xmath50 can be made to vanish at infinity .    when the operator @xmath93 acts on the fourier series each individual terms incurs a factor @xmath99 thus , the differential equations ( [ eq : linear_sys_nu_1 ] ) and ( [ eq : linear_sys_eta_1 ] ) lead to recursion relations for the fourier coefficients @xmath100 and @xmath101 of the form @xmath102     \nonumber \\    & - &       ( { \mbox{\boldmath{$k$}}}_{\rm h } { { \mbox{\boldmath{$\cdot$}}}}{\mbox{\boldmath{$v$}}}_0 )   [ ( q_{n-1}-k^2 )    \alpha_{n-1 } - ( q_{n+1}-k^2 ) \alpha_{n+1 } ]   \ , ,     \nonumber \\    \\    0 & = & 2(s + \eta q_n)\beta_n     \nonumber \\    & - & i \ ,
( { \mbox{\boldmath{$k$}}}_{\rm h } { { \mbox{\boldmath{$\cdot$}}}}{\mbox{\boldmath{$b$}}}_0 ) ( \alpha_{n-1}+\alpha_{n+1 } )     -       ( { \mbox{\boldmath{$k$}}}_{\rm h } { { \mbox{\boldmath{$\cdot$}}}}{\mbox{\boldmath{$v$}}}_0 ) ( \beta_{n-1}-\beta_{n+1 } ) \ , .    \nonumber \\\end{aligned}\ ] ] it is convenient to use the natural scales provided by the primary mri mode and rescale the variables : @xmath103 , @xmath104 , @xmath105 , @xmath106 .
we thus obtain @xmath107     \nonumber \\    & + &     \frac{({\mbox{\boldmath{$k$}}}_{\rm h } { { \mbox{\boldmath{$\cdot$}}}}\check{{\mbox{\boldmath{$v$}}}}_0)}{2q_n } \frac{v_0}{b_0 }   [ ( q_{n-1}-1 ) \alpha_{n-1 } - ( q_{n+1}-1 ) \alpha_{n+1 } ] \ , ,     \nonumber \\
\label{eq : eigensystem_v}\end{aligned}\ ] ] @xmath108 where the versors @xmath109 and @xmath110 provide the direction of the mri velocity and magnetic fields , for the primary mri mode is not an independent variable ; it can be calculated for each set of values @xmath111 , see eq .  54 in @xcite .
it is convenient to scale the fourier coefficients using @xmath112 , rather than @xmath113 , because the ratio @xmath25 is proportional to @xmath0 for @xmath114 , see eq .
( [ eq : v0_to_b0_non_ideal ] ) . ] .    the system of coupled linear equations  ( [ eq : eigensystem_v ] ) and ( [ eq : eigensystem_b ] ) , with @xmath115 , can be written in matrix form as @xmath116 , where @xmath117 is a band - diagonal , complex , non - hermitian matrix and the components of the eigenvector @xmath118 are defined according to @xmath119 for @xmath120 and @xmath121 for @xmath122 .
boundary conditions such that @xmath123 must be imposed on the eigenvectors in order to ensure that the fourier series ( [ eq : fourier_v ] ) and ( [ eq : fourier_b ] ) converge as @xmath124 .    in practice we set @xmath125 for @xmath126 and diagonalize the @xmath127 matrix @xmath117 for a given set of values @xmath128 .
the core of the algorithm employed to diagonalize the matrix @xmath117 is based on the lapack routine zgeevx , using the option of applying a balancing transformation to improve the conditioning of the eigenvalues and eigenvectors .
we solve for the finite set of @xmath129 eigenvalues @xmath130 and eigenvectors @xmath131 ( if needed ) , taking increasingly large values of @xmath132 until convergence to the desired accuracy is reached ; @xmath133 seems to do a very good job across the parameter space that we explored .
we consider a grid defined by the parameters @xmath134 and @xmath135 and solve the system of equations  ( [ eq : eigensystem_v ] ) and ( [ eq : eigensystem_b ] ) as a function of the variables @xmath136 focusing our attention on the stability of the fastest growing primary mri mode @xmath137 .
for each pair of values @xmath138 , we solve for @xmath139 over a grid of values according to @xmath140 , @xmath141 , and @xmath142 .
we find the fastest growing mode for a fixed value of the primary mri magnetic field and iterate over @xmath112 , in steps of 0.1 , until the fastest secondary instabilities matches the growth rate of the fastest mri mode , i.e. , @xmath143 .
we denominate the corresponding value of @xmath112 by @xmath144 and we refer to it as the `` saturation '' amplitude of the mri magnetic field @xcite to @xmath145 times larger than @xmath144 . note
that in neglecting the temporal dependence of the background , it is assumed that the secondary growth rates are large compared to the primary growth rate ; thus our definition of saturation implies an extrapolation to the regime where this assumption is not strictly satisfied . ] .
the real and imaginary parts corresponding to the eigenvalues with fastest growth rates , @xmath146 , of various secondary modes that feed off the fastest primary mri mode are shown in the upper and lower panels of figure  [ fig : s_max_pi ] , respectively .
the different panels in this figure correspond to @xmath147 with @xmath148 , and thus @xmath149 , from left to right . the various curves in each panel correspond to different values of the parameter @xmath50 ; and @xmath150 in all of the cases shown .
figure  [ fig : s_max_pi ] illustrates a generic result , described in @xcite for a small range of @xmath151 .
the fastest parasitic modes are non - axisymmetric ( @xmath152 ) , have the same vertical periodicity as the primary mode ( @xmath153 ) , and have purely real growth rates .
the fastest modes have horizontal wavevectors , @xmath154 , that are nearly aligned with either the velocity or the magnetic field of the primary mode .
the first type dominate for @xmath1 and correspond to kelvin - helmholtz modes that reach their maximum growth rates for @xmath155 , i.e. , the direction of the mri velocity field , @xmath156 .
the second type dominate for @xmath2 and are related to tearing modes that grow the fastest along the direction @xmath157 , i.e. , the direction of the mri magnetic field , @xmath158 . note that because the angle @xmath159 for @xmath160 @xcite , the directions parallel to @xmath83 and perpendicular to @xmath84 do not differ by more than a few degrees .
table  [ table : mri_pi_mode_data ] provides information concerning the growth rates and geometric structure of both the fastest growing primary mri mode and the associated fastest parasitic mode for a range of elsasser numbers @xmath161 , with @xmath162 .
as shown in figure  [ fig : s_max_pi ] , there is a sharp transition in behavior around a critical elsasser number of order unity .
the low and high elsasser number regime is characterized by asymptotic behaviors that are already evident for @xmath163 .
tearing modes and kelvin - helmholtz modes present the fastest growth at low and high elsasser numbers , respectively .
the physical properties of both types of modes are insensitive , was explored in @xcite . here
, we mainly focus our attention on the regime where @xmath164 , which exhibits asymptotic behaviors that are more likely to be relevant to astrophysical and experimental environments . ] to the value of @xmath36 provided @xmath165 .
the manifest asymptotic behavior inferred from table  [ table : mri_pi_mode_data ] simplifies the characterization of the physical nature of the parasitic modes as a function of the elsasser number .
thus , we first analyze the asymptotic regimes al large and low elsasser number in section  [ sec : asymptotia ] and postpone the explicit analysis of the physical structure of the eigenmodes until section  [ sec : eigenmodes ] .
it will then be evident that it is only necessary to understand the nature of the modes with moderate elsasser number @xmath166 and that the characteristics of the modes for @xmath114 and @xmath4 can be obtained using the scaling relations presented in section  [ sec : asymptotia ] .
cccrrccccc 1.e@xmath16703 & 7.50e@xmath16704 & 8.66e@xmath16704 & 0.0 & 90.1 & 90.0 & 0.48 & 4.9 & 0.0017 & tm + 1.e@xmath16702 & 7.50e@xmath16703 & 8.66e@xmath16703 & 0.0 & 90.6 & 90.0 & 0.48 & 4.9 & 0.0173 & tm + 1.e@xmath16701 & 7.39e@xmath16702 & 8.55e@xmath16702 & 3.5 & 95.6 & 95.0 & 0.48 & 4.9 & 0.1689 &
tm + 1.e@xmath16800 & 4.28e@xmath16701 & 5.16e@xmath16701 & 22.1 & 119.7 & 30.0 & 0.59 & 4.5 & 0.7187 & kh + 1.e@xmath16801 & 6.96e@xmath16701 & 8.88e@xmath16701 & 40.6 & 132.9 & 40.0 & 0.59 & 3.8 & 0.7816 & kh + 1.e@xmath16802 & 7.44e@xmath16701 & 9.59e@xmath16701 & 44.5 & 134.7 & 45.0 & 0.59 & 3.8 & 0.7751 & kh + 1.e@xmath16803 & 7.49e@xmath16701 & 9.67e@xmath16701 & 45.0 & 135.0 & 45.0 & 0.59 & 3.8 & 0.7743 & kh + 1.e@xmath16804 & 7.50e@xmath16701 & 9.67e@xmath16701 & 45.0 & 135.0 & 45.0 & 0.59 & 3.8 & 0.7742 & kh
identifying in a systematic way the fastest growing parasitic instabilities as a function of the dissipation coefficients allows us to investigate their asymptotic behavior in various regimes . this in turn
enable us to uncover simple scaling relations for their growth rates in terms of the amplitude of the mri velocity and magnetic fields . in order to highlight the physics underlying the various scalings we work with variables with physical dimensions in this section .      in order to better appreciate the asymptotic behavior of the parasitic instabilities it is necessary to consider the asymptotic behavior of the primary mri modes .
the quantities that define the structure of these modes , as well as some relationships between various timescales and amplitudes , are of particular importance and we briefly summarize them here .      in the case of a keplerian shear profile the maximum growth achieved by the mri in the ideal mhd regime , i.e. , @xmath169 , is @xmath170 and corresponds to the mode with wavenumber @xmath171 the fastest growth rate is related to the alfvn frequency , @xmath172 , associated with the background magnetic field @xmath37 , via @xmath173 the magnetic field and velocity field amplitudes for the fastest primary mri mode are related via @xcite @xmath174 the angles characterizing the planes which contain the velocity and magnetic fields of the fastest growing mri mode are given by , see figure  [ fig : sinob_cosob ] , @xmath175      in the inviscid , resistive mhd limit , i.e. , @xmath176 and @xmath114 , the fastest mri growth rate @xmath177 corresponds to the mode with wavenumber @xmath178 in this case , the various relevant inverse timescales are related via @xmath179 where @xmath180 is the alfvn frequency and @xmath181 is the inverse of the resistive timescale across a lengthscale of the order of the mri wavelength , i.e. , @xmath182 .
the magnetic field and velocity field amplitudes for the fastest primary mri mode are related via ( see eq .
92 in * ? ? ?
* ) @xmath183 and thus @xmath184 in the limit @xmath114 .
the planes containing the fastest growing mri velocity and magnetic fields are characterized by the angles , see figure  [ fig : sinob_cosob ] , @xmath185      the dependence of the ( dimensionless ) growth rate of the fastest parasitic modes on the amplitude of the primary mri magnetic field , @xmath112 , is shown in figure [ fig : s_max_pi_scaling ] .
the various sets of data points corresponds to all the pairs of values @xmath138 that can be formed with @xmath186 and @xmath187 .
the various colors ( and sizes ) correspond to different values of @xmath36 .
the values of @xmath0 associated with each set of data points are shown on the right .
open circles identify kelvin - helmholtz modes that feed off the mri velocity field at @xmath188 , which exhibit the fastest growing rates for @xmath1 .
the growth rates for the kelvin - helmholtz modes is almost insensitive to dissipation for @xmath189 .
filled circles correspond to tearing modes , that feed off the mri currents at @xmath190 .
these are the fastest growing secondary instabilities for @xmath2 .
the following ( dimensionless ) equations , shown as dashed lines in figure  [ fig : s_max_pi_scaling ] , provide an excellent description of the fastest growth rates associated with kelvin - helmholtz modes @xmath191 while the scaling of the growth rates of the tearing - modes is well described by @xmath192    although these are not shown in figure  [ fig : s_max_pi_scaling ] , the growth rates of the ( sub - dominant ) kelvin - helmholtz modes obeys @xmath193 there does not seem to be simple scaling relations for the ( sub - dominant ) tearing modes in the regime @xmath176 and @xmath194 for @xmath190 .
the critical value of the elsasser number at which the fastest growing tearing modes ( along the direction @xmath190 ) grow faster than the fastest growing kelvin - helmholtz modes ( along the direction @xmath195 ) is very close to ( but slightly less than ) unity , see figure  [ fig : s_max_pi ] , i.e. , @xmath196 the existence of a an order unity critical elsasser number that dictates the linear evolution of the mri has already been appreciated on both numerical @xcite and analytical grounds ( see , e.g. , * ? ? ?
* ; * ? ? ?
here we posit that this critical elsasser number _ also _ distinguishes which type of parasitic instability dominates the subsequent evolution of the mri . in the remainder of this section
we argue , and demonstrate in section  [ sec : eigenmodes ] , that there are two regimes such that the fastest growing parasitic modes correspond to @xmath197    .
the data shown corresponds to all the pairs of values @xmath138 that can be formed with @xmath198 and @xmath199 .
the various colors correspond to different values of @xmath36 for fixed @xmath200 . for the sake of clarity
, the size of the data points also decreases with @xmath36 for fixed @xmath200 .
the values of @xmath0 associated with each set of data points is shown on the right .
open circles represent kelvin - helmholtz modes , which exhibit the fastest growing rates for @xmath1 .
the growth rates for the kelvin - helmholtz modes is almost insensitive to dissipation for @xmath201 .
filled circles represent tearing modes , which exhibit the fastest growing rates for @xmath2 . in all cases ,
the data points are drawn for @xmath202 .
the dashed lines , proportional to @xmath112 , represent fits to the data whose analytical expressions are provided in section  [ sec : asymptotia ] . ]      for the fastest growing kelvin - helmholtz instabilities with @xmath194 and @xmath176 , the ratio of the horizontal parasitic wavenumber to the wavenumber corresponding to the fastest growing primary mri mode depends very weakly on either @xmath36 or @xmath0 , see table  [ table : mri_pi_mode_data ] , @xmath203 therefore , restoring the physical dimensions into equation  ( [ eq : s_max_kh_ideal ] ) and using equation ( [ eq : v0_to_b0_ideal ] ) leads to @xmath204 where @xmath205 is given by equation  ( [ eq : omega_a_ideal ] ) and @xmath206 note that @xmath207 corresponds to the growth rate associated with a kelvin - helmholtz instability feeding off a velocity field discontinuity of amplitude @xmath113 .
the amplitude of the mri magnetic field needed in order for the growth rate of the parasites to match the growth of the primary mode can be obtained by equating ( [ eq : gamma_max_ideal ] ) and ( [ eq : s_max_kh_ideal ] ) , which leads to @xmath208 in agreement with the results of @xcite .      in the regime @xmath176 and @xmath2
, the fastest growing kelvin - helmholtz instabilities grow slower than the fastest growing tearing modes at the same mri field amplitude .
the ratio of the horizontal parasitic wavenumber to the wavenumber corresponding to the fastest growing primary mri mode depends very weakly on either @xmath36 or @xmath0 , see table  [ table : mri_pi_mode_data ] , @xmath209 therefore , restoring the physical dimensions into equation  ( [ eq : s_max_kh_non_ideal ] ) and using equation ( [ eq : v0_to_b0_non_ideal ] ) leads to @xmath210 where @xmath207 is the kelvin - helmholtz growth rate defined in equation  ( [ eq : kh_pi ] ) and @xmath205 is the alfvn frequency defined in equation  ( [ eq : omega_a_non_ideal ] ) . note that despite the fact that these are kelvin - helmholtz modes , the resistive timescale enters because the resistivity implicitly affects the amplitude of the mri velocity field which is the energy source for these parasitic modes .
, i.e. , @xmath211 . the data shown corresponds to all the pairs of values @xmath138 that can be formed with @xmath212 and @xmath213 .
the various colors correspond to different values of @xmath36 .
there are two clear asymptotic regimes that correspond to @xmath214 and @xmath215 , which are associated with kelvin - helmholtz and tearing modes respectively , see section  [ sec : asymptotia ] . for @xmath216 ,
the viscous quenching of the kelvin - helmholtz modes leads to larger saturation amplitudes for the primary mri modes ; note that these numbers correspond to a highly viscous regime . ]      for the fastest growing tearing modes with @xmath176 and @xmath2 , the ratio of the horizontal parasitic wavenumber to the wavenumber corresponding to the fastest growing primary mri mode depends very weakly on either @xmath36 or @xmath0 , see table  [ table : mri_pi_mode_data ] , @xmath217 therefore , restoring the physical dimensions into equation  ( [ eq : s_max_tm ] ) leads to @xmath218 where we have used equation  ( [ eq : gamma_max_non_ideal_relations ] ) and defined @xmath219 as the alfvn frequency associated with the horizontal mri magnetic field @xmath112 .
the amplitude of the magnetic field that the mri needs to grow in order for the growth rate of the parasites to be as large as that corresponding to the primary mode can be obtained by equating ( [ eq : gamma_max_non_ideal ] ) and ( [ eq : s_max_tm ] ) .
we obtain @xmath220    we can formulate a heuristic argument for the existence of an asymptotic limit for the amplitude @xmath112 in this limit .
the analysis of the classical tearing mode instability , where @xmath221 , with @xmath17 fixed , leads to growth rates @xmath222 with @xmath223 or @xmath224 depending on the details of how the problem is solved in the resistive layer . for @xmath225
, there are only two characteristic scales in the problem , the resistive timescale @xmath226 and the alfvn timescale @xmath227 .
therefore , the growth rate of the parasites must be @xmath228 , with @xmath229 . at a fixed scale @xmath17 ,
the growth rate of the secondaries is proportional to a positive power of the resistivity , @xmath230 .
however , if the scale @xmath17 is given by the fastest growing mri mode then @xmath231 and thus @xmath232 . in this case , because @xmath233 , the amplitude @xmath112 at which @xmath234 must be independent of @xmath9 , i.e. , @xmath235 therefore , @xmath236 const . for @xmath114 .    because the growth rate of the parasitic modes and the mri decay at the same rate for @xmath114 , the amplitude to which the primary mri magnetic field needs to grow to in order for both growth rates to become comparable is independent of the elsasser number
this means that , as resistivity increases the magnetic field that can be generated by the mri before the tearing modes become dynamically important reaches an asymptotic value . in this respect
, one could argue that the ability of the mri to tap energy from the differential rotation does not inherently diminish as resistivity increases .
the instability can achieve this goal via unstable modes with increasingly longer wavelengths that operate at increasingly slower rates .
the situation is quite different for the stresses , however .    , where @xmath84 corresponds to the angle subtended by the magnetic field component of the fastest primary mri mode and the radial direction .
the data shown corresponds to all the pairs of values @xmath138 that can be formed with @xmath212 and @xmath213 . the various symbols ( colors ) correspond to different values of @xmath36 . in the ideal limit , @xmath237 , @xmath238 , while in the inviscid , resistive limit , i.e. , @xmath176 and @xmath114 , @xmath239 .
]      we calculate the dimensionless stress at saturation as @xmath240 , where @xmath241 , is the sum of the reynolds and maxwell stresses @xmath242 integrating these expressions we obtain the dimensionless stress @xmath243 in terms of the parameter and @xmath244 , with @xmath245 , in order to compare results with @xcite .
the definitions adopted here do not affect the expression or the numerical value for the product @xmath246 . ]
@xmath247 , @xmath248 \ , .
\nonumber \\
\label{eq : alpha_beta}\end{aligned}\ ] ] using equations ( [ eq : v0_to_b0_ideal ] ) , ( [ eq : angles_ideal ] ) , and ( [ eq : b0_sat_ideal ] ) for the ideal mhd limit and equations ( [ eq : v0_to_b0_non_ideal ] ) , ( [ eq : angles_non_ideal ] ) , and ( [ eq : b0_sat_non_ideal ] ) for the inviscid , resistive mhd limit , we obtain the asymptotic limit of the dimensionless stress at saturation in equation  [ eq : alpha_beta ] as @xmath249 for @xmath4 and @xmath250 for @xmath114 .
we can relate the initial value of the @xmath251-parameter , @xmath252 , with an estimate for the corresponding value at saturation via @xmath253 . using equation ( [ eq : b0_sat_ideal ] ) and ( [ eq : b0_sat_non_ideal ] )
this leads to @xmath254    it is worth mentioning that numerical simulations of mri driven turbulence ( with no explicit dissipation ) carried out over a wide range of physical conditions lead to saturation values for the parameters @xmath243 and @xmath255 that vary by several orders of magnitude ( see , e.g. * ? ? ?
* and references therein ) .
however , their product remains roughly constant with @xmath256 @xcite .
we emphasize that , in our calculations , it is the behavior of the angle @xmath257 , see figure  [ fig : sinob_cosob ] , and not of the saturation amplitude @xmath258 const .
, see figure  [ fig : b0_sat ] , as the elsasser number decreases that brings about the asymptotic decay of the stress @xmath259 for @xmath260 , see figure  [ fig : alpha_beta_sat ] .
( because @xmath261 , the contribution of the reynolds stress to the total stress is negligible in the very resistive limit . )     if saturation occurs when the fastest parasitic and primary mri growth rates match . the data shown corresponds to all the pairs of values @xmath138 that can be formed with @xmath212 and @xmath135 .
the various symbols ( colors ) correspond to different values of @xmath36 . in the ideal limit , @xmath237 , @xmath262 , while in the inviscid , resistive limit , i.e. , @xmath176 and @xmath263 , @xmath264 . despite the fact that the magnetic field at saturation asymptotes to a constant value , see fig .
[ fig : b0_sat ] , due to the behavior of the function @xmath265 for @xmath114 , see fig .
[ fig : sinob_cosob ] , the dimensionless stress decreases at low elsasser number . ]
throughout this paper we have stated that the fastest parasites are related to kelvin - helmholtz and tearing mode instabilities in the regimes @xmath266 and @xmath2 , respectively .
the reason to delay presenting the rigorous evidence supporting these claims up to this point is based on the existence of the asymptotic regimes and scaling laws presented in section  [ sec : asymptotia ] .
we can now focus on the region of parameter space with `` moderate '' elsasser number , @xmath3 , where changes in @xmath0 produce non - trivial modifications to the parasitic mode structure . with this knowledge , and
the insight gained in section  [ sec : asymptotia ] , it is straightforward to describe the structure of these modes in the limits of large and small elsasser numbers .    as a prelude to the calculations and description of the eigenmodes associated with equations ( [ eq : eigensystem_v ] ) and ( [ eq : eigensystem_b ] ) , we briefly describe some of the key physical signatures that can be used to recognize kelvin - helmholtz ( shear driven ) and tearing mode ( resistive ) instabilities .
the simplest example of a physical system subject to the kelvin - helmholtz instability is that of a uniform , incompressible , shear flow with a discontinuous velocity profile @xmath267 for @xmath268 and @xmath269 for @xmath270 .
this flow can be envisaged as a vortex sheet with an associated vorticity field @xmath271 for @xmath272 and @xmath273 for @xmath274 .
the ensuing kelvin - helmholtz instability emerges as a perturbation in the vorticity field . for a given wavenumber
, the instability can be visualized as an alternating set of perturbations , @xmath275 which tend to increase / decrease the local background vorticity ( see fig .
7.1.2 in * ? ? ?
the kelvin - helmholtz instability reaches the largest growth rates in the absence of dissipation . in the ideal limit
, the instability grows at a rate proportional to the magnitude of the wavenumber parallel to the direction of the velocity field and the magnitude of the velocity discontinuity .
finite viscosity modifies this picture in a thin boundary layer but the unstable modes with longer wavelengths are unaffected .
the tearing - mode instability belongs to the class of spontaneous ( as opposed to driven ) resistive instabilities .
these are instabilities for which non - zero resistivity allows the plasma to access configurations with lower energies by enabling modifications to the field structure ( i.e. , magnetic reconnection ) that would be topologically inaccessible due to ideal constraints .
perhaps the simplest physical system subject to the tearing mode instability is that of a uniform , incompressible , plasma threaded by a magnetic field that reverses direction in a narrow region i.e. , a current sheet , and is uniform far away . the magnetic field given by @xmath276
is usually employed as a convenient model to understand the physics involved in the stability of a thin current sheet @xmath277 of characteristic width @xmath278 .
the detailed solution of this problem depends on the careful analysis of the physical processes that take place in the reconnection region @xcite .
however , asymptotic analytical treatments that match approximate solutions inside and outside the narrow region where the magnetic field reverses provide important insight into the generic features of the solutions ( see , e.g. , * ? ? ?
* ; * ? ? ?
a robust result is that the timescale for the growth of the instability is intermediate between the alfvn timescale and the ( usually much slower ) resistive timescale across the reconnection region .
the tearing mode instability is driven by the lorentz force that acts toward the resonant ( or null ) surface where the magnetic field goes through zero , i.e. , @xmath272 .
non - zero resistivity allows the field lines to slip through the plasma enabling the reconnection of field lines of opposite polarity .
this process takes place in various points along the null surface ( determining the wavelength of the unstable mode ) and gives rise to the formation of magnetic islands .
these are a set of nested magnetic surfaces whose magnetic axis lie along the direction parallel to the background current sheet and are centered in the so - called `` @xmath279-points '' .
the magnetic islands are confined by a separatrix whose intersections with the null surface defines the so - called `` @xmath280-points '' .
finite resistivity allows lines to reconnect at the nodal @xmath280-points and then contract toward the axes of the magnetic islands passing through the @xmath279-points . for a given wavenumber
the instability can thus be visualized as an alternating set of perturbations in the current density , @xmath281 , that tend to increase / decrease the local background current density ( see , e.g. , * ? ? ?
* ; * ? ? ?
although the hallmark of the tearing mode is its current density pattern , the associated velocity field displays a distinct vorticity field pattern ( see , figure 5.7 in * ? ? ?
in order to understand the physical nature of the secondary modes that are likely to play a role in the evolution of the primary mri modes , it is of special interest to understand the fluid dynamics along the directions associated with the fastest growing parasites , @xmath282 . as discussed in section  [ sec : parasitic instabilities ] , for fixed values of the dissipation coefficients , the growth rates of the secondary instabilities peak around specific directions which are almost aligned with either the velocity or magnetic fields of the primary mri mode , i.e. , @xmath155 for @xmath4 and @xmath157 for @xmath114 , see figure  [ fig : s_max_pi ] .
thus , motivated by the characteristics of the kelvin - helmholtz and tearing modes described above , it is of particular interest to understand the structure of the vorticity and current densities associated with the parasitic modes in the directions that are perpendicular to the planes defined by @xmath85 , simply as @xmath85 ] and @xmath284 , when @xmath285 .    in order to best
illustrate the underlying physics of the fastest parasitic modes , it is useful to define a new coordinate system @xmath286 that is rotated with respect to the coordinate system defined by @xmath287 with a rotation angle @xmath288 about the @xmath284 direction .
we define the coordinates @xmath289 and @xmath290 in such a way that they increase along the direction of the versors @xmath85 and @xmath86 , respectively , see figure  [ fig : mri_representation_2d ] , with @xmath291 in this coordinate system , the velocity and magnetic fields of the parasitic modes have components @xmath292 since in these coordinates @xmath293 , the secondary perturbations @xmath59 and @xmath60 in equations ( [ eq : fourier_v ] ) and ( [ eq : fourier_b ] ) become @xmath294        because we are considering an incompressible fluid , we can calculate the projection of the velocity and magnetic field of a given parasitic mode onto the plane defined by @xmath295 , i.e. , @xmath296 just in terms of the fourier coefficients @xmath297 and @xmath298 obtained by solving the eigenvalue problem defined by equations ( [ eq : eigensystem_v ] ) and ( [ eq : eigensystem_b ] ) . using
that @xmath299 and @xmath300 we obtain @xmath301 @xmath302 note that the components of the three - dimensional velocity and magnetic fields that are orthogonal to these planes are independent of the coordinate @xmath290 .
thus , using the fact that the divergence is invariant under rotations , we conclude that the divergence of the two - dimensional vector fields that lie on the plane @xmath303 should vanish , i.e. , @xmath304 .
the components of the vorticity and current density perpendicular to the plane defined by @xmath303 , i.e. , @xmath305 and @xmath306 , are given by @xmath307 where the action of the curl operators in the rightmost equalities is given by @xmath308 therefore , using equations ( [ eq : v_parallel ] ) and ( [ eq : b_parallel ] ) , we can obtain expressions for @xmath309 and @xmath310 in terms of the fourier coefficients @xmath297 and @xmath298 , @xmath311 \alpha_n     e^{-i(n+k_z)z }   e^{-ik_{\rm h } h } \,\,\
, , \nonumber \\
\\    \delta { \mbox{\boldmath{$j$}}}_\perp & = &    \frac{i \check{{\mbox{\boldmath{$k$}}}}_{\rm p}}{k_{\rm h } }    \sum_{n=-\infty}^{\infty } [ k_{\rm h}^2+(n+k_z)^2 ] \beta_n     e^{-i(n+k_z)z }   e^{-ik_{\rm h } h } \,\,\ , .
\nonumber \\ \end{aligned}\ ] ]      the vertical lagrangian displacement @xmath312 provides useful complementary information to help us identify the physical nature of the eigenmodes .
let us then briefly describe how to calculate @xmath312 in terms of the fourier coefficients @xmath297 that result from solving the eigenvalue problem .
the rate of change in the lagrangian displacement @xmath313 with respect to a point moving with velocity @xmath314 is given by @xmath315 in the periodic background provided by the primary mri mode , the lagrangian displacement is of the form @xmath316}$ ] with @xmath317 .
therefore , considering the mri velocity field , @xmath318 , as the background velocity , we can obtain @xmath312 in terms of the vertical eulerian velocity @xmath59 as @xmath319          in section  [ sec : asymptotia ] we showed that the eigenvalues corresponding to the growth rates of the parasitic modes reached well defined asymptotic regimes for @xmath114 and @xmath320 in the limit @xmath176 .
this must also be true for the eigenmodes .
it is then only necessary to explore in detail the behavior of the modes for moderate values of @xmath0 .
we thus focus our attention in the region of parameter space spanned by @xmath321 .
consider as an example if we were to chose @xmath322 or @xmath323 we would obtain the same results for the range of magnetic prandtl numbers @xmath324 or @xmath325 .
] the set of parameters @xmath148 and @xmath326 .
the structure of the corresponding fastest parasitic modes are shown in figure  [ fig : vorticity_current_no_background ] , from left to tight .
the projections of the velocity and magnetic fields onto the planes defined by @xmath303 , i.e. , @xmath327 and @xmath328 , are shown with white arrows in the upper and lower set of panels , respectively .
the color contours show the projection of the vorticity and current density along the direction @xmath86 ( perpendicular to the page ) , i.e. , @xmath329 and @xmath310 .
the versor characterizing the direction of fastest growth , @xmath85 , for the elsasser number @xmath330 points in the direction @xmath331 , see fig .
( [ fig : s_max_pi ] ) and table  [ table : mri_pi_mode_data ] .
this mode feeds off the current density of the primary mri mode .
the corresponding mode structure resembles closely the perturbations in the current density and induced vorticity patterns expected in the analysis of the stability of a set of equidistant current sheets distributed along the @xmath284 direction and alternating sense according to @xmath332 , ( see , figure 5.7 in * ? ? ?
the current density of the secondary modes presents maxima and minima along the planes @xmath333 where the magnetic field of the primary mode , @xmath334 , reverses sign .
thus , the fluctuations induced by these fastest resistive secondary modes tend to promote reconnection of the mri field ; these are the locations where the currents associated with the mri magnetic field vanish .
thus , the presence of these current density perturbations seem to be needed to satisfy the periodic constraints ( on the scale of the unstable mode ) rather than due to the unstable configuration presented by the mri currents themselves . ]
this can be better appreciated by analyzing the vertical lagrangian displacement defined in equation  ( [ eq : xi_z ] ) , which is shown in the leftmost panel of figure  [ fig : lagrangian_z ] .
the observed mode structure is qualitatively insensitive to the value of the elsasser number as long as @xmath2 and @xmath335 , the growth rates and lengthscales associated with each value of @xmath0 change , of course , as discussed in section  [ sec : asymptotia ] .
we thus conclude that the fastest parasitic modes correspond to tearing modes for @xmath2 .
these parasitic modes are enabled by non - zero resistivity and are thus absent in the ideal mhd regime studied by @xcite .
the versors @xmath85 characterizing the direction of fastest growth for the elsasser numbers @xmath336 point in the direction @xmath155 for the three rightmost panels in fig .
( [ fig : s_max_pi ] ) , see also table  [ table : mri_pi_mode_data ] .
these modes feed off the shear in the velocity field of the corresponding primary mri modes .
the velocity and vorticity fields show a periodic structure similar to what is expected from the stability analysis of a periodic set of equidistant vortex sheets distributed along the @xmath284 direction and alternating sense according to @xmath332 .
the mode structure in the velocity and vorticity fields , as well as the growth rates and lengthscales , associated with each value of @xmath0 are quantitatively insensitive to the value of the elsasser number as long as @xmath1 and @xmath337 , as discussed in section  [ sec : asymptotia ] .
the white arrows in the lower panels show the magnetic field projected onto the plane defined by @xmath303 .
the color contours correspond to the associated current densities along the directions @xmath86 .
the structures seen in these panels show some evolution as a function of @xmath194 .
however , these modes are not intrinsically modified from pure kelvin - helmholtz modes since the amplitude of the fluctuations in the magnetic and current density fields is much smaller than the fluctuations in the velocity and vorticity fields .
we thus conclude that the fastest parasitic modes correspond to kelvin - helmholtz modes for @xmath194 . in the limit @xmath4 , these correspond of course to the kelvin - helmholtz modes alluded to in @xcite .
note that the current density of the secondary modes vanishes along , and also in the vicinity of , the planes @xmath338 where the magnetic field of the primary mode , @xmath339 , changes sign . the corresponding vertical lagrangian displacement associated with these modes for @xmath336
are shown in the three rightmost panels in figure  [ fig : lagrangian_z ] , respectively .
thus , the fluctuations induced by the fastest growing kelvin - helmholtz secondary instabilities tend to bend the horizontal mri magnetic field without directly promoting their reconnection .
the upper and lower panels of figure  [ fig : vorticity_current_background ] show the projections of the total ( primary plus secondary ) velocity and magnetic fields onto the planes defined by @xmath303 , i.e. , @xmath340 and @xmath341 , where @xmath342 the four panels correspond to @xmath321 , from left to right .
the color contours show the projection of the total vorticity and current density along the direction @xmath86 ( perpendicular to the page ) , i.e. , @xmath343 and @xmath344 .
the contributions to the vorticity and current density associated with the background velocity and magnetic fields of the primary mri mode are given by @xmath345 in all cases , the amplitude of the primary mode @xmath113 and @xmath112 is such that the fastest secondary modes grow as fast as the primary mri mode , i.e. , @xmath346 .
the addition of the mri fields facilitates the identification of the velocity and magnetic fields that result from the influence of the secondary instabilities with the more familiar structures that are expected from kelvin - helmholtz and tearing mode instabilities in periodic backgrounds .
in particular , it highlights the presence of `` @xmath279 '' and `` @xmath280 '' points in the case of the tearing mode and the wave - like structure of the velocity field in the cases associated with kelvin - helmholtz modes .
note that the contributions to the total velocity and vorticity field from the mri mode is relatively small in the case of @xmath347 ( i.e. , the tearing mode in the leftmost panel ) .
this is because the ratio @xmath25 is small and @xmath348 .
analogously , the contributions to the total magnetic and current density field from the mri mode is relatively small in the case of @xmath349 ( i.e. , the kelvin - helmholtz mode in the rightmost panel ) . in this case ,
we have solved for the spectrum of parasitic instabilities that feed off exact mri modes for a broad range of dissipation coefficients . we parameterized the effects of viscosity and resistivity via the elsasser number number @xmath0 and its viscous counterpart @xmath36 .
our approach to this problem enabled us to explore dissipative regimes that are relevant to astrophysical and laboratory conditions that lie beyond the regime accessible to current numerical simulations .
we explored in a systematic way the physical properties of parasitic modes for a grid of parameters defined by @xmath351 and magnetic prandtl numbers @xmath352 .
this allowed us to identify the existence of a critical elsasser number of order unity and two distinct asymptotic regimes for quasi - ideal , @xmath4 , and inviscid , resistive mhd , @xmath114 .
the dynamics associated with arbitrary combinations of the various parameters involved ( six after choosing the fastest mri mode to set the growth rate and geometry of the primary ) , can be in general rather complex @xcite .
we have thus chosen to focus our attention on the fastest growing parasitic modes .
axisymmetric modes , i.e. , @xmath353 , present growth rates that can be rather small compared to the fastest parasites .
however , they might deserve some attention as they seem to play an important role in the saturation of 2d ( axisymmetric ) simulations @xcite .
we showed that important properties of the fastest secondary instabilities found in @xcite for @xmath354 and @xmath355 are generic .
the fastest parasitic modes are non - axisymmetric , have purely real growth rates , have the same vertical periodicity as the primary mode upon which the feed , and horizontal wavelengths that are roughly twice as long .
their wavevectors @xmath64 are nearly aligned with either the velocity or the magnetic field of the primary mode .
the first type dominate for @xmath1 and correspond to kelvin - helmholtz modes that feed off the vertical shear induced by the mri and reach their maximum growth rates along the direction of the mri velocity field .
the second type dominate for @xmath2 and are related to tearing modes that feed off the mri currents and grow the fastest along the direction of the mri magnetic field .
we provided a set of scaling laws that describe the growth rates of the fastest growing parasitic kelvin - helmholtz and tearing instabilities in two important asymptotic regimes @xmath356    following @xcite , we consider the `` saturation amplitude '' @xmath144 as the amplitude that the magnetic field produced by the mri must have grown to in order for the instantaneous growth rate of the fastest parasitic mode , @xmath357 , to match that of the fastest primary mode , @xmath30 . because both the growth rate of the parasitic modes and the mri decay proportionally to @xmath0 for @xmath114 ,
the amplitude to which the mri magnetic field needs to grow to in order for these growth rates to become comparable is independent of the elsasser number .
as resistivity increases the magnetic field that can be generated by the mri before the tearing modes become dynamically important reaches an asymptotic value .
however , due to the alignment of the mri magnetic field with the azimuthal direction at low elsasser number , the situation is quite different for the stresses .
the stresses produced by the mri primary modes , at the time when the fastest tearing modes have growth rates similar to their own , decay linearly with decreasing elsasser number , i.e. , @xmath358 for @xmath114 .
this behavior seems consistent with numerical simulations of resistive mhd shearing boxes , where the stresses at saturation decrease linearly with elsasser number ( see table  1 in @xcite ; note that the associated elsasser numbers in the more resistive runs are somewhat larger than unity ) .    in the region of parameter space where our analysis overlaps with the regime accessible to numerical simulations , the values of @xmath359 are within factors of a few of the values obtained in the turbulent regime .
as pointed out in @xcite , the dependence on the stresses at saturation on the magnetic prandtl number is less pronounced and the predicted value of @xmath243 is smaller , by a factor of 6 at @xmath360 , than that found in the non - linear simulations of @xcite .
part of these differences can be accounted for with a more sensible operational definition of saturation .
@xcite provide a way to estimate the amplitude of the fields at saturation when primary and secondary modes reach comparable amplitudes based on the values obtained when they reach equal growth rates .
they conclude that the overshoot factor , which depends logarithmically on the initial amplitude of the parasite , is likely between @xmath361 and @xmath145 .
these arguments should also be applicable at low elsasser numbers since they mainly rely on the fact that the growth rate of the secondary modes is linear in the amplitude of the primary mri mode .
it is then conceivable that the asymptotic analysis carried out in this paper could provide , within factors of a few , a reasonable description of the saturation level in a wide region of parameter space .
we explored in detail the physics of the eigenmodes with fastest growth rates in the regime where variations in the elsasser number lead to non - trivial modifications to the structure of the parasitic modes . to this end
we emphasized the importance of understanding the structure of the vorticity and current density patterns associated with the secondary instabilities and used these as a mean to confirm their association with kelvin - helmholtz and tearing mode instabilities .
our analysis suggests a strategy to recognize the presence of parasitic modes and attempt to identify their nature : ( i ) evolve the simulation until the breakdown of the initial exponential growth or subsequent peaks in stress or magnetic energy .
( ii ) dissect the simulation domain in planes perpendicular to the mid - plane that contain the @xmath56-axis .
( iii ) project the velocity and magnetic field in these planes and take the curls in order to obtain the corresponding vorticity and current density .
( iv ) determine whether these resemble what is expected from the kelvin - helmholtz or tearing mode instabilities discussed in section  [ sec : eigenmodes ] .
understanding from first principles the saturation of the mhd turbulence , and the associated angular momentum transport , driven by the mri under realistic physical conditions ( including , ambipolar diffusion , compressibility , stratification , radiation , etc . ) is a challenging endeavor .
the results presented throughout this study provide guidance toward understanding some of the key physical processes that are likely to play a role in halting the exponential growth of the mri and , more speculatively , in the subsequent turbulent state ( @xcite , see also @xcite ) .
for instance , there is numerical evidence suggesting that ohmic heating due to the reconnection of mri field lines is an important source of energy in shearing box simulations of resistive mhd @xcite . in allowing for non - zero resistivity and calculating the currents associated with the secondary instabilities , we have gone one step forward in establishing the chain of processes that enable the conversion of gravitational energy into thermal energy in differentially rotating , magnetized , non - ideal plasmas .    as a final comment , we mention two limitations that result from our working assumptions that are worth emphasizing due to their potential significance ( see also the related discussion in section 2.4 in @xcite ) : ( i ) we have neglected the effects of shear on the dynamics of the secondary modes .
non - axisymmetric parasitic modes will shear linearly in time @xcite ; therefore , assuming fixed horizontal versors to characterize a parasitic mode is an approximation .
( ii ) we have not considered the explicit coupling between the evolution of the mri - modes and the secondary instabilities . as
the secondary modes grow they drain energy from the primary modes ; therefore the secondary growth rates , that rely on the amplitude of the mri modes , could be affected .
both of these effects are likely to diminish the ability of the parasitic modes discussed here to slow the growth of the mri , or the subsequent `` channel flows '' .
thus the extrapolation of the results presented here to the non - linear regime should be complemented with the pertinent quota of skepticism .
having said this , at present , the properties of parasitic modes provide valuable analytical guidance and a basic framework to design and interpret tailored numerical experiments in order to shed light into the non - linear saturation of the mri .