The Voyager 2 crossing of the termination shock indicated that most of the upstream energy from the thermal solar wind ions was transferred to pickup ions (PUIs) and other energetic particles downstream of the shock. We use hybrid simulations at the termination shock for the Voyager 2, flank, and tail directions to evaluate the distributions of different ion species downstream of the shock over the energy range of 0.52–55 keV. Here, we extend the work of Gkioulidou et al., which showed an energy-dependent discrepancy between modeled and energetic neutral atom (ENA) observations, and fit distributions to a hybrid model to show that a population of PUIs accelerated via diffusive shock acceleration (DSA) to become low-energy anomalous cosmic rays (ACRs) can bridge the gap between modeled and observed ENA fluxes. Our results with the inclusion of DSA via hybrid fitting give entirely new and novel evidence that DSA at the termination shock is likely to be an important physical process. These ACRs carry a significant fraction of the energy density at the termination shock (22%, 13%, and 19% in the Voyager 2, flank, and tail directions, respectively). Using these ACRs in global ENA modeling of the heliosphere from 0.52 to 55 keV, we find that scaling factors as large as 1.8–2.5 are no longer required to match ENA observations at energies of ∼1–4 keV. Large discrepancies between modeled and observed ENAs only remain over energies of 4–20 keV, indicating that there may be a further acceleration mechanism in the heliosheath at these energies.

The shape of the heliosphere is currently under active debate. Energetic neutral atoms (ENAs) offer the best method for investigating the global structure of the heliosphere. To date, the Interstellar Boundary Explorer (IBEX) and the Ion and Neutral Camera (INCA) that was on board Cassini provide the only global ENA observations of the heliosphere. While extensive modeling has been done at IBEX-Hi energies (0.52–6 keV), no global ENA modeling has been conducted for INCA energies (5.2–55 keV). Here, we use an ENA model of the heliosphere based on hybrid results that capture the heating and acceleration of pickup ions (PUIs) at the termination shock to compare modeled global ENA results with IBEX-Hi and INCA observations using both a long- and short-tail model of the heliosphere. We find that the modeled ENA results for the two heliotail configurations produce similar results from the IBEX-Hi through the INCA energies. We conclude from our modeled ENAs, which only include PUI acceleration at the termination shock, that ENA observations in currently available energy ranges are insufficient for probing the shape and length of the heliotail. However, as a prediction for the future IMAP-Ultra mission (3–300 keV) we present modeled ENA maps at 80 keV, where the cooling length (∼600 au) is greater than the distance where the long- and short-heliotail models differ (∼400 au), and find that IMAP-Ultra should be able to identify the shape of the heliotail, predicting differences in the north lobe to downwind flux ratio between the models at 48%.

Heliospheric energetic neutral atoms (ENAs) originate from energetic ions that are neutralized by charge exchange with neutral atoms in the heliosheath and very local interstellar medium (VLISM). Since neutral atoms are unaffected by electromagnetic fields, they propagate ballistically with the same speeds as parent particles. Consequently, measurements of ENA distributions allow one to remotely image the energetic ion distributions in the heliosheath and VLISM. The origin of the energetic ions that spawn ENAs is still debated, particularly at energies higher than ∼keV. In this work, we summarize five possible sources of energetic ions in the heliosheath that cover the ENA energy from a few keV to hundreds of keV. Three sources of the energetic ions are related to pickup ions (PUIs): those PUIs transmitted across the heliospheric termination shock (HTS), those reflected once or multiple times at the HTS, i.e., reflected PUIs, and those PUIs multiply reflected and further accelerated by the HTS. Two other kinds of ions that can be considered are ions transmitted from the suprathermal tail of the PUI distribution and other particles accelerated at the HTS. By way of illustration, we use these energetic particle distributions, taking account of their evolution in the heliosheath, to calculate the ENA intensities and to analyze the characteristics of ENA spectra observed at 1 au.

The structure of shocks and turbulence are strongly modified during the acceleration of cosmic rays (CRs) at a shock wave. The pressure and the collisionless viscous stress decelerate the incoming thermal gas and thus modify the shock structure. A CR streaming instability ahead of the shock generates the turbulence on which CRs scatter. The turbulent magnetic field in turn determines the CR diffusion coefficient and further affects the CR energy spectrum and pressure distribution. The dissipation of turbulence contributes to heating the thermal gas. Within a multicomponent fluid framework, CRs and thermal gas are treated as fluids and are closely coupled to the turbulence. The system equations comprise the gas dynamic equations, the CR pressure evolution equation, and the turbulence transport equations, and we adopt typical parameters for the hot ionized interstellar medium. It is shown that the shock has no discontinuity but possesses a narrow but smooth transition. The self-generated turbulent magnetic field is much stronger than both the large-scale magnetic field and the preexisting turbulent magnetic field. The resulting CR diffusion coefficient is substantially suppressed and is more than three orders smaller near the shock than it is far upstream. The results are qualitatively consistent with certain observations.

No abstract