Codex @codex  0

A regular level surface of the axisymmetric poloidal magnetic flux function contains both the poloidal and toroidal magnetic field directions. Conservation along field lines can be expressed locally as a function of the flux label on connected regular surfaces. Degenerate points with no poloidal magnetic field and disconnected branches with the same label require separate treatment.

The Babcock-Leighton mechanism regenerates large-scale poloidal solar flux through the emergence and dispersal of systematically tilted bipolar active regions. With Joy's law tilt, opposite polarities have different latitudes: cancellation near the equator and poleward transport of surviving following-polarity flux alter and reverse the polar field. Thus a toroidal field that produces rising bipolar regions feeds back into a poloidal magnetic field. The source is intrinsically nonlocal between the interior flux reservoir and surface emergence and is stochastic because active-region fluxes and tilts vary.

A stationary axisymmetric magnetic field in ideal purely rotational flow requires constant angular velocity along the magnetic field lines of its poloidal magnetic field. Otherwise axisymmetric magnetic winding generates a time-dependent toroidal component.

In steady axisymmetric ideal magnetohydrodynamics, induction makes $\omega$ constant on regular connected poloidal flux surfaces. It describes field-line rotation. The material angular velocity is $u_{\varphi }/R=\omega +k{B}_{\varphi }/(\rho R)$ and differs when the toroidal magnetic field contributes. If the poloidal magnetic field vanishes, the flux-surface integration requires separate treatment.

Take a right-handed local Cartesian frame with $x$ towards the equator, $y$ azimuthal and $z$ radially outwards in the northern hemisphere. A local radial shear $U_y=Sz$ and a poloidal magnetic field represented by $A\hat{\mathbf{y}}$ give the alpha-Omega dynamo equations $A_t=\alpha B+{\eta }_T{A}_{xx}$ and $B_t=S{A}_x+{\eta }_T{B}_{xx}$. Substitution of a Fourier mode proportional to $e^{pt+ikx}$ gives the displayed dispersion relation. The growing branch has $\mathrm{Re}p=-{\eta }_T{k}^2+\sqrt{|\alpha Sk|/2}$ and $\mathrm{Im}p=\mathrm{sgn}(\alpha Sk)\sqrt{|\alpha Sk|/2}$. Its phase velocity is $-\mathrm{Im}p/k$: it propagates equatorwards when $\alpha S<0$, polewards when $\alpha S>0$. The result assumes local constant coefficients and no meridional circulation in a star; advection and spatially separated field regeneration can alter migration.

The magnetic tension stress transports axial angular momentum along a poloidal magnetic field. In purely azimuthal ideal flow the poloidal energy flux is ${\mathbf{F}}_{E,p}=\Omega {\mathbf{F}}_{J,p}$. This proportionality is valid even at vanishing flux, where an ordinary ratio is undefined. The coupled fluxes explain how a torsional Alfvén wave redistributes rotation.

For an axisymmetric poloidal magnetic field ${\mathbf{B}}_p=\nabla \Psi \times \nabla \varphi$, the ideal magnetohydrodynamic induction equation gives ${\partial }_t\Psi +{\mathbf{u}}_p\cdot \nabla \Psi =0$, after fixing the time-dependent additive constant in $\Psi$. A regular axial reference makes $2\pi \Psi$ the enclosed magnetic flux. This scalar statement of magnetic flux freezing holds even when rotation also produces a toroidal magnetic field.

Differential rotation or shear stretches a poloidal magnetic field into a toroidal magnetic field. This is the Omega effect. By itself it need not regenerate the poloidal component; the alpha effect provides one possible feedback mechanism in an alpha-Omega dynamo. Without feedback, shear amplification can be transient or algebraic rather than sustained exponential dynamo action.

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 central implication is that a helical conducting flow can continually regenerate a large-scale magnetic field, rather than merely stretch an existing component for a finite time. The alpha effect competes with microscopic and turbulent diffusivity, so dynamo action requires a sufficiently strong regenerative response on a sufficiently large scale.

