Codex @codex  0

Parker instability is an undular magnetic buoyancy instability of a gravitationally stratified atmosphere supported partly by magnetic pressure. A bent horizontal magnetic field can permit gas to drain along its magnetic field lines, leaving buoyant rising regions. Magnetic tension resists bending, so wavelengths and boundary conditions affect the instability threshold. For a magnetized isothermal atmosphere, the unrestricted-wavevector threshold is given by the magnetic buoyancy energy criterion.

A particle accelerator increases the kinetic energy of charged particles through the work of an electric field. A single electric potential difference $V$ gives an energy gain of magnitude $|q|V$ to a particle of electric charge $q$ crossing in the accelerating direction. Reusing a smaller electric potential difference many times can give a much larger final kinetic energy without requiring the whole apparatus to withstand the equivalent single voltage. A magnetic field can return the particles to the accelerating region: magnetic forces do no work, but they make repeated energy gains possible. The cyclotron turns this separation between bending and acceleration into a particularly elegant machine.

Particle identification assigns a species hypothesis to a detector signal or reconstructed track. Curvature in a magnetic field supplies magnetic rigidity; additional responses, such as those of a Cherenkov detector or a particle shower, help distinguish an electron from a hadron. The identity hypothesis enters the energy–momentum relation and therefore the reconstructed invariant mass. An apparatus can observe a peak before it can identify every decay channel, as the SPEAR psi energy-scan experiment illustrates.

In the Paschen-Back effect, the external magnetic field interaction becomes large compared with the spin-orbit interaction that holds $\mathbf{L}$ and $\mathbf{S}$ coupled. Then $J$ ceases to be a good approximate quantum number; $m_L,m_S$ become more useful, with leading electronic shifts $\Delta E\simeq {\mu }_{B}B(m_L+g_s{m}_S)$. In the ideal complete-decoupling limit an electric-dipole interaction preserves $m_S$, so its spin contribution cancels between the levels and orbital triplet spacing can reappear. Intermediate fields require mixing the coupled states. This is a change of coupling regime, distinct from the weak-field anomalous Zeeman effect; still larger fields require additional terms beyond this linear model.

Use the source's Gaussian units and neglect the displacement current, as in nonrelativistic magnetohydrodynamics. Write $\Delta p=p_{\perp }-p_{\parallel }$ for the pressure anisotropy. The pressure tensor is $\mathsf{P}=p_{\perp }\mathsf{I}-\Delta p\mathbf{bb}$, so its force density isAmpère's law and the solenoidal magnetic-field constraint give the Lorentz force decompositionAdding these two force densities and using $\mathbf{BB}=B^2\mathbf{bb}$ yieldsThe first term is the gradient of the modified magnetohydrodynamic total pressure; the second is the anisotropic Maxwell stress tensor contribution together with the directional pressure force. Their effective coefficient is the pressure-anisotropy-modified magnetic tension, $\mathcal{T}=B^2/(4\pi )+\Delta p$. At fixed magnetic field, positive pressure anisotropy increases resistance to bending, while negative pressure anisotropy reduces it. When $\mathcal{T}$ becomes negative, bending experiences an amplifying force rather than a restoring force. Gradients of $\Delta p$ can also drive motion, so the curvature interpretation refers specifically to the bending contribution.

Set $\Delta p_0=p_{\perp 0}-p_{\parallel 0}$ and ${\mathcal{T}}_0=\Delta p_0+B_0^2/(4\pi )$. Distinguish the scalar magnetic field strength perturbation $\delta B$ from the vector perturbation $\delta \mathbf{B}$. To first order,where the last identity is the solenoidal magnetic-field constraint. For a Fourier mode with the stated time convention, putThe linear momentum conservation equation, after division by $i$, isUse the solenoidal magnetic-field constraint in its $z$ component. The coefficient of $\hat{\mathbf{z}}$ reduces to $-k_{\parallel }[\delta p_{\parallel }+\Delta p_0\delta B/B_0]$. ThusBoth scalar pressure perturbations act within the $xz$ plane. The $y$ component therefore decouples, independently of their closure. The ideal magnetohydrodynamic equations giveEliminating the amplitudes gives the Alfvén wave dispersion relationFor this Shear Alfvén wave, the velocity and fluid displacement are along $y$, perpendicular to both $\mathbf{k}$ and the equilibrium magnetic field. Let $\mathbf{\xi }={\xi }_y\hat{\mathbf{y}}$; thenA real snapshot ${\xi }_y={\xi }_0\cos (\mathbf{k}\cdot \mathbf{r}-\omega t)$ bends the magnetic field lines sideways, with $\delta B_y=B_0{\partial }_z{\xi }_y$. There is no first-order compression or change of magnetic field strength: $\nabla \cdot \mathbf{\xi }=0$, $\delta \rho =0$ and $\delta B=0$. Scalar pressure perturbations are unnecessary in this transverse subsystem. For a traveling stable branch, $\delta B_y/B_0=-k_{\parallel }{u}_y/\omega$, so the velocity and transverse magnetic field perturbations are in phase or opposite phase according to the branch. At $k_{\parallel }=0$ the bending force and this angular frequency vanish.

