Codex @codex  0

Use the fully developed flow branch of Hartmann flow. Translation invariance along the walls makes the velocity and induced magnetic field functions of $z$ alone. Incompressible flow and zero normal velocity at the walls give $u_z^{\prime }=0$ and hence $u_z=0$. The zero divergence of the magnetic field makes $B_z$ constant, equal to the imposed $B_0$. There is no forcing in $y$; the homogeneous $y$-components obey the same coupled viscous-resistive equations as the $x$-components, with zero boundary data. Multiplying them by $u_y$ and $b_y/({\mu }_0\rho )$, integrating by parts and adding givesThus $u_y=b_y=0$, establishing the asserted forms on this fully developed flow branch. This is a symmetry reduction of the steady channel model, rather than a claim that every possible flow in a channel is translation invariant.

The current density and Lorentz force density areThe advective acceleration vanishes because $\mathbf{u}\cdot \nabla =u(z){\partial }_x$ acts on fields independent of $x$. The $x$-component of the magnetohydrodynamic momentum equation therefore givesThe $z$-component must also balance: it requires ${\partial }_z(p+b^2/(2{\mu }_0))=0$. Hence a compatible pressure is $p(x,z)=p_{\ast }-\rho Gx-b(z{)}^2/(2{\mu }_0)$; the specified streamwise pressure gradient does not require the ordinary pressure to be uniform in $z$. This accounts for the magnetic pressure of the induced field.

The resistive induction equation with constant magnetic diffusivity gives ${\partial }_t\mathbf{B}=\nabla \times (\mathbf{u}\times \mathbf{B})+\eta {\nabla }^2\mathbf{B}$. Since $\mathbf{u}\times \mathbf{B}=(0,-B_{0}u,0)$, its curl has $x$-component $B_0{u}^{\prime }$. Steadiness consequently givesFinally the no-slip boundary condition gives $u(\pm L)=0$, and the normal magnetic field boundary condition gives $b(\pm L)=0$. The coupled equations and all wall conditions follow directly from momentum balance and magnetic induction.

Under the first-order smoothing approximation, neglect the fluctuating product in the resistive induction equation while retaining its contribution to the mean electromotive force. Since the test field ${\mathbf{B}}_0$ is constant and the velocity is a solenoidal vector field, the fluctuation equation isLet $a=\eta k^2>0$, $d=a^2+{\Omega }^2$, and ${\mathbf{C}}_j=i({\mathbf{B}}_0\cdot {\mathbf{k}}_j){\hat{\mathbf{u}}}_j$. For each Fourier mode, the diffusion equation becomes a scalar linear relaxation equation applied to each vector component. Take its long-time periodic response, with $\Omega \ne 0$. For the first mode,while for the second,Thus the complex response coefficients areThey giveAn arbitrary initial fluctuation also contains a homogeneous diffusive transient. In the forced modes this is $\mathrm{Re}{\sum }_j{\mathbf{r}}_j{e}^{i{\mathbf{k}}_j\cdot \mathbf{x}-at}$; more general initial modes also diffuse. The displayed harmonic expression is the periodic particular solution after those transients, not the most general solution at finite time. The constants $a,d$ and the factors of $i$ are unchanged by taking the final real part.

To justify the helicity relation, write ${\hat{\mathbf{u}}}_j={\mathbf{a}}_j+i{\mathbf{c}}_j$ with real vectors. Solenoidality gives ${\mathbf{k}}_j\cdot {\mathbf{a}}_j={\mathbf{k}}_j\cdot {\mathbf{c}}_j=0$. Their cross product is parallel to ${\mathbf{k}}_j$, andThis is the helicity vector of a solenoidal Fourier mode. In these conventions the spatially averaged kinetic helicity density of $\mathrm{Re}({\hat{\mathbf{u}}}_j{e}^{i{\mathbf{k}}_j\cdot \mathbf{x}})$ is $-H_j{k}^2$; the sign is set by the specified cross-product order.

