Codex @codex  0

Current helicity correlates a magnetic field with its curl, which is proportional to electric current density. It differs from magnetic helicity, which uses a magnetic vector potential. For a solenoidal monochromatic field of wavenumber $k$ in a periodic domain, choose ${\mathbf{a}}_b=\nabla \times \mathbf{b}/k^2$. Then $\nabla \times {\mathbf{a}}_b=\mathbf{b}$ and $H_c=k^2⟨{\mathbf{a}}_b\cdot \mathbf{b}⟩$.

Dynamo quenching is the reduction or modification of magnetic-field generation by magnetic feedback on the generating flow and its correlations. The Lorentz force can change shear or turbulence, while magnetic helicity affects mean-field regeneration. A phenomenological scalar example is $\alpha (B)={\alpha }_0/(1+B^2/B_{\ast }^2)$. In the illustrative amplitude equation $\dot{B}=[c\alpha (B)-d]B$ with $c,d>0$, a growing field saturates at $B^2=B_{\ast }^2(c{\alpha }_0/d-1)$ when $c{\alpha }_0>d$. This demonstrates the negative feedback mechanism without claiming that algebraic quenching is a complete solar closure.

An open flux tube carries magnetic flux through its end caps. Its ordinary magnetic helicity can therefore change under a single-valued gauge transformation. For $\mathbf{B}=aR{\mathbf{e}}_{\varphi }+b{\mathbf{e}}_z$ inside $R<R_0$ and zero outside, a regular family of magnetic vector potentials is $A_{\varphi }=bR/2$, $A_z=C-a{R}^2/2$ inside and $A_{\varphi }=b{R}_0^2/(2R)$, $A_z=C-a{R}_0^2/2$ outside. Their tangential components agree at $R_0$. Inside, $\mathbf{A}\cdot \mathbf{B}=bC$, so the helicity per unit length is $\pi bC{R}_0^2$. Changing $C$ is the gauge transformation $\Lambda =\Delta Cz$, and changes this result whenever $b\ne 0$. The convention $A_z=0$ outside gives $C=a{R}_0^2/2$ and helicity per unit length $\pi ab{R}_0^4/2$, while $C=0$ gives zero. Without an end-cap or gauge convention the ordinary helicity is not unique.

In a bounded regular domain, fix the global extension or circulation convention for the magnetic vector potential. A bounded inverse of curl on tangent solenoidal vector fields gives $\Vert \mathbf{A}{\Vert }_2\le C_D\Vert \mathbf{B}{\Vert }_2$. The Cauchy-Schwarz inequality then gives $|H_M|\le C_D\Vert \mathbf{B}{\Vert }_2^2=2{\mu }_0{C}_D{E}_B$. One construction extends the tangent field by zero outside the domain and uses the whole-space Biot-Savart law; its kernel of size $|\mathbf{x}-\mathbf{y}{|}^{-2}$ is integrable over a bounded region, giving the required $L^2$ bound. Nonzero conserved magnetic helicity therefore prevents complete decay of magnetic energy.

The magnetic helicity integrates the dot product of a magnetic vector potential and its magnetic field. It measures linkage and twisting of magnetic flux. For a solenoidal vector field $\mathbf{B}$ tangent to the boundary, it is unchanged by the gauge transformation $\mathbf{A}\mapsto \mathbf{A}+\nabla \chi$, since the change is ${\int }_V\nabla \cdot (\chi \mathbf{B})dV={\int }_{\partial V}\chi \mathbf{B}\cdot d\mathbf{S}=0$. If its density satisfies ${\partial }_{t}h=\nabla \cdot (f\mathbf{B}-h\mathbf{u})$ with both vector fields tangent to the boundary of a fixed volume, the divergence theorem gives conservation of its integral. This is an instance of helicity conservation by tangent boundary conditions.

