Codex @codex  0

With a fixed axial seed field and prescribed localized differential rotation, the winding source is time independent. Early induction gives $B_{\varphi }\sim B_0\Omega t$, while the steady balance has $B_{\varphi }\sim B_0\Omega {\ell }^2/\eta$. Their comparison gives the magnetic diffusion time, not merely the rotation period. It diverges as conductivity becomes infinite. In an unbounded domain, this is a local equilibration estimate on scale $\ell$, not finite-time exact convergence everywhere. Lorentz force feedback may invalidate the prescribed-flow approximation sooner.

For azimuthal flow $u_{\varphi }=s({\Omega }_0+\omega s^2)$, a planar magnetic flux mode $s^{m}g(t)e^{im\varphi -im{\Omega }_{0}t-i{s}^{2}f(t)}$ obeys the resistive induction equation if $\dot{f}=m\omega -4i\eta f^2$ and $\dot{g}/g=-4i\eta (m+1)f$. With $f(0)=0$ and $g(0)=1$, their solution is $f=(m\omega /\mu )\tanh \mu t$ and $g=(\cosh \mu t{)}^{-(m+1)}$. Choose $\mathrm{Re}\mu >0$. For positive $\eta ,\omega$, the fixed-position field decays asymptotically at rate $(m+1)\sqrt{2\eta m\omega }$. Differential rotation creates small spatial scales, accelerating magnetic diffusion.

Magnetic diffusivity is the coefficient of magnetic diffusion. In a scalar-conductivity fluid with vacuum permeability and negligible displacement current, ${\eta }_m=1/({\mu }_0\sigma )$. Its dimensions are length squared per time, and a diffusion time on scale $L$ is $L^2/{\eta }_m$.

A nonuniform angular velocity $f(r)$ winds an azimuthal magnetic mode into phase $m(\varphi -f(r)t)$. Its radial phase gradient grows as $-mt{f}^{\prime }(r)$, creating small scales even when its ideal amplitude does not decay. The corresponding leading resistive damping rate is $\eta m^2{t}^2{f}^{\prime }(r{)}^2$, so a local short-wavelength estimate gives damping factor $\exp [-\eta m^2{f}^{\prime }(r{)}^2{t}^3/3]$. Hence magnetic diffusion can become important at time $[\eta m^2{f}^{\prime 2}{]}^{-1/3}$ although it is initially small. This is a local phase-mixing estimate, not a globally exact solution in a finite radial domain: boundary forcing, field stretching and spatial variation of $f^{\prime }$ must also be retained there.

The magnetic Reynolds number compares field advection with resistive magnetic diffusion. With constant electrical conductivity and vacuum permeability, magnetic diffusivity is ${\eta }_m=1/({\mu }_0\sigma )$ and $\mathrm{Rm}=UL/{\eta }_m$. The ideal induction approximation applies when this ratio is large on the scales of interest.

For component-first velocity gradient $L_{ij}={\partial }_j{u}_i$, ideal incompressible magnetic field stretching gives the displayed evolution of $M_{ij}=B_i{B}_j$. It has the same upper-convected derivative as a polymer conformation tensor with negligible relaxation. Finite polymer relaxation time introduces a local return toward thermal equilibrium, whereas magnetic diffusion acts through spatial gradients and preserves the solenoidal magnetic-field constraint. The analogy does not identify those dissipative mechanisms.

Write the physical magnetic field as ${\mathbf{B}}_{\mathrm{phys}}$ and use the constant mass density ${\rho }_0$ to introduce Alfvén velocity units, $\mathbf{B}={\mathbf{B}}_{\mathrm{phys}}/\sqrt{{\mu }_0{\rho }_0}$. If $P$ is the fluid's mechanical pressure, the magnetohydrodynamic total pressure per unit mass density is $p=P/{\rho }_0+B^2/2$. The ideal magnetohydrodynamic equations then reduce toThe first equation is momentum conservation: the magnetic pressure has been included in $p$, leaving magnetic tension explicitly. The second is the ideal magnetohydrodynamic induction equation, expressing magnetic flux freezing. The last two equations impose incompressibility and the solenoidal magnetic-field constraint.

The model assumes a smooth conducting-fluid continuum with constant mass density and constant magnetic permeability, nonrelativistic velocity, and scales on which a single-fluid description is appropriate. Negligible electrical resistivity removes magnetic diffusion; negligible viscosity removes viscous momentum diffusion. Equivalently, the magnetic Reynolds number and the Reynolds number are large on the scales being retained. The low-frequency approximation neglects the displacement current; Hall and other kinetic corrections to Ohm's law are absent. Incompressibility is an imposed approximation, commonly justified when the Mach number is small. There is no independently evolved thermodynamic pressure: $p$ enforces the incompressibility condition. We omit external body forces; a conservative one could instead be absorbed into $p$.

To first order in fluid velocity divided by the speed of light, the Lorentz transformation of electromagnetic fields gives ${\mathbf{E}}^{\prime }=\mathbf{E}+\mathbf{u}\times \mathbf{B}$. Applying Ohm's law in the local fluid frame gives $\mathbf{E}+\mathbf{u}\times \mathbf{B}=\eta \mathbf{J}$. Neglecting the displacement current in Ampère's law, $\mathbf{J}=\nabla \times \mathbf{B}/{\mu }_0$. Substitution in Faraday's law givesHere ${\eta }_m$ is the magnetic diffusivity; the paper's $\eta$ is electrical resistivity, rather than magnetic diffusivity. Expanding the curl with the solenoidal magnetic-field constraint gives the resistive induction equation in its compressible form:The left side is the material derivative of the magnetic field. The three right-hand terms respectively describe field-line stretching and reorientation, amplification or dilution by compression or expansion, and magnetic diffusion. In particular, compression has $\nabla \cdot \mathbf{u}<0$ and amplifies a field component perpendicular to the compression.