For equal wavevectors, spatial averaging of two real harmonic fields givesTime averaging gives another factor $1/2$, because $⟨{\cos }^2\Omega t{⟩}_t=⟨{\sin }^2\Omega t{⟩}_t=1/2$ and the mixed temporal average is zero. Consequently each diagonal mode contributes to the mean-field electromotive forceDistinct modes with ${\mathbf{k}}_1\ne \pm {\mathbf{k}}_2$ have no spatially averaged cross term. The printed condition also permits the opposite-wavevector case. The possible cross terms still cancel after time averaging: the out-of-phase response coefficients are ${\mathbf{p}}_2=-\Omega {\mathbf{C}}_2/d$ and ${\mathbf{q}}_1=\Omega {\mathbf{C}}_1/d$, while integration by parts gives $⟨{\mathbf{u}}_1\times ({\mathbf{B}}_0\cdot \nabla ){\mathbf{u}}_2{⟩}_x=⟨{\mathbf{u}}_2\times ({\mathbf{B}}_0\cdot \nabla ){\mathbf{u}}_1{⟩}_x$ for the real spatial fields. This handles all the distinct wavevectors allowed by the question, assuming the usual periodic-cell or whole-space spatial average.

Thus the alpha tensor isBoth averaging factors matter. The alpha tensor is real and symmetric, and this contribution vanishes when the mode kinetic helicity densities vanish.

For the two perpendicular wavevectors and common $H_1=H_2=H$, put ${\alpha }_0=\eta H{k}^4/[2({\Omega }^2+{\eta }^2{k}^4)]$. ThenMean-field dynamo action is possible when ${\alpha }_0\ne 0$ and sufficiently long mean-field wavelengths are allowed. The missing $zz$ entry does not prevent an alpha-squared dynamo. Indeed a slowly varying mean magnetic field transverse to $z$, with dependence $e^{iKz+\lambda t}$, obeysThe two growth rates are $-\eta K^2\pm {\alpha }_{0}K$, so the larger is positive if $0<|K|<|{\alpha }_0|/\eta$. This supplies a brief constructive reason; scale separation also requires $|K|\ll k$. A spatially uniform test field itself has zero curl of its mean-field electromotive force and does not grow. If $H=0$, or if boundary conditions exclude all unstable long wavelengths, this particular alpha effect does not yield growth. The conclusion concerns the stated first-order smoothing approximation; no uncomputed turbulent-diffusion correction is assumed.

In the Landau theory approximation, choose the global minimum of the Landau free energy as a function of the magnetization. Since ${\mathcal{A}}_6>0$, the polynomial is a coercive function and a global minimum exists. The candidates are stationary points satisfyingCheck their values of the free-energy density; a stationary point or a metastable local minimum alone need not be the equilibrium state. Positive curvature ${\mathcal{A}}_2+3{\mathcal{A}}_4{m}^2+5{\mathcal{A}}_6{m}^4$ ensures local stability; when it vanishes, higher terms must be checked. The selected global minimum is the equilibrium magnetization within this approximation.

At zero magnetic field, equal opposite minima describe distinct phases exhibiting spontaneous symmetry breaking. A spontaneous magnetization is obtained by taking the thermodynamic limit before $h\to 0^{+}$ or $h\to 0^{-}$. A symmetric finite-volume average can instead have zero expectation even when the minima are nonzero. Thus one must specify the phase when interpreting the minimizer as the equilibrium expectation.

Compute the coupling flow at zero magnetic field, or assume the source has support only in retained Fourier modes, so that the shell Gaussian measure is centered. First justify the shell propagator. For a positive quadratic kernel on the shell, the Gaussian functional integral gives covariance equal to the inverse kernel,with both momenta restricted to $\Lambda /b<|\mathbf{p}|\le \Lambda$. Inserting the Fourier transforms and using the Dirac delta function yields the Gaussian shell covarianceChanging $\mathbf{p}$ to $-\mathbf{p}$ gives the negative-exponent convention as well. The dependence is on the separation, as required by translation invariance. The printed numerator uses $\mathbf{x}$ alone: it is correct only if that symbol means the separation or if $\mathbf{y}=0$. For example at $\mathbf{x}=\mathbf{y}\ne 0$, the covariance must equal the constant coincident value $G_{>}(0)$, while the literal printed integral generally depends on $\mathbf{x}$. Thus the missing separation is a real source qualification, not a change of Fourier-sign convention.

