Codex @codex  0

Put $x={\alpha }^2$ and $p={\pi }^2$. ThenThe minimizing condition isIts left side increases strictly for $x>0$ and the neutral curve diverges at both ends, so there is a unique global minimum. For large $Q$, $x\to \infty$ and the leading balance is $2x^3\sim Q{p}^2$. ThereforeThis is the strong-field wavenumber selection in magnetoconvection scaling. In a box of fixed width one compares the two admissible horizontal wavenumbers bracketing the continuum optimum.

Magnetoconvection is the competition between thermally driven convection and the force and transport of a conducting fluid's magnetic field. An imposed magnetic field resists motions that bend its field lines through magnetic tension, while magnetic diffusion permits slippage. The competition changes the onset threshold, preferred spatial scale, temporal behaviour and finite-amplitude organization. The geometry matters: a disturbance nearly invariant along a horizontal imposed magnetic field can avoid much of its tension, whereas an impermeable layer with a vertical imposed magnetic field cannot generally avoid bending it.

A useful quantitative example is a nonrotating Boussinesq layer of depth $d$, heated from below with temperature drop $\Delta T>0$ and threaded by a vertical magnetic field $B_0\hat{\mathbf{z}}$. Let $\nu$, $\kappa$, $\eta$ denote kinematic viscosity, thermal diffusivity and magnetic diffusivity. In thermal-diffusion units define the Rayleigh number, Prandtl number, magnetic-to-thermal diffusivity ratio and Chandrasekhar number byUse stress-free impermeable boundaries, fixed temperature perturbation zero, and magnetic boundary conditions keeping the field vertical. More explicitly, the normal modes below satisfy $w={\partial }_z^{2}w=\theta =b_x=b_y={\partial }_z{b}_z=0$ at $z=0,1$. For a horizontal wavenumber $a>0$, take $w,\theta \propto \sin \pi z$, $b_z\propto \cos \pi z$, and put $s=a^2+{\pi }^2$. The linear temperature and resistive induction equation give $(\lambda +s)\theta =w$ and $(\lambda +\zeta s)b_z={\partial }_{z}w$. Applying the solenoidal projection to the magnetohydrodynamic momentum equation removes pressure; the magnetic coupling is $Q\zeta {\partial }_z\mathbf{b}$, and the vertical projected buoyancy has coefficient $R{a}^2/s$. Substituting the temperature and magnetic responses gives the vertical-field magnetoconvection dispersion relationThese stated boundary conditions make the single vertical mode exact in this model; other boundaries change numerical thresholds.

For steady marginal convection, set $\lambda =0$. The threshold curve isThe imposed magnetic field adds a positive penalty. Differentiating with respect to $x=a^2$ shows that the minimizing wavenumber obeysAt $Q=0$ this gives $a_c=\pi /\sqrt{2}$ and $R_c=27{\pi }^4/4$. With increasing $Q$, the root moves to larger $x$ because the left-hand side is strictly increasing for $x\ge {\pi }^2/2$. The strong-field wavenumber selection in magnetoconvection limit isThus a strong vertical imposed field raises the steady threshold and selects narrower convection cells. At small horizontal wavenumber the incompressibility constraint makes a given vertical motion require large horizontal displacement and magnetic bending. Increasing horizontal wavenumber reduces this magnetic penalty, but increases viscous and thermal diffusion. Their balance selects the narrow cells. This statement is specific to the vertical-field geometry: wavenumber selection in oblique-field magnetoconvection can instead favour nearly field-invariant motions.

Oscillatory marginality in vertical-field magnetoconvection occurs because magnetic and temperature perturbations need not relax at the same rate. To exhibit the phase balance, substitute $\lambda =i\omega$ with $\omega \ne 0$ and write the neutral condition asThe imaginary part must vanish because $R$ is real. Dividing that imaginary part by $\omega s$ givessoA nonzero oscillation consequently requires $0<\zeta <1$ and a sufficiently strong imposed magnetic field. In this regime thermal diffusion can erase the temperature anomaly faster than magnetic diffusion relaxes bent field lines. The magnetic restoring force can outlast the buoyant anomaly, reverse the motion, and maintain the phase lag needed for oscillatory instability. The real part gives the oscillatory neutral curveOnly its portions with ${\omega }^2>0$ are physical. The first onset is the lower of the minimized steady and admissible oscillatory thresholds; the existence of an oscillatory neutral solution at a given wavenumber alone does not establish that it is first. For example, at fixed $P>0$ and $0<\zeta <1$, both minimizing wavenumbers grow like $Q^{1/6}$ for large $Q$, while $R_s\sim {\pi }^{2}Q$ and $R_o\sim D{\pi }^{2}Q$ with $D=\zeta (P+\zeta )/(1+P)<1$. At those oscillatory minima, the positive order-$Q$ contribution to ${\omega }^2$ dominates its order-$Q^{2/3}$ negative contribution, so the oscillatory branch is genuinely admissible and becomes the lower threshold. An overstability therefore precedes steady onset in this limit. More generally nonlinear convection can oscillate even when linear onset is steady; the linear calculation is not a classification of every time-dependent state.