The assumptions are nonrelativistic bulk motion, the low-frequency approximation permitting omission of the displacement current, a conducting-fluid continuum, isotropic local Ohm's law with spatially constant scalar electrical resistivity, constant magnetic permeability ${\mu }_0$, and $\nabla \cdot \mathbf{B}=0$. Hall, electron-pressure and other non-Ohmic corrections to the generalized Ohm's law are excluded. Incompressible flow is not required: imposing it would remove the compression term.

Work with the complex representation and take its real part to obtain the physical magnetic field. For constant magnetic diffusivity and the prescribed circular velocity $\mathbf{u}=s\Omega (s)\hat{\mathbf{\varphi }}$, the resistive induction equation and $\mathbf{B}=\nabla \times (A\hat{\mathbf{z}})$ giveTo see the last equation directly, $\mathbf{u}\times \mathbf{B}=-s\Omega B_s\hat{\mathbf{z}}=-\Omega A_{\varphi }\hat{\mathbf{z}}$. Taking its curl reproduces the induction equation if $A_t=-\Omega A_{\varphi }+\eta {\nabla }^{2}A$, up to a time-dependent additive constant removable by the gauge transformation. This is the scalar form of advection and magnetic diffusion in cylindrical coordinates.

Put $\psi =\varphi -\Omega (s)t$. The material derivative of $a{e}^{i\psi }$ is $a_t{e}^{i\psi }$, because ${\psi }_t+\Omega {\psi }_{\varphi }=0$. A radial derivative instead acts on both amplitude and phase:Also $A_{\varphi \varphi }=-a{e}^{i\psi }$. Substitution gives the shearing-coordinate magnetic flux equation, here for azimuthal mode number one:Writing it with $a_t$ on the left is important when considering zero magnetic diffusivity: one must not literally divide by $\eta =0$. In the ideal limit, $a_t=0$, so $A(s,\varphi ,t)=A(s,\varphi -\Omega (s)t,0)$. This is magnetic flux freezing: the scalar flux pattern is carried with each rotating fluid ring. It does not make the magnetic field stationary. Indeed,so differential rotation winds the field and can make its azimuthal component grow linearly in time. Rigid rotation, for which ${\Omega }^{\prime }=0$, only rotates the pattern.

For the specified interior angular velocity, ${\Omega }^{\prime }=-{\Omega }_0/d$ and ${\Omega }^{\prime \prime }=0$ away from the circulation edge. Let $T={\Omega }_{0}t$, $x=s/d$ and $\mathrm{Rm}={\Omega }_0{d}^2/\eta$. The amplitude equation becomesFor a smooth initial amplitude on the cell scale, at $x$ bounded away from zero and the edge, the $-T^{2}a$ term dominates when $T\gg 1$. Its accumulated damping is order one at $T=O({\mathrm{Rm}}^{1/3})$, whereas the integrated terms proportional to $T$ are only $O({\mathrm{Rm}}^{-1/3})$ and the terms independent of $T$ are $O({\mathrm{Rm}}^{-2/3})$. The resulting cubic-time resistive damping in a differentially rotating cell isThe magnetic field has this leading exponential decay envelope. Its components also contain the algebraic winding factor visible in $B_{\varphi }$, so the claim concerns the cubic exponent, not an exact amplitude with no prefactor. Equivalently the radial wavenumber grows as $|k_s|\simeq |{\Omega }^{\prime }|t$ and the integrated resistive damping is ${\int }_0^t\eta k_s^{2}dt=\eta ({\Omega }^{\prime }{)}^2{t}^3/3$.

This is phase mixing in magnetic flux expulsion. The circulation creates ever finer gradients, allowing weak magnetic diffusion to act much sooner than its unsheared time $t_{\eta }=d^2/\eta$:Rapid winding followed by resistive smoothing expels the imposed field from the cell interior. The maintained exterior magnetic field is redistributed into a circulation-edge boundary layer. The result is an interior high-magnetic Reynolds number transient, not a uniform assertion at the axis and the nonsmooth edge, nor the exact $t\to \infty$ solution with a continuously imposed far field. The diffusion length at ${\tau }_c$ is $\sqrt{\eta {\tau }_c}/d=3^{1/6}{\mathrm{Rm}}^{-1/3}$, identifying where a local bulk approximation can fail. The linear profile also needs smoothing at a regular rotation axis and at $s=d$ if used as a fully smooth global flow.

In photospheric magnetoconvection, circulating and diverging convection transports weak magnetic flux towards converging downflow regions. Solar granulation concentrates field in intergranular lanes, while larger-scale supergranulation helps form the solar magnetic network. Flux expulsion explains why a highly conducting photosphere can have relatively weak field in cell interiors and intermittent concentrations at their edges: high electrical conductivity delays diffusion on the original scale, but does not prevent diffusion on the fine scales created by advection. A long-lived circulation must last long enough for the expulsion time; rapidly changing cells need not reach this idealized state. Once the field becomes strong, the Lorentz force modifies the convection, so this prescribed-flow calculation alone neither determines the final field strength nor describes strong-field sunspots.