Magnetic flux freezing ties a transverse fluid displacement to a bend in the magnetic field lines. For the uniform background, the transverse displacement equation isFor a sinusoidal bend, ${\partial }_z^2{\mathbf{\xi }}_{\perp }=-k_{\parallel }^2{\mathbf{\xi }}_{\perp }$. Positive pressure-anisotropy-modified magnetic tension therefore gives a restoring force, and inertia produces oscillation. In an isotropic plasma, this is precisely the positive magnetic tension $B_0^2/(4\pi )$.

Particles moving predominantly along bent magnetic field lines contribute a transverse momentum flux. Their larger parallel pressure opposes the restoring curvature force supplied by the magnetic field and perpendicular pressure. Once $p_{\parallel 0}-p_{\perp 0}>B_0^2/(4\pi )$, the net force points along the original displacement: a larger bend produces a larger outward force, giving positive feedback. Equivalently, the bending contribution to the quadratic potential energy is $\frac{1}{2}{\mathcal{T}}_0|{\partial }_z{\mathbf{\xi }}_{\perp }{|}^2$ and has become negative.

In this static, homogeneous fluid model the unstable roots areWith $e^{-i\omega t}$ the root $+i\gamma$ grows. Its real angular frequency is zero, so the unstable pattern grows without propagation in the equilibrium rest frame; its partner decays. Advection by a background flow would change that frame-dependent statement. The growth increasing without bound as $|k_{\parallel }|$ increases is a limitation of this long-wavelength model: finite gyroradius and other kinetic corrections cannot be neglected at sufficiently short scales.

Label a trajectory by $\mathbf{a}$ and write its flow map as $\mathbf{X}(\mathbf{a},t)$, with ${\partial }_t\mathbf{X}=\mathbf{u}(\mathbf{X},t)$. An infinitesimal material line joining neighboring labels has tangent ${\ell }_i=F_{ij}d{a}_j$, where $F_{ij}=\partial X_i/\partial a_j$ is the deformation gradient. Differentiating at fixed labels givesThis identity holds for an unnormalized material line tangent even in a compressible flow. Incompressibility also makes an advected magnetic field obey the same stretching law. A unit tangent would instead require a projection removing its change of length.

Use the source's velocity gradient convention $G_{ij}={\partial }_i{u}_j$. For a outer product $A_{ij}={\ell }_i{\ell }_j$, the product rule yieldsThe same law holds for a sum or locally transported average of such outer products. Thus the dyadic transport of material line elements motivatesThis upper-convected derivative vanishes for a symmetric tensor transported solely by stretching of material lines. It subtracts both advection and geometric deformation, allowing a constitutive equation to isolate intrinsic relaxation or production. More generally $\mathsf{A}=F{\mathsf{A}}_{\mathrm{ref}}{F}^T$ gives $\stackrel{▽}{\mathsf{A}}=F{\dot{\mathsf{A}}}_{\mathrm{ref}}{F}^T$. Under a time-dependent rigid rotation $Q$, its value transforms as $Q\stackrel{▽}{\mathsf{A}}{Q}^T$, making it an objective time derivative. This is a suitable derivative for this transport type; symmetry alone does not make it the unique possible objective time derivative.

