Globally Distributed Energetic Neutral Atom Maps for the "Croissant" Heliosphere

A recent study by Opher et al. (2015) suggested the heliosphere has a "croissant" shape, where the heliosheath plasma is confined by the toroidal solar magnetic field. The "croissant" heliosphere is in contrast to the classically accepted view of a comet-like tail. We investigate the effect of the "croissant" heliosphere model on energetic neutral atom (ENA) maps. Regardless of the existence of a split tail, the confinement of the heliosheath plasma should appear in ENA maps. ENA maps from the Interstellar Boundary Explorer (IBEX) have shown two high latitude lobes with excess ENA flux at higher energies in the tail of the heliosphere. These lobes could be a signature of the confinement of the heliosheath plasma, while some have argued they are caused by the fast/slow solar wind profile. Here we present ENA maps of the "croissant" heliosphere, focusing on understanding the effect of the heliosheath plasma collimation by the solar magnetic field while using a uniform solar wind. We incorporate pick-up ions (PUIs) into our model based on Malama et al. (2006) and Zank et al. (2010). We use the neutral solution from our MHD model to determine the angular variation of the PUIs, and include the extinction of PUIs in the heliosheath. In the presence of a uniform solar wind, we find that the collimation in the "croissant" heliosphere does manifest itself into two high latitude lobes of increased ENA flux in the downwind direction.Comment: 14 pages, 1 table, 7 figures, Accepted for publication in Ap

Using hydrogen energetic neutral atoms to study the heliosphere

The interaction between the solar wind and the partially ionized gas of the local interstellar medium (ISM) creates a bubble known as the heliosphere. Classically, the shape of the heliosphere has been regarded as comet-like, with a long tail pointed in the direction opposite the Sun’s motion through the ISM. In this view, the solar magnetic field was assumed to have a negligible effect on the global structure of the heliosphere. Recent advances in numerical modeling have revealed the importance of the solar magnetic field in its ability to confine and collimate the solar wind plasma, and the shape of the heliosphere has been called into question. Energetic neutral atoms (ENAs) are created throughout the heliosphere via charge exchange. The separate contributions of the solar magnetic field topology and the solar wind structure to ENA observations is largely unexplored. The Interstellar Boundary Explorer (IBEX) has been providing a global perspective of the heliosphere through ENA maps with energies ranging from 0.2 to 6 keV. In this dissertation, three-dimensional magnetohydrodynamic simulations of the heliosphere are used as input to an ENA model designed to produce synthetic ENA maps. I compare modeled ENA maps with IBEX observations to investigate how different heliospheric conditions and properties affect ENAs created in the heliosphere, and therefore how ENA observations can be used to understand the heliosphere. First, I investigate the effect of the solar wind collimation by the solar magnetic field on ENA maps in the case of a solar wind without latitudinal variation. I find that even in the absence of variations of the solar wind, two lobes of strong ENA flux form at high latitudes, similar to what is observed by IBEX at high energies. Second, I test the effect of a latitudinally-varying solar wind on ENAs both with and without the inclusion of the solar magnetic field. I show that the latitudinal variations of the solar wind during solar minimum creates a structured ENA profile with latitude, corresponding to the profile observed at 1 AU, but that the solar magnetic field significantly enhances ENA flux in the region where the solar wind is confined. Lastly, I investigate the effect of the solar cycle on ENAs and how changing solar wind conditions (e.g. density, temperature, velocity) affect the heliosphere over time. I demonstrate that, given changes in the solar cycle, there is a significant evolution in the modeled ENA flux due to the changes in the solar wind profile and the solar magnetic field, which is also seen by ENA observations

The Structure of the Large-Scale Heliosphere as Seen by Current Models.

This review summarizes the current state of research aiming at a description of the global heliosphere using both analytical and numerical modeling efforts, particularly in view of the overall plasma/neutral flow and magnetic field structure, and its relation to energetic neutral atoms. Being part of a larger volume on current heliospheric research, it also lays out a number of key concepts and describes several classic, though still relevant early works on the topic. Regarding numerical simulations, emphasis is put on magnetohydrodynamic (MHD), multi-fluid, kinetic-MHD, and hybrid modeling frameworks. Finally, open issues relating to the physical relevance of so-called "croissant" models of the heliosphere, as well as the general (dis)agreement of model predictions with observations are highlighted and critically discussed

The Heliosphere and Local Interstellar Medium from Neutral Atom Observations at Energies Below 10 keV.

