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

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

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

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

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

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

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

Use incompressible ideal MHD at scales above the ion skin depth, a strong uniform mean field, and a statistically steady local energy cascade with mean flux $ϵ$ per unit mass. Assume balanced counterpropagating wave populations, comparable velocity and magnetic amplitudes in velocity units, and no additional scale-dependent alignment or intermittency correction. Write $\mathbf{b}=\delta \mathbf{B}/\sqrt{4\pi \rho }$ and introduce the fluctuation Elsässer variables ${\mathbf{z}}^{\pm }=\mathbf{u}\pm \mathbf{b}$. Their equations areOnly opposite populations distort one another. A single population has no leading nonlinear self-interaction, so the balance assumption is essential to a one-amplitude cascade estimate.

Let $\delta z_l$ denote the typical balanced amplitude and let $l_{\perp },l_{\parallel }$ be coherence lengths relative to the local mean field. The Alfvén wave time and eddy turnover time areFor $\chi \ll 1$, one counterpropagating encounter changes the packet fractionally by order $\chi$. Assume phases decorrelate between encounters. After $N$ encounters, random accumulation is of order $\sqrt{N}\chi$; an order-one change requires $N\sim {\chi }^{-2}$. The weak wave cascade time is thereforeWeakness by itself does not prove this random walk estimate; phase decorrelation and the interaction geometry are additional assumptions.

For the Iroshnikov-Kraichnan spectrum, impose isotropy, $l_{\parallel }\sim l_{\perp }\sim l$, in addition to $\delta z_l\ll v_A$. Then ${\tau }_A=l/v_A\ll {\tau }_{\mathrm{nl}}=l/\delta z_l$ andHence $\delta z_l\sim (ϵv_{A}l{)}^{1/4}$. A local shell-integrated turbulent energy spectrum has $\delta z_l^2\sim kE(k)$, $k\sim l^{-1}$, soThe isotropy is imposed phenomenologically; a strong guide field actually distinguishes directions and does not force this isotropic picture. The spectrum is not an automatically valid asymptotic result for every strongly magnetized cascade.

For the anisotropic weak Alfvénic cascade, keep $l_{\parallel }$ approximately fixed while cascading to smaller $l_{\perp }$, and require $\chi \ll 1$. This geometry follows from resonant three-wave Alfvén interactions: frequency conservation together with parallel-wavenumber conservation makes one member of an oppositely propagating resonant triad have zero parallel wavenumber, so the propagating modes retain their parallel wavenumber in the leading resonant transfer. Thus there is no leading parallel cascade in this ideal weak-wave ordering. Using the fixed parallel scale,Therefore $\delta z_l\sim (ϵv_A/l_{\parallel }{)}^{1/4}{l}_{\perp }^{1/2}$. For the one-dimensional perpendicular spectrum, $\delta z_l^2\sim k_{\perp }{E}_{\perp }(k_{\perp })$ givesUnlike the isotropic estimate, the parallel length does not decrease with the perpendicular length. Since $\chi =\delta z_l{l}_{\parallel }/(v_A{l}_{\perp })\propto l_{\perp }^{-1/2}$, the interactions become stronger down the cascade. If ${\chi }_0\ll 1$ at outer perpendicular scale $L_{\perp }$, the weak approximation breaks at $l_{\perp ,\ast }\sim L_{\perp }{\chi }_0^2$. This is the weak-to-strong transition of an Alfvénic cascade.

For Goldreich–Sridhar turbulence, abandon the weak-collision estimate and impose critical balance, ${\tau }_A\sim {\tau }_{\mathrm{nl}}$. The local cascade is strong, with transfer time of order one eddy turnover time. ThenThe parallel scale is set dynamically, not held fixed:At smaller scales $l_{\parallel }/l_{\perp }$ grows, so the eddies are increasingly elongated along the local magnetic field. Strong turbulence does not mean $\delta z_l\sim v_A$: small transverse amplitudes can still be strong when the eddies are sufficiently anisotropic to satisfy critical balance. Imbalance, alignment or nonlocal interactions would require different cascade reasoning.

Finally, in the usual isotropic hydrodynamic Kolmogorov 1941 theory inertial-range argument, the only relevant dimensional quantities are $ϵ$ and $k$. With $[ϵ]=L^2{T}^{-3}$ and $[E]=L^3{T}^{-2}$, writing $E\propto {ϵ}^a{k}^b$ forces $a=2/3$, $b=-5/3$. This uses the physical assumptions that outer-scale, viscous and intermittency effects do not supply another leading dependence; dimensions alone are not a proof of those assumptions.

A magnetic guide field supplies the independent speed $v_A$, and anisotropy adds independent parallel lengths. Even before adding anisotropy, dimensions permitwith arbitrary dimensionless $\Phi$. For example $\Phi (x)\propto x^{1/2}$ reproduces the isotropic Iroshnikov-Kraichnan spectrum estimate. This dimensional freedom of an MHD spectrum explains why propagation, decorrelation, anisotropy and cascade-time assumptions must select among the spectra: purely dimensional considerations do not select a unique MHD exponent.