Codex @codex  0

For a Cartesian change of basis represented by an orthogonal matrix $Q$, vectors transform as $E_i^{\prime }=Q_{ia}{E}_a$ and $B_i^{\prime }=Q_{ia}{B}_a$. Their squared norms are invariant and $Q_{ia}{Q}_{jb}{\delta }_{ab}={\delta }_{ij}$. Substitution into the given expression therefore givesthe transformation law of a Cartesian second-rank tensor. The magnetic field's extra axial sign under an improper physical reflection, if included, appears twice and cancels in its quadratic contribution.

Write the Maxwell stress tensor as $T=\mathbf{EE}+\mathbf{BB}-\frac{1}{2}(E^2+B^2)I$. Its divergence isThe identity $(\mathbf{A}\cdot \nabla )\mathbf{A}-\frac{1}{2}\nabla A^2=-\mathbf{A}\times (\nabla \times \mathbf{A})$ follows by contracting two Levi-Civita symbols, or directly by differentiating components. Using Maxwell's equations consequently yieldsHence the local conservation of electromagnetic momentum isThe two subtracted terms are the Lorentz force density, while $\mathbf{E}\times \mathbf{B}$ is the electromagnetic momentum density in the units used here.

Write $\mathbf{u}={\partial }_t\mathbf{\xi }$. In the uniform equilibrium, the linearized ideal magnetohydrodynamic equations reduce toIntegrating from the undisplaced reference state, and using the uniform background magnetic field, givesThese are the perturbations induced by the fluid displacement; independent time-independent changes to the reference state are excluded. The linear Lorentz force density is $(\nabla \times \delta \mathbf{B})\times \mathbf{B}/{\mu }_0$. Its magnetic pressure and magnetic tension parts giveFor a Fourier mode, put $s=\mathbf{k}\cdot \mathbf{\xi }$, $q=\mathbf{B}\cdot \mathbf{\xi }$ and $K=\mathbf{k}\cdot \mathbf{B}$. Then $\delta p=-i\gamma ps$ and $\delta \mathbf{B}=i(K\mathbf{\xi }-\mathbf{B}s)$. Substitution givesIf $s=q=0$, the fluid displacement is perpendicular to both the wave vector and the magnetic field. Only magnetic tension restores it, andThis is the Alfvén wave. The Alfvén velocity is the field-directed propagation vector: $\omega =\pm \mathbf{k}\cdot {\mathbf{v}}_a$. The signed phase velocity normal to the wavefront is $\pm v_a\cos \theta \hat{\mathbf{k}}$; the group velocity is $\pm {\mathbf{v}}_a$. The distinction matters for oblique propagation. For generic directions the Alfvén wave has one transverse polarization, with $\delta \rho =\delta p=0$.

To obtain the other magnetohydrodynamic waves, take the scalar products with $\mathbf{k}$ and $\mathbf{B}$:The determinant condition for a nonzero pair $(s,q)$ isThus the fast magnetosonic wave and slow magnetosonic wave haveTheir fluid displacements lie in the plane of $\mathbf{k}$ and $\mathbf{B}$ and are generally compressive. In the fast magnetosonic wave, gas and magnetic pressure provide the stronger restoring combination; in the slow magnetosonic wave, their perturbations oppose one another. The phase speeds depend on direction. They satisfy $v_{\mathrm{sl}}\le \min (c_s,v_a)$ and $v_{\mathrm{f}}\ge \max (c_s,v_a)$. For parallel propagation the two speeds are $c_s$ and $v_a$, with a degeneracy between a transverse branch and the Alfvén wave. For perpendicular propagation $v_{\mathrm{f}}^2=c_s^2+v_a^2$ and $v_{\mathrm{sl}}=0$; the latter is a nonpropagating limiting disturbance. These special directions require interpreting the polarizations by continuity rather than assuming three distinct nonzero frequencies.