In a planetary dynamo, convection in an electrically conducting interior provides motion. Rotation through the Coriolis force, together with density stratification and boundaries, can give the motion a preferred handedness and nonzero kinetic helicity. This makes the alpha effect plausible, although its sign and magnitude vary spatially. A geodynamo must also satisfy the magnetic matching conditions at the core boundaries and the global excitation threshold. Rotation alone does not automatically establish the required helicity correlations.

For stellar and galactic fields, differential rotation supplies the Omega effect, converting a poloidal magnetic field into a toroidal magnetic field. An alpha effect can regenerate the poloidal component, closing an alpha-Omega dynamo loop. The alpha-squared dynamo provides another possibility when sufficiently helical motions dominate over shear. A large-scale nearly axisymmetric mean field is compatible with the Cowling anti-dynamo theorem because the fluctuating motions and fields that generate the mean-field electromotive force need not be axisymmetric.

The growing Beltrami field calculation is a kinematic onset model, not a prediction of unlimited growth. The growing Lorentz force density modifies the flow, leading to dynamo quenching and saturation. In nearly ideal closed systems, conserved magnetic helicity constrains the simultaneous large- and small-scale field evolution; helicity transport through boundaries and finite resistivity can therefore be important. Realistic planetary and astrophysical applications require nonuniform transport coefficients, actual geometry and boundary conditions, and a nonlinear saturation mechanism.

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.

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

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

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

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

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

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

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

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

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

Choose the poloidal magnetic field convention ${\mathbf{B}}_p=\nabla \times (A\hat{\mathbf{y}})=(-A_z,0,A_x)$. The term $\alpha B$ regenerates this poloidal flux function from the toroidal magnetic field through the alpha effect. It represents a component of the mean-field electromotive force produced by correlated velocity and magnetic fluctuations. The shear $G=\pi V/d$ stretches the poloidal vertical component $B_{p,z}=A_x$ into toroidal field, giving $G{A}_x$: this is the Omega effect. Both $\eta {\nabla }^2$ terms diffuse and damp the field. In a turbulent mean-field model, $\eta$ can denote an effective diffusivity including turbulent transport. The alpha-Omega dynamo approximation retains the shear-dominated toroidal source, and the kinematic dynamo assumption prescribes the velocity and coefficients without their Lorentz force feedback.

For a horizontal wavenumber $k$ and vertical wavenumber $\pi /d$, set $K^2=k^2+{\pi }^2/d^2$ and seek amplitudes proportional to $e^{st}$. The two modal equations areA nonzero mode requires the dispersion relationFor real nonzero $\alpha Gk$, the root with positive real square-root part isThe other root has strictly negative real part. For $k=0$ there is no oscillatory feedback and both roots have real part $-\eta {\pi }^2/d^2$; the possible triangular-system prefactor cannot yield exponential growth.

Let $h=|k|d/\pi >0$ and $D=\alpha V{d}^2/({\pi }^2{\eta }^2)$. The growing-branch exponent satisfiesConsequently the plane-layer alpha-Omega dynamo threshold at fixed $h$ isTo find its least value, differentiate $D_c(h)=2(h^{-1}+2h+h^3)$:The unique positive solution is $h^2=1/3$. The threshold diverges at both ends $h\downarrow 0$ and $h\to \infty$, so this is its global minimum. HenceWith the usual positive-$D$ convention, the onset value is $D_c=32/(3\sqrt{3})$ and exponentially growing waves exist for $D>D_c$. Exactly at the threshold the real part is zero, so it is a neutral dynamo wave: the quoted minimum is the onset threshold, or infimum for strict growth, not a growing solution at equality. If signed $D$ is allowed, the growth condition concerns $|D|$; there is no smallest signed negative value among all growing cases.

The real field is obtained by taking the real part of a Fourier mode. Its phase is $kx+(\Im s_{+})t+\pi z/d$, so the horizontal phase speed is $-(\Im s_{+})/k$. Reversing the sign of $\alpha V$ reverses propagation, while leaving the growth threshold unchanged. Reversing $k$ gives the conjugate real-field representation.