For regular globally defined electromagnetic potentials, use $\mathbf{E}=-\nabla \varphi -{\partial }_t\mathbf{A}$, $\mathbf{B}=\nabla \times \mathbf{A}$ and the homogeneous Maxwell equations to prove the displayed identity. With vanishing spatial surface flux, its spacetime integral is minus one half the change of magnetic helicity. Zero integral therefore needs matching temporal helicity boundaries, not merely localized electric fields and magnetic fields. The Gaussian potential $a(t)e^{-r^2}(-y,x,\lambda )$ has helicity $2\lambda (\pi /2{)}^{3/2}a(t{)}^2$ and supplies an explicit spatially localized counterexample when $a$ changes from zero to one.

Magnetic relaxation converts magnetic energy into material motion through the Lorentz force density, and dissipates the resulting kinetic energy through viscosity. In a perfect conductor, magnetic flux freezing preserves field-line connectivity throughout every smooth stage. A regular limiting state has zero velocity and is a magnetostatic equilibrium, but its magnetic field can contain current sheets. Conservation of one global magnetic helicity is weaker than preservation of the complete frozen field-line structure.

For smooth ideal magnetohydrodynamics, write $\mathbf{B}=\nabla \times \mathbf{A}$ and ${\partial }_t\mathbf{A}=\mathbf{u}\times \mathbf{B}-\nabla \phi$. The continuity equation and induction equation give $D_t(\mathbf{B}/\rho )=(\mathbf{B}/\rho )\cdot \nabla \mathbf{u}$. The vector-potential equation becomes $D_t\mathbf{A}=\nabla (\mathbf{u}\cdot \mathbf{A}-\phi )-(\nabla \mathbf{u}{)}^T\mathbf{A}$. Contracting these two equations cancels their velocity-gradient terms and provesConsequently the magnetic helicity in a material volume satisfies $d{H}_m/dt={\int }_{\partial V}(\mathbf{u}\cdot \mathbf{A}-\phi )\mathbf{B}\cdot \mathbf{n}dS$. It is conserved when the magnetic field is tangent to the boundary. A single-valued gauge transformation changes $H_m$ by ${\int }_{\partial V}\Lambda \mathbf{B}\cdot \mathbf{n}dS$, which also vanishes under that condition.

Choose a magnetic vector potential with $\mathbf{B}=\nabla \times \mathbf{A}$. The resistive induction equation permits the evolution equationA gauge transformation changes $H={\int }_V\mathbf{A}\cdot \mathbf{B}dV$ by a surface integral of the normal magnetic field, which is zero here. Thus the integrated magnetic helicity has gauge invariance under the stated boundary conditions.

Differentiate both factors in $H$. Use the stated vanishing-boundary assumption to remove the helicity surface flux; this can be implemented in the usual closed-boundary gauge $\mathbf{n}\times \mathbf{A}=0$, or by extending the localized field and integrating over all space with sufficient decay. The curl integration by parts identity then givesThe motional term is zero pointwise because $(\mathbf{u}\times \mathbf{B})\cdot \mathbf{B}=0$. The scalar-potential term integrates to zero by $\nabla \cdot \mathbf{B}=0$ and the boundary condition. The resistive magnetic-helicity balance is thereforeUnlike the resistive magnetic energy contribution, this signed integral need not have a fixed sign.

A magnetic field varying over a length $\ell$ has $|\nabla \times \mathbf{B}|$ of order $B/\ell$. Consequently, the resistive magnetic energy loss contains two inverse powers of $\ell$, while the resistive magnetic-helicity balance contains only one. Small-scale gradients can therefore dissipate substantial magnetic energy while retaining an approximately conserved large-scale magnetic helicity.

A dimensionally consistent version is the selective decay bound for magnetic helicity. With $D_B={\eta }_m{\int }_V|\nabla \times \mathbf{B}{|}^{2}dV/{\mu }_0$, the Cauchy-Schwarz inequality givesFor bounded $E_B$ and bounded, possibly nonvanishing $D_B$ as ${\eta }_m\to 0$, the resistive change in magnetic helicity tends to zero, even though resistive magnetic energy loss may remain finite. Magnetic energy is preferentially dissipated, while magnetic helicity is approximately conserved in this bounded-energy, bounded-dissipation limit.

