Use the modulus of continuity on the circle, $\omega (f,\delta )={\sup }_{|h|\le \delta }\Vert f(\cdot +h)-f{\Vert }_{\infty }$. Splitting a displacement into steps of length at most $\delta$ and applying the triangle inequality givesThe Fejér kernel is nonnegative and has integral one, soWe therefore need its first absolute moment. The inequalities $|\sin (nt/2)|\le n|\sin (t/2)|$ and $\sin (|t|/2)\ge |t|/\pi$ on $[-\pi ,\pi ]$ implySplitting at $1/n$ givesTaking $\delta ={\delta }_n=(\ln n)/n$ in the preceding error bound proves the Fejér first-moment approximation bound:The constants do not depend on $f$ or $n$. By the Heine-Cantor theorem, a continuous function on the compact circle is uniformly continuous, hence $\omega (f,{\delta }_n)\to 0$ as ${\delta }_n\to 0$. Thus the Fejér sums converge uniformly to $f$. The index restriction is necessary: at $n=1$ the printed scale is zero, while ${\sigma }_{1}f$ is the mean of $f$ and generally differs from it.

In a thin disk dominated by a central stellar mass, the vertical gravitational potential is harmonic to leading order, so its vertical acceleration is $-{\Omega }_0^{2}z$. For a locally isothermal background, $p=c_s^2\rho$ with $c_s$ independent of height. Vertical hydrostatic equilibrium therefore givesThis is the isothermal vertical structure in stellar gravity. On the upper half of the disk the downward gravity magnitude is $g={\Omega }_0^{2}z>0$; the density profile extends evenly through the midplane.

For an adiabatically displaced parcel kept in pressure equilibrium, the parcel's logarithmic density gradient is ${\gamma }^{-1}d\log p/dz$, whereas the environmental density gradient is $d\log \rho /dz$. A small upward displacement $\eta$ consequently produces a density contrast $\delta \rho /\rho =[{\gamma }^{-1}d\log p/dz-d\log \rho /dz]\eta$. Its buoyancy acceleration is $-g\delta \rho /\rho =-N^2\eta$. Thus the Brunt–Väisälä frequency isBoth background gradients equal $-z/H^2$ in the isothermal atmosphere, soAn isothermal background is therefore stably stratified for adiabatic perturbations with $\gamma >1$, but not buoyant for isothermal perturbations. In applying the adiabatic equation of state to a stratified atmosphere, its proportionality relates Lagrangian parcel changes: $\Delta p=(\gamma p/\rho )\Delta \rho$. Eulerian perturbations generally obey $p^{\prime }+\eta p_z=(\gamma p/\rho )({\rho }^{\prime }+\eta {\rho }_z)$ instead. Thus the printed simple proportionality must be interpreted as a parcel relation for the buoyancy calculation. For $\gamma =1$ in this isothermal background, the gradient terms cancel and the Eulerian proportionality holds too.

Now take isothermal perturbations, so $\gamma =1$ and $N^2=0$. The two supplied wave equations reduce toFor the nonsingular oscillatory case $\omega \ne 0$ and ${\omega }^2\ne {\Omega }_0^2$, substitute $u_z={\xi }_z/{\omega }^2$. With $x=z/H$, the resulting Hermite differential equation isThe exceptional frequencies must be treated through the original velocity-pressure equations rather than division by these factors.

The physical vertical condition is finite perturbation energy when integrated with the Gaussian density. In particular the pressure contribution is proportional to $\int e^{-x^2/2}|\xi {|}^{2}dx$ for nonzero frequency. Generic nonterminating solutions have an exponentially growing component at at least one vertical end and fail this condition; acceptable modes are the terminating Hermite polynomials. This is a finite-energy condition, not pointwise boundedness of the velocity in the vanishing-density atmosphere.