In a steady flow the ideal magnetohydrodynamic induction equation gives $\nabla \times (\mathbf{u}\times \mathbf{B})=0$. Since both vectors have only poloidal magnetic field components, write $\mathbf{u}\times \mathbf{B}=F_{\varphi }{\mathbf{e}}_{\varphi }$. Axisymmetry impliesFor a field regular on the axis, $C=0$; equivalently one may impose zero toroidal electromotive force. Under this physical condition, $\mathbf{u}$ is parallel to $\mathbf{B}$. Write $\rho \mathbf{u}=k\mathbf{B}$. Continuity and $\nabla \cdot \mathbf{B}=0$ then give $\mathbf{B}\cdot \nabla k=0$. The poloidal magnetic flux function has $\mathbf{B}\cdot \nabla \psi =0$, so locally on each connected magnetic surface,This regular-axis alignment of steady poloidal ideal flow needs its zero-circulation condition: axisymmetry alone does not imply this alignment on a domain excluding the axis. An explicit counterexample is available in a cylindrical annulus. Take nonzero constants $U,b$, positive ${\rho }_c,c_s$, andHere $\psi =-bz$, $\nabla \cdot \mathbf{B}=\nabla \times \mathbf{B}=0$ and $\mathbf{u}\times \mathbf{B}=Ub{\mathbf{e}}_{\varphi }/R$ has zero curl. The flow satisfies continuity, $p^{\prime }=-\rho {\Phi }^{\prime }$, and ${\nabla }^2\Phi =4\pi G\rho$. Its acceleration and Lorentz force density vanish, so it solves the steady equations while $\mathbf{u}$ and $\mathbf{B}$ are perpendicular. The nonzero toroidal circulation and the axial singularity explain why this counterexample is excluded from a regular polar accretion column. The remaining derivation uses the regular, aligned branch.

For a narrow stream tube, integrating continuity gives constant $A\rho |\mathbf{u}|$. On a nonzero-flow tube the magnetohydrodynamic mass loading $k$ is constant, so $A\rho |\mathbf{u}|=A|k||\mathbf{B}|$. Therefore $A|\mathbf{B}|$ is constant and $A\propto |\mathbf{B}{|}^{-1}$. This also follows directly from conserved magnetic flux through a flux tube.

For constant isothermal sound speed, $\nabla p/\rho =c_s^2\nabla \ln (\rho /{\rho }_0)$, with arbitrary reference density ${\rho }_0$. Project momentum along $\mathbf{B}$. The Lorentz force density has no component in that direction, and alignment implies $\mathbf{B}\cdot (\mathbf{u}\cdot \nabla )\mathbf{u}=\mathbf{B}\cdot \nabla (|\mathbf{u}{|}^2/2)$. Thus the isothermal magnetic Bernoulli integral isChanging ${\rho }_0$ only shifts the Bernoulli function by a constant. Since $k$ is constant along the magnetic field,Differentiating the isothermal magnetic Bernoulli integral therefore gives

In the slender polar accretion column, the dipolar flux-tube area is $A(z)=A_{\ast }(z/R_{\ast }{)}^3$. The decreasing axial field is a leading approximation, not an exactly solenoidal field $(0,0,B_z(z))$ throughout a cylinder. Indeed $\nabla \cdot \mathbf{B}={\partial }_z{B}_z\ne 0$ for that literal field. Near the axis a small radial component $B_R\simeq -(R/2){\partial }_z{B}_z=3R{B}_z/(2z)$ supplies the required radial divergence. Its magnitude is smaller by $R/z$, even though its divergence is leading order. Keeping this expanding flux tube while neglecting transverse forces yields the intended one-dimensional model.

Use the positive inward speed $v=-u_z$; the signed sonic velocity is $u_z=-c_s$. Conservation of mass and the Bernoulli equation giveEliminating ${\rho }^{\prime }$ yields the isothermal dipolar accretion equation,A smooth transonic branch requires both sides to vanish at its sonic point, soThis is the critical point of the flow equation; a generic subsonic solution need not cross it. A crossing in the exterior column requires $z_s\ge R_{\ast }$, and a sonic point strictly outside the star requires $z_s>R_{\ast }$. Differentiating the equation at the crossing gives $2(v_s^{\prime }{)}^2=3c_s^2/z_s^2$. For inward accretion that accelerates toward the star, $v_s^{\prime }=-\sqrt{3/2}{c}_s/z_s$.

