Codex @codex  0

For zero mean velocity and magnetic field in a shearing box, the perturbation energy obeys $\dot{E}\le 2(|A|-4{\pi }^2\eta /L^2)E$ when kinematic viscosity and magnetic diffusivity agree. Bound the work of Reynolds stress and Maxwell stress tensor by $2|A|E$, and use the spectral gap of a shearing-periodic box to bound dissipation below. Thus the displayed bound on magnetic Reynolds number guarantees exponential decay and excludes sustained turbulence.

In a steady axisymmetric star, a torque due to Reynolds stress, Maxwell stress tensor or viscosity must balance advection of specific angular momentum by stellar meridional circulation. Writing the nonadvective angular momentum transport flux as ${\mathbf{F}}_J$ gives the displayed balance. Thus stresses that maintain differential rotation also drive a meridional response. The stellar thermal-wind balance supplies a separate leading force constraint, not a replacement for this torque equation.

With the momentum-flux sign convention used for the Maxwell stress tensor, Maxwell's equations implywhere $\mathbf{g}$ is the electromagnetic momentum density and ${\mathbf{f}}_{\mathrm{L}}$ is the Lorentz force density exerted by the field on matter. Thus total field-plus-matter momentum is locally conserved.

In a planar magnetohydrodynamic shock with both velocity and magnetic field normal to the front, the normal magnetic field is continuous. Its normal Maxwell stress tensor contribution is consequently the same on both sides, and the Poynting vector vanishes because $\mathbf{u}\times \mathbf{B}=0$. The remaining jump conditions are the ordinary Rankine-Hugoniot conditions for a perfect gas. This is the purely parallel branch; an oblique or switch-on shock has additional transverse fields and is not covered by this reduction.

Use the source's Gaussian units and neglect the displacement current, as in nonrelativistic magnetohydrodynamics. Write $\Delta p=p_{\perp }-p_{\parallel }$ for the pressure anisotropy. The pressure tensor is $\mathsf{P}=p_{\perp }\mathsf{I}-\Delta p\mathbf{bb}$, so its force density isAmpère's law and the solenoidal magnetic-field constraint give the Lorentz force decompositionAdding these two force densities and using $\mathbf{BB}=B^2\mathbf{bb}$ yieldsThe first term is the gradient of the modified magnetohydrodynamic total pressure; the second is the anisotropic Maxwell stress tensor contribution together with the directional pressure force. Their effective coefficient is the pressure-anisotropy-modified magnetic tension, $\mathcal{T}=B^2/(4\pi )+\Delta p$. At fixed magnetic field, positive pressure anisotropy increases resistance to bending, while negative pressure anisotropy reduces it. When $\mathcal{T}$ becomes negative, bending experiences an amplifying force rather than a restoring force. Gradients of $\Delta p$ can also drive motion, so the curvature interpretation refers specifically to the bending contribution.

The magnetic second-moment tensor is $M_{ij}=B_i{B}_j$. This is the magnetic dyadic called the Maxwell tensor in the paper; the isotropic part of the usual Maxwell stress tensor is already contained in the magnetohydrodynamic total pressure. The solenoidal magnetic-field constraint givesConsequently the magnetohydrodynamic momentum equation in part a becomes $D_t{u}_i=-{\partial }_{i}p+{\partial }_j{M}_{ij}$, where $D_t={\partial }_t+u_n{\partial }_n$ is the material derivative. Applying its product rule and the ideal magnetohydrodynamic induction equation givesThus the evolution uses only $\mathbf{u},M,p$, together with ${\partial }_i{u}_i=0$. To determine $p$, take the divergence of the magnetohydrodynamic momentum equation. Since ${\partial }_t{\partial }_i{u}_i=0$ and $u_j{\partial }_j{\partial }_i{u}_i=0$,This is the pressure Poisson equation for incompressible flow, with boundary conditions supplied by the imposed normal acceleration and normal stress, or by periodicity or decay at infinity as appropriate. For example, the normal component of the magnetohydrodynamic momentum equation requires $\mathbf{n}\cdot \nabla p=\mathbf{n}\cdot (\nabla \cdot M-D_t\mathbf{u})$ at a prescribed boundary. The pressure is determined up to a spatially constant function of time.