To prove the allowed values of $\alpha$, let the nonzero polynomial have degree $n$ and leading coefficient $a_n$. The highest power in ${\xi }_{xx}-x{\xi }_x+\alpha \xi$ has coefficient $(\alpha -n)a_n$. Hence $\alpha =n\in \{0,1,2,\dots \}$. Conversely, writing $\xi ={\sum }_j{a}_j{x}^j$ gives the recurrenceFor $\alpha =n$, choose the parity of $n$ and iterate until the series terminates at degree $n$. Thus all allowed Hermite vertical modes of an isothermal disk havewhere ${\mathrm{He}}_n$ is the Probabilists' Hermite polynomial. In particular, ${\mathrm{He}}_0=1$ and ${\mathrm{He}}_1=x$.

Let $K=kH$ and $s={\omega }^2/{\Omega }_0^2$. Eliminating the denominator for regular modes givesFor $n=0$, the nonzero-frequency solution isIts pressure variable is independent of height and $u_z=0$. These are inertial-acoustic waves, restored by epicyclic motion and pressure; at short radial wavelength they become ordinary acoustic waves. The multiplied dispersion polynomial also has $s=0$. This is not an extra oscillatory solution represented by nonzero $\xi$, since $\xi =i\omega p^{\prime }/\rho$ vanishes at zero frequency for finite pressure perturbation. The full horizontal equations have a zero-frequency balanced mode of an axisymmetric disk, with no radial velocity and azimuthal Coriolis acceleration balancing a radial pressure gradient. It must be recovered without dividing by $\omega$, rather than counted as a propagating wave. At $k=0$, the propagating branch limits to an epicyclic oscillation at $|\omega |={\Omega }_0$, another value at which the eliminated form must be interpreted by a limit.

For $n=1$, solve $(s-1{)}^2=K^{2}s$. The two branches areBoth signs of frequency occur on each squared-frequency branch. The upper branch is acoustic in the large-$|K|$ limit, $|{\omega }_{+}|\sim |k|c_s$; the lower branch is an inertial wave with $|{\omega }_{-}|\sim {\Omega }_0/|K|$. At $|K|\ll 1$, both frequencies approach ${\Omega }_0$, reflecting the equality of the radial and vertical epicyclic frequencies in a Keplerian disk. Here $\xi \propto z$, so $u_z={\xi }_z/{\omega }^2$ is independent of height: these modes have a coherent vertical displacement of the disk column, rather than a vertical breathing motion.

For a local non-axisymmetric disturbance proportional to $e^{i(kR+m\varphi -\omega t)}$, the material frequency is $\hat{\omega }=\omega -m{\Omega }_0$. In the same local radial-wave approximation, replace $\omega$ by $\hat{\omega }$. For $n=1,m=1$, this givesThe four frequencies can equivalently be written $\omega /{\Omega }_0=1\pm (\sqrt{K^2+4}\pm |K|)/2$, with independent choices of the two signs. This is a local dispersion relation, not a claim that Doppler substitution supplies an exact global non-axisymmetric solution in a differentially rotating disk. The radial-wave ordering neglects azimuthal gradients compared with radial gradients; $|k|R\gg 1$ can coexist with $|k|H\ll 1$ in a thin disk.

Write $w=\omega /{\Omega }_0$. For the two low inertial-frame frequency branches, $|w|\ll 1$, while the material frequency remains near $-{\Omega }_0$. Expanding the dispersion relation givesThereforeFor signed $k$, the two leading branches are conventionally labelled by $\pm kH/2$. These are bending waves of a Keplerian disk: an $m=1$ vertical displacement varies once around an annulus and corresponds to a small tilt of its orbital plane. A radial variation of the tilt is a warp. The waves transport the warp radially at group speeds $c_s/2$ in opposite directions. At zero radial wavenumber their low-frequency limit is a stationary rigid tilt; the apparent slow pattern does not mean that the orbiting gas has slow vertical motion in its own frame.

