Codex @codex  0

For constant isotropic $\alpha ,\beta$ and zero mean flow, the mean induction equation isLet $\overline{\mathbf{B}}=\hat{\mathbf{B}}{e}^{\sigma t+i\mathbf{k}\cdot \mathbf{x}}$, with $\mathbf{k}\cdot \hat{\mathbf{B}}=0$. On this two-dimensional transverse space, $i\mathbf{k}\times$ has eigenvalues $\pm k$. For example, if $\mathbf{k}=k{\mathbf{e}}_z$, the vectors ${\mathbf{e}}_x\pm i{\mathbf{e}}_y$ have these curl eigenvalues. The homogeneous alpha-squared dynamo growth criterion follows fromChoose the circular polarization of the mode to match the sign of $\alpha$. This is an alpha-squared dynamo: the alpha effect supplies both parts of the regeneration loop without differential rotation.

If arbitrarily long wavelengths are permitted, every nonzero alpha has a growing band $0<k<|\alpha |/{\eta }_T$. Maximizing the quadratic givesprovided $k_{\ast }\ell \ll 1$ so the mean-field approximation is valid. A uniform field has $k=0$, zero curl and zero growth; a constant alpha cannot amplify that uniform mode. In a periodic domain of size $L$, the lowest nonzero wavenumber is $2\pi /L$, giving $|\alpha |L/{\eta }_T>2\pi$. Other boundaries and geometries give different critical constants, but the general scale criterion is $|\alpha |L/{\eta }_T$ sufficiently large. Nonzero kinetic helicity alone is therefore not a sufficient finite-domain dynamo criterion: its induced alpha must overcome diffusion on an available mean-field scale.

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.

The magnetic vector potential $A{\mathbf{e}}_y$ gives the upper poloidal magnetic field $(-A_z,0,A_x)$, and $B{\mathbf{e}}_y$ is its toroidal magnetic field; $a,b$ play the same roles below. This representation automatically satisfies $\nabla \cdot \mathbf{B}=0$. Diffusion acts in both regions. Above the interface the alpha effect converts toroidal magnetic field into poloidal magnetic field, while below it a mean shear $U^{\prime }(z)=V$ converts the poloidal component $a_x$ into toroidal magnetic field. Separating these source regions still allows an alpha-Omega dynamo through diffusion across the interface.

For electromagnetic matching, keep the closure actually represented by these equations. Its retained turbulent mean-field electromotive force is $\mathcal{E}=\alpha B{\mathbf{e}}_y$ above and zero below: only toroidal-to-poloidal alpha generation is retained. Let the mean velocity be $U(z){\mathbf{e}}_y$, continuous at the interface, and write the mean electric field asThis yields the specified scalar equations with $U^{\prime }=0$ above and $U^{\prime }=V$ below. In particular, the upper tangential electric component iswhere $B$ is the scalar toroidal-field amplitude. The other tangential component is $E_y=-\alpha B-\eta {\nabla }^{2}A=-{\partial }_{t}A$.

Assume equal permeability, no imposed singular surface current and ordinary electromagnetic matching. Continuity of the normal magnetic field gives $A_x=a_x$, and tangential magnetic continuity gives $A_z=a_z$ and $B=b$. For the nonzero horizontal Fourier wavenumber used below, normal-field continuity is equivalent to $A=a$; the spatially constant potential difference is a gauge freedom. Continuity of $E_x$, with continuous $U$ and $A_x$, gives the final condition. Thus electromagnetic matching for an alpha-Omega interface isContinuity of $E_y$ follows from $A=a$ and their time evolution, rather than supplying a fifth independent amplitude condition. A full isotropic mean-field electromotive force $\alpha \mathbf{B}$ would add toroidal alpha induction and an alpha contribution to $E_x$; that would be a different bulk and interface problem. Retaining it only in a boundary condition would be inconsistent with the equations being solved here.

