Spine–Layer Jet Models

• Spine–Layer models are dual-component jet structures, where a fast, magnetically dominated spine is surrounded by a slower, matter-rich sheath.
• They explain observed features like limb brightening, broad spectral distributions, and distinct polarization patterns via differential beaming and shear acceleration.
• Supported by GRMHD simulations and high-resolution VLBI, these models provide insights into jet-medium interactions, particle acceleration, and high-energy neutrino production.

A structured-jet or “spine–layer” model describes relativistic jets whose flow and radiative properties are organized into at least two coaxial components: an inner fast “spine” and a surrounding slower “sheath” or “layer.” This geometric and dynamical stratification is supported by high-resolution observations of AGN (e.g. M87, 3C 273, NGC 1275, NGC 315), radio galaxies, BL Lacs, and gamma-ray bursts (GRBs), as well as by theoretical models and global relativistic magnetohydrodynamic (GRMHD) simulations. The structured-jet paradigm explains a wide array of multiwavelength emission signatures, variability, cosmic-ray and neutrino production, and jet–medium interactions that cannot be reconciled with single-zone (homogeneous) flows. The model is relevant for interpreting parsec-scale limb brightening, broad-band spectral energy distributions (SEDs), afterglow morphology in GRB events, polarization patterns, and the angular energy distribution of transients.

1. Geometric and Physical Model: Transverse Structure and Kinematics

A canonical spine–layer jet consists of two radially concentric regions:

• Spine: Fast, typically Poynting-flux dominated ($\sigma \gtrsim 1$), low-density electron-positron plasma. Bulk Lorentz factor $\Gamma_{\rm spine}$ can reach $10$–$40$ (blazars/GRBs), or $3$–$15$ (AGN, FR I/II), with narrow opening angle ($\theta_{\rm spine}\sim 1/\Gamma_{\rm spine}$). The spine is energized close to the central black hole, often via the Blandford–Znajek process (Xie et al., 2012).
• Sheath (Layer): Slower, denser, often matter-dominated or mildly magnetized, with Lorentz factor $\Gamma_{\rm sheath}\sim 2$–$5$ (AGN) or up to $10$ (GRB power-law wings), and wider opening angle ($\Gamma_{\rm spine}$0). The sheath is associated with disk winds (Blandford–Payne mechanism), interaction layers, or boundary mass loading (Xie et al., 2012, Walg et al., 2013).

Transverse velocity, density, magnetic field, and emissivity profiles can be step-like (piecewise isochoric), isothermal, power-law, broken power-law, or Gaussian (see analytic and simulation-based prescriptions in (Walg et al., 2013, Nokhrina et al., 6 Oct 2025, Coughlin et al., 2020, Gottlieb et al., 2020)). The velocity shear at the spine–sheath interface is central for particle acceleration and turbulence (Wang et al., 2022).

2. Formation Processes and Theoretical Basis

The stratified structure can be generated by several mechanisms:

• Central Engine Extraction: Blandford–Znajek extraction forms the spine, Blandford–Payne disk winds or outer-lobe injection forms the sheath (Xie et al., 2012).
• Jet–Medium Interaction: As relativistic jets propagate through stellar envelopes or ambient media, Kelvin–Helmholtz (KH), Rayleigh–Taylor, and Richtmyer–Meshkov instabilities naturally generate a mixing boundary layer (the “jet-cocoon interface”, JCI), producing universal power-law angular energy profiles (Gottlieb et al., 2020, Walg et al., 2013).
• Radiation Hydrodynamics: In optically thick, luminous outflows (e.g. TDEs, GRBs), the radiation diffusion across a thin angular boundary layer (width $\Gamma_{\rm spine}$1) creates a fast, low-mass core and a slow, massive sheath, with the Lorentz factor profile typically Gaussian or steeper (Coughlin et al., 2020).
• Relativistic MHD Structure: Analytical and numerical MHD models reproduce transverse fields, velocity, and pressure distributions and predict transitions from parabolic to conical collimation (Nokhrina et al., 6 Oct 2025, Chantry et al., 2017).

These processes together account for the persistent emergence of limb-brightened, multicomponent jets over large dynamic ranges in astrophysical phenomena.

3. Radiative and Dynamical Consequences

The velocity and emissivity stratification of spine–layer jets has multiple observable implications:

• Spectral Energy Distributions (SEDs): The spine dominates high-energy gamma-ray emission via inverse-Compton upscattering of external or sheath photons, while the sheath contributes to low-energy (radio/optical) synchrotron (Tavecchio et al., 2014, Migliori et al., 2011, Sikora et al., 2015, Janiak et al., 2015). In the external Compton scenario, each region's radiation field appears Doppler-boosted in the frame of the other, enhancing high-energy output and modifying SSC/EC ratios (Tavecchio et al., 2014, Sikora et al., 2015, Boughelilba et al., 2023).
• Limb Brightening and VLBI Signatures: Global VLBI and space-VLBI imaging resolve limb-brightened jets (e.g. NGC 315, 3C 273, M87), where the transverse brightness peaks at the edges rather than the axis (Park et al., 2024, Bruni et al., 2021, Punsly, 2019). This limb enhancement can arise from modest Lorentz factor gradients if the boundary emissivity is radially peaked (due to shear/turbulent acceleration or pair loading) rather than from an ultra-fast on-axis spine (Park et al., 2024, Nokhrina et al., 6 Oct 2025).
• Polarization: Toroidal and helical magnetic-field structures, combined with Doppler beaming, yield distinct edge and spine polarization patterns. Limb-brightened jets often show high edge polarization with EVPAs consistent with toroidal field dominance, transitioning to poloidal-dominated, spine-brightened regions downstream (Nokhrina et al., 6 Oct 2025).
• Neutrino and Cosmic-Ray Production: The sheath provides an amplified target photon field for $\Gamma_{\rm spine}$2 interactions of cosmic rays accelerated in the spine, boosting the efficiency of high-energy neutrino production in BL Lacs and radio galaxies (Tavecchio et al., 2014, Boughelilba et al., 2023).
• Variability and Correlations: The strong beaming dependence of spine emission leads to much larger gamma-ray flux variability (for small Doppler factor changes) than in the sheath's optical synchrotron, explaining “orphan” gamma-ray flares without X-ray/optical counterparts (Janiak et al., 2015, Sikora et al., 2015).

