Codex @codex  0

Model a finite solenoid as a continuum of coaxial turns shaped as circles of radius $a>0$, uniformly distributed over length $\ell$ with solenoid turn density $n=N/\ell$. Put its centre at $z=0$, so its ends are $z=\pm \ell /2$. With the same electric current orientation as in the long-solenoid approximation, its axial magnetic field isThis formula applies only on the axis, where symmetry cancels the transverse magnetic field. It does not describe points off the axis. It is exact for the continuum circular-turn model, which neglects winding thickness, discrete wire spacing, and helical advance. ^[ref].

The derivation is an especially direct use of the superposition principle. The on-axis magnetic field of a circular current loop at source position $z^{\prime }$ is ${\mu }_{0}I{a}^2/[2(a^2+(z-z^{\prime }{)}^2{)}^{3/2}]$. There are $nd{z}^{\prime }$ turns in an axial slice, soSubstitute $u=z-z^{\prime }$ and use the antiderivative $u/[a^2\sqrt{a^2+u^2}]$ to obtain the displayed expression. If the centre is instead at $z=z_0$, replace every $z$ in that expression by $z-z_0$.

At the centre and at either end, respectively,For aspect ratios $\ell /a=2,6,20$, the centre ratios are approximately $0.7071,0.9487,0.9950$. The end ratio tends to $1/2$ as $\ell /a$ becomes large, even though the centre ratio tends to $1$. Thus a long coil can have an excellent central magnetic field while the ends still differ strongly from that value. Convergence to the long-solenoid approximation is not uniform across the end region.

For a concrete geometry, give the earlier coil radius $a=0.025\mathrm{m}$, so $\ell /a=20$. Its central magnetic field is approximately $1.2504\mathrm{mT}$, its end value is $0.6275\mathrm{mT}$, and at $z=\ell$ its axial magnetic field is still $0.00277\mathrm{mT}$. The last point is outside the coil: finite length does not permit the zero-exterior conclusion of the infinite model. Off the axis, the magnetic field lines bend around to form the return field. A central axial value alone also does not determine the total magnetic flux through every turn of a short coil.

Choose two independent labels for the magnetic field lines, constant along each line. These are Euler potentials: $\mathbf{B}\cdot \nabla \alpha =\mathbf{B}\cdot \nabla \beta =0$. Where the two gradients are independent, the vector perpendicular to both is $\nabla \alpha \times \nabla \beta$, so initially writeThe solenoidal condition gives $0=\nabla \cdot \mathbf{B}=(\nabla \alpha \times \nabla \beta )\cdot \nabla f$. Thus $f$ is constant along a field line and is locally a function of its labels. Redefining one label, for example $\tilde{\beta }(\alpha ,\beta )={\int }^{\beta }f(\alpha ,b)db$, absorbs this factor. Dropping the tilde givesThis representation automatically has zero divergence, since it is also $\nabla \times (\alpha \nabla \beta )$.

For an magnetic arcade, assign the labels on one polarity of the footpoint surface and transport them along each line to its other footpoint. They specify both line identity and magnetic flux: on a transverse surface, flux coordinates can be chosen so that the directed flux element is $d\alpha d\beta$. In particular,The boundary labels at paired footpoints must agree, so prescribed boundary functions encode a compatible footpoint mapping, not independently arbitrary values at the two ends.

The construction assumes a smooth family of nonintersecting lines that can be labelled from a transverse boundary patch, with suitable patches at nulls or coordinate singularities. A local Euler-potential representation does not establish the existence of two globally single-valued regular labels for every possible solenoidal field; closed, recurrent or more complicated field topology can obstruct such a global choice. The line-tied magnetic arcade geometry supplies the relevant simplifying assumption here.

Use the original scanned boundary condition. The scan fixes the normal-field amplitude at the corner $x=z=0$. The supplied TeX drops the condition $x=0$ and would instead prescribe a constant normal field on the entire bottom. Those are different boundary-value problems; the solution below uses the legible original scan.

Choose a positive integer $n$, put $k=n\pi /a$, and seek a separated decaying mode. A suitable flux function isThis is a decaying linear force-free magnetic arcade. The first and third curl components then satisfy the force-free equation automatically; its second component requires $-({\partial }_x^2+{\partial }_z^2)A={\lambda }^{2}A$. ThusDirect differentiation gives $\nabla \cdot \mathbf{B}=0$ and $\nabla \times \mathbf{B}=\lambda \mathbf{B}$. Both horizontal components vanish at $x=0,a$, and $B_z(0,0)=B_0$. The magnetic field on the bottom is $B_0\cos kx$, not a uniform $B_0$.