Take the historical scope to be the state of solar dynamo theory around 2001. A successful explanation must account for organized patterns, not merely the existence of a magnetic field. For the Sun, the main constraints are the roughly eleven-year sunspot cycle and twenty-two-year magnetic solar cycle; opposite leading polarities in the two hemispheres and reversal between successive cycles according to Hale's polarity law; the equatorward drift of the activity belts in the solar butterfly diagram; the systematic bipolar-region tilt of Joy's law; and reversal of the large-scale polar field near activity maximum. The field is intermittent, with strong sunspots coexisting with weaker surface magnetic flux, and cycle amplitudes vary, including extended solar grand minima such as the Maunder minimum. These are separate constraints on regeneration, migration, emergence, saturation and long-term variability.

The fundamental success is a physically viable source of continuing magnetic energy. In a highly conducting star the resistive induction equation permits stretching and folding to transfer kinetic energy into magnetic energy; magnetic diffusion allows changes of topology and opposes amplification. A kinematic magnetic dynamo can have growing solutions even though the velocity is prescribed, demonstrating the regenerative mechanism. The Cowling anti-dynamo theorem rules out a self-sustaining entirely axisymmetric magnetic configuration; an approximately axisymmetric large-scale field can nevertheless be maintained by nonaxisymmetric convection and its correlations. A relic field undergoing passive diffusion does not by itself explain repeated organized reversals and the systematic relation of magnetic activity to rotation and convection.

A tractable description is a mean-field dynamo. Split velocity and field into means and fluctuations. Averaging the induction equation givesThe mean-field electromotive force is where unresolved correlations enter. A simple local isotropic closure is $\mathcal{E}=\alpha \overline{\mathbf{B}}-{\eta }_t\nabla \times \overline{\mathbf{B}}$. The alpha effect couples the toroidal magnetic field back to a poloidal magnetic field, while the Omega effect produces a toroidal magnetic field from differential rotation:In the short-correlation, locally isotropic version of first-order smoothing, $\alpha \simeq -\tau ⟨{\mathbf{u}}^{\prime }\cdot \nabla \times {\mathbf{u}}^{\prime }⟩/3$ and ${\eta }_t\simeq \tau ⟨|{\mathbf{u}}^{\prime }{|}^2⟩/3$. Rotation and stratification can produce kinetic helicity with opposite signs in the two hemispheres, allowing a hemispherically organized alpha effect. These formulae explain a possible mechanism, but their controlled assumptions are not automatically satisfied by stellar turbulence at very large magnetic Reynolds number.

An alpha-Omega dynamo naturally supports oscillation and migration rather than just monotone amplification. For example, in the northern hemisphere choose local right-handed directions $x$ equatorwards, $y$ azimuthally and $z$ radially outwards. With $U_y=Sz$ and poloidal potential $A\hat{\mathbf{y}}$, the local equations areFor $e^{pt+ikx}$ the local alpha-Omega dynamo wave dispersion isThe growing Parker dynamo wave has phase velocity $-\mathrm{Im}p/k$, equatorward for $\alpha S<0$. Thus a migration direction, an excitation threshold and a finite oscillation period arise from the feedback between the alpha effect and shear, rather than being independently imposed. The local example illustrates the mechanism; it does not determine the global stellar period, parity or latitude range. Dipole and quadrupole parity in a mean-field dynamo depend on the spatial coefficients, coupling between hemispheres and boundary conditions.

A major difficulty is matching that mechanism to the measured internal rotation. Helioseismology constrains differential rotation and the tachocline, giving a plausible strong-shear region near the base of the convection zone. At low northern latitudes a positive radial shear combined with the usual positive convective $\alpha$ would give poleward propagation in the simple local wave model, rather than the equatorward solar butterfly diagram. A different sign or location of the alpha effect, an interface solar dynamo, and magnetic transport must therefore be examined rather than assuming that any alpha-Omega dynamo predicts the observed migration. The stable region below the convection zone may store strong toroidal magnetic flux; magnetic buoyancy can then produce rising tubes and bipolar active regions. Flux-tube emergence also has to account for their low latitudes and Joy's law tilts, not simply for a large interior field.

The Babcock-Leighton mechanism supplies a more directly observable route to regeneration. Tilted bipolar regions disperse; cancellation of opposite polarities and poleward transport of surviving flux alter the polar poloidal magnetic field. Differential rotation subsequently winds that field into a toroidal magnetic field. A flux-transport solar dynamo combines this surface source with meridional circulation in a star and turbulent magnetic diffusivity. Models with poleward surface transport and an equatorward deep return flow can reproduce an equatorward activity belt together with poleward surface-flux migration and an appropriate phase relation of the polar field. This was already an explicit model construction in the 1999 flux-transport calculation. Its success remains conditional on the assumed circulation, transport coefficients and emergence prescription; the internal return flow was not then measured well enough to make these quantities unique. The source requires emerged active regions and therefore also leaves the recovery from a state with almost no spots as a nontrivial problem.

Saturation is another limitation of kinematic models. Exponential growth cannot persist indefinitely: the Lorentz force alters the shear and convection, and changes the correlations that generate the alpha effect. Algebraic dynamo quenching can produce bounded cycles but fitting a quenching coefficient is not a first-principles amplitude prediction. Near conservation of magnetic helicity in a highly conducting fluid adds constraints on the growth and saturation of a large-scale field; boundaries and helicity transport matter. Magnetic buoyancy, intermittency and the interaction of mean and fluctuating fields likewise make a local isotropic closure incomplete. Fluctuating regeneration and flow, or nonlinear interactions, can produce cycle irregularity and solar grand minima, but reproducing one irregular time series does not establish a unique explanation of the Maunder minimum.