For equivalence to one resolved magnetic field, the initial $M$ must actually be $\mathbf{B}{\mathbf{B}}^T$ for a divergence-free vector field. In particular it is a symmetric positive semidefinite matrix of matrix rank at most one. The tensor equations also admit more general stresses; the averaged $M$ in part d need not have matrix rank one.

Stellar rotation is a problem of coupled force balance, heat transport and angular momentum transport. A star need not rotate as a rigid body: in the Sun, the convective envelope rotates faster near the equator than near the poles, while much of the radiative interior rotates more nearly uniformly. The tachocline is the thin transition near the base of the convection zone. Surface tracking and Doppler effect measurements constrain the exterior flow; solar rotational mode splitting in helioseismology constrains interior stellar rotation through integrals of the form $\delta {\omega }_{n\ell m}\simeq m\int K_{n\ell m}\Omega dV$. Multiple kernels are needed, so surface rotation is not a measurement of every depth.

A useful mean velocity decomposition is $\mathbf{U}=\varpi \Omega {\mathbf{e}}_{\varphi }+{\mathbf{v}}_m$, where $\varpi =r\sin \theta$ is distance from the axis and ${\mathbf{v}}_m$ is stellar meridional circulation. In a frame of constant angular velocity ${\mathbf{\Omega }}_0$, use the relative velocity $\mathbf{V}=\mathbf{U}-{\mathbf{\Omega }}_0\times \mathbf{r}$. Its momentum equation includesThe Coriolis force changes the direction of a parcel's motion but does no work. Thermal driving, release of differential rotational kinetic energy, contraction, or external torques supply the energy of the flows. Solid-body stellar rotation in a suitable hydrostatic equilibrium does not by itself require a meridional flow.

For slowly circulating, approximately axisymmetric material, the leading mechanical balance in an inertial frame is $\nabla p/\rho =-\nabla \Phi +\varpi {\Omega }^2{\mathbf{e}}_{\varpi }$. Taking its curl gives the stellar thermal-wind balanceHere $z$ is the axial coordinate, $S$ is specific entropy and $c_p$ is the specific heat at fixed pressure. The last form uses approximately radial hydrostatic equilibrium and uniform composition. If the material is barotropic, the right-hand side vanishes and angular velocity is constant on cylinders. Small latitudinal entropy differences can instead support the observed noncylindrical differential rotation. Stresses and inertial terms modify this leading balance. The equation relates shear to baroclinicity; it does not give the circulation speed or explain how the shear is maintained.

That second question requires an angular-momentum budget. With $j={\varpi }^2\Omega$ and density-weighted mean quantities, conservation of angular momentum can be written schematically asThe three contributions are Reynolds stress, magnetic Maxwell stress tensor and viscosity transport. Appropriate density weighting replaces the simple correlations in a strongly compressible average. Since $\nabla \cdot (\rho {\mathbf{v}}_m)=0$ in a steady star, the stationary equation gives gyroscopic pumping in a star:A torque that maintains differential rotation therefore drives a meridional response unless another stress cancels it. Conversely the circulation transports angular momentum and feeds back on the rotation. An assumed meridional pattern without its driving torque is not a closed dynamical model.