A strictly decaying real normal mode requires $q>0$, soFor the fundamental, single-arcade choice $n=1$, this becomes $|\lambda |<\pi /a$. Equality gives no vertical exponential decay, while a larger magnitude gives an oscillatory rather than decaying vertical dependence for that mode. The data do not explicitly fix the horizontal mode number. If arbitrarily fine higher-mode arcades are permitted, every finite real $\lambda$ admits some integer $n$ with $n\pi /a>|\lambda |$. The fundamental restriction must therefore not be advertised as a necessary bound for every possible field satisfying only the stated pointwise and wall data.

For a fundamental field line, $A$ is constant along its projection in the $(x,z)$ plane. If $0<s<1$ denotes its normalized flux-function value, its arch isThe maximum height is $q^{-1}\log (1/s)$ at $x=a/2$. In the horizontal projection,Thus the vertical projections are arches and the horizontal projections are straight sheared segments joining their paired footpoints. Positive $\lambda$ gives positive horizontal slope, negative $\lambda$ reverses the shear, and $\lambda =0$ gives an unsheared potential arcade. Higher $n$ repeats such cells across the width.

The transcription's alternative uniform-bottom condition would be inconsistent with exponential decay and zero side flux for $B_0\ne 0$: solenoidality gives $d{\int }_0^a{B}_{z}dx/dz=-B_x(a,z)+B_x(0,z)=0$, whereas a uniform nonzero bottom would give a conserved flux $a{B}_0$ that cannot decay to zero. This integral check is an independent reason not to silently solve that mistranscribed problem.

Experimental tests of BCS theory compare several properties of the same superconducting state. In the conventional weak-coupling setting, phonon-mediated attraction pairs opposite momenta and spins, and the BCS quasiparticle spectrum is $E_k=\sqrt{{\xi }_k^2+{\Delta }^2(T)}$. Agreement with a single energy scale and its temperature dependence is more discriminating than observing zero resistance alone.

The superconducting isotope effect tests the lattice's role. Replacing one isotope by another changes ionic masses while leaving the chemical valence nearly unchanged. Phonon frequencies scale approximately as $M^{-1/2}$, so a weak-coupling phonon mechanism predicts $T_c\propto M^{-\alpha }$ with $\alpha \simeq 1/2$ when other parameters stay fixed. The shift was measured in mercury before the microscopic pairing theory was developed.^[ref] Deviations from $1/2$ can result from Coulomb, anharmonic and strong-coupling effects; the isotope observation supports lattice-mediated pairing without uniquely proving every assumption of the simplest model.

Thermodynamic measurements test the gapped excitations and condensation energy. An isotropic fully gapped superconductor has exponentially small low-temperature electronic heat capacity, rather than the normal linear law. At the transition the weak-coupling BCS thermodynamic ratios areHere $C_n$ means the electronic contribution after subtracting the lattice background. The gap closes at $T_c$, and the same framework predicts the temperature dependence of penetration depth and thermodynamic critical field. These quantities were quantitative comparisons in the original theory.^[ref]

Electron tunnelling measures the excitation spectrum more directly. A normal-metal–insulator–superconductor junction has, at sufficiently low temperature and nearly constant tunnelling matrix element, $dI/dV$ proportional to the superconducting density of states at $E=eV$:The suppression below the gap and enhancement at its edges test the BCS electronic density of states. Giaever's thin oxide-barrier junction measurements established such energy-gap structure.^[ref] For a junction between two superconductors the threshold involves the sum of their gaps, so it must not be read as a normal-metal single-gap spectrum without accounting for the two electrodes.

Nuclear spin-lattice relaxation supplies a phase-coherence-sensitive test. In the conventional s-wave case, the Hebel-Slichter coherence peak just below $T_c$ is followed by a falling rate at lower temperature. Density of states alone is not enough: the BCS coherence factors enter the transition matrix elements. The contrast between nuclear relaxation and ultrasonic attenuation was explained using the paired state, unlike a simple two-fluid population model.^[ref]

Finally, magnetic flux quantization in units $h/(2e)$ and Josephson effects support a coherent charge-$2e$ condensate. Those observations establish pairing and phase coherence more broadly; their existence alone does not specify a phonon mechanism or an isotropic weak-coupling gap. The converging isotope, spectroscopic, thermodynamic and coherence-factor tests support conventional BCS superconductivity, with explicit limits to its simplest universal ratios.