A passive disk is heated predominantly by intercepted stellar radiation, rather than by its own viscous dissipation. A rough luminosity criterion, ignoring geometric factors, isThe intrinsic accretion disk luminosity is of order $GM\dot{M}/(2R_{\ast })$, with additional accretion energy potentially released at the stellar surface. The local test is more precise: absorbed irradiation must dominate the viscous flux at the radius being considered. For an ordinary steady Keplerian accretion disk, sufficiently far outside its inner edge,where $\alpha$ is the irradiation grazing angle of a disk. Reflection, very small incidence angles and shadowing can invalidate a simple total-luminosity comparison. The following calculation assumes an illuminated, absorbing surface.

At a point $(R,z_s)$ on the upper surface, the outward unit normal in the meridional plane is proportional to $(-z_s^{\prime },1)$, and a light ray from the point source is directed along $(R,z_s)$. Its flux through the upper face is the incident inverse-square flux times minus their scalar product:For a thin disk with small surface slope this becomesThe tangent slope and ray slope differ by $\alpha \simeq z_s^{\prime }-z_s/R$. Positive incidence requires disk flaring, $d(z_s/R)/dR>0$, and an unobstructed ray from the star. A negative value of the formula would mean that direct illumination of this face is absent, not that it receives a negative heating flux. The hierarchy $z_s\gg R_{\ast }$ makes the finite angular size of the star negligible here.

Let $a=\mathcal{R}/\mu$, so the ideal gas equation of state gives the isothermal sound speed $c_s^2=a{T}_c$. In the stellar vertical gravitational potential, hydrostatic equilibrium has vertical acceleration $-{\Omega }_K^{2}z$, where ${\Omega }_K^2=GM/R^3$. Integration of $c_s^2{\partial }_z\rho =-\rho {\Omega }_K^{2}z$ gives a Gaussian atmosphere with disk scale heightConsequently $q=z_s/R=fH/R$ satisfies $T_c=GM{q}^2/(a{f}^{2}R)$. Substitute this into the stipulated Stefan–Boltzmann law $F_d=\sigma T_c^4$ and the irradiation flux:Here $f$ specifies surface height in units of disk scale height and is unrelated to the wind-energy fraction in the preceding solution. This separable differential equation integrates exactly toIf the thin-disc solution is formally matched to $q(R_M)=1$, then $B=1-R_0^2/R_M^2$, givingIt follows that the flaring passive disk asymptotic isThis applies in the thin intermediate region $R\ll \min (R_0,R_M)$, in addition to the point-source and stellar-gravity assumptions. The corrections are of relative order $R^2/R_0^2+R^2/R_M^2$. In particular, the printed hierarchy $R\ll R_M$ must still be accompanied by the initial thin-disc condition $q\ll 1$: if $R_M\gg R_0$, radii larger than $R_0$ cannot satisfy that approximation. The matching at $q=1$ is an extrapolation used to set an integration constant; the derivation is not a valid treatment of a thick surface at $R_M$ itself.

Identify the stellar surface temperature with its effective radiating temperature, so that luminosity of a spherical blackbody gives $L_{\ast }=4\pi R_{\ast }^2\sigma T_{\ast }^4$. The given stellar surface scale height has $H_{\ast }^2=a{T}_{\ast }{R}_{\ast }^3/(GM)$, hence $GM/(a{T}_{\ast }{R}_{\ast })=(R_{\ast }/H_{\ast }{)}^2$. Substituting into $R_0$ yieldsThus a small stellar surface scale height leaves a wide interval above $R_{\ast }$ and below $R_0$ in which the thin disk flaring law can apply. The point-source approximation imposes the stronger lower cutoff obtained from $z_s\gg R_{\ast }$:This cutoff is still much smaller than $R_0$ when $R_0\gg R_{\ast }$. The full interval is therefore $R_{\ast }^{7/9}{R}_0^{2/9}\ll R\ll \min (R_0,R_M)$, whenever these bounds overlap. An assumed thickening radius must lie above the lower cutoff for all of the original approximations to hold simultaneously. The stellar luminosity-temperature relation is essential to the final comparison; independent arbitrary values of $L_{\ast }$ and $T_{\ast }$ would not imply it.