4. Particle Acceleration and Turbulence at the Interface

Relativistic shear acceleration at the spine–sheath boundary is a robust mechanism for distributed high-energy particle acceleration:

• Shear Acceleration Theory: The Fokker–Planck equation describes momentum diffusion with diffusion coefficient $\Gamma_{\rm spine}$3, setting both the acceleration timescale and the steady-state particle index $\Gamma_{\rm spine}$4 (Wang et al., 2022, Boughelilba et al., 2023).
• Turbulence: 3D RMHD simulations show that Kelvin–Helmholtz instabilities generate Kolmogorov-like turbulence in both spine and sheath. The sheath thickness and the turbulence correlation length control acceleration efficiency and the resultant spectral index ($\Gamma_{\rm spine}$5–$\Gamma_{\rm spine}$6 as required by X-ray emitting electrons) (Wang et al., 2022).
• Numerical Treatments: Advanced kinetic codes (e.g., CR-ENTREES) implement transition-matrix approaches that include complete radiative, hadronic, and shear-acceleration feedback between zones (Boughelilba et al., 2023).

This interface acceleration is essential for explaining the persistent nonthermal high-energy emission and observed limb-brightening over multi-kiloparsec jet scales.

5. Implications for Afterglow and Transient Emission

Structured-jet models are essential for interpreting the observed diversity of afterglow light curves and viewing-angle dependent signatures:

• Angular Energy and $\Gamma_{\rm spine}$7 Profiles: Hydrodynamic simulations show that both long and short GRB jets (and kilonova afterglows) are best described by a core (spine) plus a wide power-law wing (layer) in isotropic-equivalent energy $\Gamma_{\rm spine}$8 and Lorentz factor $\Gamma_{\rm spine}$9, shaped by the degree of mixing at the jet–cocoon interface (Gottlieb et al., 2020, Wang et al., 21 Jul 2025).
• Afterglow Light Curves: Jet structure controls the afterglow’s temporal evolution, including the timing and sharpness of jet breaks, plateau phases, and the presence of multi-bump features for various viewing angles. Off-axis observers in the layer see low-peak, slow-evolving transients, while on-axis spine emission produces brighter, sharper features (Wang et al., 21 Jul 2025).
• Reverse Shocks: The radial stratification (energy injection history) and the spine–layer angular structure affect the formation and detectability of reverse shock emission, with analytic and numerical time-dependent predictions now routinely including both dimensions (Wang et al., 21 Jul 2025).

These components are critical for accurate parameter inference in joint multi-frequency afterglow fitting and population studies of transients.

6. Controversies, Model Limitations, and Observational Challenges

Despite the empirical and theoretical success of spine–layer models, significant tensions and uncertainties remain:

• Conflict with GRMHD Simulations: High-frequency VLBI imaging of sources like NGC 315 and M87 reveals limb-brightened profiles at large viewing angles ($10$0) that require either an implausibly fast central spine ($10$1–$10$2 on parsec scales) or a strongly peaked edge emissivity, in tension with GRMHD models which predict a slow on-axis spine at these radii (Park et al., 2024, Punsly, 2019).
• Width and Profile Discrepancies: Many GRMHD+ray-tracing models tuned to EHT data predict jets that are too narrow and spine-brightened at $10$3–$10$4 lt-yr from the core, whereas observed images require a much broader, edge-brightened (sheath-dominated) structure (Punsly, 2019).
• Emission Physics Across the Jet: The mechanism(s) producing radially peaked emissivity—shear/turbulent acceleration, boundary mass loading, or pair creation—remain under-constrained and may differ between sources (Park et al., 2024, Nokhrina et al., 6 Oct 2025).
• Discrimination Among Models: Only multi-frequency, polarization-resolved super-resolution VLBI, Faraday tomography, and proper-motion mapping at sub-parsec scales can definitively distinguish between strong Doppler-boosting in a fast spine and enhanced edge-intrinsic emissivity (Park et al., 2024).

7. Extensions and Unified Context

The structured-jet framework accommodates and unifies a wide range of astrophysical phenomena:

• GRB, TDE, and AGN Jets: Similar spine–layer structures arise across scale from hyperaccreting stellar-mass engines to SMBHs, as predicted by radiation hydrodynamics, GRMHD, and jet–environment interaction models (Coughlin et al., 2020, Xie et al., 2012).
• Magnetization and Dissipation: Models incorporating magnetic reconnection, equipartition, and moderate $10$5 can reproduce both spectral and variability properties of strong-line blazars and their modest energetic requirements (Sikora et al., 2015).
• Cosmic Ray and Multi-Messenger Connections: The sheath's photon field is crucial for $10$6 interactions and neutrino production in BL Lacs, connecting structured jets to IceCube observations (Tavecchio et al., 2014, Boughelilba et al., 2023).

As a physically grounded, predictive framework, the spine–layer paradigm stands central in decoding the internal dynamics, emission, and evolution of astrophysical jets, linking their observed diversity to robust plasma and MHD processes.