Codex @codex  0

For an axisymmetric magnetic field, the solenoidal condition becomes ${\partial }_R(R{B}_R)+{\partial }_z(R{B}_z)=0$. Locally introduce a poloidal magnetic flux function $\Psi$ byThe divergence condition is exactly the integrability condition for these two derivatives. The additive constant in $\Psi$ fixes the reference flux. For a field regular on the axis choose $\Psi (0,z,t)=0$. Then the magnetic flux through a horizontal disk of radius $R$ is $2\pi \Psi (R,z,t)$, and the flux between two such surfaces is $2\pi \Delta \Psi$.

Since ${\mathbf{B}}_p\cdot \nabla \Psi =0$ and ${\partial }_{\varphi }\Psi =0$, the level surfaces of $\Psi$ contain the complete magnetic field, including its toroidal magnetic field component. Thus regular levels are poloidal flux surfaces. Their meridional contours trace the projected field lines. The representation is local on regular regions; disconnected flux surfaces with the same label need not share additional physical data.

The azimuthal component of $\mathbf{u}\times \mathbf{B}$ is $u_z{B}_R-u_R{B}_z=-{\mathbf{u}}_p\cdot \nabla \Psi /R$. The ideal magnetohydrodynamic induction equation, equivalently Faraday's law with $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$, gives ${\Psi }_t=R(\mathbf{u}\times \mathbf{B}{)}_{\varphi }$ up to a time-dependent reference constant. Fixing that reference gives the material advection of an axisymmetric magnetic flux function:Axisymmetry removes the azimuthal contribution to the material derivative. This is magnetic flux freezing expressed using one scalar flux label.

The azimuthal component itself is generally not an advected scalar. Including the cylindrical basis-vector derivatives in the induction equation givesUsing mass conservation, $D\rho /Dt=-\rho \nabla \cdot \mathbf{u}$, and $DR/Dt=u_R$, this becomes the toroidal magnetic-field winding equation:Thus $B_{\varphi }/(\rho R)$ is conserved following the flow if the angular velocity is constant along the poloidal magnetic field, including the case of no rotation. Otherwise differential rotation winds poloidal field into toroidal field.

There is also an explicit flux-conservation form,Integrating this over a meridional region converts its change of toroidal flux to boundary terms. Hence such a flux is conserved when the indicated boundary transport vanishes. More generally magnetic flux freezing conserves flux through every material surface, but a meridional coordinate surface need not remain material when the fluid rotates. Toroidal flux conservation therefore does not imply that $B_{\varphi }$ itself is materially constant.

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.

Let $L$ be a core-scale length, $B_{\ast }$ a representative strength of the prescribed poloidal magnetic field, and $v_A=B_{\ast }/\sqrt{{\mu }_0\rho }$ the Alfvén speed. A consistent weakly damped, small-amplitude ordering isThe first ratio is the Lehnert number. It makes the wave slow relative to rotation and allows almost geostrophic flow, with the azimuthal velocity approximately constant along a rotation-axis cylinder. The small wave amplitude gives $|B_{\varphi }|\sim \sqrt{{\mu }_0\rho }|U|\ll B_{\ast }$. Poloidal advection then changes neither the background field nor the wave substantially over one wave period. If there are no-slip boundary layers, also require their damping rate $\sqrt{\nu \Omega }/L\ll \sigma$, or otherwise retain their torque; a small bulk viscosity alone does not remove that boundary torque.

In the rotating magnetohydrodynamic momentum equation, axisymmetry removes the azimuthal pressure derivative. The neglected nonlinear azimuthal advection has size $|{\mathbf{u}}_p|U/L$, and the neglected viscous term has size $\nu U/L^2$. The retained azimuthal Coriolis force is $2\Omega u_r$. A typical compatible local balance has $|u_r|\sim \sigma U/\Omega$, small compared with $U$ when the Lehnert number is small. The azimuthal Lorentz force density can be written using $\nabla \cdot {\mathbf{B}}_p=0$:The physical prefactor is $1/({\mu }_0\rho )$: the extra $i$ printed inside the first reduced momentum equation is a typographical error, also contradicted by the subsequent real-valued system.

In the resistive induction equation, neglecting poloidal advection of the perturbation and magnetic diffusion leaves the winding of the prescribed poloidal magnetic field by differential rotation. ThusThis is the small-amplitude torsional Alfvén wave approximation, with a time-independent background and no external azimuthal forcing.