Take $\dot{M}>0$ to mean inward accretion rate, and let $S$ count loss through both faces per unit projected area of the accretion disk. An annulus of width $dR$ loses $2\pi RSdR$ to the wind. In steady state, its incoming mass flux at its outer edge equals the inward flux at its inner edge plus this loss. Hence mass conservation givesIf a loss rate is instead defined per individual face, it must be doubled before substituting for $S$ here.

Use Keplerian rotation ${\Omega }^2=GM/R^3$. The vertically integrated viscous dissipation, summed over both disk faces, isHere $\nu$ is the kinematic viscosity, and $\Sigma$ is the surface density of a disk; for height-dependent viscosity the product means $\int \rho \nu dz$. Material already in a circular orbit has specific orbital energy $v_K^2/2-GM/R=-GM/(2R)$. Giving it zero total binding energy therefore costs an additional $GM/(2R)$ per unit mass. The total launch speed is the escape velocity $\sqrt{2GM/R}$, but one must not charge the wind for the orbital kinetic energy it already possesses. Equating wind power to $f{Q}^{+}$ givesCombining this with mass conservation yields the wind-powered mass loss from a steady accretion disk relation

Define the outward viscous torque in an accretion disk asThe net inward angular momentum flux is $\dot{M}h-\mathcal{G}$. Because each unit mass lost to the wind carries the local specific angular momentum $h$, annular conservation of angular momentum readsIn the wind region, $\mathcal{G}=Rh{\dot{M}}^{\prime }/(3f)$. For constant $f$, differentiation and cancellation of $h$ givesThus the requested Cauchy-Euler differential equation is

For $f=1/3$, a power-law trial $\dot{M}\propto R^p$ gives $p^2+p/2-1/2=0$, with roots $p=1/2,-1$. Therefore $\dot{M}=A\sqrt{R}+B/R$ in $R_{\ast }\le R\le R_0$. A zero inner viscous torque in an accretion disk implies ${\dot{M}}^{\prime }(R_{\ast })=0$, so $B=A{R}_{\ast }^{3/2}/2$. There is no point mass sink at $R_0$, so the accretion rate is continuous there and equals ${\dot{M}}_0$. The complete result isIts derivative is nonnegative throughout the wind region, as required for a nonnegative loss rate $S$. The derivative can jump to zero outside $R_0$: the energy relation applies only inside the wind zone and does not impose ${\dot{M}}^{\prime }(R_0)=0$ on its inner limit. The viscous torque in an accretion disk remains continuous and extends outside via ${\mathcal{G}}^{\prime }={\dot{M}}_0{h}^{\prime }$.

At the inner boundary the exact accreted fraction isThe last expression is a leading-order asymptotic equality, not the exact finite-radius fraction. Most of the supplied mass is expelled before reaching the star when $R_{\ast }/R_0$ is small.

A planetary magnetic field is evidence about electrical conductivity, fluid motion and thermal evolution, but its presence alone does not identify a present planetary dynamo. One must separate internal global fields, crustal remanent magnetization and currents driven by the solar wind. The following discussion uses the observational picture available by 2009, including the early MESSENGER flybys and the Voyager, Galileo and Cassini results.

Outside a planet's internal electric currents, the internal component is approximately a potential magnetic field. Its spherical harmonic of degree $\ell$ decreases as $r^{-(\ell +2)}$. Measurements at different positions and altitudes can therefore distinguish a deep, low-degree source from localized crustal anomalies and external magnetospheric currents. Repeated observations, radio emissions, aurorae and particle measurements supply complementary evidence for a magnetosphere, while slow secular variation and records of past reversals test the internal field's evolution.