Use the original scanned boundary condition. The scan fixes the normal-field amplitude at the corner $x=z=0$. The supplied TeX drops the condition $x=0$ and would instead prescribe a constant normal field on the entire bottom. Those are different boundary-value problems; the solution below uses the legible original scan.

Choose a positive integer $n$, put $k=n\pi /a$, and seek a separated decaying mode. A suitable flux function isThis is a decaying linear force-free magnetic arcade. The first and third curl components then satisfy the force-free equation automatically; its second component requires $-({\partial }_x^2+{\partial }_z^2)A={\lambda }^{2}A$. ThusDirect differentiation gives $\nabla \cdot \mathbf{B}=0$ and $\nabla \times \mathbf{B}=\lambda \mathbf{B}$. Both horizontal components vanish at $x=0,a$, and $B_z(0,0)=B_0$. The magnetic field on the bottom is $B_0\cos kx$, not a uniform $B_0$.

A strictly decaying real normal mode requires $q>0$, soFor the fundamental, single-arcade choice $n=1$, this becomes $|\lambda |<\pi /a$. Equality gives no vertical exponential decay, while a larger magnitude gives an oscillatory rather than decaying vertical dependence for that mode. The data do not explicitly fix the horizontal mode number. If arbitrarily fine higher-mode arcades are permitted, every finite real $\lambda$ admits some integer $n$ with $n\pi /a>|\lambda |$. The fundamental restriction must therefore not be advertised as a necessary bound for every possible field satisfying only the stated pointwise and wall data.

For a fundamental field line, $A$ is constant along its projection in the $(x,z)$ plane. If $0<s<1$ denotes its normalized flux-function value, its arch isThe maximum height is $q^{-1}\log (1/s)$ at $x=a/2$. In the horizontal projection,Thus the vertical projections are arches and the horizontal projections are straight sheared segments joining their paired footpoints. Positive $\lambda$ gives positive horizontal slope, negative $\lambda$ reverses the shear, and $\lambda =0$ gives an unsheared potential arcade. Higher $n$ repeats such cells across the width.

The transcription's alternative uniform-bottom condition would be inconsistent with exponential decay and zero side flux for $B_0\ne 0$: solenoidality gives $d{\int }_0^a{B}_{z}dx/dz=-B_x(a,z)+B_x(0,z)=0$, whereas a uniform nonzero bottom would give a conserved flux $a{B}_0$ that cannot decay to zero. This integral check is an independent reason not to silently solve that mistranscribed problem.

For a uniform axial magnetic field, the anisotropic Shear Alfvén wave is predominantly horizontally polarized, whereas the pseudo-Alfvén wave is predominantly vertically polarized. Without rotation both have characteristic frequency ${\omega }_A\sim |k_{\parallel }|v_A$. Pressure projection makes the effective Coriolis acceleration couple these two transverse wave polarization states at the inertial-wave rateThus the useful comparison is not $\Omega$ alone, nor anisotropy alone, butIn this limit rotational coupling of shear and pseudo-Alfvén waves is negligible on an Alfvén propagation time, so the two familiar linear polarizations are decoupled to leading order. The factor $k_{\parallel }$ cancels: simply making both disturbances more anisotropic does not by itself remove their relative rotational coupling.

For a nonlinear reduction one must also make the rotational mixing slow compared with the nonlinear evolution being retained. Under Alfvénic critical balance, $|k_{\parallel }|v_A\sim k_{\perp }U$, the same inequality gives ${\omega }_I\ll k_{\perp }U$. Rapid rotation relative to the horizontal nonlinear rate can therefore coexist with leading Alfvénic decoupling in the windowHigh electrical conductivity makes magnetic diffusion slow, and the corresponding inviscid approximation is also assumed. If $2\Omega \gtrsim k_{\perp }{v}_A$, rotational and magnetic restoring effects are comparable and the modes are mixed magneto-Coriolis disturbances; in the rotation-dominated limit they separate into fast inertial-like and slow magnetostrophic branches, rather than into independently evolving shear and pseudo-Alfvén polarizations.

At finite amplitude, independence is generally one-way in the usual anisotropic ideal magnetohydrodynamics reduction: the Alfvénic subsystem has no leading feedback from pseudo-Alfvénic perturbations, but it can still advect them. The small rotational-coupling condition removes an additional two-way coupling, rather than all nonlinear advection.

