Codex @codex  0

Magnetic subcriticality refers to a mass-to-flux ratio below its critical mass-to-flux ratio. For an idealized thin layer threaded by a uniform frozen-in perpendicular magnetic field, define in SI unitsHere $\Sigma$ is the surface density and $B_{z0}$ the perpendicular magnetic field. The magnetically subcritical case has $\lambda <1$, or $B_{z0}^2>\pi G{\mu }_0{\Sigma }^2$. In the specified ideal layer, magnetic support prevents the linear gravitational collapse studied by Takenori Nakano and Takashi Nakamura. Their original Gaussian-unit threshold is $\Sigma /|B_{z0}|<1/(2\pi \sqrt{G})$. The numerical coefficient is tied to geometry and field assumptions.^[ref]

A complementary sufficient-contraction result applies to an isolated cold initially resting razor-thin disk. If $B_{z0}^2<\pi G{\mu }_0{\Sigma }^2$ pointwise, the absolute gravity-equivalent magnetic surface density is smaller than the actual surface density. The positive-kernel comparison of thin-disk fields then makes the integrand in the horizontal virial balance of a cold magnetized fluid negative above and below the disk. Its radial second moment therefore begins to decrease. This weak-field inequality is on the magnetically supercritical disk side of the mass-to-flux convention.

The second-moment conclusion is global. It does not establish inward acceleration of every fluid element, and it does not prove unlimited subsequent collapse. Pressure, rotation, magnetic flux loss and different boundary-field configurations require their own analysis.

Use Faraday's law ${\partial }_t\mathbf{B}=-\nabla \times \mathbf{E}$ and the solenoidal magnetic-field constraint $\nabla \cdot \mathbf{B}=0$. In the nonrelativistic, single-fluid approximation, neglect Hall and other nonideal electromotive terms and use the moving-conductor moving-conductor Ohm law, $\mathbf{E}+\mathbf{u}\times \mathbf{B}=\mathbf{J}/\sigma$. Infinite electrical conductivity with finite current gives $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$, henceThis is the ideal magnetohydrodynamic induction equation. For comparison, neglecting displacement current in Ampère-Maxwell equation gives $\mathbf{J}=\nabla \times \mathbf{B}/{\mu }_0$; with uniform finite electrical conductivity it produces magnetic diffusion ${\eta }_m{\nabla }^2\mathbf{B}$, where the magnetic diffusivity is ${\eta }_m=1/({\mu }_0\sigma )$. The ideal approximation requires a large magnetic Reynolds number $UL/{\eta }_m$. Dropping displacement current is useful for this finite-conductivity comparison, but Faraday's law and the ideal Ohm relation already suffice for the ideal induction equation.

Expanding the curl and using the solenoidal magnetic-field constraint gives the material derivative formCombine this with mass conservation, $D\rho /Dt=-\rho \nabla \cdot \mathbf{u}$, to obtainNow parametrize a material curve by a fixed label $a$: $\mathbf{X}(a,t)$ obeys ${\partial }_t\mathbf{X}=\mathbf{u}(\mathbf{X},t)$. Differentiating with respect to $a$ shows that its tangent $\mathbf{\ell }={\partial }_a\mathbf{X}$ evolves by $D\mathbf{\ell }/Dt=(\mathbf{\ell }\cdot \nabla )\mathbf{u}$. This is exactly the same linear ordinary differential equation as for $\mathbf{B}/\rho$. Initially parallel tangents remain parallel by uniqueness, with a label-dependent proportionality factor constant along each particle trajectory. Thus magnetic field lines are transported as material curves, wherever the field and fluid flow are smooth and the field is nonzero. This is the field-line form of magnetic flux freezing.

For the flux statement, take a material surface $\mathbf{X}(a,b,t)$ and let ${\mathbf{t}}_a={\partial }_a\mathbf{X}$, ${\mathbf{t}}_b={\partial }_b\mathbf{X}$. Both tangents obey $D{t}_i/Dt=({\partial }_j{u}_i)t_j$. Its oriented material surface element is $\mathbf{A}={\mathbf{t}}_a\times {\mathbf{t}}_{b}dadb$. Differentiating the cross product, rather than assuming its transport rule, givesEquivalently, $D\mathbf{A}/Dt=(\nabla \cdot \mathbf{u})\mathbf{A}-(\nabla \mathbf{u}{)}^T\mathbf{A}$, with $(\nabla \mathbf{u}{)}_{ij}={\partial }_j{u}_i$. Contracting this derived rule with the induction equation gives a pointwise cancellation:Integrating over the fixed material labels therefore proves conservation of flux through an open material surface:The surface need not be closed; its boundary is carried with the fluid. The result follows from material transport, rather than from the zero flux through a closed surface.

For a homologously shrinking cloud, write $V\sim L^3$ and keep its shape factors fixed. Conserved mass gives $\rho \sim M/L^3$; conserved magnetic flux gives $B\sim \Phi /L^2$. Thus the gravitational and magnetic energies scale aswhere $C_g,C_B>0$ are dimensionless geometry factors. Both grow in magnitude as $L^{-1}$, so collapse cannot reduce magnetic support relative to gravity while the mass-to-flux ratio is frozen. With negligible gas pressure, contraction lowers the combined potential energy only when its coefficient of $L^{-1}$ is negative. Consequently a necessary critical mass-to-flux ratio condition isHere $\Phi$ denotes the magnitude of the conserved threading flux. The numerical coefficient depends on geometry and boundary conditions; the scaling argument does not determine it or make the condition sufficient in the presence of other support.

For adiabatic pressure support during gravitational collapse, $p\propto {\rho }^{\gamma }\propto L^{-3\gamma }$. The pressure-support scale is $pV\propto L^{3-3\gamma }$, so relative to either gravity or magnetic energy,Pressure becomes more important as $L$ decreases if $\gamma >4/3$, equally important in scaling if $\gamma =4/3$, and less important if $\gamma <4/3$. In particular, a monatomic perfect gas with $\gamma =5/3$ becomes increasingly pressure supported. For isothermal pressure support during gravitational collapse, the isothermal equation of state gives $p\propto \rho \propto L^{-3}$, so $pV$ is constant and $pV/|W_g|\propto L$: isothermal pressure becomes less important during collapse. The same comparisons hold against magnetic support because its energy has the same $L^{-1}$ scaling as gravity.