For Earth, ground observatories and satellite measurements reveal a predominantly dipolar field of several tens of microtesla at the surface, tilted relative to the rotation axis and changing with time. Magnetic minerals in rocks record a much older field and repeated geomagnetic reversals. These observations, together with seismic evidence for a liquid outer core, support an operating planetary dynamo in conducting iron alloy. Core temperatures exceed the relevant Curie temperatures, so ordinary permanent remanent magnetization cannot explain the main field. Thermal convection, cooling and compositional convection associated with inner-core growth supply buoyancy; the Coriolis force organizes the motion. The observed nonaxisymmetric components and time dependence are consistent with a dynamic source, rather than a fixed bar magnet.

Mariner 10 and the first two MESSENGER flybys established a weak but coherent intrinsic field at Mercury, with roughly dipolar structure closely aligned with its rotation axis. It is much weaker than Earth's surface field, but its global character favors a deep source over an explanation solely by localized crustal patches. A planetary dynamo in a fluid part of Mercury's large iron-rich core is the natural interpretation, supported by evidence for a molten core component. Light elements can depress the alloy's freezing temperature, and thermal or compositional buoyancy can maintain motion. Slow rotation is not enough to rule out dynamo action. The limited pre-orbital measurements do not uniquely determine the detailed core structure or the mechanism responsible for the field's weakness.

For Venus, spacecraft including Pioneer Venus found no appreciable internally generated global field. Nevertheless its ionosphere interacts with the solar wind, producing magnetic draping, a bow shock and an induced magnetosphere. Those external signatures are not evidence for an internal planetary dynamo. Venus is comparable to Earth in size and may possess a substantial conducting core; the absence of a detectable dynamo therefore need not mean that all of its core is solid. Insufficient core heat loss or stable stratification can prevent the required convection. Its slow rotation can affect flow organization, but rotation rate by itself is not a sufficient explanation of the contrast with Earth or Mercury.

For Mars, Mars Global Surveyor found strong, spatially localized crustal anomalies, especially in ancient southern terrain, rather than an appreciable present global dipole. This is remanent magnetization, indicating that an ancient field magnetized the crust while it formed and cooled. The distribution of anomalies and weakly magnetized major impact basins constrain when that ancient planetary dynamo operated; crustal magnetization is a record of past activity, not a present global dynamo. Modern interaction with the solar wind also produces an induced magnetosphere. Cooling and changes in core heat extraction can explain the cessation of dynamo-driving convection; the observations alone do not prove complete solidification of the core.

For Jupiter, radio emissions, its intense radiation belts and measurements by Pioneer, Voyager and Galileo demonstrate a strong global intrinsic field and an enormous magnetosphere. The field is broadly dipolar but tilted and contains higher multipoles. It is generated in a deep conducting fluid rather than by a permanent solid magnet: the principal conducting region is metallic hydrogen in the high-pressure interior. Cooling and gravitational energy release maintain convection, while rapid rotation organizes the flow. Plasma supplied by Io and currents in the magnetosphere strongly modify the observed external field, but cannot account for the internally sourced planetary component. Jupiter's large magnetic moment also reflects its large size; a moment ratio should not be confused with the ratio of local surface field strengths.

For Saturn, Pioneer, Voyager and Cassini measurements likewise show an intrinsic field and a magnetosphere, with a dipole strikingly close to the rotation axis. A planetary dynamo in conducting metallic hydrogen is supported by the same broad interior argument as for Jupiter. Interior cooling and differentiation can provide power, including energy associated with helium separation. The near alignment does not invalidate the Cowling anti-dynamo theorem: that theorem excludes an exactly axisymmetric isolated dynamo field, not a predominantly axisymmetric exterior with undetected nonaxisymmetric structure in its source region. A stably stratified conducting layer above the active region is one possible way to attenuate the observable nonaxisymmetric components; it is a model, not a uniquely established consequence of the observations.