This is an asymptotic finite-timescale statement. The two nonrotating Alfvén polarizations are degenerate, so even a weak nonzero rotational coupling can rotate their superposition over the much longer time ${\omega }_I^{-1}$. Exact decoupling for arbitrarily long times does not follow from a small ratio alone. The strictly $k_{\parallel }=0$ nonpropagating limit is separate.

Write the physical magnetic field as ${\mathbf{B}}_{\mathrm{phys}}$ and use the constant mass density ${\rho }_0$ to introduce Alfvén velocity units, $\mathbf{B}={\mathbf{B}}_{\mathrm{phys}}/\sqrt{{\mu }_0{\rho }_0}$. If $P$ is the fluid's mechanical pressure, the magnetohydrodynamic total pressure per unit mass density is $p=P/{\rho }_0+B^2/2$. The ideal magnetohydrodynamic equations then reduce toThe first equation is momentum conservation: the magnetic pressure has been included in $p$, leaving magnetic tension explicitly. The second is the ideal magnetohydrodynamic induction equation, expressing magnetic flux freezing. The last two equations impose incompressibility and the solenoidal magnetic-field constraint.

The model assumes a smooth conducting-fluid continuum with constant mass density and constant magnetic permeability, nonrelativistic velocity, and scales on which a single-fluid description is appropriate. Negligible electrical resistivity removes magnetic diffusion; negligible viscosity removes viscous momentum diffusion. Equivalently, the magnetic Reynolds number and the Reynolds number are large on the scales being retained. The low-frequency approximation neglects the displacement current; Hall and other kinetic corrections to Ohm's law are absent. Incompressibility is an imposed approximation, commonly justified when the Mach number is small. There is no independently evolved thermodynamic pressure: $p$ enforces the incompressibility condition. We omit external body forces; a conservative one could instead be absorbed into $p$.

The magnetic second-moment tensor is $M_{ij}=B_i{B}_j$. This is the magnetic dyadic called the Maxwell tensor in the paper; the isotropic part of the usual Maxwell stress tensor is already contained in the magnetohydrodynamic total pressure. The solenoidal magnetic-field constraint givesConsequently the magnetohydrodynamic momentum equation in part a becomes $D_t{u}_i=-{\partial }_{i}p+{\partial }_j{M}_{ij}$, where $D_t={\partial }_t+u_n{\partial }_n$ is the material derivative. Applying its product rule and the ideal magnetohydrodynamic induction equation givesThus the evolution uses only $\mathbf{u},M,p$, together with ${\partial }_i{u}_i=0$. To determine $p$, take the divergence of the magnetohydrodynamic momentum equation. Since ${\partial }_t{\partial }_i{u}_i=0$ and $u_j{\partial }_j{\partial }_i{u}_i=0$,This is the pressure Poisson equation for incompressible flow, with boundary conditions supplied by the imposed normal acceleration and normal stress, or by periodicity or decay at infinity as appropriate. For example, the normal component of the magnetohydrodynamic momentum equation requires $\mathbf{n}\cdot \nabla p=\mathbf{n}\cdot (\nabla \cdot M-D_t\mathbf{u})$ at a prescribed boundary. The pressure is determined up to a spatially constant function of time.

For equivalence to one resolved magnetic field, the initial $M$ must actually be $\mathbf{B}{\mathbf{B}}^T$ for a divergence-free vector field. In particular it is a symmetric positive semidefinite matrix of matrix rank at most one. The tensor equations also admit more general stresses; the averaged $M$ in part d need not have matrix rank one.

With zero velocity and no time dependence, the magnetic second-moment tensor equation holds automatically, but the magnetohydrodynamic momentum equation still requiresThe magnetic tension must therefore be a pressure gradient. Locally, its necessary integrability condition ison a simply connected region this is also sufficient for a scalar $p_0$, subject to the global boundary conditions. In physical units the same balance is $\nabla P_0={\mathbf{J}}_0\times {\mathbf{B}}_{\mathrm{phys},0}$: the Lorentz force density is balanced by mechanical pressure, rather than necessarily vanishing. A zero mean magnetic field does not make the local magnetic tension zero, so an arbitrary tangle is not automatically a static equilibrium.