In the convection zone, outward heat transport drives convection. Rotation deflects the rising and falling parcels and makes their correlations anisotropic. Nondiffusive convective angular momentum transport can drive shear even from an initially uniform rotation rate; treating every Reynolds stress as a positive eddy viscosity would only smooth shear and miss its source. A limiting parcel that conserves $j$ when it moves outwards slows relative to uniform rotation, whereas rapid rotational constraint and correlated horizontal motions change the sign and direction of the net transport. Thus the solar equator-fast pattern requires the actual convective correlations, not a universal parcel argument. The Rossby number $u_c/(2\Omega {\ell }_c)$ compares convective inertia with rotational deflection: deep, slower, larger-scale convection is more rotationally constrained than the fastest photospheric motions. Turnover times ${\ell }_c/u_c$ range from minutes near the visible surface to weeks or longer at greater depth. On these times the Reynolds stresses can redistribute angular momentum and drive circulation, while the star's total spin changes much more slowly.

In a radiative region, rotation distorts equipotential surfaces. Under uniform rotation and uniform composition, leading hydrostatic equilibrium makes pressure, density and temperature constant on each such surface, but effective gravity varies along it. Radiative diffusion then produces a flux $\mathbf{F}=f({\Phi }_{\mathrm{eff}})\nabla {\Phi }_{\mathrm{eff}}$. Its divergence contains $f^{\prime }|\nabla {\Phi }_{\mathrm{eff}}{|}^2$, which varies over a surface and generally cannot match a nuclear source constant there. This is the radiative-equilibrium obstruction in a rotating barotropic star. Thermal imbalance drives Eddington-Sweet circulation, with heat advection enteringIts global thermal estimate isA rotational thermal imbalance of order ${ϵ}_{\Omega }L$ processes the thermal reservoir on this longer time. For the slowly rotating Sun, $t_{\mathrm{KH}}$ is of order $3\times 10^7$ years and ${ϵ}_{\Omega }$ of order $2\times 10^{-5}$, giving $t_{\mathrm{ES}}\sim 10^{12}$ years, much longer than the solar age. Stable composition gradients impede the flow further. Thus classical radiative Eddington-Sweet circulation cannot by itself establish the Sun's interior rotation rapidly. In a rapidly rotating star the same estimate is shorter and rotational mixing can compete with its evolutionary lifetime. Torque-driven convection-zone circulation and radiative thermal circulation have different controlling balances.

Magnetic fields introduce a further coupling. A poloidal magnetic field linking different layers is wound by differential rotation into a toroidal magnetic field, roughly ${\partial }_t{B}_{\varphi }\simeq \varpi {\mathbf{B}}_p\cdot \nabla \Omega$. The resulting Maxwell stress tensor transport angular momentum. Communication along a field occurs on an Alfvén speed crossing time $t_A\sim \ell \sqrt{{\mu }_0\rho }/B$, which can be far shorter than molecular diffusion if a coherent connecting field exists. In steady ideal axisymmetric induction, absence of continuing winding requires ${\mathbf{B}}_p\cdot \nabla \Omega =0$; this is constancy along the field, not automatically solid-body rotation throughout the star. The topology and strength of an interior field matter. Solar dynamo action, especially the interaction of shear with poloidal-field regeneration, couples stellar rotation to the magnetic cycle. Magnetic feedback can change shear and circulation on dynamical or cycle times; an approximately eleven-year activity cycle and twenty-two-year magnetic polarity cycle are distinct from secular spin-down.

A magnetized stellar wind also supplies an external torque. It enforces approximate corotation out to an effective Alfvén radius $R_A$, giving wind-driven magnetic braking of a solar-type star $\dot{J}\sim -{\dot{M}}_w\Omega R_A^2$. Its spin-down time is $t_J\sim I/({\dot{M}}_w{R}_A^2)$, much longer than a convective turnover time. For slowly changing moment of inertia and an unsaturated field scaling that gives $\dot{J}=-K{\Omega }^3$, integration gives ${\Omega }^{-2}={\Omega }_0^{-2}+2Kt/I$, the Skumanich rotation law. A growing or contracting star also changes $I$ and can spin up or down even without torque. Differential spin-down can create shear between layers unless internal transport couples them.