This is not a comparison of quantities with different physical units or a universal ordering of fractional decay rates. For a single Beltrami field satisfying $\nabla \times \mathbf{B}=k\mathbf{B}$, both integrated quantities have fractional decay rate $2{\eta }_m{k}^2$ without fluid motion. The boundedness assumptions and distribution of magnetic helicity across spatial scales matter.

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.

Choose a magnetic vector potential $\mathbf{A}$ with $\nabla \times \mathbf{A}=\mathbf{B}$, and define the total magnetic helicity byThe value is unchanged by a single-valued gauge transformation $\mathbf{A}\mapsto \mathbf{A}+\nabla \chi$, becauseHere both the solenoidal magnetic-field constraint and the tangent boundary field are essential. In a multiply connected domain, choose the vector potential by extending the tangent field by zero to all space and fixing its decay there, or fix the equivalent circulation convention. This avoids silently varying the non-gradient harmonic ambiguity of a vector potential on the restricted domain.

The ideal magnetohydrodynamic induction equation permitsDifferentiating $H_M$ and using integration by parts for curl,The volume product vanishes algebraically, the gauge flux vanishes since $\mathbf{B}\cdot \mathbf{n}=0$, and the last boundary term vanishes since $\mathbf{v}=0$. Thus magnetic helicity is conserved even though viscosity dissipates mechanical energy.

Through magnetic flux freezing, magnetic field lines move with the fluid and retain their linkage during every smooth ideal motion. For two thin closed flux tubes carrying fluxes ${\Phi }_1,{\Phi }_2$, their mutual contribution is $2L{\Phi }_1{\Phi }_2$, where $L$ is their linking number; self-contributions measure twisting and writhing of the tubes. This explains the topological content of magnetic helicity. Conservation of this one scalar is a consequence of the frozen structure, not a complete description of all that structure.

Let $K$ be the kinetic energy, $M$ the magnetic energy and $E=K+M$:Dot the magnetohydrodynamic momentum equation with $\mathbf{v}$ and integrate. Incompressibility and the zero boundary velocity remove the advective and pressure work. The viscous term integrates to $-\rho \nu {\int }_D|\nabla \mathbf{v}{|}^2$. ThusThe ideal magnetohydrodynamic induction equation gives, with no boundary contribution,The Lorentz force density therefore transfers energy between the field and the fluid without creating total energy. Since $\nabla \cdot \mathbf{v}=0$ and $\mathbf{v}=0$ on the boundary, another integration by parts gives $\int |\nabla \mathbf{v}{|}^2=\int |{\mathbf{\omega }}_v{|}^2$, where ${\mathbf{\omega }}_v=\nabla \times \mathbf{v}$ is vorticity. The viscous magnetic-relaxation energy identity isThe conserved nonzero magnetic helicity prevents all the field energy from disappearing. Fix a bounded inverse-curl convention for the magnetic vector potential, so that $\Vert \mathbf{A}{\Vert }_2\le C_D\Vert \mathbf{B}{\Vert }_2$. The Cauchy-Schwarz inequality gives the magnetic-energy lower bound from helicity,Consequently $E$ decreases to a finite positive limit andUnder the intended regular long-time interpretation of smooth vorticity, the dissipation cannot keep producing peaks of fixed height and arbitrarily short duration. More precisely, if $d(t)=\Vert {\mathbf{\omega }}_v(t){\Vert }_2^2$ is uniformly continuous, uniformly continuous integrable dissipation tends to zero proves $d(t)\to 0$. The Poincaré inequality and the identity above then imply $\Vert \mathbf{v}(t){\Vert }_2\to 0$.