For Uranus, Voyager 2 in 1986 observed an intrinsic field with a dipole tilted by about $59^{\circ }$ to the rotation axis and significantly displaced from the centre. Higher multipoles are substantial, so a centred dipole is an inadequate description. The field and associated magnetosphere indicate an active planetary dynamo, despite the planet's small intrinsic heat flux. That low luminosity constrains its thermal evolution and may indicate inhibited deep mixing; it does not demonstrate an absence of all conducting-fluid motion. A conducting, volatile-rich fluid shell is a plausible source region. The term ice giant describes bulk composition, not a claim that the dynamo region is ordinary frozen ice.

For Neptune, Voyager 2 in 1989 also found a highly non-dipolar, displaced intrinsic field, with a dipole tilt near $47^{\circ }$. Its field, radio emissions and magnetosphere imply another active planetary dynamo. Its larger intrinsic heat flux than Uranus supplies evidence for substantial continuing interior energy transport. Electrically conducting water-ammonia-rich mixtures at high pressure, with ionic or otherwise mobile charges, can supply the conducting region without requiring the metallic hydrogen layer characteristic of gas giants. Numerical dynamos in convecting shells surrounding more stably stratified interiors can produce strongly nonaxisymmetric, multipolar fields resembling those of the two ice giants. Such models show that their unusual field geometry is compatible with dynamo action, but field measurements do not uniquely identify the shell's composition or thickness.

The maintenance mechanism is governed by the resistive induction equation,A conducting fluid needs a sufficiently large magnetic Reynolds number and a flow geometry that regenerates magnetic flux rather than merely amplifying one component temporarily. Three-dimensional stretching and twisting can couple poloidal and toroidal fields; strictly planar regeneration and a strictly axisymmetric magnetic field are restricted by the relevant anti-dynamo theorems. Rotation influences this regeneration through the Coriolis force, but does not provide an inexhaustible power source merely by being present.

Without regeneration, the characteristic magnetic diffusion time is ${\tau }_{\eta }\sim L^2/\eta$. For a terrestrial core scale of order $10^6$ metres and diffusivity of order one square metre per second, this is of order $10^4$–$10^5$ years, far shorter than geological ages. Long-lived global fields therefore require continual energy input. A schematic magnetic energy balance iswhen boundary energy fluxes vanish. Mechanical work against the Lorentz force density supplies the energy lost as Ohmic heating. Thermal cooling, latent heat release, compositional convection and differentiation supply that mechanical power. Magnetic feedback eventually changes the flow and saturates the field. Strong conductivity without an adequate energy source is not sufficient for a dynamo.

The 2009 observations support active global dynamos in Mercury, Earth, Jupiter, Saturn, Uranus and Neptune; Venus has no detected appreciable intrinsic dynamo field, and Mars preserves evidence of an extinct earlier dynamo in its crust. Induced fields and surviving remanence explain why planets without a current global dynamo can still display magnetic phenomena. The diversity of field geometries and strengths constrains, rather than uniquely solves, the diversity of planetary conducting regions and their convective histories.

Let the layer depth be $d$, its imposed temperature drop $\Delta T$, and its dimensional vertical field $B_0\hat{\mathbf{z}}$. Use length $d$, thermal time $d^2/\kappa$, velocity $\kappa /d$, temperature $\Delta T$ and field $B_0$ as units. Write the total dimensionless field as $\hat{\mathbf{z}}+\mathbf{b}$, and temperature as the conductive profile plus $\theta$. The Boussinesq equations coupled to the resistive induction equation areThe dimensionless parameters in this normalization areHere ${\alpha }_T$ is the thermal expansion coefficient, $\nu$ is kinematic viscosity, $\kappa$ is thermal diffusivity, and $\eta$ is magnetic diffusivity. The first parameter is the Rayleigh number and $\sigma$ is the Prandtl number.

The displayed dispersion relation fixes the magnetic normalization: its $Q$ must be the thermal-diffusion magnetic-field parameter above. The conventional Chandrasekhar number is instead $Q_{\mathrm{Ch}}=B_0^2{d}^2/({\mu }_0\rho \nu \eta )$, so $Q=\zeta Q_{\mathrm{Ch}}$. If $Q$ is reserved for that conventional definition, the magnetic term in the printed polynomial requires an extra factor $\zeta$. We keep $Q$ for the printed polynomial throughout the following calculation.

