Codex @codex  0

If the horizontal magnetic field perturbations vanish at $z=\pm H$, their allowed vertical wavenumbers are $k_m=m\pi /(2H)$, $m=1,2,\dots$. The lowest wavenumber is unstable exactly whenIf that magnetic mode is stable, all higher magnetic modes are stable as well. A weak field admits longer unstable wavelengths that fit inside the slab; sufficiently strong magnetic tension stabilizes every allowed wavelength.

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

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

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

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

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

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

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

The horizontally invariant magnetized shearing-sheet equations are linear in the horizontal fields, so perturbations about the equilibrium satisfy the same equations. The fixed surface boundary conditions require $\delta B_x=\delta B_y=0$ at $z=\pm H$. Choose a normal mode withwhere the vertical wavenumber is $k=m\pi /(2H)$, $m=1,2,\dots$. With the Alfvén frequencythe four amplitude equations becomeEliminating the velocities leavesA nonzero amplitude requires the determinant of this system to vanish, yielding the ideal magnetorotational dispersion relationEquivalently, with ${\kappa }_r^2=2(2-q){\Omega }^2$,As a quadratic equation for $s^2$, its discriminant is ${\kappa }_r^4+16{\Omega }^2{\omega }_a^2>0$. If $0<{\omega }_a^2<2q{\Omega }^2$, its constant term is negative, so one root $s^2$ is positive and there is an exponentially growing mode. If ${\omega }_a^2>2q{\Omega }^2$, both the constant term and the coefficient of $s^2$ are positive, giving two negative roots and only oscillatory modes. Equality is marginal.

The lowest allowed vertical wavenumber, $k=\pi /(2H)$, is the last to be stabilized as $|B_z|$ increases. Therefore the finite-thickness magnetorotational instability criterion isThere is also a vertically uniform velocity mode with zero magnetic perturbation; for the usual orbitally stable regime $q<2$, it is just stable epicyclic motion. For $q>2$, that uniform mode is already hydrodynamically unstable, independently of the magnetic criterion. At $q=2$ it is marginal. The criterion above concerns the magnetic modes.

Putand let $q_1=nk\cos {\theta }_T$ be the downward vertical wavenumber in the lower half-space. The given form of Snell's law saysWith the $e^{-i\omega t}$ convention, write the incident, reflected, and transmitted plane waves asThe two interface conditions at $z=0$ giveSolving this linear system gives the reflection and transmission coefficients at a scalar-wave interfaceThus the exact fields are