We now average over the full height of a rotation-axis cylinder. The sphere is impenetrable. Integrating incompressible flow over the volume inside any such cylinder gives zero flux through its cylindrical side, since there is no flux through its spherical caps. ThereforeThis is why the Coriolis force disappears from the cylinder average; it need not vanish locally. An axisymmetric field in the insulating exterior is a potential field with $B_{\varphi }=0$, so continuity gives $B_{\varphi }=0$ on the sphere. Consequently both the end flux terms and the moving-endpoint terms in the magnetic-stress integral vanish. Expanding the divergence givesSince $U$ is independent of $z$, its integrated time derivative is $2z_b{U}_t$. Defining $T={\int }_{-z_b}^{z_b}{B}_r{B}_{\varphi }dz$ givesThe stationary background and the fact that $U/r$ depends only on $r,t$ also givesoTogether these form the cylinder-averaged torsional Alfvén wave system.

Assume $H>0$ on the cylinders carrying the wave, with regular endpoint limits. If $H$ vanishes on an interior cylinder, then $B_r$ and $T$ vanish there and that cylinder has no radial magnetic coupling; the displayed $H^{-1}$ energy term must be interpreted on the coupled region rather than literally dividing by zero. Since $H$ is time independent, differentiation of the proposed conserved energy of a cylinder-averaged torsional wave givesRegularity at the rotation axis makes the lower endpoint vanish, while vanishing magnetic torque at the collapsed outer cylinder makes the upper endpoint vanish. This provesThe magnetic term has a direct interpretation. Decompose $B_{\varphi }=(T/H)B_r+B_{\varphi }^{\perp }$, with $\int B_r{B}_{\varphi }^{\perp }dz=0$. Then $\int B_{\varphi }^{2}dz=T^2/H+\int (B_{\varphi }^{\perp }{)}^{2}dz$, and the induction equation gives $(B_{\varphi }^{\perp }{)}_t=0$. Thus, apart from a common dimensional factor and a stationary orthogonal magnetic component, the conserved quantity is the physical kinetic plus magnetic energy.

For the frequency estimate, put $V=U/r$ and eliminate $T$:Locally this is a wave equation with speedIt is the cylinder-averaged radial field component that matters, not simply the dipole measured at the Earth's surface. Taking a cylinder height of order $6\times 10^6\mathrm{m}$, $\rho \sim 10^4\mathrm{kg}{\mathrm{m}}^{-3}$, and an assumed internal cylindrical radial field $B_{r,\mathrm{rms}}\sim 0.2$--$2\mathrm{mT}$ givesFor a standing wave of radial length $L_c\sim 2$--$3.5\times 10^6\mathrm{m}$, the cyclic frequency is of order $c_T/(2L_c)$ and the angular frequency is of order $\pi c_T/L_c$. HenceGeometry and the assumed internal field can shift the estimate substantially. As a concrete calibration, the 2010 torsional-wave study inferred a cylindrical radial rms field near $2\mathrm{mT}$ and a roughly four-year crossing time, with an observed six-year recurrence. A much weaker assumed interior field instead gives the older multi-decadal estimate. The paper supplies no numerical interior field strength, so a unique numerical frequency cannot be deduced without such an assumption.

An anti-dynamo theorem identifies a restriction that prevents a conducting fluid from regenerating its magnetic field against resistive losses. The hypotheses matter: throughout the following arguments take positive uniform magnetic diffusivity, smooth fields, no imposed magnetic source, and the stated isolated or energy-closed boundary conditions. An anti-dynamo theorem does not rule out transient magnetic amplification, winding of an imposed field, or an externally supplied magnetic flux. The resistive induction equation is the common starting point:

First consider Backus' necessary condition for dynamo action. For a conductor contained in a sphere of radius $a$, include the potential field in the insulating exterior when defining the total magnetic energy:For an incompressible flow, with boundary conditions eliminating the mechanical surface term, such as a no-slip boundary condition for the fluid, integration by parts giveswhere $e_{ij}=({\partial }_i{u}_j+{\partial }_j{u}_i)/2$ is the rate-of-strain tensor. Its antisymmetric part does no stretching, since $B_i{B}_j$ is symmetric. If $S$ bounds the largest eigenvalue of $e$ throughout the flow, the first term is at most $2S{E}_B$.