As the heliosphere moves through the surrounding interstellar medium, a fraction of the interstellar neutral helium, hydrogen, and heavier species crossing the heliopause make it to the inner heliosphere as neutral atoms with energies ranging from few eV to several hundred eV. In addition, energetic neutral hydrogen atoms originating from solar wind protons and from pick-up ions are created through charge-exchange with interstellar atoms. This review summarizes all observations of heliospheric energetic neutral atoms and interstellar neutrals at energies below 10 keV. Most of these data were acquired with the Interstellar Boundary Explorer launched in 2008. Among many other IBEX breakthroughs, it provided the first ever all-sky maps of energetic neutral atoms from the heliosphere and enabled the science community to measure in-situ interstellar neutral hydrogen, oxygen, and neon for the first time. These observations have revolutionized and keep challenging our understanding of the heliosphere shaped by the combined forces of the local interstellar flow, the local interstellar magnetic field, and the time-dependent solar wind

Relating Energetic Ion Spectra to Energetic Neutral Atoms

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

A kinetic-magnetohydrodynamic model with adaptive mesh refinement for modeling heliosphere neutral-plasma interaction

The charge exchange between the interstellar medium (ISM) and the solar wind plasma is crucial for determining the structures of the heliosphere. Since both the neutral-ion and neutral-neutral collision mean free paths are either comparable to or larger than the size of the heliosphere, the neutral phase space distribution can deviate far away from the Maxwellian distribution. A kinetic description for the neutrals is crucial for accurately modeling the heliosphere. It is computationally challenging to run three-dimensional (3D) time-dependent kinetic simulations due to the large number of macro-particles. In this paper, we present the new highly efficient SHIELD-2 model with a kinetic model of neutrals and a magnetohydrodynamic (MHD) model for the ions and electrons. To improve the simulation efficiency, we implement adaptive mesh refinement (AMR) and particle splitting and merging algorithms for the neutral particles to reduce the particle number that is required for an accurate simulation. We present several tests to verify and demonstrate the capabilities of the model

Measuring Neutral Hydrogen Properties in the Interplanetary Medium: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033 white paper e-id. 275

International audienceThis supports development and deployment-to-space of a high-resolution spectrograph to distinguish different populations of H atoms that directly interact at the interface of the heliosheath, the region where the solar wind is subsonic. Resulting science would support Voyagers, IBEX, IMAP, New Horizons, and Interstellar Probe mission observations, and improve characterization of the solar wind

The Discrepancy between Observed and Predicted Heliospheric Energetic Neutral Atoms below Solar Wind Energy

Measuring energetic neutral atoms (ENAs) allows for the remote observation of ion populations from the frontiers of our heliosphere. In this study, we compare the ENAs observed with the IBEX-Lo instrument onboard the Interstellar Boundary Explorer with ENA predictions from two heliosphere models. In contrast to previous studies, this paper presents model-data comparisons for the energy range 50 eV–2 keV over one full solar cycle not only in the upwind direction (Voyager 1 and Voyager 2 sky directions), but also for the north pole, south pole, port tail lobe, and downwind directions. The two heliosphere models produce the same basic result: there is a large gap (1 to 2 orders of magnitude in ENA intensity at 100 eV) between ENA data and model predictions between 100 and 500 eV for all sky directions. The reason for this gap is not understood yet. While some explanations are plausible and will be investigated in future studies, other explanations are excluded

Measuring Neutral Hydrogen Properties in the Interplanetary Medium: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033 white paper e-id. 275

International audienceThis supports development and deployment-to-space of a high-resolution spectrograph to distinguish different populations of H atoms that directly interact at the interface of the heliosheath, the region where the solar wind is subsonic. Resulting science would support Voyagers, IBEX, IMAP, New Horizons, and Interstellar Probe mission observations, and improve characterization of the solar wind

The Passage of the Solar System through the Edge of the Local Bubble

The Sun moves through the interstellar medium (ISM) at a velocity of similar to 19 pc Myr-1, making the conditions outside the solar system vary with time over millions of years. Today's solar system is protected from interstellar particles by the heliosphere, the bubble formed by the solar wind as the Sun moves through the ISM, which engulfs the planets. There is geological evidence from 60Fe that Earth was in direct contact with the ISM 2-3 and 5-7 million years ago (MYA). Recent work argues that the Sun encountered a massive cold cloud 2 MYA as part of the Local Ribbon of Cold Clouds that shrunk the heliosphere and exposed Earth to the ISM. Here, we consider the effects of the passage of the Sun through the edge of the Local Bubble occurring at 6.8-0.4+0.5 MYA assuming that the Sun encountered a cloud with a density of 900 cm-3. If we consider additional turbulent motion within the cloud due to shocks, the density encountered can be as low as 283 cm-3. Clouds of this density cover a small but nonzero (less than or similar to 4.6%) fraction of the surface of the Local Bubble, making an encounter plausible. Using a state-of-the art magnetohydrodynamic model, we show that the heliosphere shrank to a scale smaller than Earth's orbit, thereby exposing Earth to cold dense ISM, consistent with 60Fe evidence. The timing of the event matches perturbations observed in the paleoclimate record recovered from deep-sea sediments. The passage through the Local Bubble's surface and contraction of the heliosphere therefore may have impacted the climate and biosphere significantly, suggesting a new driver of major events in Earth's history