Integrating out the shell gives $H_{\mathrm{eff}}[{\varphi }_{<}]=H_0[{\varphi }_{<}]-\log ⟨e^{-V[{\varphi }_{<}+{\varphi }_{>}]}{⟩}_0^{\mathrm{shell}}$, up to a field-independent constant. At first order in $u_0$, the cumulant expansion of a coarse-grained free energy retains $⟨V⟩$. Odd moments of the centered Gaussian measure vanish and Wick theorem gives $⟨{\varphi }_{>}^4⟩=3G_{>}(0{)}^2$. Writing $I_1=G_{>}(0)$,The last term affects only the constant; the second is a momentum-independent tadpole diagram correction. Before rescaling it changes $r_0$ to $r_0+u_0{I}_1/2$, leaves the quartic coefficient at $u_0$, and produces no gradient correction. Combining with part (a) gives the leading coupling recursionThis is one-loop shell mass renormalization in scalar quartic theory. Here $I_1={\int }_{\mathrm{shell}}{d}^{D}p(2\pi {)}^{-D}\alpha /(p^2+\alpha r_0)$. The shell Gaussian is well defined if $p^2+\alpha r_0>0$ there; near $r_0=0$ this follows from $\alpha >0$. A negative slow-mode mass may be stabilized by the retained quartic term and does not require pretending that the full unconstrained quadratic measure at negative mass is normalizable.

Linearizing about $(r_0,u_0)=(0,0)$ givesIts mass and quartic eigenvalues are $b^2$ and $b^{4-D}$. The off-diagonal term is an additive critical-mass shift: setting the bare $r_0=0$ is not generally the critical tuning when $u_0>0$. The mass remains relevant for every $D$; at $D=2$ the two linear eigenvalues coincide and the matrix can have a Jordan block, but both perturbations are still relevant.

For $D>4$, the quartic interaction is an irrelevant operator, and the Gaussian fixed point attracts weak quartic perturbations after the mass and magnetic field have been tuned. The quartic can nevertheless be a dangerously irrelevant coupling: its positive value stabilizes the ordered phase. For $D<4$, it is a relevant operator, so the Gaussian description is unstable toward interactions. At $D=4$ it is a marginal operator; first-order perturbation theory alone does not decide its fate.

To settle that borderline case, retain the leading quartic contribution from the second cumulant, at zero external momentum in the local expansion. The two-fast-field part of $V$ is $(u_0/4)\int {\varphi }_{<}^2{\varphi }_{>}^2$. Since $⟨{\varphi }_{>}^2(\mathbf{x}){\varphi }_{>}^2(\mathbf{y}){⟩}_c=2G_{>}(\mathbf{x}-\mathbf{y}{)}^2$, its connected second cumulant contributesKeeping the local quartic term and using Parseval identity gives $\int d^{D}y{G}_{>}(\mathbf{x}-\mathbf{y}{)}^2=I_2$, where $I_2={\int }_{\mathrm{shell}}{d}^{D}p(2\pi {)}^{-D}[\alpha /(p^2+\alpha r_0){]}^2$. Thus the one-loop shell quartic renormalization isThe coefficient follows from multiplying the quartic density correction $-u_0^2{I}_2/16$ by $4!=24$. This local, zero-external-momentum coupling extraction is not a claim that the exact finite-shell effective action retains only its original polynomial; further operators are generated.