The magnetic free-decay spectral bound for the enclosing insulating sphere isThe exterior energy is essential to this variational estimate. Its constant comes from the slowest freely decaying dipolar poloidal mode: the spherical radial matching condition reduces to $j_0(\kappa a)=0$, whose first positive root is $\kappa a=\pi$. Expanding an admissible field in the orthogonal free-decay modes, or using their variational characterization, gives the bound. A smaller conductor can only increase the minimum dissipation quotient for the enclosing-sphere comparison. Combining the two estimates yieldsThusBelow this threshold the Gronwall inequality gives exponential energy decay. Above it the estimate merely permits growth; it does not prove that the geometry regenerates the field. A time-dependent version bounds the field-amplitude exponent by the time average of the maximum strain minus the diffusive rate.

The geometric obstruction in Cowling anti-dynamo theorem is different: an isolated exactly axisymmetric magnetic field cannot be self-maintained. To see the missing regeneration mechanism, write its poloidal magnetic field in terms of the poloidal magnetic flux function $\psi$:The poloidal part of the resistive induction equation becomesThere is no term containing $B_{\varphi }$ that creates new poloidal flux. Fix the gauge by setting $\psi =0$ on the regular symmetry axis, and match to a decaying exterior field with no imposed flux. At a positive interior maximum away from the axis, the advective term vanishes and the diffusive term is nonpositive. The parabolic comparison principle therefore prevents increase of the maximum; it likewise prevents decrease of a negative minimum. In the steady problem the strong maximum principle, regularity at the axis and decay outside exclude a nonzero maximum altogether. Matching to the insulating exterior does not create a maximum at the interface: the exterior flux is harmonic in the corresponding axisymmetric operator, and the two normal-derivative conditions exclude such an interface extremum. The same homogeneous scalar diffusion problem gives decay in time on a bounded conductor with no replenishing boundary data.

Once the poloidal field has decayed, let $b=B_{\varphi }/r$ and let ${\Omega }_f=u_{\varphi }/r$ denote the fluid's angular velocity. Its equation isDifferential rotation can make toroidal field from poloidal field, but it cannot close the reverse loop. When ${\mathbf{B}}_p=0$, this too is a homogeneous advection-diffusion equation, with a regular axis and zero toroidal field at an insulating exterior. Hence it cannot sustain the remaining field. The exactly axisymmetric field, rather than merely an axisymmetric velocity, is excluded. If an exactly axisymmetric field were maintained by a nonaxisymmetric velocity, averaging the linear induction equation in azimuth would replace that velocity by its azimuthal mean and give the same contradiction. An axisymmetric flow can, however, excite a nonaxisymmetric field. Also an axisymmetric averaged field in mean-field dynamo theory has an additional fluctuating electromotive force, so this unmodified scalar proof does not apply to it.

The toroidal-velocity anti-dynamo theorem restricts the flow instead. In a conducting sphere, suppose $\mathbf{u}\cdot \mathbf{x}=0$, so the incompressible flow is tangent to every concentric sphere. No assumption of axisymmetry is made. Define $P=\mathbf{x}\cdot \mathbf{B}$. A direct Cartesian calculation is particularly useful:The last equality uses $\nabla \cdot \mathbf{B}=0$, so ${\nabla }^2(\mathbf{x}\cdot \mathbf{B})=\mathbf{x}\cdot {\nabla }^2\mathbf{B}$. Thus the radial magnetic scalar is a passive diffusing scalar, not a regenerated field. Multiplying by $P$, integrating over the sphere and using tangential incompressible advection givesIn an insulating exterior, a spherical harmonic component of degree $l\ge 1$ has $P\propto r^{-l-1}$, so its boundary contribution is negative: ${\partial }_{n}P=-(l+1)P/a$. The monopole component is absent because the magnetic field is solenoidal and regular. The resulting dissipative scalar problem makes $P$, and hence the poloidal field, decay.