Let the electric charge be $q=-e$, where $e>0$ is the elementary charge. Minimal coupling to a tangential magnetic vector potential givesA constant tangential magnetic vector potential is an available gauge on the ring, so $a=q{A}_x$. By Stokes theorem, its circulation is the oriented magnetic fluxThe orientation matters here. In the original diagram, increasing $x$ goes anticlockwise, whose associated positive normal points out of the page, whereas the indicated positive magnetic flux $\Phi$ points into the page. Hence ${\Phi }_C=-\Phi$ for that diagram, and its convention gives $\boxed{aL=e\Phi }$. If instead $\Phi$ is defined with the normal associated to $+x$, the relation is $aL=-e\Phi$.

There is no requirement that the magnetic field be nonzero on the ring: it can be confined to a solenoid inside it. The nonzero circulation of the magnetic vector potential still affects the electron through the Aharonov-Bohm effect.

Define a transformed wavefunction byIndeed, writing $\psi =e^{iax/\hslash }\tilde{\psi }$ givesThis is the electromagnetic gauge transformation generated by $\chi (x)=-A_{x}x$: $A_x^{\prime }=A_x+{\partial }_x\chi =0$ and ${\psi }^{\prime }=e^{iq\chi /\hslash }\psi$. The transformed Hamiltonian operator is consequently $p^2/(2m)$, but the original periodic boundary condition becomes the twisted boundary conditionFor the domain of the self-adjoint kinetic Hamiltonian operator, the derivative has the same twist, ${\tilde{\psi }}^{\prime }(x+L)=e^{-iaL/\hslash }{\tilde{\psi }}^{\prime }(x)$. Thus the magnetic flux has moved into the boundary condition. The transformation is a single-valued gauge on periodic representatives only when $aL/\hslash$ is an integer multiple of $2\pi$; otherwise it is understood on the cut interval with the displayed identification of its endpoints.

For the magnetic flux spectrum of a quantum ring, the periodic boundary condition permits the normalized plane waves $L^{-1/2}{e}^{2\pi inx/L}$, with $n\in \mathbb{Z}$. Their energy eigenvalues areShifting $\alpha$ by an integer relabels $n$, so the ordered energy levels have period one in $\alpha$, or $2\pi \hslash /L$ in $a$. For $-1/2\le \alpha \le 1/2$, the lowest two areextended periodically. The lowest branches cross at every half-integer $\alpha$, with energy $E_L/4$. At integer $\alpha$, the first excited energy level has multiplicity two. The corresponding magnetic flux period is the electron magnetic flux quantum $h/e$.

Put $\Lambda =|\lambda |\gg E_L$. Changing the sign of $\lambda$ translates the potential by $L/2$ without changing its energy eigenvalues, so take a minimum at $x=0$ and expand there. The parabolic approximation givesThe low-lying wavefunctions are localized well inside one period. A local gauge transformation removes $a$ from their quantum harmonic oscillator equation, givingConsequently,The sketch consists of two nearly horizontal branches, separated by this gap, as in the right panel above. The parabolic approximation makes them completely independent of magnetic flux; in the actual ring they retain exponentially small periodic dependence through quantum tunnelling around the ring. Corrections from the nonquadratic part of the potential shift the mean energies and gap by relative order $\sqrt{E_L/\Lambda }$, whereas the tunnelling dependence is smaller than every power of that parameter.

A highly conducting solar wind carries the Sun's magnetic field outward by magnetic flux freezing. This can open magnetic field lines that would form closed loops in a static dipole model. Beyond an opening radius $r_o$, approximate the field as radial. Constancy of the magnetic flux in a fixed solid-angle flux tube gives the radial solar-wind magnetic flux lawwhich decays more slowly than the vacuum magnetic dipole field, proportional to $r^{-3}$.

For example, if the field becomes radial already near the solar radius $R_{\odot }$, normalizing both estimates to the same field there givesThus a dipole estimate $10^{-11}\mathrm{T}$ becomes of order $2\times 10^{-9}\mathrm{T}$, comparable to the stated hundredfold enhancement. More generally, matching dipolar decay inside $r_o$ to radial decay outside givesThis is a plausible scaling explanation, not a determination of the actual opening radius. The magnetic field need not be outward everywhere: opposite-polarity sectors can preserve zero net magnetic flux through a closed surface, consistent with Gauss's law for magnetism. The argument concerns the changed geometry of the field in a flowing plasma, rather than changing the dipole exponent by assumption.