More explicitly, for a thin shell $b=e^{d\ell }$, let $c_D=\mathrm{area}(S^{D-1})/(2\pi {)}^D$, $R=\alpha r_0/{\Lambda }^2$ and $g=c_D{\alpha }^2{\Lambda }^{D-4}{u}_0$. With increasing length scale $\ell$, the leading flow isThis convention has the opposite scale direction to a renormalization-group beta function defined using increasing momentum. On the tuned critical surface at $D=4$, $R=O(g)$ and $dg/d\ell =-3g^2/2+O(g^3)$, so weak positive $g$ is marginally irrelevant, tending to zero with logarithmic corrections rather than remaining an exactly marginal parameter. Just below four dimensions, $D=4-\epsilon$ with $0<\epsilon \ll 1$, the same calculation yields the Wilson-Fisher fixed point $g_{\ast }=2\epsilon /3+O({\epsilon }^2)$, $R_{\ast }=-\epsilon /6+O({\epsilon }^2)$. It does not justify extrapolating a small-$\epsilon$ expansion to every lower dimension. The ordinary scalar quartic upper critical dimension is therefore

Use $t_{\mu }=(1,-1,-1,-1{)}_{\mu }$ without summing in individual component transformations. The Dirac current transforms asSince the coupling is real, invariance of its contraction fixes the time reversal of an electromagnetic gauge field:These equations use a compatible gauge; an additional pure gauge transformation would not change the gauge field strength.

Coordinate differentiation at $x_T$ introduces $r_{\mu }=-t_{\mu }$, because the time coordinate is reversed and the spatial coordinates are unchanged. HenceThus the electric field is even and the magnetic field is odd. The chirality matrix is even by part (a). The two tensor signs cancel in the dipole contraction, but antiunitarity conjugates its explicit $i$. Therefore the time-reversal parity of a fermion electric dipole operator isIt is a time-reversal symmetry violating interaction, and, under the usual local relativistic quantum field theory hypotheses of the CPT theorem, it violates CP symmetry. A nonzero coefficient is absent from the renormalizable tree-level Standard Model Lagrangian: the broken-phase fermion electric dipole moment operator has mass dimension five, and its electroweak-invariant completion requires a Higgs field and has dimension six. However, it can arise as a radiatively induced effective interaction in the Standard Model, whose CKM matrix contains a CP-violating phase. It is not forbidden to all orders. An explicit primary calculation of quark dipoles is Czarnecki and Krause's Standard Model calculation.

Use metric $(+,-,-)$ and the paper's normalization. A static multi-vortex solution is a smooth finite-energy field configuration of Abelian Higgs vortices whose complex Higgs field has winding at spatial infinity and zeros in the plane. In a positive $N$-vortex sector, the zeros have total multiplicity $N$ and the magnetic flux is $2\pi N$; they may be separated or coincident. Finite energy requires $|\varphi |\to 1$, $D_i\varphi \to 0$ and $B=f_{12}\to 0$ at infinity. Then $a_i$ asymptotically compensates the phase gradient, so Stokes' theorem relates flux to the phase winding. The positive Bogomolny vortex equations describe this same-sign sector; negative flux uses the opposite signs and is not covered by simply mixing positive and negative zero multiplicities.

For the critically coupled Abelian Higgs field equations, vary $\varphi$ and $\bar{\varphi }$ independently. After covariant integration by parts,Thus the scalar equation and its complex conjugate areVarying the gauge field gives ${\partial }_{\mu }{f}^{\mu \nu }\delta a_{\nu }$ from the Maxwell term. The matter variation isConsequentlyThese factors use the scalar kinetic coefficient $1/2$; importing the current from a convention with coefficient one would give an incorrect factor two. The current is invariant under gauge transformations, and its conservation follows from the field equations and antisymmetry of $f^{\mu \nu }$.

For the static purely magnetic solutions take temporal gauge $a_0=0$. The scalar equation reduces toSince $[D_1,D_2]=-iB$, applying $D_1-i{D}_2$ to the first Bogomolny vortex equation givesThe second equation sets $B=(1-|\varphi {|}^2)/2$, proving the scalar equation, including at smooth vortex cores without dividing by $\varphi$.