Other solar-type stars provide an important independent test. Long-term chromosphere observations already showed both cyclic and irregular activity, with rapidly rotating young stars tending to be more active and older slower rotators more often showing smooth cycles; see the original stellar activity survey. These qualitative trends support a rotational-convective origin and stellar rotation-activity feedback through wind-driven magnetic braking of a solar-type star. The relevant control need not be rotation alone: the stellar activity Rossby number, comparing rotation and the convection time, expresses why stellar structure matters. A predictive theory should explain the diversity of cycle periods and geometries, not assign a separate arbitrary $\alpha$, diffusivity and circulation to every star. Stellar activity measurements also trace magnetic heating or spots, not automatically the full internal magnetic geometry.

By 2001, dynamo theory explained credible mechanisms for regeneration, polarity cycling and migrating activity, and suitably constructed models reproduced important solar patterns. It had not supplied a unique, self-consistent quantitative explanation of the full solar and stellar activity phenomenology. The gap lay in determining transport and regeneration from the actual stellar flows, nonlinear saturation and emergence, and robust predictions across different stars. Three-dimensional calculations supported magnetic generation, but numerical diffusivities and the restricted range of resolved scales limited direct extrapolation to stellar conditions. Agreement of a prescribed-coefficient mean-field dynamo with selected observations is meaningful evidence of feasibility, not a demonstration that its detailed mechanism or parameters have been uniquely established.

Choose the planar magnetic field convention $\mathbf{B}=\nabla \times (A{\mathbf{e}}_z)$, so $B_s=s^{-1}{A}_{\varphi }$ and $B_{\varphi }=-A_s$. The resistive induction equation, with constant magnetic diffusivity and a suitable time-dependent gauge constant removed, reduces to advection and diffusion of this Cartesian magnetic flux function:For the specified differential rotation, its mth azimuthal Fourier mode obeysThe physical field is the real part of the complex mode. Substitute $a=s^{m}g(t)\exp [-im{\Omega }_{0}t-i{s}^{2}f(t)]$. A direct radial differentiation givesOn the left, rigid rotation cancels, leaving $\dot{g}/g-i{s}^2\dot{f}+im\omega s^2$. Equating the constant and quadratic coefficients yieldsThe initial condition fixes $f(0)=0$ and $g(0)=1$, since the constant C was already factored out. Set ${\mu }^2=4i\eta m\omega$. Then the exact quadratic-shear magnetic flux solution follows by differentiation:Indeed $\dot{f}=m\omega {\mathrm{sech}}^2(\mu t)$ and $4i\eta f^2=m\omega {\tanh }^2(\mu t)$. Also $\dot{g}/g=-(m+1)\mu \tanh (\mu t)=-4i\eta (m+1)f$. For an unspecified initial normalization the second formula has a multiplicative constant $C^{\prime }$, but here $C^{\prime }=1$.

For $\eta >0$ and $\omega >0$, choose $\mu =(1+i)\sqrt{2\eta m\omega }$. The asymptotic condition is $\mathrm{Re}(\mu )t\gg 1$: then $f$ tends to the finite constant $m\omega /\mu$ and $g\sim 2^{m+1}\exp [-(m+1)\mu t]$. At any fixed s, differentiating A to obtain the magnetic field multiplies it by factors involving $m/s$ and $sf$, which have no additional exponential time dependence. Therefore its envelope decays asUsing the paper's magnetic Reynolds number definition gives the precise scalingThus the decay time scales as ${\mathrm{Rm}}^{1/2}/{\Omega }_0$ for positive ${\Omega }_0$ and fixed m. The mechanism is magnetic phase mixing under differential rotation: winding produces progressively fine radial structure on which magnetic diffusion acts efficiently. For negative $\omega$, choose the square root with positive real part, replacing $\omega$ by $|\omega |$ in the decay time. If $\omega =0$, the initial harmonic mode simply rotates, with $f=0$, $g=1$; there is no shear-enhanced decay. The printed positive-root and Reynolds-number expressions implicitly use $\omega >0$.

A solar-type main sequence star has a radiative interior and an outer convection zone. Rotating convection and differential rotation continually regenerate its magnetic field by dynamo action. The field is therefore not simply an initially stored field undergoing passive magnetic diffusion. In an alpha-Omega dynamo, the Omega effect winds a poloidal magnetic field into a toroidal magnetic field, while the alpha effect associated with helical motions can regenerate poloidal flux. In a solar dynamo, the Babcock-Leighton mechanism provides another poloidal source through the emergence and dispersal of tilted active regions. Shear near the tachocline and transport through the convection zone help organize the cycle; no single dynamo model is needed to claim that rotation and convection are essential ingredients.

Rotation acts on convective motions through the Coriolis force, altering their correlations and helicity. A useful activity measure is the stellar activity Rossby number, ${\mathrm{Ro}}_{\ast }=P_{\mathrm{rot}}/{\tau }_c$, where ${\tau }_c$ is a convective turnover time. At a given stellar structure, faster rotation ordinarily increases rotational influence and produces greater magnetic activity in the unsaturated regime. Young solar-type main sequence stars commonly rotate rapidly and show strong sunspots, flares, and chromospheric and coronal emission, although their initial rotation rates have substantial scatter. Very rapid rotators enter a saturated regime, so activity does not increase without limit as rotation increases. The saturation mechanism is not uniquely fixed by the correlation itself. The observed unsaturated correlation and saturated regime are documented in primary rotation–activity measurements.