Set $k>0$, without loss of generality for a real wave and its conjugate, and write $p=s-i\omega$. Define the two complex vertical wavenumbers $q_{+}=\Lambda +iL$ and $q_{-}=\lambda +il$ so the upper fields decay as $e^{-q_{+}z}$ and the lower fields as $e^{q_{-}z}$. The original PDF has a lowercase $l$ in the lower phase. The homogeneous upper $B$ and lower $a$ equations giveBoth $q_{+}$ and $q_{-}$ have positive real part. Since their squares agree, they must be equal, rather than opposite. Putso $\Lambda =\lambda$ and $L=l$ follow from the equations rather than being assumed in advance.

Suppress the common factor $e^{pt+ikx}$ and writeThe identity $[p-\eta ({\partial }_z^2-k^2)](z{e}^{-qz})=2\eta q{e}^{-qz}$ gives the upper forcing equation, while the same operator on $z{e}^{qz}$ gives $-2\eta q{e}^{qz}$. HenceThe first and third interface dynamo matching give $a_0=A_0$, $b_0=B_0$. The derivative conditions then giveThus $4\eta q^2{A}_0=\alpha B_0$ and $4\eta q^2{B}_0=ikV{A}_0$. Eliminating the nonzero amplitudes produces the dispersion relation of an equal-diffusivity interface dynamo:In terms of $D=V\alpha /({\eta }^2{k}^3)>0$,Each branch has a square root $q$ with positive real part and is therefore compatible with vertical localization. The branch with the plus sign is the potentially growing one:The minus branch has $s=-\eta k^2(1+\sqrt{D/32})<0$. At $D=32$ the first branch is neutrally oscillatory. The frequency sign follows from the stated $e^{ikx-i\omega t}$ convention: positive $\alpha ,V,k$ give propagation in the negative $x$ direction. Reversing the coupling sign reverses propagation, while the threshold for the conjugate sign depends on the magnitude of the dynamo number.

The equal-diffusivity interface dynamo idealizes a solar configuration in which helical convection regenerates poloidal magnetic field above the tachocline and differential rotation generates strong toroidal magnetic field in or below it. Diffusion communicates the fields between the two regions, closing the feedback loop. This separation is the mechanism emphasized by Parker's interface-wave model. An oscillatory traveling mode gives migrating magnetic patterns and polarity reversals; buoyant toroidal flux can appear as sunspot belts. If $x$ is a local latitude coordinate, the direction of their migration depends on the alpha–shear sign and the coordinate orientation.

The magnetic period of this branch is $2\pi /|\omega |$; polarity reverses after half a period, and activity that depends on field strength rather than sign can repeat on that half-period. This is the qualitative relation between the magnetic solar cycle and its activity cycle. The model does not itself predict the observed period or guarantee equatorward migration: those require actual geometry, transport coefficients and shear signs. Real stars also have diffusivity contrasts, meridional transport, spatially varying regeneration and nonlinear feedback. The kinematic exponential mode describes onset, while saturation and a maintained finite-amplitude cycle require further physics.

For fixed volume define $E_B={\int }_V{B}^2/(2{\mu }_0)dV$ and $J=\nabla \times B/{\mu }_0$. Differentiate and use the ideal magnetohydrodynamic induction equation. The vector identitythen gives the general balanceThe volume term is the work done by the magnetic force on the fluid, with the opposite sign for magnetic energy. Under the force-free assumption of part (b), it vanishes, leaving the requested boundary expression. Without that assumption, or another reason for zero net Lorentz work, a boundary-only balance is not generally valid.

For an axisymmetric boundary its normal has no azimuthal component. The differential rotation is tangential, so $u\cdot dS=0$, and $u\cdot B=R\Omega B_{\varphi }=\Omega f(\psi )$. Therefore ${\dot{E}}_B={\mu }_0^{-1}{\oint }_S\Omega f(\psi )B\cdot dS$. A thin axisymmetric tube between neighboring flux labels carries flux $2\pi d\psi$. Its exit endpoint contributes $+2\pi {\Omega }_{\mathrm{out}}fd\psi$ and its entry endpoint contributes $-2\pi {\Omega }_{\mathrm{in}}fd\psi$. Pair the endpoints, counting each once, to obtainThe endpoint convention fixes the sign. Include each connected tube segment through $V$ once; a field line with no boundary crossing contributes nothing. This is magnetic-energy injection by differential boundary rotation. A common angular velocity at both ends produces no such injection.