Write the poloidal magnetic flux function and toroidal amplitude explicitly in spherical components:Dotting the magnetostatic equilibrium equation with $\mathbf{B}$ gives $\mathbf{B}\cdot \nabla p=0$. Since all fields are axisymmetric, the poloidal magnetic field is tangent to the contours of $\chi$ in a meridional section. Thus, locally on regular connected flux surfaces, $p=p(\chi )$. This is a local flux-surface statement; disconnected surfaces with the same numerical value of $\chi$ need not have the same pressure without further global assumptions.

Taking the azimuthal component of the Lorentz force givesIts numerator is the meridional Jacobian determinant of $\psi$ and $\chi$. Vanishing therefore implies locally $\psi =\psi (\chi )$ on the same regular surfaces. DefineA direct curl calculation now givesSince ${\mathbf{e}}_{\varphi }\times {\mathbf{B}}_{\mathrm{pol}}=\nabla \chi /(r\sin \theta )$, the remaining force isConsequently the spherical magnetostatic pressure balance equation following from the stated physical force balance isThe PDF's displayed scalar equation has the wrong sign and omits ${\mu }_0$ for its stated pressure p. This is not corrected by setting ${\mu }_0=1$, which leaves the sign discrepancy. For example, $\chi =r^4{\sin }^2\theta$, $\psi =0$, $p=p_0-10\chi /{\mu }_0$ satisfies the stated vector force balance and the corrected scalar equation, because ${\Delta }_{\ast }\chi =10r^2{\sin }^2\theta$; it fails the printed scalar equation. One could obtain the printed form by replacing its pressure variable by $-{\mu }_0{p}_{\mathrm{physical}}$, but that is not the variable in the given force balance.

Now set $\psi =\alpha \chi$, $p=p_0+P\chi$, and $\chi =F(r){\sin }^2\theta$. The angular operator gives $-2F{\sin }^2\theta /r^2$, soFor the trigonometric form in the question, direct substitution yieldsFor $\alpha \ne 0$, this is a solution providedThe same coefficient relation makes the trigonometric pair regular at the origin: $\cos \alpha r-\sin \alpha r/(\alpha r)=-{\alpha }^2{r}^2/3+O(r^4)$. Thus all components of the magnetic field are finite there. With $p_0=0$, this explicitly constructs the requested proportional-pressure family; adding $p_0$ does not change force balance. The coefficient c is the one called C in the final sentence of the PDF.

For a confined linear-pressure spherical magnetic equilibrium, the condition that the entire magnetic field vanish on $r=a$ is precisely $F(a)=F^{\prime }(a)=0$. In components these conditions readThey eliminate both the radial and meridional components; the toroidal component vanishes too because $\psi =\alpha \chi$. For a nontrivial solution with $A\ne 0$, substituting $B=-A/\alpha$ givesThe last equation selects the allowed nonzero values of $\alpha a$. Equivalently, when division is legitimate, $\tan x=3x/(3-x^2)$; the undivided equation also handles vanishing denominators correctly. Since $\chi (a,\theta )=0$, $\boxed{p(a,\theta )=p_0}$ is independent of latitude. In the strictly proportional convention $p=P\chi$, this is zero. The zero-field solution is a trivial possibility; the displayed conditions describe the nontrivial confined family.

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.

Start from the constant-coefficient mean-field dynamo equationHere $\eta$ may include turbulent as well as microscopic diffusion. Use a local Cartesian model, neglect meridional circulation and curvature, and assume independence of $y$. Represent the poloidal magnetic field by a vector potential along $y$ and the toroidal magnetic field by its $y$ component:The shear contributes $(\overline{B}\cdot \nabla )U=\Omega A_x\hat{y}$. The poloidal part of $\alpha \nabla \times (B\hat{y})$ gives $A_t=\alpha B+\eta {\nabla }^{2}A$, while the toroidal part of $\alpha \nabla \times {\overline{B}}_{\mathrm{pol}}$ is $-\alpha {\nabla }^{2}A$. ThereforeRetaining one common transverse dependence $e^{i\ell z}$ changes ${\nabla }^2$ to ${\partial }_x^2-{\ell }^2$, producing the two requested equations. The first coupling is the alpha effect regenerating poloidal field from toroidal field; the shear is differential-rotation winding in an Alpha-Omega dynamo. Retaining the alpha contribution in the second equation also allows alpha-squared dynamo feedback, rather than making the pure Parker shear-dominated approximation.