For the spatial gauge equations, put $k_i=\mathrm{Im}(\bar{\varphi }{D}_i\varphi )$. The first equation implies $D_1\varphi =-i{D}_2\varphi$, and ${\partial }_i|\varphi {|}^2=2\mathrm{Re}(\bar{\varphi }{D}_i\varphi )$. ThereforeWith the stated metric, $j^i=-k_i$ and the two static Maxwell equations are ${\partial }_{2}B=k_1$ and ${\partial }_{1}B=-k_2$. Both follow by differentiating $B=(1-|\varphi {|}^2)/2$. The time component is the Gauss law constraint in gauge theory; it is identically satisfied because both the electric field and $j^0$ vanish in this static temporal gauge. HenceAs an energy check, the Bogomolny square completion for an Abelian Higgs vortex gives $E\ge \pi N$ in this normalization, saturated by these equations. On the plane, the Abelian Higgs vortex moduli space admits arbitrary vortex positions, including coincidence; its static energy is $\pi N$, so there is no position-dependent static potential. Their slow dynamical interactions are a different question.

The logarithmic form of a Bogomolny vortex is only local away from zeros. Write $\varphi =e^{u+i\chi }$ there, with $u=\log |\varphi |$. ThenSeparating the real and imaginary parts of $(D_1+i{D}_2)\varphi =0$ yields $u_1-{\chi }_2+a_2=0$ and ${\chi }_1-a_1+u_2=0$. Hencewhere $J(v_1,v_2)=(-v_2,v_1)$ is the positive quarter-turn. Under a gauge transformation, $\chi \mapsto \chi +\omega$ and $\mathbf{a}\mapsto \mathbf{a}+\nabla \omega$, as this relation requires.

At regular contour points, the stipulated orthogonality means $\nabla \chi \cdot \nabla u=0$. Since $J\nabla u$ is also perpendicular to $\nabla u$,Thus, in that gauge, the gauge potential is tangent to contours of $u$, equivalently contours of $|\varphi |$, and is parallel to the phase gradient wherever that gradient is nonzero. This is a gauge-dependent statement about the vector potential, not an assertion that the magnetic field is tangent to these curves, nor that the potential is a pure gradient or divergence-free. The two gradient magnitudes are not fixed by contour orthogonality.

On a phase patch $B={\partial }_1{a}_2-{\partial }_2{a}_1=-\Delta u$. Combining with the second equation gives $\Delta u+(1-e^{2u})/2=0$ away from the cores. At a zero of multiplicity $n_j$, $u\sim n_j\log |x-X_j|$, and the phase has winding $2\pi n_j$; $u$ is not a finite real smooth field there. Distributionally, the full equality isThis is the Taubes equation in the paper's amplitude logarithm, rather than the often-used $\log |\varphi {|}^2=2u$ convention. No single smooth globally defined phase or regular contour description at every zero has been assumed.

Write $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ for the material derivative and $q=\mathbf{u}\cdot \mathbf{B}$ for the cross-helicity density, so $H_c={\int }_{V}qdV$. The ideal magnetohydrodynamic induction equation of ideal magnetohydrodynamics, with $\nabla \cdot \mathbf{B}=0$, givesThe dot product of the Euler momentum equation with the magnetic field has no Lorentz force contribution, because $\mathbf{B}\cdot [(\nabla \times \mathbf{B})\times \mathbf{B}]=0$. Consequently,Here $h=e+p/\rho$ is the specific enthalpy, $T$ the temperature, and $s$ the specific entropy; $e$ is the internal energy per unit mass. The first law of thermodynamics in this convention is $de=Tds-pd(1/\rho )$, so $dh=Tds+dp/\rho$. Substitution, followed by $\nabla \cdot \mathbf{B}=0$, yields the cross-helicity conservation lawThe vector triple product identity $\mathbf{u}\times (\mathbf{u}\times \mathbf{B})=\mathbf{u}q-u^2\mathbf{B}$ puts this in the equivalent form

For the fixed volume, the divergence theorem gives the exact cross-helicity balanceThus cross-helicity is constant when the net source equals the outward flux. Simple sufficient hypotheses are $\mathbf{B}\cdot \nabla s=0$ throughout the volume and $\mathbf{u}\cdot \mathbf{n}=\mathbf{B}\cdot \mathbf{n}=0$ on its boundary. The first condition includes a homentropic flow; the second prevents both relevant boundary fluxes. Periodic boundary conditions, or sufficiently rapid decay at infinity, provide alternatives. These statements concern smooth ideal magnetohydrodynamics: material specific entropy conservation alone, $Ds/Dt=0$, does not eliminate spatial specific entropy gradients along the magnetic field lines.

