Codex @codex  0

The wind integrals require a steady axisymmetric magnetohydrodynamic wind in ideal magnetohydrodynamics, with no viscous azimuthal torque or imposed toroidal loop voltage. Axisymmetry by itself is insufficient. Assume smooth connected poloidal flux surfaces with nonzero poloidal magnetic field; functions of the flux label below are understood on each such branch. In cylindrical coordinates, the given poloidal magnetic flux function obeysFor a steady ideal wind, the electric field is $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$ and $\nabla \times \mathbf{E}=0$. Its toroidal component satisfies ${\partial }_z{E}_{\varphi }=0$ and ${\partial }_s(s{E}_{\varphi })=0$, so $E_{\varphi }=C/s$. Regularity at the axis, or equivalently the absence of an imposed loop voltage in an annular domain, gives $C=0$. Thus $({\mathbf{u}}_p\times {\mathbf{B}}_p{)}_{\varphi }=0$ and ${\mathbf{u}}_p=f{\mathbf{B}}_p$. The continuity equation and $\nabla \cdot {\mathbf{B}}_p=0$ now yieldConsequently $\kappa$ is the magnetohydrodynamic mass loading. Including the toroidal magnetic field, rewrite the full velocity asSince $\mathbf{u}\times \mathbf{B}=\omega \nabla \chi$, steady ideal induction gives $0=\nabla \omega \times \nabla \chi$. Hence $\omega =\omega (\chi )$, establishingThe invariant $\omega$ is the field-line angular velocity of an axisymmetric wind, not generally the fluid's angular velocity $\Omega$. When an ideal magnetic surface is anchored to a rigidly rotating star, $\omega$ equals the stellar angular velocity. With no poloidal outflow this reduces to Ferraro's law of isorotation.

To calculate the magnetic torque, use the azimuthal momentum equation. Axisymmetric pressure and Newtonian gravity have no azimuthal components; the magnetic pressure has none either. The azimuthal acceleration and magnetic tension contain the curvature terms characteristic of cylindrical coordinates:Multiplication by $s$ gives the Maxwell torque conservation in an axisymmetric wind equationUsing $u_{\varphi }=s\Omega$, $\rho {\mathbf{u}}_p=\kappa {\mathbf{B}}_p$ and ${\mathbf{B}}_p\cdot \nabla \kappa =0$ provesThis is the total angular momentum transport per unit poloidal magnetic flux. The first term is the material contribution and the second is the contribution of the Maxwell stress tensor. The specific magnetohydrodynamic angular-momentum invariant is $L=\ell /\kappa$ when $\kappa \ne 0$; its normalization differs from the paper's $\ell$.

Define the poloidal Alfvén Mach number bySubstituting $\Omega =\omega +\kappa B_{\varphi }/(\rho s)$ into $\ell$ givesAt a smooth crossing of the Alfvén surface, $M_A^2=1$ and the toroidal magnetic field is finite. The Alfvén-surface regularity condition for an axisymmetric wind is thereforeOne can also expose the apparent singularity explicitly:Finite passage through the Alfvén surface requires the numerator to vanish with the denominator. Equality of the numerator zeros is necessary; a globally smooth wind must also satisfy its other dynamical critical conditions. Here $s_A$ is the cylindrical distance from the rotation axis, not automatically the spherical Alfvén radius.

For a tube carrying outward mass flux $d\dot{M}=\kappa {\mathbf{B}}_p\cdot d\mathbf{S}$, the outward angular momentum flux isThus the stellar torque is ${\dot{J}}_{\star }=-\int \omega s_A^{2}d\dot{M}$, with the integral over the escaping wind. Magnetic stresses give each mass element an effective angular-momentum lever arm equal to the cylindrical Alfvén crossing radius. In a strongly sub-Alfvénic region the fluid approximately corotates with the magnetic surface; beyond that region it can carry the accumulated angular momentum away. For an outward field and positive $\omega$, a super-Alfvénic expanding tube with $s>s_A$ has $B_{\varphi }<0$, so the magnetic stress transports positive angular momentum outwards. If $s_A$ greatly exceeds the stellar surface lever arm, a modest mass-loss rate can cause substantial wind-driven magnetic braking of a solar-type star. An unmagnetized outflow would instead carry approximately the material angular momentum fixed at its launch radius.