When that radial field is absent, the remaining tangential solenoidal field can be written $\mathbf{B}=\nabla \times (T\mathbf{x})$. Its induction equation reduces to a second passive scalar equation. Indeedand the curl of the last gradient term is zero. Therefore $T_t+\mathbf{u}\cdot \nabla T=\eta {\nabla }^{2}T$, up to an irrelevant radial gauge function, which can be removed by taking the angular mean of $T$ to be zero. Its angular variations vanish at the insulating spherical boundary, so this scalar also decays by its quadratic energy identity. Equivalently, a putative growing magnetic normal mode must first have $P=0$, after which the scalar $T$ equation rules it out. Purely toroidal spherical flow cannot close a dynamo regeneration loop, even if it is nonaxisymmetric. This does not assert that purely poloidal flow is always incapable of a dynamo, or that any azimuthally directed cylindrical flow is governed by the same spherical theorem.

Other restrictions illuminate the same distinction between stretching and regeneration. The planar anti-dynamo theorem has a similar scalar proof in the usual two-dimensional geometry. For planar velocity and magnetic field independent of the third coordinate, write the in-plane field as $\nabla \times (A\hat{\mathbf{z}})$. Its vector potential obeys $A_t+\mathbf{u}\cdot \nabla A=\eta \Delta A$ after fixing the gauge, so $\frac{1}{2}\frac{d}{dt}\int A^2=-\eta \int |\nabla A{|}^2$ under periodic or homogeneous energy-closed boundary conditions. A third field component independent of that coordinate is also a passive scalar. A spatially uniform imposed field may survive in a periodic box but is not regenerated dynamo field. Three-component motions depending on only two coordinates are not covered by this strictly planar formulation.

There are also necessary bounds in terms of velocity amplitude rather than strain. The magnetic work term is bounded by $\Vert \mathbf{u}{\Vert }_{\infty }\Vert \mathbf{B}{\Vert }_2\Vert \nabla \times \mathbf{B}{\Vert }_2$; combining this with the same free-decay inequality gives a velocity-based lower threshold on the magnetic Reynolds number. These bounds, like Backus' necessary condition for dynamo action, are necessary rather than sufficient. Pure motion that only rearranges a scalar magnetic potential cannot evade diffusion merely by becoming faster. Conversely the stagnation-flow energy amplification in Question 2 does not contradict the bounded-body theorems: it uses an unbounded three-dimensional strain field and, for its growing finite-energy modes, a slowly decaying field at infinity. Similarly, the torsional Alfvén waves in Question 3 transfer energy between velocity and a toroidal perturbation about a prescribed background; they are not a construction of that background by dynamo action. The anti-dynamo results identify failed regeneration mechanisms and the assumptions under which those failures are rigorous.

The flux representation gives $B_R=-{\psi }_z/R$ and $B_z={\psi }_R/R$. Therefore $\nabla \cdot \mathbf{B}=0$ and $\mathbf{B}\cdot \nabla \psi =0$. With magnetohydrodynamic mass loading $k=k(\psi )$,Thus the steady continuity equation is satisfied.

The cross product in the ideal magnetohydrodynamic induction equation simplifies toHence steady induction requires $\mathbf{B}\cdot \nabla (v/R)=0$. On a regular connected poloidal flux surface, this integrates locally toThis is the field-line angular velocity of an axisymmetric wind. It need not equal the gas angular velocity when $B_{\varphi }\ne 0$, because $u_{\varphi }=R\omega +k{B}_{\varphi }/\rho$. The regular-surface qualification matters: if $\psi$ is constant and the field is purely toroidal, induction alone cannot turn an arbitrary $v(R,z)/R$ into a function of $\psi$.

Put $\mathbf{F}=(\nabla \times \mathbf{B})\times \mathbf{B}$ and $E=|\mathbf{u}{|}^2/2+\Phi$. Dot the pressure-free momentum equation with $\mathbf{u}$. Since $\mathbf{B}\cdot \mathbf{F}=0$, axisymmetry givesThe azimuthal momentum equation in cylindrical coordinates similarly givesSubtracting $\omega$ times the second relation from the first, and using $\mathbf{B}\cdot \nabla \omega =0$, proves that $E-\omega R{u}_{\varphi }$ is constant along a flowing regular poloidal flux surface. For nonzero mass loading its corotating energy invariant of an axisymmetric magnetic wind isThis is the invariant labeled $ϵ$ in the question. It is a corotating mechanical energy, distinct from the total magnetohydrodynamic Bernoulli invariant, which includes electromagnetic energy flux. In standard notation it equals total specific energy minus $\omega$ times the total specific angular-momentum invariant.