To identify the asymptotic state, suppose the regular relaxation has a limit, or work on a compact invariant limiting set on which the energy is continuous. On such a set the energy is constant, so the displayed identity forces ${\mathbf{\omega }}_v=0$, and the zero boundary velocity forces $\mathbf{v}=0$. The induction equation is then stationary, and the momentum equation becomesThis is magnetostatic equilibrium, with nonzero field allowed and required here by conserved magnetic helicity. It is not necessarily a force-free magnetic field: a pressure gradient can balance the magnetic force. If the field develops a current sheet, equilibrium is understood with the corresponding interface force balance.

There is a mathematical qualification to the smoothness premise. Smoothness at each finite time alone does not prove uniform long-time bounds, uniform continuity of dissipation or convergence to one limiting field. The energy calculation proves finite total dissipation unconditionally for a smooth solution; the settling conclusion is the regular-relaxation argument just given. A rigorous global convergence theorem requires that additional long-time control. In particular, smooth vorticity does not require smooth electric current density in the limiting state.

The relevant competition is between energy reduction and magnetic flux freezing. Ideal evolution transports the entire connectivity of the magnetic field lines, including the separation and linkage of flux regions. Viscous dissipation can remove kinetic energy but cannot permit magnetic reconnection between such regions. As each region relaxes toward lower magnetic energy, neighboring flux regions can approach with different limiting tangential fields. A smooth transition can become increasingly thin; its field remains bounded while its curl, and therefore its electric current density, becomes large.

The limiting object is a current sheet. If $\mathbf{n}$ points from its minus to its plus side, its surface electric current density isFor a tangential discontinuity, ${\mathbf{B}}_{\pm }\cdot \mathbf{n}=0$. The normal balance in magnetostatic equilibrium requiresThus a finite jump of tangential field can coexist with equilibrium: the jump of magnetic pressure is balanced by the ordinary pressure jump.

This need not happen for every initial field; a field already in smooth magnetostatic equilibrium is an immediate exception. The expectation concerns topology-constrained relaxation for which a smooth limiting arrangement is obstructed. Current-sheet formation is consistent with the relaxation framework developed in Moffatt's primary analysis. Retaining complete frozen connectivity is essential; minimizing energy subject only to the value of total magnetic helicity would generally permit field rearrangements forbidden by this evolution.

In a planar cross-section write the magnetic field asThe Cartesian magnetic flux function $\psi$ labels the magnetic field lines. Ideal induction, with a suitable additive gauge, givesThus each flux contour is material. Incompressibility also preserves the area inside any material contour. Neither a magnetic island's flux distribution nor the connections of a separatrix can be freely altered during relaxation.

A configuration with two or more magnetic islands separated by an X-type saddle illustrates the mechanism. As the magnetic islands reshape to lower magnetic energy, the arms of the separatrix can be pressed together. Magnetic reconnection would change the topology, so it is excluded in a perfect conductor. The X-type structure can instead flatten into an extended interface, with oppositely directed tangential fields on its two sides. The transition shrinks and the electric current density concentrates into a current sheet. Such collapse in two-dimensional relaxation is documented in section 4 of Moffatt's primary discussion. It is a possible asymptotic outcome, not a finite-time singularity assertion or a property of every planar field.

The local equations show precisely how the interface can support equilibrium. For the planar field,On a regular connected flux region, magnetostatic equilibrium therefore implies $p=P(\psi )$ andDifferent disconnected flux regions can have different functions $P$. Their limiting fields need not join with continuous first derivatives of $\psi$ across a common separatrix; continuity of total pressure allows the tangential-field jump.

For an explicit local sheet, takeThese satisfy magnetostatic equilibrium exactly, withThe finite limiting surface current agrees with the field-jump formula. This is a local equilibrium demonstration of the sheet balance, not a claim that this particular family is the full ideal relaxation trajectory from the given initial data.

