Bell Instability and Cosmic-Ray Acceleration in AGN Ultrafast Outflow Shocks

TL;DR

This study examines magnetic-field amplification by the nonresonant hybrid (Bell) instability at reverse shocks of AGN ultrafast outflows and its consequences for cosmic-ray acceleration. Using a self-consistent 1D MHD–CR framework that evolves CR diffusion–convection alongside NRH-driven magnetic growth, the authors map how $E_{\max}$ depends on the background field $B_0$, injection efficiency $\eta$, and initial turbulence $\xi_{B,\mathrm{ini}}$, including $\,p\gamma$ cooling. A key result is a transition: for weak $B_0$ ($\lesssim 10^{-4}$ G) the NRH instability amplifies upstream turbulence and drives $E_{\max}$ to a self-regulated value largely independent of initial turbulence, while for stronger fields ($\gtrsim 10^{-3}$ G) the escaping CR current is too weak to sustain NRH, causing $E_{\max}$ to follow initial conditions and possibly fall short of the EeV regime. The work also shows that higher ISM densities can enhance NRH growth via larger shock velocities, and $\,p\gamma$ cooling can cap $E_{\max}$ at high $B_0$, implying that UFOs can reach PeV–EeV energies only under a narrow set of environmental and spectral conditions.

Abstract

We investigate magnetic-field amplification driven by the nonresonant hybrid (NRH or Bell) instability and its impact on cosmic-ray (CR) acceleration at reverse shocks of ultrafast outflows (UFOs) from active galactic nuclei (AGN). Previous kinetic studies by particle-in-cell simulations have demonstrated that when maximum CR energy is near the injection scale, NRH instability efficiently amplifies magnetic field up to the saturation level. However, the efficiency of NRH instability goes down as maximum energy increase since CR current is carried by escaping CRs near the maximum energy. We employ a one-dimensional MHD--CR framework solving telegraph-type diffusion--convection equations to trace the coupled evolution of CRs, magnetic fields, and shock dynamics under realistic parameters. We find a distinct transition with magnetic field strength: for weak background fields ($B_{0}\!\lesssim\!10^{-4}\,\mathrm{G}$), NRH instability efficiently amplifies upstream turbulence, driving a self-regulated state where $E_{\max}$ becomes independent of initial strength of magnetic turbulence. In contrast, for stronger background fields ($B_{0}\!\gtrsim\!10^{-3}\,\mathrm{G}$), the escaping CR current is too weak to drive NRH instability, and magnetic turbulence further decays through parametric instabilities, potentially reducing the acceleration efficiency. We give the physical interpretation for the transition and discuss conditions for PeV--EeV acceleration at UFO reverse shocks.

Figures (11)

• Figure 1: Comparison of $\mathrm{p}\gamma$ cooling and acceleration timescales as a function of proton energy. The blue solid line shows the $\mathrm{p}\gamma$ cooling timescale, the orange dashed line indicates the simulation runtime (the propagation timescale of the reverse shock across pc scales), and the red dot–dashed line represents the acceleration timescale described in Sec. \ref{['subsec:pgamma_cooling_UFO']}, calculated from Eq. \ref{['eq:acceleration_time_UFO']} with $\xi_{B}=0.1$, $B_0=10^{-2}~\mathrm{G}$, and $v_{\text{sh}}=5.0\times10^9~\mathrm{cm~s^{-1}}$ (see Tab. \ref{['tab:Model_Parameters_UFO']} for the fiducial model).
• Figure 2: Left: Density distribution of the global wind–ISM interaction reproduced by one-dimensional MHD simulations. Right: Extracted density profile around the reverse shock highlighted as grey region in Left panel, focusing on the local region used for particle acceleration analysis.
• Figure 3: Energy spectra of the isotropic component of CRs in the Fiducial model with $(B_0,~\eta,~\xi_{B,\mathrm{ini}})=(10^{-5}~\mathrm{G},~10^{-4},~0.1)$ listed in Tab. \ref{['tab:Model_Parameters_UFO']}. The left panel shows the case without NRH instability, while the right panel includes NRH instability. Blue, orange, green, red, and purple solid lines represent the spectra at 20, 40, 60, 80, and 100 yr after the start of the simulation, respectively. The gray dashed line denotes the power-law slope predicted by the test-particle DSA solution. Red dots indicate the maximum acceleration energy, defined as the momentum at which the spectral index of $f_0p^3$ falls below $-2.1$.
• Figure 4: Spatial profile of magnetic fluctuations $\delta B/B_0$ in the fiducial model with $(B_0,~\eta,~\xi_{B,\mathrm{ini}})=(10^{-5}~\mathrm{G},~10^{-4},~0.1)$ at $t=100~\mathrm{yr}$. Left: without NRH instability. Right: with NRH instability. The horizontal axis denotes the distance from the shock front, normalized so that the shock position is $x=0$ (in pc). The vertical axis shows $\delta B/B_0$. The orange dashed line marks the shock position, while the green dot-dashed line represents the mean amplitude of the initial fluctuations ($\delta B/B_0 \sim 0.3$).
• Figure 5: Upstream profiles for the fiducial model with $(B_0,~\eta,~\xi_{B,\mathrm{ini}})=(10^{-5}~\mathrm{G},~10^{-4},~0.1)$ at $t=100~\mathrm{yr}$. The horizontal axis indicates the distance from the shock position (in pc), with the orange dashed line marking the shock. Left: current density carried by the CR anisotropic component $j_{\text{CR}}$ (esu s$^{-1}$ cm$^{-2}$). Right: e-folding number of NRH instability $t_{\mathrm{adv}}/t_{\mathrm{NRH}}$ (dimensionless), evaluated using Eq. \ref{['eq:NRH_linear_timescale']} and Eq. \ref{['eq:definition_advection_time_Bell']}. See also Tab. \ref{['tab:NRH_timescale_compare']}.