An official website of the United States government
Abstract 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.
more »
« less
This content will become publicly available on February 3, 2026
Hybrid Simulations of Interstellar Pickup Ions at the Solar Wind Termination Shock Revisited
Abstract We revisit previous hybrid simulations of the heating and acceleration of interstellar pickup ions (PUIs) at the solar wind termination shock. In previous simulations, a relatively cold initial distribution of PUIs was assumed; and while the resulting shock-heated distribution was consistent with Voyager 2 LECP measurements at about 30 keV, the intensity of the distribution downstream of the shock in the ~1–10 keV energy range was lower than predictions based on analysis of energetic neutral atoms (ENAs) from the Interstellar Boundary Explorer-Hi and Cassini's Ion and Neutral Camera. Here we perform new simulations with more realistic initial PUI distributions. We assume the distribution is a partially filled spherical shell in velocity space with a radius that varies from 320 to 640 km s−1. We then use the distributions downstream of the shock from these new simulations to estimate the ENA flux spectrum and compare with observations. We find that the predicted ENA spectrum from the new simulations much better matches the observations over a broad range of energies. We conclude that the hybrid simulations provide reasonable predictions for the distribution of charged particles in the energy range from ~0.5 to 50 keV.
more »
« less
- Award ID(s):
- 2148653
- PAR ID:
- 10573717
- Publisher / Repository:
- The Astrophysical Journal
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 980
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 29
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract 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.more » « less
-
Abstract 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%.more » « less
-
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.more » « less
-
Abstract Interstellar pickup ions are an ubiquitous and thermodynamically important component of the solar wind plasma in the heliosphere. These PUIs are born from the ionization of the interstellar neutral gas, consisting of hydrogen, helium, and trace amounts of heavier elements, in the solar wind as the heliosphere moves through the local interstellar medium. As cold interstellar neutral atoms become ionized, they form an energetic ring beam distribution comoving with the solar wind. Subsequent scattering in pitch angle by intrinsic and self-generated turbulence and their advection with the radially expanding solar wind leads to the formation of a filled-shell PUI distribution, whose density and pressure relative to the thermal solar wind ions grows with distance from the Sun. This paper reviews the history of in situ measurements of interstellar PUIs in the heliosphere. Starting with the first detection in the 1980s, interstellar PUIs were identified by their highly nonthermal distribution with a cutoff at twice the solar wind speed. Measurements of the PUI distribution shell cutoff and the He focusing cone, a downwind region of increased density formed by the solar gravity, have helped characterize the properties of the interstellar gas from near-Earth vantage points. The preferential heating of interstellar PUIs compared to the core solar wind has become evident in the existence of suprathermal PUI tails, the nonadiabatic cooling index of the PUI distribution, and PUIs’ mediation of interplanetary shocks. Unlike the Voyager and Pioneer spacecraft, New Horizon’s Solar Wind Around Pluto (SWAP) instrument is taking the only direct measurements of interstellar PUIs in the outer heliosphere, currently out to $$\sim47~\text{au}$$ ∼ 47 au from the Sun or halfway to the heliospheric termination shock.more » « less