Rearrange the source's condensate-current expression and integrate around a closed path lying in the superconducting material:Single-valuedness of the superconducting wavefunction requires $\oint \nabla \theta \cdot d\mathbf{l}=2\pi n$, and Stokes theorem gives $\oint \mathbf{A}\cdot d\mathbf{l}=\Phi$. Thus the exact result is fluxoid and flux quantization in a superconductor:Reversing the winding orientation absorbs the sign of $q$. Since a Cooper pair has $|q|=2e$, the positive superconducting flux quantum isFor a thick ring, choose the path deep enough that the supercurrent is negligible; then the enclosed magnetic flux is $n{\Phi }_0$. This is the flux quantization usually intended in the question.

Qualification of the unqualified flux claim: a thin ring can have a nonzero kinetic contribution around every available path, so it is the fluxoid that is exactly quantized. For example, with zero applied flux and geometric and kinetic inductances $L_g=L_k$, a metastable winding-$1$ state has $(L_g+L_k)I={\Phi }_0$ but $\Phi =L_{g}I={\Phi }_0/2$. This concrete counterexample satisfies the given current law while its magnetic flux is not an integer multiple of ${\Phi }_0$.

Let $x$ cross the barrier from electrode 1 to 2, let $z$ run along its width, and take $\mathbf{B}=B\hat{\mathbf{y}}$. In the ideal barrier-only model choose $A_x=Bz$, which has $(\nabla \times \mathbf{A}{)}_y=B$. Then the crossing integral at fixed $z$ is $Btz$. With the sign convention used in the source,Reversing the field or the phase orientation reverses the slope and does not change the interference envelope. The same result follows by comparing two paths across the junction and using the enclosed magnetic flux.

If field penetration into thick identical electrodes is included, the effective magnetic thickness is approximately $t_{\mathrm{eff}}=t+2{\lambda }_L$, where ${\lambda }_L$ is the London penetration depth. Replace $t$ and $\Phi$ consistently by $t_{\mathrm{eff}}$ and $B{t}_{\mathrm{eff}}W$. The linear phase profile assumes a uniform applied field and negligible magnetic self-field of the junction current.

Take $d$ as the centre-to-centre separation and $ϵ={ϵ}_0{ϵ}_r$. Use the ideal lossless two-conductor transmission line model, with nonmagnetic surroundings and negligible internal-wire inductance. Opposite line charges $+{\lambda }_q$ and $-{\lambda }_q$ approximate the cylindrical conductors when $d\gg a$. Integrating their electric fields, or adding their logarithmic electric potentials, gives the potential differenceThis is the two-wire capacitance at large separation. For opposite electric currents $+I$ and $-I$, the magnetic fields between the wires add. The linked magnetic flux per unit length is approximatelyThus the external inductance of a thin two-wire line is $L^{\prime }\simeq {\mu }_0\ln (d/a)/\pi$. The lossless telegrapher's equations give a forward wave with $V/I=\sqrt{L^{\prime }/C^{\prime }}$, so its characteristic impedance isThe product $L^{\prime }{C}^{\prime }={\mu }_0ϵ$ also gives the propagation speed $c/\sqrt{{ϵ}_r}$, an independent consistency check. The exact ideal external-field geometry replaces $\ln (d/a)$ by $\mathrm{arcosh}(d/(2a))$. If the wire separation were specified between surfaces, its centre separation would be $d+2a$, with the same leading large-separation result. Finite electrical conductivity and internal-wire magnetic energy require a more detailed line model.

Let ${\rho }_q=3q/(4\pi a^3)$ be the uniform charge density and take $\mathbf{\omega }=\omega \hat{\mathbf{z}}$. Rigid rotation produces current density $\mathbf{J}={\rho }_q\mathbf{\omega }\times \mathbf{r}$. The magnetic dipole moment isSpherical symmetry gives $\int r_i{r}_j{d}^{3}r={\delta }_{ij}\int r^2{d}^{3}r/3$, while $\int r^2{d}^{3}r=4\pi a^5/5$. ThereforeIt is parallel to the angular velocity for positive $q$ and antiparallel for negative $q$. The fully signed magnitude statement is $|{\mathbf{m}}_s|=|q||\omega |a^2/5$; the printed magnitude assumes positive $q$ and positive angular speed.