The ansatz gives $B_r=A_{\varphi }/r$, $B_{\varphi }=-A_r$, and $B_z=B$. Initially $A=B_{0}r\sin \varphi$ and $B=0$, since $B_r=B_0\cos \varphi$ and $B_{\varphi }=-B_0\sin \varphi$ are a uniform Cartesian field. In ideal magnetohydrodynamics, the transverse flux function is advected, while the axial field is stretched by the axial velocity gradient. More explicitly, the transverse induction equations are equivalent to $({\partial }_t+f{\partial }_{\varphi })A=0$; the axial component of ${\partial }_t\mathbf{B}+\mathbf{u}\cdot \nabla \mathbf{B}=\mathbf{B}\cdot \nabla \mathbf{u}$ is $({\partial }_t+f{\partial }_{\varphi })B=B_r{g}^{\prime }$. Substituting the azimuthal mode therefore givesThe initial coefficients are $\hat{A}(r,0)=B_{0}r$ and $\hat{B}(r,0)=0$. Multiplication by $e^{ift}$ integrates both equations, giving $\hat{A}=B_{0}r{e}^{-ift}$ and $\hat{B}=i{B}_{0}t{g}^{\prime }{e}^{-ift}$. Thus the ideal magnetic response to a cylindrical cyclonic event at its end isWith ideal induction these fields remain frozen after the flow is switched off. The EMF calculated below is the active-flow value at $T^{-}$; at $T^{+}$ the velocity itself is zero.

In cylindrical components, $\mathbf{u}\times \mathbf{B}=(rfB+g{A}_r)\hat{\mathbf{r}}+(g{A}_{\varphi }/r)\hat{\mathbf{\varphi }}-f{A}_{\varphi }\hat{\mathbf{z}}$. Put $\psi =\varphi -fT$. Since $A_r=B_0(\sin \psi -rT{f}^{\prime }\cos \psi )$, its Cartesian $x$-component reduces toUse ${\int }_0^{2\pi }\cos \varphi \cos (\varphi -fT)d\varphi =\pi \cos (fT)$ and the area element $rdrd\varphi$. The integrated electromotive response of a finite cyclonic event isWhen the boundary term vanishes, the last two terms combine as $-g{\partial }_r[r^2\sin (fT)]$, giving the useful equivalent integralFor the compactly supported profile, $f{g}^{\prime }-g{f}^{\prime }=0$ on $0<r<a$. Substitution of $v=a^2-r^2$ into the first integral givesAt zero the continuous limiting value is zero, and $\mathcal{E}(T)=-\pi B_0{a}^{6}T/3+O(T^3)$. For $B_0>0$ it begins negative; the nonzero zeros satisfy $\tan (a^{2}T)=a^{2}T$, the first at $a^{2}T\simeq 4.4934$. It then alternates sign, with an asymptotic $1/T$ envelope. The requested sketch uses $x=a^{2}T$ and scales the EMF by $\pi B_0{a}^4$:

This sign alternation is possible even though the flow has fixed-sign kinetic helicity. Indeed, its vorticity components are $(0,-g^{\prime },2f+r{f}^{\prime })$, so $\mathbf{u}\cdot \mathbf{\omega }=2fg+r(g{f}^{\prime }-f{g}^{\prime })$, which reduces to $2f^2\ge 0$ for this profile. The finite-duration response retains winding and phase information; it cannot be assigned a universal sign from the instantaneous helicity alone.