The reservoir boundary condition fixes $ϵ=c_s^2\ln ({\rho }_{\infty }/{\rho }_0)$. At the sonic point, ${\Phi }_s=-3c_s^2$ and $v_s^2/2=c_s^2/2$, givingFinally use the sonic mass flux through the dipolar flux-tube area:This rate belongs to the smooth transonic branch of the idealized column, with $A_{\ast }$ the supplied total loaded area. The reservoir condition alone does not force every steady solution onto this branch. A real narrow dipolar column must match an outer flow where the slender approximation ceases to apply.

Use the ideal long-solenoid approximation: negligible pitch and end effects, vacuum permeability ${\mu }_0$, and negligible displacement current. By cylindrical symmetry the magnetic field is axial. An Ampère's law rectangular loop with one axial side inside and the other outside gives $B_{\mathrm{in}}-B_{\mathrm{out}}={\mu }_{0}NI$. A loop wholly outside shows the outside axial field is constant with radius; requiring it to vanish far away givesThe axial sign is determined by the winding direction and the right-hand rule. The small axial component of current from the spiral pitch is neglected in this idealization.

There are $Nl$ turns, and the magnetic flux through each is $\pi a^2{\mu }_{0}NI$. The flux linkage is $LI$, so the self-inductance isThe induced electromotive force opposes increasing current by Faraday's law. Together with Ohm's law, the circuit equation is $L\dot{I}+RI={\mathcal{E}}_0$. Solving this first-order ordinary differential equation with $I(0)=0$ givesThe time constant is $L/R$, and the limiting current is ${\mathcal{E}}_0/R$.

The magnetic field satisfies $\nabla \cdot \mathbf{B}=0$. Expanding the divergence of its Maxwell stress tensor givesFor the Newtonian gravitational field $\mathbf{g}=-\nabla \Phi$, Poisson equation for Newtonian gravity gives ${\partial }_j{g}_j=-4\pi G\rho$, and the gradient representation gives ${\partial }_j{g}_i={\partial }_i{g}_j$. ConsequentlyThe negative of the Newtonian gravitational stress tensor, in this force-stress convention, therefore supplies $\rho g_i$. Adding the pressure stress supplies $-{\partial }_{i}p$, soEach summand is symmetric. Absence of external gravitational sources is needed to represent the full gravitational force by this self-gravitating stress.

Integrate Poisson equation for Newtonian gravity through a narrow slab around the disk. The horizontal derivative contributions vanish as its thickness tends to zero, leaving the normal-derivative jumpThe even function symmetry of $\Phi$ makes the derivatives opposite, soIn the current-free simply connected upper half-space, Ampère's circuital law gives $\nabla \times \mathbf{B}=0$ and permits a magnetic scalar potential. Rescale it so that $\sqrt{4\pi G/{\mu }_0}\mathbf{B}=-\nabla {\Phi }_M$. The divergence-free condition makes ${\Phi }_M$ satisfy Laplace's equation, withCompare with the gravitational jump condition. Subject to the same isolated-field boundary condition at infinity, ${\Phi }_M$ is the harmonic potential of the effective surface densityThis gravity-equivalent magnetic surface density can have either sign; it is a mathematical representation of the exterior magnetic field, not physical negative mass. An imposed nondecaying field would require additional boundary data and would not be fixed by the disk surface density alone.

The Lorentz force in the uniform magnetic field giveswhere the last step is the vector triple product identity. HenceOn the other hand, $|\mathbf{r}\times \mathbf{B}{|}^2=B^2{r}^2-(\mathbf{r}\cdot \mathbf{B}{)}^2$, soThe derivatives cancel, proving the magnetic angular-momentum invariantThe magnetic Lorentz force is perpendicular to velocity, so its power is zero:Thus the kinetic energy is constant.

For the alternative kinetic energy expression, work on an interval where $r>0$ so that the radial unit vector $\mathbf{u}=\mathbf{r}/r$ is defined. Differentiating $\mathbf{u}\cdot \mathbf{u}=1$ once and twice givesSince $\mathbf{v}=\dot{r}\mathbf{u}+r\dot{\mathbf{u}}$, it follows that $\mathbf{u}\cdot \mathbf{v}=\dot{r}$ and $v^2={\dot{r}}^2+r^2|\dot{\mathbf{u}}{|}^2$. Therefore the unit-direction kinetic-energy identity isThis identity holds for any twice differentiable trajectory away from the origin, not only for motion in a magnetic field.