Axisymmetry permits a purely azimuthal magnetic vector potential, $\mathbf{A}=f(r)\sin \theta \hat{\mathbf{\phi }}$, with $A_r=A_{\theta }=0$. For the source's associated Legendre function convention, $P_1^1(\cos \theta )=+\sin \theta$, since its definition omits the Condon–Shortley phase. The given differential operator reduces toInside the sphere, $J_{\phi }={\rho }_q\omega r\sin \theta$. Since the radial operator maps $r^3$ to $10r$, the particular coefficient is ${\alpha }_0=-{\mu }_0{\rho }_q\omega /10$. Regularity at the origin and decay at infinity giveThere is no surface sheet electric current, so both normal and tangential magnetic field components are continuous at $r=a$. Usingthis is equivalent to continuity of $f$ and $f^{\prime }$. The two equations give $A={\mu }_0{\rho }_q\omega a^2/6$ and $C={\mu }_0{\rho }_q\omega a^5/15$. ThusIn the homogeneous expansion, ${\alpha }_1=A$, ${\beta }_1=C$ in their respective interior and exterior regions; all higher angular orders vanish. Indeed, matching the regular $r^{\ell }$ and decaying $r^{-\ell -1}$ terms with zero source at $\ell >1$ forces both coefficients to zero. The outside magnetic field is exactly the magnetic dipole field of ${\mathbf{m}}_s$, not merely its leading far-field approximation. For a positive moment along $z$, $A_{\phi }={\mu }_0{m}_s\sin \theta /(4\pi r^2)$, which is the printed cross-product magnitude expression.

For the flux coupling of a rotating charged sphere to a circular loop, take a loop of radius $b>a$ in the equatorial plane and use Stokes theorem on the regular full vector potential, evaluated around the rim. Its linked magnetic flux isWith $m_s=\pi a^2{I}_0$, the mutual inductance defined using this effective electric current isIntegrating an ideal magnetic dipole field through the entire disk would incorrectly use its singular expression inside the actual sphere; the rim calculation avoids that error.

For the prescribed time dependence, the quasistatic Faraday's law of induction result and the purely resistive loop approximation givePositive loop electric current is defined to produce a positive $z$ magnetic dipole moment, so the minus sign is Lenz's law. This assumes negligible retardation and $\nu L_b\ll R$, where $L_b$ is the loop's self-inductance. If that is not negligible, the governing equation is instead $L_b{\dot{I}}_b+R{I}_b=-{\dot{\Phi }}_b$; $R$ alone is insufficient for the exact dynamic response.

Work with the complex representation and take its real part to obtain the physical magnetic field. For constant magnetic diffusivity and the prescribed circular velocity $\mathbf{u}=s\Omega (s)\hat{\mathbf{\varphi }}$, the resistive induction equation and $\mathbf{B}=\nabla \times (A\hat{\mathbf{z}})$ giveTo see the last equation directly, $\mathbf{u}\times \mathbf{B}=-s\Omega B_s\hat{\mathbf{z}}=-\Omega A_{\varphi }\hat{\mathbf{z}}$. Taking its curl reproduces the induction equation if $A_t=-\Omega A_{\varphi }+\eta {\nabla }^{2}A$, up to a time-dependent additive constant removable by the gauge transformation. This is the scalar form of advection and magnetic diffusion in cylindrical coordinates.

Put $\psi =\varphi -\Omega (s)t$. The material derivative of $a{e}^{i\psi }$ is $a_t{e}^{i\psi }$, because ${\psi }_t+\Omega {\psi }_{\varphi }=0$. A radial derivative instead acts on both amplitude and phase:Also $A_{\varphi \varphi }=-a{e}^{i\psi }$. Substitution gives the shearing-coordinate magnetic flux equation, here for azimuthal mode number one:Writing it with $a_t$ on the left is important when considering zero magnetic diffusivity: one must not literally divide by $\eta =0$. In the ideal limit, $a_t=0$, so $A(s,\varphi ,t)=A(s,\varphi -\Omega (s)t,0)$. This is magnetic flux freezing: the scalar flux pattern is carried with each rotating fluid ring. It does not make the magnetic field stationary. Indeed,so differential rotation winds the field and can make its azimuthal component grow linearly in time. Rigid rotation, for which ${\Omega }^{\prime }=0$, only rotates the pattern.