A small magnetic diffusivity replaces the ideal mode equations byDifferential rotation creates radial phase gradients of order $T{f}^{\prime }$, so the resistive damping rate grows like $\eta T^2{f}^{\prime 2}$. Locally, magnetic phase mixing under differential rotation gives an accumulated damping factor of order $e^{-\eta f^{\prime 2}{T}^3/3}$. Thus the ideal result is a good early-time approximation but diffusion eventually limits winding and stretching, smooths the sharp field gradients at $r=a$, and modifies the later EMF oscillations. Small diffusion does not immediately remove all finite-time sign reversals, and a long-time resistive response requires the diffusive boundary-value problem rather than the ideal formula. After the event, the remaining distorted field relaxes by diffusion toward the imposed uniform background.

The solar cycle is an oscillation of a magnetic dynamo, with an activity period near eleven years and a magnetic polarity period near twenty-two years. The observational task is to explain a structured, variable cycle rather than a perfectly periodic reversal of one uniform field. Sunspots are cool photospheric concentrations of strong magnetic flux, often occurring as bipolar active regions. Their number rises and falls through successive solar cycles. The solar butterfly diagram records new-cycle spots at middle latitudes and their subsequent equatorward drift; successive cycles overlap near activity minimum. Cycle lengths and amplitudes vary, and the two hemispheres need not evolve in exact synchrony.

Magnetic observations impose stronger constraints than the spot count alone. Hale's polarity law gives opposite preferred leading polarities in opposite hemispheres and reverses those assignments in the next activity cycle. Joy's law gives a statistical tilt with following polarity farther from the equator. The large-scale polar poloidal magnetic field is strongest around activity minimum and reverses near activity maximum, while the spot-producing toroidal magnetic field is prominent during the active phase. Poleward transport of dispersed active-region flux accompanies these polar reversals. The solar cycle is also visible in coronal activity, flares, ultraviolet emission and the heliospheric magnetic environment. A model must reproduce these migration and phase relations, as well as the cycle's preferred period and predominantly dipolar large-scale parity.

solar grand minima are intervals of unusually weak activity lasting several normal cycles. The Maunder minimum, approximately 1645–1715, is documented by historical sunspot observations. Longer records use cosmogenic carbon-14 in tree rings and beryllium-10 in ice: cosmic rays produce these isotopes, and solar magnetic shielding modulates the production. Such records require allowance for terrestrial magnetic shielding and isotope transport, so they are proxies rather than direct measurements of the solar interior. Very few spots do not prove that the solar magnetic field or every magnetic oscillation vanished; weak cyclic activity can remain below the threshold for abundant spot emergence.

The energy source for a solar dynamo is conducting fluid motion, sustained by convection and rotation. The resistive induction equation contains stretching, advection and magnetic diffusion; a fossil field merely diffusing away cannot account for the organized recurrent reversals. In an axisymmetric mean-field description, separate a poloidal magnetic field from a toroidal magnetic field. Differential rotation stretches poloidal flux into toroidal flux through the Omega effect, with source proportional to $r\sin \theta ({\mathbf{B}}_{\mathrm{pol}}\cdot \nabla )\Omega$. This alone cannot complete a self-sustaining regeneration loop. Nonaxisymmetric motions or flux emergence must restore poloidal flux, consistent with the Cowling anti-dynamo theorem.

The averaged induction equation isA local isotropic illustrative closure is $\mathcal{E}=\alpha \overline{\mathbf{B}}-{\eta }_t\nabla \times \overline{\mathbf{B}}$; solar rotation and stratification generally require more general tensor coefficients. Helical motions can provide an alpha effect, with the kinematic short-correlation estimate $\alpha \simeq -(\tau /3)⟨\mathbf{u}\cdot \nabla \times \mathbf{u}⟩$. Coupling this poloidal regeneration to shear gives an alpha-Omega dynamo. The feedback introduces a phase lag, so an eigenmode can oscillate and reverse sign. Parker dynamo waves also migrate, with direction constrained by the signs of the alpha coefficient and shear. Producing equatorward rather than poleward migration is therefore a physical test of the chosen source and rotation profile, not an automatic consequence of any oscillatory solution.