Strictly planar fields of the form above have zero magnetic helicity: choose $\mathbf{A}=\psi {\mathbf{e}}_z$, giving $\mathbf{A}\cdot \mathbf{B}=0$. Thus this strictly planar example interprets the last part as a separate illustration of topology-constrained relaxation. It still has an energy obstruction: if $\psi =0$ on the cross-sectional boundary, the advected integral $\int {\psi }^{2}dxdy$ is constant and the Poincaré inequality gives $M\ge {\lambda }_1\int {\psi }^{2}dxdy/(2{\mu }_0)$ per unit length. If the earlier nonzero-helicity condition is retained instead, a possible extension is a field depending on only two coordinates but permitting an axial component, with periodic boundary conditions in the invariant direction. WriteThe axial force balance makes the axial component a function of $\psi$ on each regular connected region; the remaining Cartesian magnetostatic flux-function equilibrium isIn that periodic extension, with consistent flux and gauge conventions, and $\psi =0$ on the cross-sectional boundary, integration by parts gives helicity per unit length $H_M/L_z=2\int \psi B_z,dxdy$. It can be nonzero. The same separatrix-collapse mechanism still concentrates current in a sheet. Three-dimensional linkage is therefore unnecessary for sheet formation; planar frozen flux-contour topology already supplies the constraint.

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.

Take zero mean flow locally and incompressible fluctuations. Subtracting the averaged induction equation givesThe first-order smoothing approximation omits the second term on the right when solving for $\mathbf{b}$, but retains $⟨\mathbf{u}\times \mathbf{b}⟩$ in the mean equation. For a locally uniform, slowly changing mean magnetic field the forcing becomes $(\overline{\mathbf{B}}\cdot \nabla )\mathbf{u}$.

There are two controlled ways to make the omitted fluctuation term small. At small eddy magnetic Reynolds number $R_m=u_{\mathrm{rms}}\ell /\eta$, diffusion gives $b/B\sim R_m$, so the omitted $ub/\ell$ is smaller than the retained $uB/\ell$ by $R_m$. Alternatively, if $u_{\mathrm{rms}}\tau /\ell \ll 1$ for the fluctuation response time $\tau$, then $b/B\sim u_{\mathrm{rms}}\tau /\ell$ and the short-memory response is perturbative. Mean-field scale separation is additionally required. Merely calling the turbulence weak or assuming that the mean magnetic field is small does not justify the omission.

With homogeneous initial transients removed, the causal linear solution isInsert this into the electromotive correlation. Isotropy gives ${\alpha }_{ij}=\alpha {\delta }_{ij}$, so taking the trace yields the helicity formula for isotropic first-order smoothing:Using the alternating-tensor sign givesThe minus sign follows from ${ϵ}_{imn}=-{ϵ}_{min}$. This derivation specifies the diffusion and correlation weighting, rather than asserting that alpha is proportional to instantaneous helicity for arbitrary turbulence.

If $\eta \tau /{\ell }^2\ll 1$, diffusion over a correlation time is negligible. Define the helicity correlation time byFor an ordinary decaying correlation with nonzero mean kinetic helicity, this gives the familiar estimateIt is the kinetic helicity, not magnetic helicity. Opposite flow handedness reverses the alpha effect. The time-integral formula remains meaningful if the equal-time helicity vanishes or its correlation changes sign, when replacing it by a single positive correlation time may fail. In the same weak-diffusion short-memory approximation the gradient forcing $-(\mathbf{u}\cdot \nabla )\overline{\mathbf{B}}$ gives $\beta =\frac{1}{3}{\int }_0^{\infty }⟨\mathbf{u}(t)\cdot \mathbf{u}(t-\tau )⟩d\tau \simeq {\tau }_U{u}_{\mathrm{rms}}^2/3$; the two correlation times need not be equal.

In the small-$R_m$, slowly varying quasistatic regime, diffusion is instead essential. The response givesFor velocity concentrated at wavenumber $k_0$, this becomes $\alpha =-⟨\mathbf{u}\cdot \mathbf{\varpi }⟩/(3\eta k_0^2)$. Its effective memory is the magnetic diffusion time on the eddy scale. The controlled low-$R_m$ setting is developed in Moffatt's analysis of turbulent dynamo action; the causal-response calculation above also covers short-memory turbulence without assuming that diffusion is always negligible.