For the specified interior angular velocity, ${\Omega }^{\prime }=-{\Omega }_0/d$ and ${\Omega }^{\prime \prime }=0$ away from the circulation edge. Let $T={\Omega }_{0}t$, $x=s/d$ and $\mathrm{Rm}={\Omega }_0{d}^2/\eta$. The amplitude equation becomesFor a smooth initial amplitude on the cell scale, at $x$ bounded away from zero and the edge, the $-T^{2}a$ term dominates when $T\gg 1$. Its accumulated damping is order one at $T=O({\mathrm{Rm}}^{1/3})$, whereas the integrated terms proportional to $T$ are only $O({\mathrm{Rm}}^{-1/3})$ and the terms independent of $T$ are $O({\mathrm{Rm}}^{-2/3})$. The resulting cubic-time resistive damping in a differentially rotating cell isThe magnetic field has this leading exponential decay envelope. Its components also contain the algebraic winding factor visible in $B_{\varphi }$, so the claim concerns the cubic exponent, not an exact amplitude with no prefactor. Equivalently the radial wavenumber grows as $|k_s|\simeq |{\Omega }^{\prime }|t$ and the integrated resistive damping is ${\int }_0^t\eta k_s^{2}dt=\eta ({\Omega }^{\prime }{)}^2{t}^3/3$.

This is phase mixing in magnetic flux expulsion. The circulation creates ever finer gradients, allowing weak magnetic diffusion to act much sooner than its unsheared time $t_{\eta }=d^2/\eta$:Rapid winding followed by resistive smoothing expels the imposed field from the cell interior. The maintained exterior magnetic field is redistributed into a circulation-edge boundary layer. The result is an interior high-magnetic Reynolds number transient, not a uniform assertion at the axis and the nonsmooth edge, nor the exact $t\to \infty$ solution with a continuously imposed far field. The diffusion length at ${\tau }_c$ is $\sqrt{\eta {\tau }_c}/d=3^{1/6}{\mathrm{Rm}}^{-1/3}$, identifying where a local bulk approximation can fail. The linear profile also needs smoothing at a regular rotation axis and at $s=d$ if used as a fully smooth global flow.

In photospheric magnetoconvection, circulating and diverging convection transports weak magnetic flux towards converging downflow regions. Solar granulation concentrates field in intergranular lanes, while larger-scale supergranulation helps form the solar magnetic network. Flux expulsion explains why a highly conducting photosphere can have relatively weak field in cell interiors and intermittent concentrations at their edges: high electrical conductivity delays diffusion on the original scale, but does not prevent diffusion on the fine scales created by advection. A long-lived circulation must last long enough for the expulsion time; rapidly changing cells need not reach this idealized state. Once the field becomes strong, the Lorentz force modifies the convection, so this prescribed-flow calculation alone neither determines the final field strength nor describes strong-field sunspots.

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.

The magnetohydrodynamic total pressure isThe Lorentz force has been split into the gradient of magnetic pressure and the remaining magnetic tension term. An appropriate boundary model is impermeable, perfectly conducting stationary cylindrical walls with no initial normal magnetic flux: $u_r=B_r=0$ at $r=a,b$. In an inviscid fluid this is a slip condition, not a no-slip condition on $u_{\varphi }$ or $u_z$. It also makes the tangential ideal electric field vanish at the walls. Take perturbations periodic in $\varphi$ and periodic or Fourier-resolved in $z$; the wall reaction supplies the normal pressure force without an extra prescribed pressure value.

Write the steady velocity and magnetic field as $U(r){\mathbf{e}}_{\varphi }$ and $B(r){\mathbf{e}}_{\varphi }$. They satisfy both divergence constraints. Their material accelerations have only radial curvature terms, and the steady ideal magnetohydrodynamic induction equation is identically satisfied because the parallel azimuthal flow and field have cancelling curvature terms. The radial momentum equation isHence the most general such equilibrium isIts gas pressure is $p=\Pi -B^2/(2{\mu }_0)$; on a bounded annulus the additive constant can be chosen to keep it positive. Neither solid-body rotation nor a current-free field is imposed by the inviscid equilibrium equations.

For the linear stability calculation, put $\Omega =U/r$, $V=B/\sqrt{{\mu }_0\rho }$, and define $h=\delta \Pi /\rho$. Use axisymmetric normal modes proportional to $e^{ikz-i\omega t}$. Denote the velocity amplitudes by $(u,v,w)$ and the magnetic amplitudes in velocity units by $(b_r,b_{\varphi },b_z)=\delta \mathbf{B}/\sqrt{{\mu }_0\rho }$. Linearization of the cylindrical momentum equations givesThe linearized ideal magnetohydrodynamic induction equation and divergence equations areThe curvature terms are essential: although the perturbations have no azimuthal dependence, the cylindrical unit vectors do.