Finally, shear can excite Kelvin-Helmholtz instability, and convection can launch internal gravity waves into the stable radiative interior. Wave damping deposits angular momentum; the net transport depends on excitation and selective absorption and is not fixed by a wave-crossing time alone. Such stresses, magnetic coupling and circulation are candidate agents for the nearly uniform radiative rotation and the narrow tachocline. Molecular kinematic viscosity has a diffusion time ${\ell }^2/\nu$ generally too long to couple the entire solar interior over its age, whereas a turbulent effective kinematic viscosity can act much faster but is not appropriate to every stable layer. A circulation has an advection time $\ell /U_m$, thermal diffusion has time ${\ell }^2/\chi$, and magnetic diffusion has time ${\ell }^2/{\eta }_m$; these must be compared with the local stellar rotation and evolutionary times rather than conflated with them. The solar dynamical time $\sqrt{R^3/(GM)}$ is about half an hour, its rotation period is of order a month, and secular magnetic braking takes an evolutionary time. Rapid force adjustment, convective maintenance of shear, cyclic magnetic feedback and slow global spin evolution are separate processes. Their coupled conservation equations and timescale hierarchy are the basis of a consistent account of solar rotation.

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.

Under a Cartesian frame change $B_i^{\prime }=R_{ik}{B}_k$, where $R$ is an orthogonal matrix. Orthogonality gives $B_k^{\prime }{B}_k^{\prime }=B_k{B}_k$ and $R_{ik}{R}_{j\ell }{\delta }_{k\ell }={\delta }_{ij}$. Thereforewhich is the tensor component transformation law. Also $T_{ij}=T_{ji}$ directly, proving that it is a symmetric Cartesian second-rank tensor for every choice of the two scalars.

For $\mathbf{B}\ne 0$, the direction of $\mathbf{B}$ is an eigenvector with eigenvalue $(\alpha +\beta )|\mathbf{B}{|}^2$. Every vector perpendicular to $\mathbf{B}$ is an eigenvector with eigenvalue $\beta |\mathbf{B}{|}^2$, a two-dimensional eigenspace. If $\alpha =0$, these eigenvalues coincide and every nonzero vector is an eigenvector. If $\mathbf{B}=0$, the tensor is zero and again every nonzero vector is an eigenvector, now with eigenvalue zero.

For the position-dependent field with zero divergence, the product rule givesContraction of the two Levi-Civita symbols in the curl and cross product givesConsequentlyThese constants give the magnetic part of the Maxwell stress tensor, with its physical prefactor suppressed. They are forced if the identity is to hold for every divergence-free field: $\mathbf{B}=(y,0,0)$ forces $2\beta =-1$; at the origin $\mathbf{B}=(y,1,0)$ then forces $\alpha =1$.

Apply the divergence theorem separately to each component:Since the whole vector field vanishes on $S$, all components of $T$ vanish there, and hence the vector volume integral is zero. The argument assumes a differentiable field and a boundary for which the divergence theorem applies.

Use the adiabatic form of ideal magnetohydrodynamics: heat conduction and external heating are absent, so $\rho De/Dt=-p\nabla \cdot \mathbf{u}$. Infinite electrical conductivity and zero viscosity alone would not remove the thermal source terms on the cover page. If those terms are retained, the energy equation below must acquire $\nabla \cdot (\lambda \nabla T)+ϵ$ on its right-hand side.

Expanding the continuity equation gives its conservative form ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$. Using it to expand the conservative momentum derivative givesFor $\nabla \cdot \mathbf{B}=0$, the Lorentz force density is the divergence of the magnetic Maxwell stress tensor:Combining these identities with momentum balance provesThe dyadic stress contains both magnetic tension and magnetic pressure; the sign in the flux is fixed by moving the magnetic force to the left.