The reverse interaction is wind-driven magnetic braking of a solar-type star. An ionized stellar wind is coupled to the star's open magnetic field. Matter carries angular momentum, and the Maxwell stress tensor also transmits a torque. Magnetic coupling provides an effective lever arm comparable to the Alfvén radius, rather than only the stellar radius. Writing ${\dot{M}}_w>0$ for the outward mass-loss rate, a schematic angular-momentum equation iswhere a geometry factor can be absorbed into an effective $R_A$. This lever-arm interpretation follows from the angular-momentum invariant of a magnetized wind; it is the central result of the original magnetized solar-wind calculation. Magnetic braking can be effective despite a small mass-loss rate, because $R_A$ can greatly exceed the stellar radius.

One can derive the familiar rotation-age scaling in a simple unsaturated model. For a radial open field, take $B_r(R_A)\simeq B_{\ast }(R_{\ast }/R_A{)}^2$, ${\dot{M}}_w\simeq 4\pi R_A^2{\rho }_A{v}_w$, and $v_w^2\simeq B_r(R_A{)}^2/({\mu }_0{\rho }_A)$. Eliminating ${\rho }_A$ givesIf the unsaturated large-scale field satisfies $B_{\ast }\propto \Omega$ and the other factors vary slowly, this gives $\dot{J}=-K{\Omega }^3$ with positive K. Approximating I as constant on this part of the main sequence, integration yieldsThus $\boxed{\Omega \propto t^{-1/2},P_{\mathrm{rot}}\propto t^{1/2}}$ in the late-time part of this approximation. This Skumanich rotation law captures the classic observed secular decline of rotation and magnetic activity; its observational basis is the original rotation and chromospheric-emission analysis. It is a scaling model, not an exact trajectory for every star. With a saturated field, the same simplified wind model instead gives a torque approximately proportional to $-\Omega$, so initially rapid rotation can decay more nearly exponentially until the unsaturated regime is reached.

Over the main sequence lifetime, wind loss therefore generally makes the star rotate more slowly, increases its stellar activity Rossby number, and reduces its unsaturated mean activity. Later changes of moment of inertia, wind properties and open-field topology alter the simple power law. Internal transport of angular momentum can also feed the braked envelope from the radiative interior, so surface rotation is not automatically a measure of the entire star's angular momentum. In a star with weak braking, this secular history differs from the ideal cubic-torque model; an age law alone is not a fundamental dynamo equation.

The stellar rotation-activity feedback also acts on shorter time scales. A growing magnetic field exerts a Lorentz force that changes differential rotation and convective motions, extracting kinetic energy and limiting further growth through dynamo quenching. Conversely, varying shear and poloidal regeneration change the magnetic cycle. The present solar-like regime can sustain repeated solar cycles, with polarity reversals and activity modulation superposed on the much slower rotational evolution. Rotation promotes dynamo activity; dynamo-generated fields brake rotation through the wind and react back on the motions that generate them. This coupled feedback explains the broad progression from young, rapidly rotating, strongly active stars toward older, more slowly rotating, less active solar-type stars, while allowing cycles, irregular modulation and different initial histories.

Consider the complex field before taking its real part. Put $p=\hat{z}\times k=(-k_y,k_x,0)$ and $b=Rp+S\hat{z}$. Both components are perpendicular to $k$, so the solenoidal magnetic-field constraint holds. The time-dependent phase obeysbecause ${\dot{k}}_y=-k_0/T$ and $u\cdot k=k_{0}y/T$. Moreover $\dot{p}=(k_0/T,0,0)$, and magnetic stretching is $b\cdot \nabla u=(R{k}_0/T,0,0)=R\dot{p}$. Thus the changing basis exactly supplies the shear-stretching term. Finally,and ${\nabla }^2$ acts as $-k^2$. Equating coefficients verifies

Take $k_0,T>0$ and first $\alpha \ne 0$. Remove magnetic diffusion by writing $R=e^{-D}G$, $S=e^{-D}H$, whereThen $\dot{G}=-i\alpha H$, $\dot{H}=i\alpha k^{2}G$, and $\ddot{G}={\alpha }^2{k}^{2}G$. Equivalently, rescaling time by $|\alpha |$ and replacing $H$ by $\mathrm{sgn}(\alpha )H$ gives the pair used in the question. Its growing WKB approximation isThe complex relative phase follows directly from $H=i\dot{G}/\alpha$; the printed $H\sim kG$ is sufficient for magnitudes but omits that phase.

The phase-averaged magnetic energy is $(k^2|R{|}^2+|S{|}^2)/(4{\mu }_0)$, hence for a nonzero growing-mode componentAt $t\gg T$, $k\sim k_{0}t/T$, so its leading exponential isThe positive quadratic term permits strong amplification before the negative cubic diffusion term dominates. The leading maximum occurs near $t_{\mathrm{peak}}=|\alpha |T/(\eta k_0)$. Put $\epsilon =\eta k_0/|\alpha |$ and $b_0=|\alpha |k_{0}T$. A well-separated late-time growing interval requiresThe first ensures $t_{\mathrm{peak}}\gg T$; the second ensures the WKB condition at the peak, since $|\dot{k}|/(|\alpha |k^2)\sim {\epsilon }^2/b_0$. For fixed nonzero $b_0$ both follow by taking sufficiently small positive $\eta$. The peak's leading exponent is $b_0/(3{\epsilon }^2)$, whereas the eventual exponent tends to minus infinity. This is mean-field shearing-wave transient amplification, not sustained dynamo action.