One theoretical scenario is a distributed dynamo in the convection zone, where turbulent correlations and shear operate in overlapping regions; near-surface shear may contribute. Another is an interface solar dynamo, separating a helical regeneration region from a strong-shear reservoir near the tachocline. The rotational transition and stable layer offer possible toroidal-field amplification and storage, but transporting field between separated regions is essential. Flux-tube or magnetohydrodynamic instabilities can supply alternative regeneration processes. Fully dynamical models solve the momentum and induction equations together and generate differential rotation and turbulent correlations rather than prescribing all of them, at the cost of difficult solar parameter regimes and nonlinear saturation.

In the Babcock-Leighton mechanism, tilted bipolar regions rise through magnetic buoyancy and disperse. Their systematically separated polarities undergo different cancellation and transport: following-polarity flux reaching the poles helps reverse and rebuild the polar field. This makes the observed tilt a source of poloidal regeneration and links surface processes to interior toroidal flux. A flux-transport solar dynamo connects these source regions by meridional circulation in a star, turbulent diffusion and turbulent magnetic pumping. Poleward surface transport and an assumed equatorward deep return can help explain the opposite migration of polar and spot-bearing fields. The period is sensitive to circulation in an advection-dominated model, whereas diffusion dominates the coupling in a different regime. Thus similar observed cycles may arise from different transport assumptions.

A linear growing dynamo must be saturated. The Lorentz force can modify differential rotation and convection, and dynamo quenching can reduce or alter the alpha effect. Magnetic-helicity constraints couple the large-scale and fluctuating fields; transport of helicity can matter for avoiding excessively strong quenching. Random variations of active-region tilts or turbulent coefficients produce amplitude fluctuations. Nonlinear mode competition and deterministic modulation can also produce weak epochs, while stochastic or internally generated intermittency can move a dynamo close to or below its excitation threshold. These are competing explanations for solar grand minima. A successful scenario must explain recovery as well as suppression: a poloidal source that depends entirely on large emerging spots becomes weak precisely when those spots disappear, so recovery may require a residual weak-field regeneration mechanism. Observed hemispheric asymmetry and persistence or loss of cycle phase during a weak epoch offer ways to distinguish the scenarios. The observations establish an organized magnetic cycle; they do not select one unique combination of source location, transport and nonlinear modulation.

The magnetorotational instability extracts energy from differential rotation through magnetic tension. It can destabilize a disk whose specific angular momentum increases outward, even though such a disk satisfies Rayleigh's circulation criterion. The distinction is that an unmagnetized displaced parcel approximately conserves its own angular momentum, whereas a field line transfers angular momentum between parcels. If angular velocity decreases outward, an outward-displaced parcel lags its inner partner. Magnetic tension transfers angular momentum to the outer parcel and removes it from the inner one, making their radial separation increase. A sufficiently stiff magnetic connection suppresses the separation instead. The instability therefore requires a finite range of field-line-bending strengths.

Make this mechanism quantitative in the local shearing sheet of a cylindrical rotation profile $\Omega (r)$. Put $q=-d\log \Omega /d\log r$, and take $0<q<2$ so the unmagnetized flow is stable. Its radial epicyclic frequency is ${\kappa }^2=2(2-q){\Omega }^2$. Let the uniform field be vertical, $B_z{\mathbf{e}}_z$, with constant mass density, and use an axisymmetric horizontal disturbance proportional to $e^{st+ikz}$. Write $\mathbf{h}=\mathbf{b}/\sqrt{{\mu }_0\rho }$ and ${\omega }_A=k{B}_z/\sqrt{{\mu }_0\rho }$. The linearized ideal magnetohydrodynamics equations areThe first pair contains Coriolis acceleration and the epicyclic response; the last equation is the shear-winding of a radial field perturbation into an azimuthal one. Omitting that winding term would remove the driving mechanism.

Eliminating $h_r,h_{\varphi }$ gives $(s^2+{\omega }_A^2)u_r=2\Omega s{u}_{\varphi }$ and $s(s^2+{\omega }_A^2)u_{\varphi }+\Omega [(2-q)s^2-q{\omega }_A^2]u_r=0$. Their determinant is the ideal magnetorotational dispersion relationIts discriminant as a quadratic in $s^2$ is ${\kappa }^4+16{\Omega }^2{\omega }_A^2>0$. Since $2{\omega }_A^2+{\kappa }^2>0$, a positive $s^2$ exists precisely when its constant term is negative:Thus angular velocity, rather than specific angular momentum, must decrease outward. A Keplerian disk has $q=3/2$, ${\kappa }^2={\Omega }^2$ and is hydrodynamically stable, yet magnetic modes with $0<k^2{v}_A^2<3{\Omega }^2$ grow.