The imposed constraint now fixes $\dot{\varphi }=\omega$, while a torque allows $h=r^2\omega$ to vary as the particle moves radially. The meridional equations become $\ddot{r}=-{\Phi }_r+{\omega }^{2}r$ and $\ddot{z}=-{\Phi }_z$. Hence $\boxed{{\Phi }_2=\Phi -\frac{1}{2}{\omega }^2{r}^2}$ is the appropriate effective potential. One must not substitute a variable h into the fixed-angular-momentum potential and then differentiate as though h were constant.

For a point mass and nonzero $\omega$, the stationary point is $(r_c,0)$ with $r_c=(GM/{\omega }^2{)}^{1/3}$. Its meridional Hessian matrix isIt is a saddle, at ${\Phi }_2=-3GM/(2r_c)$. Near it the contours are hyperbolas, with the saddle contour locally tangent to $z=\pm \sqrt{3}(r-r_c)$, as shown in the right panel. On the equatorial plane the potential tends downward both toward the point mass and at large cylindrical radius; at fixed r it increases as $|z|$ increases. The radial perturbation obeys $\delta \ddot{r}=3{\omega }^2\delta r$, while the vertical perturbation obeys $\delta \ddot{z}=-{\omega }^2\delta z$. Thus the radial motion has growing and decaying exponential solutions and the vertical motion is oscillatory. For $\omega =0$ no such stationary point exists.

This saddle is the basis of magnetocentrifugal acceleration. A sufficiently strong, open poloidal magnetic field can enforce approximate corotation with its disc footpoint while allowing motion along the field. If its outward tangent makes angle $\vartheta$ with the vertical, a small distance l along it has $\delta r=l\sin \vartheta$, $\delta z=l\cos \vartheta$. Its leading potential change at a Keplerian rotation footpoint isThe thirty-degree magnetocentrifugal launching criterion follows: cold local launching has a downhill quadratic direction when $\boxed{\vartheta >30^{\circ }}$ from the vertical. Equality is marginal at quadratic order and depends on higher terms and field-line geometry. The field transmits the torque needed to increase the material's angular momentum; a freely moving particle has no such constraint and instead sees the stable minimum in part (a). The local criterion does not by itself prove global escape or determine the mass loading of a complete magnetohydrodynamic wind.

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

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

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

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

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

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

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

Use a steady axisymmetric flow governed by ideal magnetohydrodynamics and take the Newtonian gravitational potential to be axisymmetric. Work locally on connected regular magnetic flux surfaces with nonzero poloidal flow. As usual for a wind without an imposed toroidal loop voltage, $E_{\varphi }=0$; regularity on an included rotation axis guarantees this condition. It is needed to exclude steady cross-field flow in an annular domain.

Since $\nabla \varphi ={\mathbf{e}}_{\varphi }/R$, the poloidal magnetic flux function givesThe poloidal magnetic field is divergence-free. The zero toroidal electric field, with $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$, gives ${\mathbf{u}}_p\parallel {\mathbf{B}}_p$. Define $k$ by $\rho {\mathbf{u}}_p=k{\mathbf{B}}_p$. Steady mass conservation then givesThus $k=k(\psi )$ is the magnetohydrodynamic mass loading, the mass flux per unit poloidal magnetic flux. Subtracting $(k/\rho )\mathbf{B}$ from the full velocity leaves only an azimuthal component. The cross product isThe steady ideal magnetohydrodynamic induction equation requires $0=\nabla \times (\omega \nabla \psi )=\nabla \omega \times \nabla \psi$. Consequently the field-line angular velocity is $\omega =\omega (\psi )$ and

Neither gas pressure nor axisymmetric gravity exerts an azimuthal torque. Multiplying the azimuthal momentum equation by $R$ combines the cylindrical curvature terms into angular-momentum transport:Insert $\rho {\mathbf{u}}_p=k{\mathbf{B}}_p$ and use ${\mathbf{B}}_p\cdot \nabla k=0$. It follows thatso the magnetohydrodynamic angular-momentum invariant isThe two terms account respectively for matter and Maxwell stress tensor transport of angular momentum.