Growth from the initial instant is possible explicitly: choose $S(0)=i\mathrm{sgn}(\alpha )k_{0}R(0)$ with $R(0)\ne 0$. The initial logarithmic energy derivative is $2(|\alpha |k_0-\eta k_0^2)>0$ when $\epsilon <1$. Growth is not asserted for every initial polarization, for example an identically zero field. If $\alpha =0$, the amplitudes simply diffuse and the energy is proportional to $[k^2|R_0{|}^2+|S_0{|}^2]e^{-2D}$; shear can still give algebraic transient amplification of an in-plane field for small diffusivity, but a purely vertical field only decays. Positive diffusivity is essential to the final decay.

Consider a local shearing sheet centered at a circular orbit. Take uniform mass density, background velocity ${\mathbf{u}}_0=-2Ax{\mathbf{e}}_y$, rotation $\mathbf{\Omega }=\Omega {\mathbf{e}}_z$ and a weak uniform vertical magnetic field ${\mathbf{B}}_0=B_0{\mathbf{e}}_z$. HereUse ideal magnetohydrodynamics, neglect viscosity, magnetic diffusion and vertical stratification, and take incompressible axisymmetric disturbances proportional to $e^{ikz-i\omega t}$. They have horizontal velocity and field perturbations. Let $\mathbf{b}=\delta \mathbf{B}/\sqrt{{\mu }_0\rho }$ be the magnetic perturbation in velocity units, $v_A=B_0/\sqrt{{\mu }_0\rho }$ the Alfvén speed and ${\omega }_A=k{v}_A$ the Alfvén frequency. Magnetic pressure can be included in the perturbation's total pressure; the horizontal Lorentz force for this mode is the magnetic tension term $i{\omega }_A\mathbf{b}$.

Linearizing momentum and induction givesThe term $-2A{b}_x$ winds a radial magnetic perturbation into an azimuthal one. Omitting it would remove the shear's magnetic energy source and give the wrong stability criterion.

A useful elimination that also covers neutral modes introduces the fluid displacement $\mathbf{\xi }$. Frozen-in flux gives $\mathbf{b}=i{\omega }_A\mathbf{\xi }$, while the Eulerian perturbation velocity is $u_x={\dot{\xi }}_x$, $u_y={\dot{\xi }}_y+2A{\xi }_x$, because the displacement samples the gradient of the basic flow. The equations become the spring model of magnetorotational instability:For the stated time dependence, their determinant isThus the ideal magnetorotational dispersion relation isIt follows without dividing by $\omega$, so the neutral stability boundary is included.

To analyze it, put $a={\omega }_A^2$ and $Y={\omega }^2$. ThenFor a hydrodynamically stable rotation law, ${\kappa }_r^2>0$, the sum of the two roots is positive. Their product is negative precisely when $0<a<4\Omega A$. One squared frequency is then negative, giving an exponentially growing mode $\omega =i\gamma$ and its decaying partner. Therefore the magnetorotational instability, or MRI, has criterionAngular velocity must decrease outwards, although specific angular momentum may increase outwards and satisfy Rayleigh's circulation criterion. This is why a Keplerian disc can be stable to the axisymmetric hydrodynamic disturbances of Question 2 yet unstable once magnetic coupling is introduced. At zero magnetic frequency the unstable root becomes neutral, while sufficiently large magnetic frequency restores stability. For increasing angular velocity, $A<0$ with $\Omega >0$, this ideal vertical-field instability is absent.

The growth rate isFor $0<A<\Omega$, differentiation with respect to $a$ sets $\sqrt{{\kappa }_r^4+16{\Omega }^{2}a}=4{\Omega }^2$, yielding the maximum growth rate of ideal magnetorotational instability:For Keplerian rotation, $A=3\Omega /4$, soThus the fastest growth is on an orbital timescale. The field strength selects its wavelength, $k_{\max }\propto 1/v_A$, rather than the maximum ideal growth rate. At fixed $k$, growth does approach zero as $B_0\to 0$; retaining the finite maximum in that limit requires shorter wavelengths, which magnetic diffusion can eventually exclude. Unstable wavelengths satisfy $\lambda >2\pi v_A/\sqrt{4\Omega A}$ in this local model. A finite disc must accommodate such a wavelength, so a sufficiently strong field can suppress the available vertical modes through excessive tension.

Physically, a bent field line acts like a spring joining neighboring fluid elements. The inner element rotates faster and pulls the outer element forward, transferring angular momentum outward. The inner element loses angular momentum and moves inward; the outer element gains it and moves outward. In a rotation law with outward-decreasing angular velocity, these motions increase their angular-velocity difference and the magnetic coupling produces further separation. A weak restoring spring therefore creates a feedback instability by enabling angular-momentum exchange. Without that exchange, each displaced element keeps its angular momentum and a positive epicyclic frequency restores its orbit. Strong magnetic tension instead prevents the separation. The free energy is the differential rotation, not the initial magnetic energy.