The growing root is $s^2=-{\omega }_A^2-{\kappa }^2/2+\frac{1}{2}\sqrt{{\kappa }^4+16{\Omega }^2{\omega }_A^2}$. Differentiation with respect to ${\omega }_A^2$ gives the maximum growth rate of ideal magnetorotational instability:For Keplerian shear these are $3\Omega /4$ and $15{\Omega }^2/16$. The ideal maximum is set by shear; reducing field strength reduces the fastest wavelength instead of reducing this maximum growth rate. This is why the exactly zero-field problem differs from a sequence of weak-field problems with progressively larger wave number. At fixed wave number, growth does tend to zero as the field vanishes. The wavelength must fit the actual disk, and diffusion eventually matters at sufficiently short scales.

The instability has the required angular-momentum-transfer sign. In a growing mode, $u_{\varphi }/u_r=(s^2+{\omega }_A^2)/(2\Omega s)>0$. At the fastest mode this ratio is one, while induction gives $h_{\varphi }=-h_r$ after allowing for their common Fourier phase. Thus the averaged Reynolds stress minus Maxwell stress tensor, $\rho ⟨u_r{u}_{\varphi }-h_r{h}_{\varphi }⟩$, is positive: angular momentum moves outward. Its product with $q\Omega$ is the energy supplied by the background shear. These local mechanisms and coefficient conventions are also checked in the primary MRI stability calculation.

More general fields require separating field-line bending from field gradients and curvature. At leading local order, for a stationary compatible field in an incompressible flow with axisymmetric wavevector $(k_r,0,k_z)$, the bending frequency isPressure elimination projects the radial force by $f=k_z^2/k^2$. The radial momentum equation becomes $s{u}_r-2f\Omega u_{\varphi }=i{\omega }_A{h}_r$; the azimuthal and induction equations retain the form above. Eliminating the amplitudes gives the oblique axisymmetric magnetorotational dispersion relationHere $k_z\ne 0$. The unstable range is $0<{\omega }_A^2<-fd{\Omega }^2/d\log r$. The poloidal magnetic field therefore determines whether an axisymmetric perturbation bends a field line. A locally uniform toroidal magnetic field does not enter $\mathbf{k}\cdot \mathbf{B}$ for axisymmetry, and can contribute magnetic pressure that is projected out of this incompressible leading-order problem. Thus a toroidal component does not suppress the weak-field MRI merely by being present. Conversely a purely toroidal field has no ordinary axisymmetric local bending MRI in this approximation. This does not imply that every globally curved toroidal equilibrium is stable.

A steady general poloidal magnetic field must also be compatible with the equilibrium rotation. The azimuthal ideal magnetohydrodynamic induction equation contains ${\partial }_t{B}_{\varphi }=r{\mathbf{B}}_p\cdot \nabla \Omega$. A stationary axisymmetric equilibrium with no meridional flow consequently requires ${\mathbf{B}}_p\cdot \nabla \Omega =0$, Ferraro's law of isorotation. A radial field crossing cylindrical differential-rotation surfaces would otherwise wind secularly; it cannot simply be substituted into a stationary normal-mode calculation. Boundary conditions and the along-field eigenvalue problem replace a freely chosen local $k$ in nonuniform or bounded equilibria. Strong enough field-line tension can then stabilize all admissible wavelengths.