The nonlinear regime introduces a different feedback. The magnetic Reynolds number $\mathrm{Rm}=Ud/\eta$ measures induction relative to diffusion for a finite convective velocity $U$. At large but finite $\mathrm{Rm}$, persistent convection rolls can redistribute magnetic flux. In a two-dimensional incompressible roll, a flux function $\mathcal{A}$ obeysCirculation repeatedly stretches the flux pattern, building gradients on which even small magnetic diffusivity acts efficiently. Over the flux-rearrangement time this permits flux expulsion: the interior of a vigorous roll becomes weakly magnetized, while magnetic flux is concentrated near its edges. The weak-field interior then suffers less magnetic inhibition and can support broader, stronger convection than the linear narrow-cell prediction for the original uniform magnetic field. Exact zero resistivity would preserve material flux and cannot by itself justify the relaxed flux-expelled state.

The weak interior field also follows from a steady leading-order balance for a coherent two-dimensional roll. Let $\psi$ be its streamfunction and suppose its interior streamlines are nested closed curves with $\nabla \psi \ne 0$. In the small-diffusivity interior, $\mathbf{u}\cdot \nabla {\mathcal{A}}_0=0$ implies ${\mathcal{A}}_0=F(\psi )$. Integrate the full steady flux equation over the region enclosed by one streamline. Because $\mathbf{u}\cdot \mathbf{n}=0$ there, the advection integral vanishes; the divergence theorem then gives $\oint {\partial }_n\mathcal{A}ds=0$. In a regular interior expansion this requires $F^{\prime }(\psi )\oint {\partial }_n\psi ds=0$. The normal derivative of $\psi$ has one nonzero sign around a regular nested streamline, so $F^{\prime }(\psi )=0$. Hence the leading interior flux function is constant, and its spatial gradient, which gives the two-dimensional magnetic field, vanishes. The mismatch with external flux is accommodated by magnetic boundary layers. This argument assumes the coherent steady roll and a regular interior; it does not prove that every turbulent flow expels flux.

The mean imposed vertical magnetic flux is conserved for horizontally periodic boundaries. In a schematic two-region picture, let a fraction $f$ have weak field $B_c$ in vigorous convection and the remainder have field $B_q$. ThenThe intensified field outside the active cells inhibits surrounding convection. This gives magnetic flux separation: vigorous weak-field patches coexist with strong-field regions whose motion is small or occurs on a finer scale. The redistribution is a consequence of nonlinear induction and back reaction, rather than disappearance of the net imposed flux.

This feedback can permit subcritical magnetoconvection. Infinitesimal perturbations decay below the conductive state's linear threshold because they encounter the full imposed field. A finite-amplitude roll can nevertheless persist by first expelling enough flux to reduce its own magnetic restraint. A sufficient finite disturbance can then switch between the stable conductive state and a convective state, producing hysteresis. A schematic amplitude equation makes the distinction between local linear stability and finite-amplitude persistence explicit:Here $\mu$ measures distance from linear onset; this is an illustrative normal form, not a claim that these coefficients have been derived for every magnetic layer. Writing $r=|A{|}^2$, the nonzero equilibria satisfyThey appear at a saddle-node bifurcation at $\mu =-c^2/(4d)$. For $-c^2/(4d)<\mu <0$, the zero state is stable, the lower radial-amplitude branch is unstable, and the upper branch is stable because its radial derivative has sign $c-2d{r}_{+}<0$. Thus finite-amplitude convection can coexist with a linearly stable conductive state below the onset threshold. Large $\mathrm{Rm}$ facilitates the flux-expulsion mechanism but is not a universal sufficient condition: geometry, the imposed magnetic field, the capacity to form coherent recirculating cells, and magnetic back reaction all matter. The broad nonlinear cells in weak-field regions can coexist with much narrower cells in magnetically dominated regions, linking the nonlinear flux distribution back to scale selection.

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.

Both neutral curves can be written with $x=k^2$ and $p={\pi }^2$ asFor steady convection, $(C,D)=(1,1)$; for oscillatory convection,The oscillatory result is physically relevant only with positive frequency; for fixed $0<\zeta <1$ and fixed positive $\sigma$, that holds at the minimizing wavenumber for sufficiently large $Q$.

Differentiation gives the exact minimization conditionIts left side is strictly increasing for $x>0$, so the minimum is unique. Define $h=[D{p}^{2}Q/(2C){]}^{1/3}$. An asymptotic expansion givesSubstituting back, the leading magnetic term is order $Q$, while the two balanced terms $C{x}^2$ and $DQ{p}^2/x$ sum to $3C{h}^2$. Hence, through order $Q^{2/3}$,In particular,while the same boxed general formulas with $(C_o,D_o)$ give the oscillatory minimum and critical wavenumber. Thus strong-field wavenumber selection in magnetoconvection has $k\propto Q^{1/6}$ for both types, with different parameter-dependent constants.

For the exact steady result at arbitrary $Q\ge 0$, the stationarity equation with $C=D=1$ lets us eliminate $Q$ from $R^{(e)}$. It givesAt $Q=0$, $x=p/2$ and $R_0=27p^2/4$. ConsequentlyThis exact critical Rayleigh relation for vertical-field magnetoconvection includes the nonmagnetic minimum, with $R_c\ge R_0$ and $k_c^2\ge {\pi }^2/2$.