The origin of the very-high-energy radiation along the jet of Centaurus A

Problem Statement
Centaurus A (Cen A), the nearest active galactic nucleus (AGN), exhibits very-high-energy (VHE) gamma-ray emission (> 240 GeV) detected by the HESS observatory, indicating the presence of ultra-relativistic electrons along its jet. However, the specific acceleration mechanism responsible for these particles remains uncertain. While various authors have proposed that jet substructures, known as "knots," act as efficient particle accelerators, a quantitative model linking the hydrodynamics of these knots to the observed VHE spectrum has been lacking. Specifically, it is unclear if the collision of the jet with stellar winds (SW) can generate the necessary shock conditions to accelerate electrons to the required energies and reproduce the multi-wavelength morphology.

Methodology
The authors investigate the hypothesis that stationary knots in Cen A's jet result from interactions between the relativistic jet and the stellar winds of powerful stars (specifically Wolf-Rayet stars). The study employs a multi-step modeling approach:

1. Relativistic Hydrodynamic (RHD) Simulations: Using the PLUTO code, the authors simulate the interaction between the jet (assumed velocity $v_j = 0.6c$, power $L_j \sim 2 \times 10^{43}$ erg s$^{-1}$) and a stellar wind (characterized by $V_\infty = 3000$ km s$^{-1}$ and $\dot{M} = 10^{-4} M_\odot$ yr$^{-1}$). The simulations are performed in an axisymmetric 2D cylindrical frame to resolve the double-shock structure (bow shock and wind shock) and the contact discontinuity.
2. Particle Acceleration and Cooling: Based on diffusive shock acceleration (Fermi-I) theory, the authors assume electrons are accelerated at the shock waves with a power-law spectrum and an exponential cutoff at a maximum energy ($E_{max}$). The maximum energy is determined by balancing the acceleration timescale against synchrotron and Inverse Compton (IC) cooling timescales. The cooling depends on local magnetic fields (derived from thermal pressure via an equipartition parameter $\beta_m$) and radiation fields (core SSC, starlight, dust lane, CMB, and WR star radiation).
3. Radiative Modeling: The non-thermal Spectral Energy Distribution (SED) from radio to gamma-rays is calculated using the Naima library. The model is constrained by radio (VLA) and X-ray (Chandra) flux measurements and spectral indices for specific knots (AX1A, AX1C, AX2, BX2).
4. Morphological Comparison: Synthetic X-ray images are generated by projecting the simulated emissivity and applying a Gaussian filter to match Chandra's resolution, allowing for a direct comparison with observed knot sizes and shapes.

Key Contributions

• First Combined Analysis: This work represents the first attempt to combine fluid simulations, shock acceleration theory, and multi-wavelength observational data (radio, X-ray, and VHE gamma-ray) to test the jet-stellar wind collision hypothesis for Cen A's knots.
• Hydrodynamic Structure: The simulations successfully reproduce the double-shock structure expected from jet-SW interactions, providing a physical basis for the stationary nature of the observed knots.
• Extended Emission Modeling: The authors move beyond modeling only the knot centroids to include the extended emission from the entire perturbed/post-shock region, which is crucial for explaining the integrated VHE flux.
• Parameter Constraints: The study constrains key physical parameters, including the acceleration efficiency ($\chi$) and the magnetic field strength ($\beta_m$), by fitting the model to radio and X-ray data.

Results

• Electron Energies: The simulations suggest that electrons can be accelerated to energies of $\sim 3-4$ PeV at the AX2 and BX2 knots, and up to $\sim 6$ TeV at the inner AX1A/AX1C knots (assuming lower acceleration efficiency). These energies are sufficient to produce the observed VHE gamma-rays via IC scattering.
• Spectral Energy Distribution: The model satisfactorily reproduces the radio and X-ray fluxes for knots AX2 and BX2. However, for AX1A and AX1C, the standard model overestimates X-ray emission; a significantly lower acceleration efficiency ($\chi \approx 2 \times 10^{-3}$) is required to match the data, suggesting these inner knots may have different physical conditions (e.g., more relativistic shocks or different magnetization).
• Gamma-Ray Origin: The emission from the knot centroids alone is insufficient to explain the HESS data. However, the extended emission from the entire shocked region (downstream of the bow shock) provides a satisfactory explanation for the VHE gamma-ray flux detected by HESS.
• Fermi-LAT Excess: The model suggests that the Fermi-LAT signal (GeV range) may be explained by electrons that have cooled as they propagate along the jet (up to $\sim 3.6$ kpc) or are trapped in the giant lobes, rather than being produced directly at the knot centroids.
• Morphology: The predicted X-ray knot sizes for AX1C, AX2, and BX2 agree with Chandra measurements within 5–25%. The AX1A knot is overestimated in size, possibly indicating a weaker stellar wind or alternative formation mechanisms.
• UHECR Potential: While the knots can accelerate nuclei to $\sim 0.2 Z$ EeV in an optimistic Bohm diffusion scenario, the authors find this likely insufficient to explain the highest-energy UHECRs ($\sim 10^{19}$ eV) directly. However, the knots could inject heavier nuclei (from WR stars) that undergo subsequent acceleration in the lobes.

Significance and Claims
The paper claims that particle acceleration at jet-stellar wind interaction sites provides a plausible explanation for the radio, X-ray, and VHE gamma-ray emission observed in Cen A. The authors emphasize that the extended emission from the shocked regions, rather than just the knot centroids, is the dominant contributor to the VHE flux.

The study highlights that:

1. The stationary knots are likely formed by collisions with Wolf-Rayet stars.
2. The VHE gamma-ray spectrum is consistent with electrons accelerated at these shocks and cooling via IC scattering on dust and starlight radiation fields.
3. Future observations with the Cherenkov Telescope Array (CTAO) could spatially resolve the contributions of different knot complexes (A, B, C, E vs. F, G) in the 10–100 TeV range, offering a direct test of the acceleration efficiency and the dominant knots in the VHE band.

The authors remain modest regarding Ultra-High Energy Cosmic Rays (UHECRs), noting that while knots may inject heavy nuclei, the current jet activity likely cannot accelerate them to the highest observed energies without further re-acceleration in the lobes. Similarly, they acknowledge that diffusion effects, not fully captured in their advection-dominated model, may be necessary to fully explain the photon counts and knot radii.