Statistical isotropy and statistical reflection symmetry allow the spectral tensor to contain only ${\delta }_{ij}$ and $k_i{k}_j$:The solenoidal magnetic-field constraint is essential here. Taking its Fourier transform gives $k_i{\hat{B}}_i=0$, hence $k_i{\hat{C}}_{ij}=0$ and $A+D{k}^2=0$. Thus the spectral tensor is a scalar multiple of the transverse projector of a vector field. Its trace is $2A$, so the prescribed normalization givesThis is the solenoidal isotropic spectral tensor. Without the solenoidal condition, isotropy and reflection symmetry would leave two scalar functions. Without reflection symmetry, an antisymmetric term proportional to $i{ϵ}_{ijm}{k}_m$ could additionally encode magnetic helicity.

Put $D_t={\partial }_t+\mathbf{u}\cdot \nabla$. The continuity equation and the expanded ideal magnetohydrodynamic induction equation areso $D_t(B_i/\rho )=(B_j/\rho ){\partial }_j{u}_i$. Because $\mathbf{B}=\nabla \times \mathbf{A}$, the induction equation permits the magnetic vector potential equation ${\partial }_t\mathbf{A}=\mathbf{u}\times \mathbf{B}-\nabla \phi$, with an arbitrary scalar electric potential $\phi$ carrying the gauge choice. Using $(\mathbf{u}\times \nabla \times \mathbf{A}{)}_i=u_j{\partial }_i{A}_j-u_j{\partial }_j{A}_i$ yieldsContracting with $B_i/\rho$ and differentiating that factor givesThe last two terms cancel by exchanging their dummy indices. Thus the requested material derivative isThis is material transport of magnetic helicity density.

For a material volume, mass conservation gives $d{\int }_{V(t)}\rho hdV/dt={\int }_{V(t)}\rho D_{t}hdV$. Apply this with $h=\mathbf{A}\cdot \mathbf{B}/\rho$, and use $\nabla \cdot \mathbf{B}=0$ and the divergence theorem:If the magnetic field lines do not penetrate the boundary, $\mathbf{B}\cdot \mathbf{n}=0$, so $\boxed{d{H}_m/dt=0}$.

For uniqueness under a regular, single-valued gauge transformation $\mathbf{A}\mapsto \mathbf{A}+\nabla \Lambda$,Consequently this magnetic helicity is both conserved and gauge independent. As usual, vector potentials are understood to be globally regular potentials for the same field, rather than unrelated harmonic-circulation choices on a multiply connected restricted domain.

Write $m_0={\mu }_0\rho$ and treat the imposed mean magnetic field as uniform on the forcing scale. In first-order smoothing, retain terms linear in the fluctuating fields in their equations, while keeping their quadratic correlation when calculating the mean-field electromotive force. The inertial fluctuation, nonlinear induction term and quadratic fluctuating Lorentz force density are neglected in the steady response. The linear magnetic force isAbsorb the gradient into the total pressure. The resulting equations arewith $\nabla \cdot \mathbf{u}=\nabla \cdot \mathbf{b}=0$. Taking the divergence of the first equation gives ${\nabla }^{2}P=0$ for solenoidal forcing. In a periodic or homogeneous Fourier calculation its nonconstant part therefore vanishes.

For each Fourier mode with wavevector $\mathbf{q}$, $|\mathbf{q}|=k$, define $\ell =\mathbf{B}\cdot \mathbf{q}$ andSolving the two coupled linear equations gives the monochromatic coupled magnetic and velocity responseThe magnetic coupling suppresses modes with large projection of their wavevector along the imposed field.

