Fluctuations Driven by Multicomponent Pickup Ion Distributions in the Outer Heliosheath

Abstract The stability of a realistic multicomponent pickup ion (PUI) velocity distribution derived from a global model of neutral atoms in the heliosphere, which treats hydrogen and helium atoms self-consistently and includes equations for electrons and helium ions, is investigated using linear instability analysis and hybrid simulations. Linear instability analysis shows that the excited oblique mirror waves and the parallel/quasi-parallel Alfvén-cyclotron (AC) waves have lower growth rates than those obtained previously by A. Mousavi et al. for the PUI velocity distributions given by J. Heerikhuisen et al. The PUI scattering by each of the two modes alone is studied. In contrast to the previous investigations, our current simulations using the updated realistic distributions indicate that mirror waves alone do not effectively scatter PUIs in pitch angle. Instead, they primarily contribute to reducing the thermal spread anisotropy of the PUIs originating from the neutral solar wind. The unstable AC waves exhibit lower growth rates but higher saturation levels than the mirror waves. Two-dimensional (2D) simulation results show that when all unstable waves are present, the predominant contributor to the fluctuating magnetic field energy is the AC mode. The AC waves quickly scatter the PUIs with pitch angles away from 90 ∘ toward isotropy, while the PUIs near 90 ∘ pitch angle maintain a degree of anisotropy within our simulation timeframe. Moreover, several 1D and 2D hybrid simulations with different numbers of particles per cell are performed to examine the impact of numerical noise on PUI scattering. Finally, the implications of these results for the Interstellar Boundary Explorer energetic neutral atom ribbon are discussed.