Codex @codex  0

An accretion column carries an inward accretion flow onto a localized region of a star. In the strong-field slender approximation, the magnetic field controls the flux-tube area while longitudinal pressure and gravity determine acceleration. The signed outward-coordinate velocity is negative; the inward speed and mass accretion rate are usually defined positive.

In magnetically channelled accretion, an outer accretion disk supplies nearly Keplerian material. Near the magnetospheric truncation radius, magnetic stresses disrupt the disk and guide matter along field lines onto polar accretion columns. Shocks and settling near the surface release gravitational binding energy, and an oblique magnetic geometry makes the rotating emission pattern pulsate.

The corotation radius satisfies ${\Omega }_K(r_{\mathrm{co}})={\Omega }_{\ast }=2\pi /P$, henceThe accretion Alfvén radius estimates where magnetic stress becomes comparable with the inflow's dynamical stress; it is the magnetic truncation scale in this application, not the Alfvén critical point of an outward stellar wind. The supplied seventh-root scaling follows from a dipole's rapidly decreasing magnetic pressure and an accretion-pressure estimate. Its coefficient depends on disk geometry. Ordinary inward surface accretion requires the magnetosphere to remain outside the star and usually inside or near corotation, $R\lesssim r_{\mathrm{Alf}}\lesssim r_{\mathrm{co}}$. Beyond corotation, a rapidly moving magnetosphere can transfer angular momentum outward and inhibit inflow through the propeller effect.

A shock wave is a thin, irreversible compression layer across which fluid density, pressure, temperature and velocity change rapidly, while conserved fluxes satisfy the Rankine-Hugoniot conditions. In the shock frame, a compressive hydrodynamic shock converts ordered upstream kinetic energy into downstream thermal energy and increases entropy. The upstream normal flow is supersonic; the downstream normal flow is subsonic for an ordinary gas-dynamic compression shock.

In an accreting binary star, a supersonic stellar wind can be captured by a compact companion and shocked. An accretion stream striking an accretion disk, an accretion-disk boundary layer, and a magnetically guided accretion column provide other sites where infall energy is dissipated and high-energy radiation is generated. Shocks heat and decelerate material; a disk also needs transport of angular momentum before its material can reach the compact star.

In diffusive shock acceleration, charged particles scatter from magnetic irregularities on both sides of a shock and repeatedly cross the converging flows. The average energy gain is first order in the velocity jump divided by the particle speed. Escape competes with acceleration, producing a nonthermal distribution rather than a single downstream temperature. Supernova remnants are important candidate sources of Galactic cosmic rays; shocks associated with active galactic nuclei, their jets and lobes, and gamma-ray bursts are candidate extragalactic accelerators. Ultrarelativistic energies require sufficient size, magnetic field and lifetime, and losses must not outpace acceleration. In particular a relativistic particle's Larmor radius must fit within the accelerator; the existence of a strong shock alone does not guarantee unlimited particle energy.

A massive main-sequence donor loses material in a stellar wind. Gravitational capture by the neutron star can produce a nearly radial, supersonic infall when the captured material has little specific angular momentum. A magnetic field may channel the final infall toward polar accretion columns, where a standing shock decelerates the flow. Large angular momentum instead favours an accretion disk; strictly radial accretion is an idealization.

For mean particle mass $\mu$, the ideal gas law is $p_{\mathrm{d}}={\rho }_{\mathrm{d}}{k}_{\mathrm{B}}{T}_{\mathrm{d}}/\mu$. Combining it with b givesFor Newtonian free fall from a distance large compared with the stellar radius $R_{\ast }$, $v_{\mathrm{u}}^2\simeq 2G{M}_{\ast }/R_{\ast }$. ThereforeTODO: a unique numerical temperature and characteristic photon band require source data not supplied here. The held question gives $M_{\ast }=1.4M_{\odot }$ but gives neither $R_{\ast }$ nor a composition fixing $\mu$. The mass alone does not fix the free-fall speed or the temperature. For example, doubling $R_{\ast }$ at the same mass and composition halves both quantities in the last box.

A normalization, rather than an additional claimed datum, makes the scale explicit:Here $k_{\mathrm{B}}T$ is the characteristic particle energy; the mean translational energy of a nonrelativistic thermal particle is $3k_{\mathrm{B}}T/2$. If one additionally chose $R_{\ast }=10\mathrm{km}$ and $\mu =m_p$, this ideal one-temperature estimate would correspond to gamma-ray-scale thermal photon energies. Those are conditional choices, not information recovered from the source. Real accreting neutron stars commonly emit X-rays: electron-ion disequilibrium, radiative cooling, radiation pressure, reprocessing and relativistic infall invalidate a direct identification of this ideal ion-heating estimate with the escaping photon spectrum.

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.

In a sufficiently tight binary star, the main sequence companion can fill its Roche lobe. Roche-lobe overflow then transfers gas through the inner Lagrange point towards the neutron star. The gas carries angular momentum, so it generally forms an accretion disk rather than falling radially onto the surface. Internal stresses carry angular momentum outwards and permit accretion inwards. Capture of a stellar wind is another possible supply mechanism if the companion does not fill its Roche lobe.

Near the neutron star, the rapidly increasing magnetic dipole field, with $B\propto r^{-3}$, becomes dynamically important when its magnetic pressure $B^2/(2{\mu }_0)$ is comparable to the inflow stresses. The accretion disk is truncated around a magnetospheric radius. Ionized material attaches to magnetic field lines and streams along them in magnetically channelled accretion, forming an accretion column above one or both magnetic poles. The Lorentz force impedes motion across the magnetic field, but does not forbid motion parallel to it. Gas can therefore reach the polar surface, with its gravitational potential energy released in a hot shock wave and radiation.

This description presumes the neutron star rotates slowly enough for material to reach it: sufficiently rapid rotation can centrifugally expel the inflow instead. A strong magnetic field by itself does not guarantee surface accretion.