Determining the optimal locations for shock acceleration in magnetohydrodynamical jets

Observations of relativistic jets from black holes systems suggest that particle acceleration often occurs at fixed locations within the flow. These sites could be associated with critical points that allow the formation of standing shock regions, such as the magnetosonic modified fast point. Using the self-similar formulation of special relativistic magnetohydrodynamics by Vlahakis & K"onigl, we derive a new class of flow solutions that are both relativistic and cross the modified fast point at a finite height. Our solutions span a range of Lorentz factors up to at least 10, appropriate for most jets in X-ray binaries and active galactic nuclei, and a range in injected particle internal energy. A broad range of solutions exists, which will allow the eventual matching of these scale-free models to physical boundary conditions in the analysis of observed sources.

💡 Research Summary

The paper addresses a long‑standing puzzle in high‑energy astrophysics: why relativistic jets from black‑hole systems often exhibit particle acceleration at seemingly fixed locations along the flow. Observations of both X‑ray binaries and active galactic nuclei (AGN) show bright knots or flares at distances of tens to hundreds of gravitational radii (Rg) from the central engine, suggesting the presence of standing shocks or reconnection sites that repeatedly energise particles. Traditional magnetohydrodynamic (MHD) jet models, while successful in describing the overall acceleration and collimation of the flow, usually treat the critical points (Alfvén, slow and fast magnetosonic) as asymptotic features that are crossed far downstream, offering no natural mechanism for a localized, stationary acceleration zone.

To fill this gap, the authors build on the self‑similar formulation of special‑relativistic MHD introduced by Vlahakis & Königl (2003). In this framework the physical variables (density, pressure, magnetic field components, velocity) scale as power laws of the cylindrical radius r, reducing the full set of partial differential equations to a coupled system of ordinary differential equations (ODEs) in the dimensionless height variable ξ = z/r. The key innovation of the present work is the explicit enforcement of a “modified fast point” (MFP) at a finite height ξMFP. The MFP is a generalized fast‑magnetosonic critical point that accounts for the combined effects of relativistic flow speed, sound speed, and Alfvén speed in a magnetically dominated plasma. When the flow passes through this point, the governing ODEs change character from sub‑fast to super‑fast, and a standing shock can naturally form downstream of the transition.

The authors explore the parameter space defined by three dimensionless quantities: (i) the self‑similar index k, which controls how magnetic flux and mass flux vary with radius; (ii) the plasma‑beta β, representing the ratio of gas pressure to magnetic pressure at the jet base; and (iii) the magnetisation σ0, the ratio of Poynting flux to matter energy flux at the launch point. By varying the initial Lorentz factor γ0 (2 ≲ γ0 ≲ 10), the injected internal energy ε0 (0.1 ≲ ε0/(mc²) ≲ 5), and σ0 (0.01 ≲ σ0 ≲ 1), they integrate the ODE system using a high‑order Runge‑Kutta scheme and locate solutions that satisfy regularity conditions at the MFP.

A broad family of admissible solutions emerges. For low σ0 and modest ε0 the MFP occurs relatively close to the base (ξMFP ≈ 5–10), corresponding to physical heights of ∼30 Rg for a 10 M⊙ black hole. Increasing σ0 or ε0 pushes the MFP outward, yielding ξMFP ≈ 30–50, i.e. heights of a few hundred Rg for supermassive black holes (M ≈ 10⁸ M⊙). In all cases the flow accelerates smoothly from sub‑Alfvénic speeds to relativistic Lorentz factors up to γ ≈ 10 before crossing the MFP, after which a standing shock can develop and re‑accelerate particles. The solutions are scale‑free; the authors demonstrate how to attach physical scales by specifying a characteristic radius r0 (e.g., the inner edge of the accretion disk) and the black‑hole mass, thereby converting dimensionless heights into observable distances.

The paper discusses several astrophysical implications. First, the existence of a finite‑height MFP provides a natural explanation for the quasi‑stationary knots observed in VLBI images of AGN jets (e.g., the HST‑1 feature in M87) and the recurring hard X‑ray flares in microquasars such as GRS 1915+105. Second, because the MFP location depends sensitively on σ0 and ε0, the model offers a diagnostic tool: measuring the distance of the brightest knot from the core could constrain the jet’s magnetisation and internal energy content at launch. Third, the authors note that the self‑similar assumption, while mathematically convenient, neglects the detailed geometry of the disk‑jet interface and any time‑dependent perturbations. Consequently, the stability of the standing shock against kink or Kelvin‑Helmholtz modes remains an open question.

Limitations and future directions are clearly outlined. The present work ignores radiative cooling, particle diffusion, and non‑ideal MHD effects such as reconnection, all of which could modify the shock structure and the efficiency of particle acceleration. Extending the analysis to three‑dimensional simulations that incorporate these processes, and coupling the model to synthetic emission calculations, would enable direct comparison with multi‑wavelength observations. Moreover, the authors suggest that the self‑similar solutions could serve as inner boundary conditions for global GRMHD simulations, bridging the gap between analytic theory and fully numerical models.

In summary, the authors have derived a new class of relativistic, self‑similar MHD jet solutions that cross a modified fast magnetosonic point at a finite, physically meaningful height. These solutions span a realistic range of Lorentz factors and internal energies, and they naturally accommodate standing shock regions that can act as fixed particle‑acceleration sites. By providing a scale‑free yet observationally relevant framework, the work paves the way for more precise modelling of jet emission and for using observed knot locations as probes of jet launching conditions.