For the energy equation, first dot momentum balance with $\mathbf{u}$ and add the adiabatic internal-energy equation. The continuity equation converts their material derivatives into conservative derivatives, while $\mathbf{u}\cdot \nabla p+p\nabla \cdot \mathbf{u}=\nabla \cdot (p\mathbf{u})$. The result isThe ideal magnetohydrodynamic induction equation gives the magnetic counterpart. In fact,Thus the outward Poynting vector is $(\mathbf{B}\times \mathbf{u})\times \mathbf{B}/{\mu }_0$, and magnetic work cancels the fluid work. Finally, the continuity equation impliesAdding all three identities gives ideal magnetohydrodynamic energy conservation in the required local form:whereThe magnetic energy coefficient is $1/(2{\mu }_0)$, as in the original PDF. For a self-consistent, time-dependent gravitational field this local balance should not be confused with a global self-gravitating energy integral: the latter uses gravitational field energy, or the familiar one-half factor in $\int \rho \Phi$ after the field is eliminated.

Work in the shock frame of a thin planar discontinuity, with normal speed $u>0$. A smooth gravitational potential has no finite jump through a vanishingly thin front; its advected contribution cancels using the common mass flux. Integrating the conservation laws gives constant normal mass, momentum and energy fluxes, together with the normal-field and tangential-electric-field conditions of the ideal magnetohydrodynamic shock conditions.

For a parallel magnetohydrodynamic shock, $\mathbf{B}=B_n\mathbf{n}$ and $\mathbf{u}=u\mathbf{n}$. The magnetic solenoidal constraint makes $B_n$ the same on both sides. The normal momentum flux is $\rho u^2+p-B_n^2/(2{\mu }_0)$; its magnetic term is therefore identical upstream and downstream. The Poynting vector is zero because $\mathbf{u}\times \mathbf{B}=0$. Hence the nontrivial jumps are exactly the Rankine-Hugoniot conditions for a perfect gas:with specific enthalpy $h=\gamma p/[(\gamma -1)\rho ]$. Physically, motion along the straight field neither bends it nor compresses its transverse flux. The field consequently provides no changed stress or energy flux in this purely parallel geometry. This conclusion does not apply to a shock that develops a transverse magnetic field.

For a perpendicular magnetohydrodynamic shock, the magnetic field is tangential and the velocity normal. Steady induction makes $uB$ continuous, while the continuity equation makes $\rho u$ continuous. Dividing gives the magnetic flux freezing relationPut $x={\rho }_2/{\rho }_1$, so $u_2=u_1/x$ and $B_2=x{B}_1$. Constant normal momentum flux now reads $\rho u^2+p+B^2/(2{\mu }_0)$. Substitution givesThe magnetic correction is negative for compression: part of the incoming ram pressure builds magnetic pressure instead of thermal pressure.

Here the Poynting vector has normal component $u{B}^2/{\mu }_0$. Divide the constant energy flux by the nonzero mass flux to obtainWith $v_{A1}^2=B_1^2/({\mu }_0{\rho }_1)$, this gives the second expression for downstream pressure:In this expression $x$ is the only downstream unknown. Equating the two pressure expressions and multiplying by $2\gamma x/{\rho }_1$ yieldsup to an overall minus sign, where $c_1^2=\gamma p_1/{\rho }_1$ andThe discarded factor $x=1$ is the unchanged state, not a compressive shock. This derives the polynomial for the compression ratio of a perpendicular magnetohydrodynamic shock without omitting its trivial branch.

For the usual physical adiabatic index $1<\gamma \le 2$, $F$ is strictly increasing for positive $x$, has $F(0)<0$, and has one positive root. Sincethat root is greater than one precisely when $u_1^2>c_1^2+v_{A1}^2$. The necessity is not limited to $\gamma \le 2$: for any $\gamma >1$, eliminating $u_1^2$ using $F(x)=0$ givesFor $x>1$, positive downstream pressure requires $D>0$ when $c_1^2+v_{A1}^2>0$. Moreover,The bracket is positive throughout $1<x<(\gamma +1)/(\gamma -1)$, even when $\gamma >2$. The exactly cold unmagnetized limit gives the familiar strong-shock upper compression ratio directly. Thus the physical compressive branch requiresThe PDF's final $v_1$ denotes this same incoming speed $u_1$. The right side is the squared speed of a perpendicular fast magnetosonic wave: the upstream flow must be super-fast-magnetosonic so that this compressive signal cannot run upstream against the incoming flow. At equality the shock amplitude vanishes into a characteristic disturbance.