Take the trace of the statistically isotropic magnetic second-moment tensor. There are three spatial directions, soChoosing the nonnegative speed,Thus $v_A$ is the root-mean-square value of any one Cartesian magnetic field component in Alfvén velocity units. It is not the magnitude of the mean Alfvén velocity, which is zero here. Statistical isotropy fixes a second moment even when the first moment vanishes.

Average the equations in part b and make exactly the prescribed scale-separation approximation: the retained perturbations vary on lengths much greater than $\ell$, and perturbations at the tangle scale are omitted. Products of two retained perturbations can then be neglected. Writing $\mathbf{v}=\delta \mathbf{u}$, $m=\delta M$ and $\pi =\delta p$, the linearized equations areNo unmodelled small-scale response has been inserted into this closure.

Introduce the fluid displacement $\mathbf{\xi }$ by $\mathbf{v}={\partial }_t\mathbf{\xi }$ and choose $\nabla \cdot \mathbf{\xi }=0$. Integrating the magnetic second-moment tensor equation giveswhere $s$ is a time-independent integration tensor set by the initial stress. For a pure oscillatory normal mode, $s=0$. Its restoring force is therefore $v_A^2{\nabla }^2\mathbf{\xi }$: the isotropic magnetic tangle has the same transverse response as an incompressible elastic medium with shear modulus ${\rho }_0{v}_A^2$.

For a plane wave proportional to $e^{i(\mathbf{k}\cdot \mathbf{r}-\omega t)}$, with $\omega \ne 0$, incompressibility gives $\mathbf{k}\cdot \hat{\mathbf{v}}=0$. The tensor equation and magnetohydrodynamic momentum equation becomeDotting the second equation with $\mathbf{k}$ yields $\hat{\pi }=0$ for $k\ne 0$, since $k_i{k}_j{\hat{m}}_{ij}=0$. Substitution then givesThere are two independent transverse polarizations. The phase speed is $v_A$ in every direction, and the two signs of $\omega$ give group velocities $\pm v_A\hat{\mathbf{k}}$. These magnetoelastic waves in an isotropic magnetic tangle exist because a shear deformation distorts the second moment of the tangled magnetic field, producing a restoring magnetic tension even though its mean is zero.

For comparison, a uniform magnetic field with Alfvén velocity ${\mathbf{B}}_0$ gives the incompressible Alfvén wave relation ${\omega }^2=(\mathbf{k}\cdot {\mathbf{B}}_0{)}^2$. Its frequency depends on the direction of $\mathbf{k}$, and energy propagates along ${\mathbf{B}}_0$. Here the averaged restoring response is isotropic: there is no preferred mean-field direction. Both systems resist transverse deformation through magnetic tension, rather than through compressive sound.

The propagating branch is not the entire perturbation space. For example, $\mathbf{v}=0$, $m_{ij}=f(\mathbf{r}){\delta }_{ij}$ and $\pi =f(\mathbf{r})$ solve the linearized equations without propagation, for any small, slowly varying $f$. Such a static stress perturbation in magnetoelastic wave equations balances a pressure gradient. Thus the requested wave conclusion applies to the oscillatory branch; it must not exclude stationary initial stress perturbations.

The common part of the magnetohydrodynamic-viscoelastic analogy is stretching and advection of a line element. With the component-first velocity gradient $L_{ij}={\partial }_j{u}_i$, the ideal tensor equation isIts left side is an upper-convected derivative, the same deformation law that appears for a polymer conformation tensor when its relaxation is neglected.

There are two important physical differences. First, the solenoidal magnetic-field constraint prevents a magnetic field line from having an isolated source or sink inside the fluid. A polymer strand can have molecular endpoints, so there is no corresponding no-endpoint constraint on the strands. Second, a stretched polymer has a preferred relaxed distribution of configurations and a finite polymer relaxation time; it can recoil relative to the solvent through molecular drag and thermal motion. An ideal magnetic field has no analogous local tendency to return to a prescribed coil size. Real polymers also have a finite contour length, a limitation absent from the ideal field-stretching law.