Linearize about the conductive state. Absorbing linear magnetic pressure into $p$, the magnetic force becomes $\sigma Q{\partial }_z\mathbf{b}$. The pressure-eliminated vertical momentum equation isThe normal mode implicit in the given relation has $w,\theta \propto \sin \pi z$ and $b_z\propto \cos \pi z$, as for stress-free, fixed-temperature boundaries with the corresponding vertical-field magnetic conditions. Put $x=k^2$, $p_0={\pi }^2$ and $m=x+p_0={\beta }^2$. Taking growth proportional to $e^{st+i\mathbf{k}\cdot {\mathbf{x}}_h}$ gives the amplitude equationsTheir determinant isThis is the stated vertical-field magnetoconvection dispersion relation after accounting for its field-strength convention. The determinant derivation remains valid even at values where an intermediate division by $s+m$ or $s+\zeta m$ would not be allowed.

For steady marginality set $s=0$. One obtainsFor oscillatory marginality write $s=i\omega$ with $\omega \ne 0$. Dividing the polynomial by $m$ gives a real cubic $s^3+A{s}^2+Bs+C_0$, whereIts real and imaginary parts imply ${\omega }^2=B$ and $C_0=AB$. Solving for the oscillatory convection threshold givesandA physical oscillatory branch therefore requires $0<\zeta <1$ and positive frequency squared. In particular an algebraic continuation of $R_o$ into ${\omega }^2<0$ is not an oscillatory onset.

The steady-Hopf merger in vertical-field magnetoconvection occurs at zero frequency:Substitution verifies $R_s=R_o$. Equivalently the merged horizontal wavenumber satisfiesA finite-wavenumber intersection exists only when this expression is positive. At the intersection the growth polynomial has a double zero root: the nonzero-frequency neutral branch ends there.

For the requested merger at a magnetoconvection neutral-curve minimum, first note that a curve of the form $R=[A_{\ast }{m}^3+C_{\ast }m]/x$ has derivative zero exactly whenIts unique positive minimum has $x>p_0/2$, because the left side increases there from zero to infinity. At this minimum $R_c=2A_{\ast }{m}^3/p_0$.

Merger at the steady minimum. For $R_s$, the stationarity condition is $m_s^2(2x_s-p_0)=Q{p}_0^2/\zeta$. Combining it with the merger equation givesThe corresponding critical Rayleigh number is $R_{s,c}=2(x_s+p_0{)}^3/p_0$. For fixed positive $\sigma$ and large $Q$, set $\zeta =1-{\delta }_s$ and retain the leading powers. ThenMore generally, before taking $\zeta \to 1$, the steady minimum obeys $k_s^2\sim [Q{\pi }^4/(2\zeta ){]}^{1/3}$.

Merger at the oscillatory minimum. Define $A_o=(\sigma +\zeta )(1+\zeta )/\sigma$. The stationarity condition for $R_o$ isCombining it with the zero-frequency relation instead givesHere $R_{o,c}=2A_o(x_o+p_0{)}^3/p_0$. For fixed $\sigma >0$ and $Q\gg 1$, setting $\zeta =1-{\delta }_o$ givesThese two merger conditions are distinct parameter surfaces, not simultaneous assertions about the same minimum. Both have the strong-field wavenumber selection in magnetoconvection scaling $k=O(Q^{1/6})$.

If the first subleading Rayleigh corrections are desired, let $X_s=[Q{p}_0^2/(2\zeta ){]}^{1/3}$ and $X_o=[\sigma Q{p}_0^2/(2(1+\sigma )(1+\zeta )){]}^{1/3}$. With parameters bounded in the above regime,These retain the $O(Q^{2/3})$ corrections, including the effect of the $O(Q^{-1/3})$ diffusivity deficit at a merger.

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.