To calculate the mean-field electromotive force, group each pair of opposite wavevectors into a real mode and average its spatial phase. Put $\mathbf{n}=\mathbf{q}/k$, and define its forcing correlationsFor any solenoidal real mode $\mathbf{v}=\mathbf{a}\cos (\mathbf{q}\cdot \mathbf{x})+\mathbf{c}\sin (\mathbf{q}\cdot \mathbf{x})$, its two amplitude vectors are perpendicular to $\mathbf{q}$. Direct multiplication givesAlso, integration by parts gives $⟨(L\mathbf{g})\times (L\mathbf{f})⟩=-{\ell }^2⟨\mathbf{f}\times \mathbf{g}⟩$, where $L=\mathbf{B}\cdot \nabla$. Expanding the response products therefore yieldsSumming over independent real modes, or integrating over their angular spectrum, is the required general answer for $⟨\mathbf{u}\times \mathbf{b}⟩$.

At zero field, $\mathcal{E}(0)={\sum }_{\mathbf{q}}{\mathbf{C}}_{\mathbf{q}}/(\nu \eta k^4)$. Thus the stated zero-field condition removes the unweighted sum of the cross correlations. For anisotropic correlated forcing it does not necessarily remove their differently weighted sum at nonzero field. The displayed general formula retains that allowed contribution. It disappears if cross forcing is absent mode by mode, or under the jointly isotropic statistics assumed for the final part.

The helicity content can also be expressed directly using the actual fluctuating fields. Define their kinetic helicity $H_u(\mathbf{q})=⟨{\mathbf{u}}_{\mathbf{q}}\cdot \nabla \times {\mathbf{u}}_{\mathbf{q}}⟩$ and current helicity $H_b(\mathbf{q})=⟨{\mathbf{b}}_{\mathbf{q}}\cdot \nabla \times {\mathbf{b}}_{\mathbf{q}}⟩$. Substitute $\mathbf{f}=\nu k^2\mathbf{u}-L\mathbf{b}/m_0$ and $\mathbf{g}=\eta k^2\mathbf{b}-L\mathbf{u}$ into their cross correlation. This gives the exact identityIt shows the competing kinetic and magnetic contributions without introducing a spurious singularity when the factor on the left vanishes. For a monochromatic field, ${\mathbf{a}}_b=\nabla \times \mathbf{b}/k^2$ is a magnetic vector potential, so $H_b=k^2⟨{\mathbf{a}}_b\cdot \mathbf{b}⟩$: current helicity is $k^2$ times magnetic helicity. In the isotropic weak-field limit, with zero-field fluctuations ${\mathbf{u}}_0,{\mathbf{b}}_0$, the scalar alpha effect isThis identifies both signs and their different viscous and magnetic relaxation times.

For jointly isotropic forcing let $h_f,h_g$ now denote the total forcing helicities, and set $S_h=\nu h_g/m_0-\eta h_f$. The cross-forcing contribution is even in $\mathbf{B}$, whereas isotropy requires any vector depending only on $\mathbf{B}$ to be parallel to it and odd under reversal. Hence that contribution vanishes. A natural symmetric alpha tensor isOnly its longitudinal coefficient contributes to $\mathcal{E}={\alpha }_{\parallel }\mathbf{B}$. Introduce $h=|B|/[k\sqrt{m_0\nu \eta }]$ and $A_0=\nu \eta k^4$. Angular integration givesTherefore the strong-field quenching of an isotropic electromotive force isprovided $S_h\ne 0$. If the two helicity contributions cancel exactly, this electromotive force vanishes.

There is a tensor qualification to the printed scaling claim. The natural transverse coefficient instead contains $\frac{1}{2}{\int }_0^1(1-{\mu }^2)(1+h^2{\mu }^2{)}^{-2}d\mu \sim \pi /(8h)$ and scales as $|B{|}^{-1}$. The relation ${\mathcal{E}}_i={\alpha }_{ij}{B}_j$ by itself leaves transverse components undetermined. One may choose the effective representative ${\alpha }_{ij}={\alpha }_{\parallel }{\delta }_{ij}$, for which the requested $|B{|}^{-3}$ scaling holds. It is the observable longitudinal response that is fixed, rather than every component of a unique field-dependent tensor.