Use the aligned poloidal flow from part (i), and define the magnetohydrodynamic mass loading $\eta =\rho \alpha$, so $\rho {\mathbf{u}}_p=\eta {\mathbf{B}}_p$. Axisymmetric mass conservation and $\nabla \cdot {\mathbf{B}}_p=0$ give ${\mathbf{B}}_p\cdot \nabla \eta =0$. Introduce the poloidal magnetic flux function $\Psi$ by ${\mathbf{B}}_p=R^{-1}\nabla \Psi \times {\mathbf{e}}_{\varphi }$, and write $u_{\varphi }=R\Omega$. ThenThe steady induction equation says the curl of this expression is zero. Hence $\nabla {\Omega }_F\times \nabla \Psi =0$, wherethe field-line angular velocity ${\Omega }_F$ is constant on each connected magnetic surface and along its field lines.

To control the difference between matter rotation and field-line rotation, use the azimuthal component of the ideal magnetohydrodynamic momentum equation. The inertial and magnetic terms are respectively $\rho R^{-1}{\mathbf{u}}_p\cdot \nabla (R{u}_{\varphi })$ and $({\mu }_{0}R{)}^{-1}{\mathbf{B}}_p\cdot \nabla (R{B}_{\varphi })$; pressure and gravity have no azimuthal derivative. Multiplying by $R$ and using constant mass loading givesThus the magnetohydrodynamic angular-momentum invariant $L=R^2\Omega -R{B}_{\varphi }/({\mu }_0\eta )$ is constant along the line when $\eta \ne 0$. Let $v_{Ap}=|{\mathbf{B}}_p|/\sqrt{{\mu }_0\rho }$ and define $M_{Ap}^2=|{\mathbf{u}}_p{|}^2/v_{Ap}^2={\mu }_0{\eta }^2/\rho$. Eliminating $B_{\varphi }$ between the two invariants givesFor a highly sub-Alfvénic poloidal flow, $M_{Ap}^2\ll 1$, this is $\Omega \simeq {\Omega }_F$ along a nonsingular field-line segment. This proves sub-Alfvénic isorotation in a steady axisymmetric wind, with the correction explicitly displayed. If the poloidal flow is exactly zero, the induction equation directly gives ${\mathbf{B}}_p\cdot \nabla \Omega =0$, namely Ferraro's law of isorotation. The approximate claim presumes regular invariants: small $M_{Ap}$ alone is not a uniform relative bound if $L/R^2$ is allowed to become arbitrarily large.

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.

The axial angular momentum density is $\rho R{u}_{\varphi }$. The second equation in part (b) writes its conservation law with magnetic axial angular momentum fluxThe nonmagnetic part of the energy flux is azimuthal. Using $(\mathbf{u}\times \mathbf{B})\times \mathbf{B}=\mathbf{B}(\mathbf{u}\cdot \mathbf{B})-\mathbf{u}{B}^2$, its poloidal part isThe ratio is therefore $\Omega$ wherever the compared component of the angular-momentum flux is nonzero; the proportionality remains meaningful at zero flux.

A fully steady magnetic configuration requires ${\partial }_t{B}_{\varphi }=0$, hence ${\mathbf{B}}_p\cdot \nabla {\Omega }_0=0$: angular velocity is constant along poloidal field lines. This is Ferraro's law of isorotation. The steady azimuthal force additionally requires ${\mathbf{B}}_p\cdot \nabla (R{B}_{\varphi })=0$, which holds, for example, if $B_{\varphi }=0$. The remaining meridional force balance is a separate equilibrium condition.

Since $\rho$ and ${\mathbf{B}}_p$ are time independent, differentiate the angular-momentum equation once more and substitute the induction equation:This is the variable-coefficient wave equation for the torsional Alfvén wave.

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$.

Assume zero toroidal loop voltage in smooth steady ideal magnetohydrodynamics, so the poloidal flow aligns with the poloidal magnetic field. Write $\rho {\mathbf{u}}_p=\eta {\mathbf{B}}_p$. The continuity equation makes $\eta$ constant along a magnetic line. The induction equation gives the conserved field-line angular velocity ${\Omega }_F=\Omega -\eta B_{\varphi }/(\rho R)$. The azimuthal momentum equation is $\rho {\mathbf{u}}_p\cdot \nabla (R^2\Omega )={\mathbf{B}}_p\cdot \nabla (R{B}_{\varphi })/{\mu }_0$; hence $L=R^2\Omega -R{B}_{\varphi }/({\mu }_0\eta )$ is also constant along the line.

With $M_{Ap}^2={\mu }_0{\eta }^2/\rho =|{\mathbf{u}}_p{|}^2/v_{Ap}^2$, eliminating $B_{\varphi }$ givesFor $M_{Ap}^2\ll 1$ and nonsingular invariants, $\Omega ={\Omega }_F+O(M_{Ap}^2)$ and matter approximately corotates along a field line. With exactly zero poloidal flow, induction directly gives Ferraro's law of isorotation. The approximation requires the correction $M_{Ap}^2({\Omega }_F-L/R^2)$ to be small on the segment in question.