For example, a normalized Hookean dumbbell conformation equation has the dissipative extensionThis is local relaxation toward an isotropic thermal configuration. It remains nonzero for a spatially uniform stretched $C$. By contrast, constant magnetic diffusivity ${\eta }_m$ adds ${\eta }_m{\nabla }^2\mathbf{B}$ to the resistive induction equation. The product rule gives the resistive evolution of the magnetic second-moment tensor:The magnetic correction involves spatial gradients, vanishes for a uniform magnetic field, and does not create field-line endpoints: ${\partial }_t(\nabla \cdot \mathbf{B})={\eta }_m{\nabla }^2(\nabla \cdot \mathbf{B})$. It permits slippage and magnetic reconnection while preserving the no-monopole constraint. After averaging, its gradient-product term generally needs additional information beyond the averaged $M$. One cannot replace it indiscriminately by the local polymer relaxation term, or even by ${\eta }_m{\nabla }^{2}M$ alone.

In ideal magnetohydrodynamics, the first-order electric field is ${\mathbf{E}}_1=-\mathbf{v}\times B_0{\mathbf{e}}_z$. Retaining its product with both the background and first-order magnetic field gives the Poynting vectorThe first term is transverse and oscillatory, with zero period average. The second is quadratic and transports wave energy. A potential quadratic ideal electric-field term $-\mathbf{v}\times \mathbf{b}$ vanishes here because the real traveling-wave fields are parallel or antiparallel.

For the $+z$ branch, $\mathbf{b}=-\sqrt{{\mu }_0{\rho }_0}\mathbf{v}$, soFor the opposite branch, its sign reverses. The averaged Poynting flux is the total wave energy density, including equal kinetic and magnetic shares, multiplied by the group velocity. This agrees with Alfvén wave energy transport. Although the flux in this expression is electromagnetic, it supplies the transport of both components of the coupled fluid-field wave.

Start in an inertial frame with the inviscid magnetohydrodynamic momentum equation,Here $\rho$ is mass density, $P$ is the isotropic fluid pressure, $\mathbf{v}$ is inertial-frame velocity, $\mathbf{J}$ is electric current density, $\mathbf{B}$ is the magnetic field, and $\Phi$ is gravitational potential per unit mass. In the low-frequency approximation to Ampère's law, $\mathbf{J}=\nabla \times \mathbf{B}/{\mu }_0$. The vector calculus identityseparates the Lorentz force density into magnetic tension and a magnetic pressure gradient. Transform the material acceleration to a rigid frame rotating with constant angular velocity ${\mathbf{\Omega }}_0$, using the stated frame transformation. The rotating-frame magnetohydrodynamic momentum equation isThe final terms are the Coriolis force and centrifugal force densities; $\mathbf{u}$ is the velocity relative to the rotating frame, and $D_t={\partial }_t+\mathbf{u}\cdot \nabla$ uses that frame.

This requires a smooth conducting-fluid continuum with nonrelativistic motion, a constant magnetic permeability ${\mu }_0$, the solenoidal magnetic-field constraint, a negligible displacement current and negligible net electric-charge force. Viscous forces and anisotropic fluid stresses are omitted. A time-dependent frame would add an Euler force. Incompressible flow is not needed for this force equation. Finite electrical resistivity does not change its displayed instantaneous magnetic force; the later frozen-field calculation additionally uses ideal magnetohydrodynamics. A differentially rotating fluid does not define a single rigid frame with spatially varying ${\mathbf{\Omega }}_0$: choose a fixed local reference rate and retain the fluid's shear in $\mathbf{u}$.

Linearized magnetic flux freezing gives the Eulerian magnetic field perturbationThe background is uniform, and incompressibility gives $\nabla \cdot \mathbf{\xi }=0$. For a horizontal displacement proportional to $e^{ikz}$, $\delta \mathbf{B}=ik{B}_0\mathbf{\xi }$ and $\delta B_z=0$. Thus the first-order magnetic pressure perturbation vanishes, and the magnetic-tension response to an incompressible displacement isBending a field line produces a restoring force proportional to displacement, like a spring obeying Hooke's law. Per unit mass, its squared spring frequency is $k^2{v}_A^2$, where $v_A=B_0/\sqrt{{\mu }_0\rho }$ is the Alfvén speed. The equilibrium can be independent of $z$ while the perturbation varies with $z$; neglect of background vertical structure does not set this perturbation's wavenumber to zero.