The instability provides an important route to turbulence and angular momentum transport in sufficiently conducting accretion discs. The radial flux of azimuthal momentum includes the stressThe first term is a Reynolds stress and the second is the magnetic contribution from the Maxwell stress tensor, with the sign appropriate to momentum flux. For the growing branch, take $u_x$ real and positive. The linear equations give $u_y/u_x=({\gamma }^2+a)/(2\Omega \gamma )>0$ and $b_y/b_x=({\gamma }^2+a-4\Omega A)/(2\Omega \gamma )<0$. The last sign follows from $({\gamma }^2+a)({\gamma }^2+a-4\Omega A)+4{\Omega }^2{\gamma }^2=0$. Thus both the velocity and magnetic contributions carry angular momentum outward. The associated energy extraction is $-T_{xy}rd\Omega /dr=2A{T}_{xy}>0$. Inward accretion can then accompany outward angular-momentum transport, and dissipation of the resulting fluctuations heats the disc. This gives a physical mechanism behind an effective turbulent-viscosity description such as an alpha disk, whereas microscopic viscosity alone is generally too weak for rapid accretion.

MRI is an exponentially growing magnetic instability of outward-decreasing angular velocity; hydrodynamic shearing-wave amplification is transient and does not by itself establish sustained turbulence. Nonlinear saturation, field regeneration and transport efficiency depend on field geometry, boundaries and thermodynamics. Poor conductivity or important nonideal magnetic effects can suppress or modify the ideal instability, so the local dispersion relation is a mechanism and growth criterion rather than a universal prescription for the saturated stress.

Let $\ell$ be the characteristic radius over which the differential rotation changes, with angular velocity of order $\Omega$. The winding source has magnitude $F\sim B_0\Omega$. For times too short for appreciable magnetic diffusion, the zero initial toroidal magnetic field therefore grows as $B\sim B_0\Omega t$. In the steady balance, $F$ is opposed by $\eta B/{\ell }^2$, giving $B_{\mathrm{steady}}\sim B_0\Omega {\ell }^2/\eta$. Comparing these amplitudes gives the diffusive establishment time of a wound toroidal field:With $R_m=\Omega {\ell }^2/\eta$, this is $R_m$ rotation times and the steady toroidal amplitude is of order $R_m{B}_0$. Thus high electrical conductivity makes establishment very slow even though winding begins immediately.

More precisely, if $\mathcal{L}={\nabla }^2-s^{-2}$ on regular toroidal magnetic fields, the solution with zero initial value satisfies $B(t)=B_{\mathrm{steady}}-e^{\eta t\mathcal{L}}{B}_{\mathrm{steady}}$. The residual relaxes diffusively. In the unbounded domain, ${\ell }^2/\eta$ estimates equilibration of the rotating region; exact establishment at every radius is only an asymptotic statement, with diffusion on scale $r$ taking order $r^2/\eta$. As $\eta \to 0$, the time diverges and the prescribed-flow ideal solution continues its secular winding. In a real fluid, magnetic tension and Lorentz force feedback can modify the differential rotation before the very large kinematic steady amplitude is attained.

Take zero mean flow locally and incompressible fluctuations. Subtracting the averaged induction equation givesThe first-order smoothing approximation omits the second term on the right when solving for $\mathbf{b}$, but retains $⟨\mathbf{u}\times \mathbf{b}⟩$ in the mean equation. For a locally uniform, slowly changing mean magnetic field the forcing becomes $(\overline{\mathbf{B}}\cdot \nabla )\mathbf{u}$.

There are two controlled ways to make the omitted fluctuation term small. At small eddy magnetic Reynolds number $R_m=u_{\mathrm{rms}}\ell /\eta$, diffusion gives $b/B\sim R_m$, so the omitted $ub/\ell$ is smaller than the retained $uB/\ell$ by $R_m$. Alternatively, if $u_{\mathrm{rms}}\tau /\ell \ll 1$ for the fluctuation response time $\tau$, then $b/B\sim u_{\mathrm{rms}}\tau /\ell$ and the short-memory response is perturbative. Mean-field scale separation is additionally required. Merely calling the turbulence weak or assuming that the mean magnetic field is small does not justify the omission.

With homogeneous initial transients removed, the causal linear solution isInsert this into the electromotive correlation. Isotropy gives ${\alpha }_{ij}=\alpha {\delta }_{ij}$, so taking the trace yields the helicity formula for isotropic first-order smoothing:Using the alternating-tensor sign givesThe minus sign follows from ${ϵ}_{imn}=-{ϵ}_{min}$. This derivation specifies the diffusion and correlation weighting, rather than asserting that alpha is proportional to instantaneous helicity for arbitrary turbulence.

If $\eta \tau /{\ell }^2\ll 1$, diffusion over a correlation time is negligible. Define the helicity correlation time byFor an ordinary decaying correlation with nonzero mean kinetic helicity, this gives the familiar estimateIt is the kinetic helicity, not magnetic helicity. Opposite flow handedness reverses the alpha effect. The time-integral formula remains meaningful if the equal-time helicity vanishes or its correlation changes sign, when replacing it by a single positive correlation time may fail. In the same weak-diffusion short-memory approximation the gradient forcing $-(\mathbf{u}\cdot \nabla )\overline{\mathbf{B}}$ gives $\beta =\frac{1}{3}{\int }_0^{\infty }⟨\mathbf{u}(t)\cdot \mathbf{u}(t-\tau )⟩d\tau \simeq {\tau }_U{u}_{\mathrm{rms}}^2/3$; the two correlation times need not be equal.