Define the perturbation energy per unit mass by $E=⟨v^2/2+B^2/(2{\mu }_0\rho )⟩$. Dot the momentum equation with $\mathbf{v}$ and the induction equation with $\mathbf{B}/({\mu }_0\rho )$, then take volume averages. The Coriolis force does no work. Pressure work and nonlinear advection are divergences and integrate to zero, as does the background $x{\partial }_y$ transport. The two magnetic exchange terms add toThe velocity shear term supplies $2A⟨v_x{v}_y⟩$, while magnetic shear supplies $-2A⟨B_x{B}_y⟩/({\mu }_0\rho )$. Integration by parts gives the viscous and resistive losses as $-\nu ⟨|\nabla \mathbf{v}{|}^2⟩$ and $-\eta ⟨|\nabla \mathbf{B}{|}^2⟩/({\mu }_0\rho )$. For a solenoidal field $\mathbf{F}$, another integration by parts showsAll boundary terms cancel by the shearing-periodic identities. The resulting balance isThe source is work extracted from maintained differential rotation by Reynolds stress and the magnetic part of the Maxwell stress tensor; internal magnetic/kinetic exchange has canceled.

For the decay bound, the elementary inequality $2|ab|\le a^2+b^2$ givesA uniform spectral gap of a shearing-periodic box is needed; ordinary radial periodicity cannot simply be assumed. At each time the compatible Fourier modes areThey are orthogonal over the box. Their wavevectors areFor every nonconstant mode, $|\mathbf{k}|\ge 2\pi /L$, where $L=\max (L_x,L_y,L_z)$: if $m$ or $p$ is nonzero, its corresponding unsheared component supplies the bound; otherwise $n\ne 0$ and the radial component does. Zero mean excludes precisely the constant mode. The Parseval identity therefore proves the time-uniform Poincaré inequalityApply it componentwise to $\mathbf{v}$ and $\mathbf{B}$. With $\nu =\eta$, the energy decay criterion for a viscous-resistive shearing box becomesHenceThis excludes sustained finite-amplitude turbulence as well as infinitesimal growth. It is a sufficient stability condition, not a claim that every box above the threshold must become turbulent. The zero-mean restrictions matter: an imposed uniform magnetic field or bulk epicycle is not controlled by this spatial spectral-gap inequality.

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.

In the rapidly rotating magnetohydrodynamics of a conducting core, the magnetostrophic balance neglects inertia and, initially, viscosity:Take a spherical or suitable axisymmetric impenetrable container, with gravity/buoyancy having no azimuthal component and no applied mechanical torque. Small Rossby number is the slow-motion assumption. It does not imply that an arbitrary magnetic field is compatible with the balance.

Use cylindrical coordinates $(s,\varphi ,z)$ and integrate the azimuthal momentum equation over the fluid portion $C_s$ of a cylinder. The pressure term integrates to zero because it is a periodic $\varphi$ derivative. The Coriolis contribution is proportional to the radial flux ${\int }_{C_s}{u}_{s}sd\varphi dz$. By incompressibility and impermeability it is zero: closing the cylinder with the container's caps and applying the divergence theorem produces no net flux. There is no azimuthal buoyancy term. ConsequentlyThis is Taylor's condition. Equivalently define the torque per radial thickness by $\mathcal{T}(s)=s^2{\int }_{C_s}(J\times B{)}_{\varphi }d\varphi dz$; it must vanish. In a shell, cylinders split by the inner boundary require the corresponding connected-region conditions. These compatibility conditions concern rotating-fluid torque, distinct from force-free magnetic relaxation in a plasma. Their primary origin is Taylor's rotating-fluid analysis.