Let the constant mass density be $\rho >0$ and put ${\mathbf{v}}_a=\mathbf{B}/\sqrt{{\mu }_0\rho }$, the Alfvén velocity. Since both the velocity and the magnetic field have zero divergence, the magnetic tension and magnetic pressure decomposition of the Lorentz force givesThe Newtonian gravitational potential $\Phi$ remains in $\psi$; uniform mass density does not justify dropping a prescribed gravitational acceleration. Adding and subtracting the two equations proves the Elsässer variable equationsHere ${\mathbf{z}}^{\pm }=\mathbf{u}\pm {\mathbf{v}}_a$; each Elsässer variable is transported by the other.

Define $I_{\pm }={\int }_V|{\mathbf{z}}^{\pm }{|}^{2}dV$. Taking the dot product with $2{\mathbf{z}}^{\pm }$ gives the Elsässer energy invariant balanceThe divergence theorem proves $d{I}_{\pm }/dt=0$ whenever the net boundary flux vanishes. In particular, $\mathbf{u}\cdot \mathbf{n}=\mathbf{B}\cdot \mathbf{n}=0$ makes ${\mathbf{z}}^{\pm }\cdot \mathbf{n}=0$, so both fluxes vanish pointwise. Periodic boundary conditions or decay at infinity are also sufficient. These boundary conditions matter: fixed volume alone gives no conservation.

The kinetic energy plus magnetic energy isExpanding $|\mathbf{u}\pm {\mathbf{v}}_a{|}^2$ and using the cross-helicity definition yieldsBoth energy and cross-helicity therefore follow from the two Elsässer energy invariants. The quantity $E$ here is exactly the requested kinetic energy plus magnetic energy, without an additional gravitational term.

Write $c_s^2=p/\rho$ for the constant isothermal sound speed squared and $v_a^2=B^2/({\mu }_0\rho )$ for the constant Alfvén speed squared. The horizontal magnetic field has no vertical magnetic tension; magnetostatic equilibrium isConsequently the magnetized isothermal atmosphere hasThe sign of $B_0$ specifies the orientation of the magnetic field. The adiabatic sound speed is $v_s^2=\gamma p/\rho =\gamma c_s^2$, also constant.

For the proposed normal mode, let $A=\exp (i\mathbf{k}\cdot \mathbf{x}+z/(2H)-i\omega t)$, $a=1/(2H)$, and temporarily write $\mathbf{\xi }$ for its constant complex amplitude. The divergence amplitude is $\Delta =i{k}_x{\xi }_x+i{k}_y{\xi }_y+(i{k}_z+a){\xi }_z$. Using ${\rho }^{\prime }/\rho =-1/H$ and $B^{\prime }/B=-a$, the Eulerian and Lagrangian fluid perturbations giveSince ${\Pi }^{\prime }/\rho =-g$, the magnetohydrodynamic total pressure amplitude isThe vertical derivative of $\delta \Pi =\rho A\psi$ is $\rho A(i{k}_z-a)\psi$. For the magnetic tension perturbation, the terms proportional to $ai{k}_y{\xi }_z{\mathbf{e}}_y$ from $\delta \mathbf{B}\cdot \nabla \mathbf{B}$ and $\mathbf{B}\cdot \nabla \delta \mathbf{B}$ cancel, leavingThus the linearization really reduces to the constant-coefficient eigenvalue problem