In the small-$R_m$, slowly varying quasistatic regime, diffusion is instead essential. The response givesFor velocity concentrated at wavenumber $k_0$, this becomes $\alpha =-⟨\mathbf{u}\cdot \mathbf{\varpi }⟩/(3\eta k_0^2)$. Its effective memory is the magnetic diffusion time on the eddy scale. The controlled low-$R_m$ setting is developed in Moffatt's analysis of turbulent dynamo action; the causal-response calculation above also covers short-memory turbulence without assuming that diffusion is always negligible.

The geodynamo operates in the electrically conducting liquid of Earth's outer core. buoyancy from thermal and compositional convection supplies motion, while the Coriolis force organizes rotating flow. Rotation together with radial stratification permits kinetic helicity, typically with opposite signs in the two hemispheres. In mean-field language this allows an alpha effect that converts poloidal and toroidal components into one another. A large-scale alpha-squared dynamo is therefore a plausible regeneration mechanism even without dominant differential rotation.

The growth condition compares this regeneration with the sum of molecular and turbulent magnetic diffusion on the core scale. A spherical model must impose the appropriate magnetic matching at the core boundaries and permit spatially varying, generally anisotropic alpha and diffusivity tensors. Rotation selects an axis, and stratification and boundaries select positions, so homogeneous isotropic turbulence is only a local idealization. Such a model can favor a predominantly axial dipole, but the preferred parity and threshold must be found from its global eigenproblem rather than deduced from the plane-wave formula alone. If mean shear is significant, alpha–Omega regeneration can also contribute.

The first-order-smoothing helicity formula is a useful guide to the sign and magnitude of alpha, not an automatically controlled numerical closure for the core. The important convective motions need not have small eddy magnetic Reynolds number or short enough correlation time; scale separation can also be limited. Furthermore, the Lorentz force modifies the flow and its correlations once the field grows. Saturation therefore requires nonlinear feedback, including dynamo quenching and magnetic-helicity constraints where relevant.

The observed field varies and occasionally reverses. A constant-coefficient growing kinematic eigenmode does not itself explain these reversals or select a saturated amplitude. Spatial and temporal variations, competing modes, and nonlinear flow–field interactions must be included. mean-field dynamo theory supplies a physical regeneration mechanism and a testable threshold, while the geodynamo requires a rotating, bounded, nonlinear realization of that mechanism.

Magnetoconvection is the competition between thermally driven convection and the force and transport of a conducting fluid's magnetic field. An imposed magnetic field resists motions that bend its field lines through magnetic tension, while magnetic diffusion permits slippage. The competition changes the onset threshold, preferred spatial scale, temporal behaviour and finite-amplitude organization. The geometry matters: a disturbance nearly invariant along a horizontal imposed magnetic field can avoid much of its tension, whereas an impermeable layer with a vertical imposed magnetic field cannot generally avoid bending it.

A useful quantitative example is a nonrotating Boussinesq layer of depth $d$, heated from below with temperature drop $\Delta T>0$ and threaded by a vertical magnetic field $B_0\hat{\mathbf{z}}$. Let $\nu$, $\kappa$, $\eta$ denote kinematic viscosity, thermal diffusivity and magnetic diffusivity. In thermal-diffusion units define the Rayleigh number, Prandtl number, magnetic-to-thermal diffusivity ratio and Chandrasekhar number byUse stress-free impermeable boundaries, fixed temperature perturbation zero, and magnetic boundary conditions keeping the field vertical. More explicitly, the normal modes below satisfy $w={\partial }_z^{2}w=\theta =b_x=b_y={\partial }_z{b}_z=0$ at $z=0,1$. For a horizontal wavenumber $a>0$, take $w,\theta \propto \sin \pi z$, $b_z\propto \cos \pi z$, and put $s=a^2+{\pi }^2$. The linear temperature and resistive induction equation give $(\lambda +s)\theta =w$ and $(\lambda +\zeta s)b_z={\partial }_{z}w$. Applying the solenoidal projection to the magnetohydrodynamic momentum equation removes pressure; the magnetic coupling is $Q\zeta {\partial }_z\mathbf{b}$, and the vertical projected buoyancy has coefficient $R{a}^2/s$. Substituting the temperature and magnetic responses gives the vertical-field magnetoconvection dispersion relationThese stated boundary conditions make the single vertical mode exact in this model; other boundaries change numerical thresholds.

For steady marginal convection, set $\lambda =0$. The threshold curve isThe imposed magnetic field adds a positive penalty. Differentiating with respect to $x=a^2$ shows that the minimizing wavenumber obeysAt $Q=0$ this gives $a_c=\pi /\sqrt{2}$ and $R_c=27{\pi }^4/4$. With increasing $Q$, the root moves to larger $x$ because the left-hand side is strictly increasing for $x\ge {\pi }^2/2$. The strong-field wavenumber selection in magnetoconvection limit isThus a strong vertical imposed field raises the steady threshold and selects narrower convection cells. At small horizontal wavenumber the incompressibility constraint makes a given vertical motion require large horizontal displacement and magnetic bending. Increasing horizontal wavenumber reduces this magnetic penalty, but increases viscous and thermal diffusion. Their balance selects the narrow cells. This statement is specific to the vertical-field geometry: wavenumber selection in oblique-field magnetoconvection can instead favour nearly field-invariant motions.