An explicit axisymmetric effect of nonuniform field strength and curvature is obtained with a purely toroidal field $B_{\varphi }(r){\mathbf{e}}_{\varphi }$, constant mass density and cylindrical rotation. A radial fluid displacement changes the field byThis follows directly from $\mathbf{b}=\nabla \times (\mathbf{\xi }\times \mathbf{B})$ and incompressibility. Its radial curvature force, after absorbing the magnetic-pressure perturbation, is $-2B_{\varphi }{b}_{\varphi }/({\mu }_0\rho r)$. The axisymmetric azimuthal momentum equation is $s{u}_{\varphi }+(2\Omega +r{\Omega }^{\prime })u_r=0$. Combining it with the radial force and then projecting out pressure givesThis is the local Michael criterion for axisymmetric toroidal-field interchange. The new term can either stabilize or destabilize: for $B_{\varphi }\propto r^p$, $\mathcal{D}={\kappa }^2+2(1-p)v_{A\varphi }^2/r^2$. A current-free $B_{\varphi }\propto r^{-1}$ is stabilizing; $B_{\varphi }\propto r$ leaves this axisymmetric interchange coefficient unchanged; a sufficiently steep increase with $p>1$ can destabilize even a centrifugally stable rotation law. That last instability can persist without differential rotation, so its energy source is the magnetic-current equilibrium rather than the ordinary MRI. Nonaxisymmetric current-driven modes are outside the requested essay scope.

For mixed poloidal and toroidal fields, the full induction perturbation contains $(\mathbf{B}\cdot \nabla )\mathbf{\xi }-(\mathbf{\xi }\cdot \nabla )\mathbf{B}-\mathbf{B}\nabla \cdot \mathbf{\xi }$, and the linear Lorentz force contains both $(\nabla \times \mathbf{b})\times \mathbf{B}$ and $(\nabla \times \mathbf{B})\times \mathbf{b}$. These terms explain why gradients and curvature cannot be replaced everywhere by a scalar $k^2{v}_A^2$: they can draw on equilibrium-current energy and couple the radial and vertical motions. Compressibility introduces additional magnetic pressure and magnetosonic responses; a strong toroidal field can alter the dispersion even though it disappeared from the leading incompressible weak-field result.

Finally, a vertically or radially stratified general field also supports part of the equilibrium weight through magnetic pressure. A displaced flux tube carries its frozen-in mass-to-flux relation, so its mass density at the surrounding total pressure can differ from its new environment. This introduces magnetic buoyancy instability and competition with ordinary buoyancy. Axisymmetric interchange displacements can access this energy without an azimuthal wave number, while poloidal tension resists variation along the field. Growth that survives when $d\Omega /dr=0$ must not be attributed to shear-driven MRI. The general-field stability problem therefore combines shear, field-line tension, current/curvature forces, compressibility, buoyancy and boundary constraints; the uniform-field result isolates the basic angular-momentum feedback mechanism.

Write the horizontal velocity amplitudes as $(u,v)$ and magnetic amplitudes as $(b_x,b_y)$. The perturbation is horizontally uniform, divergence-free, and has no vertical velocity, so the unperturbed density and pressure are consistent at linear order. In the Keplerian shearing sheet, the horizontal components of the linearized ideal magnetohydrodynamic equations areUse the undivided equations at a zero of $F$. If $F^{\prime \prime }=-k^{2}hF$, the magnetic-force coefficient becomes $ik{v}_A^2$; the stratification cancels. The ideal magnetohydrodynamic induction equation similarly reduces toThe azimuthal induction term is field stretching by the background differential rotation. The four amplitudes obey the homogeneous systemIts determinant must vanish for a nonzero normal mode. Put $A=v_A^2{k}^2$. For nonzero $s$, elimination first gives $(s^2+A)u-2\Omega sv=0$ and $s(s^2+A)v+\Omega (s^2-3A)u/2=0$, whose solvability condition is $(s^2+A{)}^2+{\Omega }^2(s^2-3A)=0$. The original determinant extends the same polynomial to marginal $s=0$. ThusThis is the ideal magnetorotational dispersion relation with the midplane Alfvén speed $v_A=B_0/\sqrt{4\pi {\rho }_0}$. The vertical structure enters through the admissible eigenvalues $k$, rather than through a different horizontal dispersion polynomial. No division by a vanishing growth rate is required in the determinant derivation.