Near the footpoint, negligible $k{B}_{\varphi }/\rho$ gives $u_{\varphi }\simeq R\omega$, while the poloidal velocity is parallel to the poloidal magnetic field. The gas approximately corotates and slides along a rotating field line as a bead on a wire, with effective potentialA decreasing effective potential converts corotating potential energy into poloidal kinetic energy. This is the mechanism of magnetocentrifugal acceleration; local launching does not alone guarantee global escape.

For an attractive point mass, take $\Phi =-GM/(R^2+z^2{)}^{1/2}$ and ${\omega }^2=GM/R_0^3$. Parameterize the outward straight line by $R=R_0+s\sin \alpha$, $z=s\cos \alpha$. The first derivative of ${\Phi }_{\mathrm{eff}}$ vanishes at $s=0$ because the footpoint is in circular radial balance. The meridional Hessian there is ${\omega }^2\mathrm{diag}(-3,1)$, soThus the thirty-degree magnetocentrifugal launching criterion in its strict negative-curvature form isAbove this angle the equilibrium is unstable to an outward displacement; below it the near-footpoint potential initially rises. The source's proportionality constant must be negative to describe attraction.

The marginal straight-line launch at thirty degrees has zero quadratic curvature. For precisely the straight line and point-mass potential here, let $x=s/R_0$ and $a=\sin \alpha$. Direct expansion givesAt $\alpha =\pi /6$, this is $-7x^3/16+O(x^4)$, so the outward direction is downhill at cubic order. Consequently the strict inequality is the nondegenerate quadratic criterion; it should not be interpreted as excluding every one-sided nonlinear marginal launch at equality. An exactly resting bead remains an equilibrium until displaced.

In a steady flow the ideal magnetohydrodynamic induction equation gives $\nabla \times (\mathbf{u}\times \mathbf{B})=0$. Since both vectors have only poloidal magnetic field components, write $\mathbf{u}\times \mathbf{B}=F_{\varphi }{\mathbf{e}}_{\varphi }$. Axisymmetry impliesFor a field regular on the axis, $C=0$; equivalently one may impose zero toroidal electromotive force. Under this physical condition, $\mathbf{u}$ is parallel to $\mathbf{B}$. Write $\rho \mathbf{u}=k\mathbf{B}$. Continuity and $\nabla \cdot \mathbf{B}=0$ then give $\mathbf{B}\cdot \nabla k=0$. The poloidal magnetic flux function has $\mathbf{B}\cdot \nabla \psi =0$, so locally on each connected magnetic surface,This regular-axis alignment of steady poloidal ideal flow needs its zero-circulation condition: axisymmetry alone does not imply this alignment on a domain excluding the axis. An explicit counterexample is available in a cylindrical annulus. Take nonzero constants $U,b$, positive ${\rho }_c,c_s$, andHere $\psi =-bz$, $\nabla \cdot \mathbf{B}=\nabla \times \mathbf{B}=0$ and $\mathbf{u}\times \mathbf{B}=Ub{\mathbf{e}}_{\varphi }/R$ has zero curl. The flow satisfies continuity, $p^{\prime }=-\rho {\Phi }^{\prime }$, and ${\nabla }^2\Phi =4\pi G\rho$. Its acceleration and Lorentz force density vanish, so it solves the steady equations while $\mathbf{u}$ and $\mathbf{B}$ are perpendicular. The nonzero toroidal circulation and the axial singularity explain why this counterexample is excluded from a regular polar accretion column. The remaining derivation uses the regular, aligned branch.

For a narrow stream tube, integrating continuity gives constant $A\rho |\mathbf{u}|$. On a nonzero-flow tube the magnetohydrodynamic mass loading $k$ is constant, so $A\rho |\mathbf{u}|=A|k||\mathbf{B}|$. Therefore $A|\mathbf{B}|$ is constant and $A\propto |\mathbf{B}{|}^{-1}$. This also follows directly from conserved magnetic flux through a flux tube.