A magnetostrophic Taylor state allows construction of the flow's ageostrophic part. Taking the curl givesIntegrate along rotation-axis lines and impose impermeability, then use incompressibility and periodicity to reconstruct a particular solution $u_a$. In a sphere the boundary geometry supplies the independent endpoint conditions; a cylinder with parallel planar end normals has additional degeneracies. The homogeneous momentum equation admits any differential cylindrical rotation $u_g=s\omega (s)\hat{\varphi }$. Its Coriolis force is a radial pressure gradient, so momentum balance alone cannot determine $\omega$.

The additional equation is dynamic preservation of the constraint. LetThis is the first variation of the cylindrical Lorentz torque. The resistive induction equation supplies $\dot{B}=\nabla \times (u\times B)+\eta {\nabla }^{2}B$. Define $C_a=\nabla \times (u_a\times B)+\eta {\nabla }^{2}B$ and $L_B[\omega ]=\nabla \times (s\omega \hat{\varphi }\times B)$. Then the geostrophic flow preserving the Taylor constraint solvesFor a given field this is linear in $\omega$; it prescribes the cylindrical rotation needed to keep the field on the Taylor-state manifold. The magnetic boundary conditions must also be imposed on $\dot{B}$. Boundary terms cannot generally be discarded for insulating exteriors; their role in a correct construction is described in Hardy and collaborators' primary analysis.

For example, when magnetic surface torque vanishes, as for compatible $B\cdot n=0$ conditions, the Maxwell stress tensor givesCylindrical rotation induces component changes ${\dot{B}}_s=-\omega {\partial }_{\varphi }{B}_s$ and ${\dot{B}}_{\varphi }=s{B}_s{\omega }^{\prime }-\omega {\partial }_{\varphi }{B}_{\varphi }$. Their product variation integrates to $s{\omega }^{\prime }\int B_s^{2}d\varphi dz$, since the azimuthal derivative integrates to zero. Thus the geostrophic part of the preservation equation reduces toRegularity and compatible boundary conditions fix the nontrivial integration freedom; total angular momentum fixes the solid-body rotation null mode. If $B_s$ vanishes on an entire cylinder, the magnetic coupling coefficient degenerates and generic uniqueness must not be asserted. The simplified expression is conditional on the surface-torque assumption, rather than a general replacement for the boundary-compatible operator equation.

With small positive viscosity the momentum equation adds $\rho \nu {\nabla }^{2}u$. Its cylinder integral becomes the viscous regularization of the Taylor constraintMagnetic torque need no longer vanish exactly; it can be transmitted by viscous stress. Bulk viscosity is order $E=\nu /(2\Omega L^2)$ at fixed velocity. For ordinary no-slip rotating boundaries, thin Ekman layers have thickness of order $L\sqrt{E}$ and shear stress of order $\rho \Omega LU\sqrt{E}$, larger than the bulk order-$E$ correction. Boundary pumping and this friction help select the otherwise undetermined geostrophic flow. In the inviscid interior only impermeability is imposed; no-slip matching is supplied by the thin layer, so simply imposing all viscous boundary conditions on the zero-viscosity equations is incorrect.

For bounded velocities, the viscous torque tends to zero as $E\to 0$, recovering the leading Taylor condition. A prescribed order-one uncompensated magnetic torque would instead require a velocity growing like $E^{-1/2}$ under Ekman friction, or $E^{-1}$ under bulk friction alone. Such growth eventually invalidates the assumed small Rossby number: the limit is singular, not a regular guarantee that every imposed magnetic field has an inviscid solution. With stress-free boundaries the leading no-slip friction is absent, and bulk stresses, magnetic boundary conditions and angular-momentum normalization require their own treatment. Retaining inertia permits torque-driven torsional motion rather than an exact magnetostrophic constraint. Thus zero viscosity gives a magnetic compatibility condition and an induction-based selection of cylindrical rotation, while small viscosity replaces it by a torque balance whose bounded solutions approach compatible Taylor states.

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