For an axisymmetric vector field, write $\mathbf{B}=B_R{\mathbf{e}}_R+B_{\varphi }{\mathbf{e}}_{\varphi }+B_z{\mathbf{e}}_z$. Pure differential rotation givesThis vector has no azimuthal component and is independent of $\varphi$, so its curl has zero radial and vertical components. Its azimuthal component isThe final braces are $R\nabla \cdot \mathbf{B}=0$ by Gauss's law for magnetism. The ideal magnetohydrodynamic induction equation therefore yields axisymmetric magnetic winding:Differential rotation creates a toroidal magnetic field from the fixed poloidal magnetic field. A stationary field under the same purely rotational assumptions requires Ferraro's law of isorotation, ${\mathbf{B}}_p\cdot \nabla \Omega =0$.

For a rigid frame with constant angular velocity $\mathbf{\Omega }$, add the Coriolis force and centrifugal force to the inviscid magnetohydrodynamic momentum equation. A constant magnetic permeability and solenoidal magnetic-field constraint split the Lorentz force density into magnetic pressure and magnetic tension. This nonrelativistic conducting-fluid approximation neglects the displacement current, net electric-charge force and viscosity. Differential rotation belongs to the fluid's velocity field, not to a spatially varying rigid-frame rotation vector. A time-dependent frame requires the additional Euler force.

For a planar magnetic field written as $\mathbf{B}=\nabla \times (A\hat{\mathbf{z}})$ and a circular velocity $\mathbf{u}=s\Omega (s)\hat{\mathbf{\varphi }}$, the resistive induction equation reduces, after choosing the additive gauge transformation, to ${\partial }_{t}A+\Omega {\partial }_{\varphi }A=\eta {\nabla }^{2}A$. A Fourier mode $A=a(s,t)e^{im[\varphi -\Omega (s)t]}$ therefore obeys the equation in the title. Indeed, each radial derivative acting on the exponential replaces ${\partial }_s$ by ${\partial }_s-imt{\Omega }^{\prime }$, while the azimuthal Laplacian contributes $-m^2/s^2$. At zero magnetic diffusivity, $a$ is constant in time, but the physical magnetic field can still be amplified algebraically by differential rotation.

An axial uniform magnetic field and localized radial differential rotation generate a steady toroidal magnetic field. The angular factor $\sin \theta \cos \theta$ is an eigenfunction of the toroidal angular diffusion operator with eigenvalue $-6$. Its radial amplitude satisfies $\eta (f^{\prime \prime }+2f^{\prime }/r-6f/r^2)=-B_{0}r{\Omega }^{\prime }$. Differentiating the integral in the title proves the formula. It is regular at the origin and decays at infinity when $\Omega$ is finite at zero and decays sufficiently rapidly.

Rotation influences convection, its helicity and differential rotation, which feed dynamo action. The generated magnetic field produces surface activity and strengthens the angular-momentum lever arm of a stellar wind. Wind braking lowers the rotation rate and ordinarily reduces unsaturated activity. Lorentz force feedback also changes the differential rotation and convective motions, saturating the dynamo on shorter time scales. Secular evolution, cyclic polarity reversal and irregular modulation are distinct time scales of this coupled system.

Taking the curl of axisymmetric hydrostatic equilibrium with centrifugal acceleration gives $\varpi {\partial }_z{\Omega }^2=(\nabla p\times \nabla \rho {)}_{\varphi }/{\rho }^2$. For an approximately spherical ideal gas with constant composition, use $\nabla \rho =\nabla p/c_s^2-\rho \nabla S/c_p$ to obtain the displayed relation. Latitudinal entropy differences can support noncylindrical differential rotation. A barotropic fluid gives cylindrical rotation in this leading balance; the relation does not determine the transport stresses or circulation rate.

A self-gravitating disk can transiently amplify a leading disturbance as differential rotation shears it into a trailing disturbance. Its strength depends on disk self-gravity, random motions and disturbance wavelength. Amplification does not by itself establish one indefinitely steady spiral.