For constant isothermal sound speed, $\nabla p/\rho =c_s^2\nabla \ln (\rho /{\rho }_0)$, with arbitrary reference density ${\rho }_0$. Project momentum along $\mathbf{B}$. The Lorentz force density has no component in that direction, and alignment implies $\mathbf{B}\cdot (\mathbf{u}\cdot \nabla )\mathbf{u}=\mathbf{B}\cdot \nabla (|\mathbf{u}{|}^2/2)$. Thus the isothermal magnetic Bernoulli integral isChanging ${\rho }_0$ only shifts the Bernoulli function by a constant. Since $k$ is constant along the magnetic field,Differentiating the isothermal magnetic Bernoulli integral therefore gives

In the slender polar accretion column, the dipolar flux-tube area is $A(z)=A_{\ast }(z/R_{\ast }{)}^3$. The decreasing axial field is a leading approximation, not an exactly solenoidal field $(0,0,B_z(z))$ throughout a cylinder. Indeed $\nabla \cdot \mathbf{B}={\partial }_z{B}_z\ne 0$ for that literal field. Near the axis a small radial component $B_R\simeq -(R/2){\partial }_z{B}_z=3R{B}_z/(2z)$ supplies the required radial divergence. Its magnitude is smaller by $R/z$, even though its divergence is leading order. Keeping this expanding flux tube while neglecting transverse forces yields the intended one-dimensional model.

Use the positive inward speed $v=-u_z$; the signed sonic velocity is $u_z=-c_s$. Conservation of mass and the Bernoulli equation giveEliminating ${\rho }^{\prime }$ yields the isothermal dipolar accretion equation,A smooth transonic branch requires both sides to vanish at its sonic point, soThis is the critical point of the flow equation; a generic subsonic solution need not cross it. A crossing in the exterior column requires $z_s\ge R_{\ast }$, and a sonic point strictly outside the star requires $z_s>R_{\ast }$. Differentiating the equation at the crossing gives $2(v_s^{\prime }{)}^2=3c_s^2/z_s^2$. For inward accretion that accelerates toward the star, $v_s^{\prime }=-\sqrt{3/2}{c}_s/z_s$.

The reservoir boundary condition fixes $ϵ=c_s^2\ln ({\rho }_{\infty }/{\rho }_0)$. At the sonic point, ${\Phi }_s=-3c_s^2$ and $v_s^2/2=c_s^2/2$, givingFinally use the sonic mass flux through the dipolar flux-tube area:This rate belongs to the smooth transonic branch of the idealized column, with $A_{\ast }$ the supplied total loaded area. The reservoir condition alone does not force every steady solution onto this branch. A real narrow dipolar column must match an outer flow where the slender approximation ceases to apply.

Use axisymmetry and $\nabla \cdot \mathbf{B}=0$. Since the purely azimuthal velocity has zero divergence, the ideal magnetohydrodynamic induction equation becomes ${\partial }_t\mathbf{B}=(\mathbf{B}\cdot \nabla )\mathbf{u}-(\mathbf{u}\cdot \nabla )\mathbf{B}$. Differentiation of the cylindrical unit vectors contributes to its azimuthal component, givingThe azimuthal component of the magnetic tension force isThere is no azimuthal pressure or gravitational force and no azimuthal advective acceleration for this motion. Because the poloidal magnetic field is divergence-free,These coupled induction and tension equations describe a torsional Alfvén wave.

Use the poloidal magnetic flux function convention ${\mathbf{B}}_p=R^{-1}\nabla \Psi \times {\mathbf{e}}_{\varphi }$, so $2\pi \Psi$ is magnetic flux up to a reference constant. Its poloidal magnetic field components areThe prescribed surface flux therefore givesThe constant offset ${\Psi }_1$ does not affect either field component. If instead $\Psi$ denotes the full physical flux rather than flux divided by $2\pi$, both component formulas acquire a common $1/(2\pi )$ factor; the subsequent logarithmic derivative relation is unchanged.

The Gauss's law for magnetism constraint is $R^{-1}{\partial }_R(R{B}_R)+{\partial }_z{B}_z=0$. Since $\alpha$ is constant with radius on the surface, it gives the power-law poloidal field near a disk surface relationFor a smooth one-sided field above the surface, $B_z(R,z)=B_z(R,0)[1-(\delta -1)\tan \alpha z/R]+O(z^2)$. Thus the field initially decreases with height if $\delta >1$, increases if $0<\delta <1$, and has zero first height derivative if $\delta =1$.