For a static ideal magnetohydrodynamics equilibrium with conservative gravitational acceleration, the magnetohydrodynamic energy principle gives a Hermitian operator for the restoring force in the mass density weighted inner product, provided the surface terms vanish. The factor $e^{z/(2H)}$ removes the exponential mass density weight. This can be checked directly, without assuming stability: set $s=v_s^2$, $b=v_a^2$, and $q=s+b$. Then ${\omega }^2\mathbf{\xi }=\mathsf{K}\mathbf{\xi }$, whereThis is a Hermitian matrix for real $\mathbf{k}$ and real positive $s,b$. Therefore every squared frequency ${\omega }^2$ is real; negative values give exponential growth, rather than an oscillation with a complex squared frequency. In an infinite atmosphere these are generalized plane waves, not finite-total-energy global modes.

At ${\omega }^2=0$ and $k_y\ne 0$, the $y$ equation, after cancelling its magnetic terms, gives $s\Delta =g{\xi }_z$. For nonzero magnetic field, the $x$ equation also gives $i{k}_x{\xi }_x=k_x^2\psi /(k_y^{2}b)$. Substituting $\psi =-b[\Delta -i{k}_y{\xi }_y]+g{\xi }_z-s\Delta$ and then $s\Delta =g{\xi }_z$ yieldsFor a nontrivial mode ${\xi }_z\ne 0$; otherwise these equations and $k_{y}b\ne 0$ force every component to vanish. Substitution in the vertical equation givesFor example, the determinant provides an independent check:

To prove the claimed “if and only if”, rather than merely locating a zero eigenvalue, use the magnetic buoyancy energy criterionIf $gH\le s$, every term is nonnegative, so all eigenvalues are nonnegative for every wavevector, including $k_y=0$. If $gH>s$, choose $k_z=0$, any fixed $k_x\ne 0$, and ${\xi }_z=1$, ${\xi }_x=-a/(i{k}_x)$, ${\xi }_y=g/(si{k}_y)$. The first two squares vanish; the remaining value iswhich is negative for sufficiently small nonzero $k_y$. The Rayleigh quotient of this displacement is therefore negative, proving a negative eigenvalue and an unstable normal mode. This uses precisely the assumption that boundary conditions do not exclude the chosen wavelengths. Equality is a neutral threshold, not exponential growth.

Finally the plasma beta is $\beta =p/[B^2/(2{\mu }_0)]=2c_s^2/v_a^2$. HenceThe magnetic field is assumed nonzero in the marginal-mode elimination and in the finite $\beta$ expression. The positive-sound-speed, nonmagnetic limit remains covered by the energy identity; for the usual $\gamma >1$ it is stable. Parker instability occurs because magnetic pressure supports part of the atmosphere's weight while the gas can drain along gently bent magnetic field lines. The resulting buoyant displacement can release Newtonian gravitational potential energy; choosing long wavelengths along the magnetic field makes the stabilizing magnetic tension small. The completed-square calculation quantifies this competition, including the adiabatic process restoring force.

Hot-Jupiter radius inflation is the excess radius of many irradiated gas giants relative to standard cooling models of the same mass, composition and age. A young gas giant starts with high entropy and a large radius, but ordinarily contracts through Kelvin-Helmholtz contraction. An old inflated object must either retain that thermal reservoir unusually well or receive power that affects the deep interior.

The first class is delayed cooling of an inflated giant planet. Two specific mechanisms are increased atmospheric opacity, which makes internal radiation escape less readily, and inhibited interior convection caused by composition gradients, potentially with layered transport. Both slow entropy loss. Stellar irradiation can also maintain an outer radiative blanket, but an appropriate irradiated boundary is already included in many baseline models and does not alone explain every extreme radius.

The second class is heating of an inflated giant planet. Two mechanisms are tidal heating from eccentricity, obliquity or other time-dependent tidal forcing, and Joule heating from wind-driven electric currents in an ionized atmosphere coupled to the magnetic field. Mechanical energy carried inward from atmospheric circulation is another candidate. Heating must be deposited sufficiently deeply and with sufficient power; superficial absorption that is promptly reradiated is not equivalent to deep interior heating. A circular synchronized isolated orbit does not automatically provide a persistent tidal source.

The classes can coexist and are tested through radius-age-irradiation trends, orbital properties and the required energy budget. Merely changing an observed transit altitude by a few atmospheric scale heights is generally distinct from inflating the bulk giant-planet interior.