Interstellar neutral atoms provide a remote diagnostic of plasma in the outer heliosheath and the very local interstellar medium via charge exchange collisions that convert ions into atoms and vice versa. So far, most studies of interstellar atoms assumed that daughter hydrogen atoms directly inherit the kinetic properties of parent protons. This assumption neglects angular scattering of the interacting particles. However, for low relative velocities, as expected for charge exchanges in the outer heliosheath, this scattering is significant. In this study, we present how the parameters of daughter populations depend on the relative velocity and temperatures of parent populations. For this purpose, we numerically compute collision terms with and without this scattering. We find that the secondary population of interstellar hydrogen atoms, for the parent populations with the relative bulk velocity of 20 km s$^{-1}$ and equal temperatures of 7500 K, has ~2 km s$^{-1}$ higher bulk velocity if the scattering is taken into account. Additionally, temperatures are higher by ~2400 K and ~1200 K in parallel and perpendicular direction to the relative motion of parent populations, respectively. Moreover, a significant departure of secondary atoms from the Maxwell-Boltzmann distribution is expected for high relative velocities of parent populations. This process affects the distribution and density of interstellar atoms in the heliosphere and production of pickup ions. Thus, we show that angular scattering in charge exchange collisions is important to include in analyses of interstellar neutral atoms and pickup ions observed at 1 au and in the outer heliosphere.

Energetic neutral atoms (ENA) are an important tool for investigating the structure of the heliosphere. Recently, it was observed that fluxes of ENAs (with energy $\le$ 55 keV) coming from the upwind and downwind regions of the heliosphere are similar in strength. This led the authors of these observations to hypothesize that the heliosphere is bubble-like rather than comet-like, meaning that it has no extended tail. We investigate the directional distribution of the ENA flux for a wide energy range (3--88 keV) including the observations from IBEX (Interstellar Boundary Explorer), INCA (Ion and Neutral Camera, on board Cassini), and HSTOF (High energy Suprathermal Time Of Flight sensor, on board SOHO, Solar and Heliospheric Observatory). An essential element is the model of pickup ion acceleration at the termination shock (TS) proposed by Zank. We use state of the art models of the global heliosphere, interstellar neutral gas density, and pickup ion distributions. The results, based on the "comet-like" model of the heliosphere, are close in flux magnitude to ENA observations by IBEX, HSTOF and partly by INCA (except for the 5.2-13.5 keV energy channel). We find that the ENA flux from the tail dominates at high energy (in agreement with HSTOF, but not INCA). At low energy, our comet-like model produces the similar strengths of the ENA fluxes from the upwind and downwind directions, which, therefore, removes this as a compelling argument for a bubble-like heliosphere.

Models play an important role in our understanding of the global structure of the solar wind and its interaction with the interstellar medium. A critical ingredient in many types of models are the charge-exchange collisions between ions and neutrals. Some ambiguity exists in the charge-exchange cross-section for protons and hydrogen atoms, depending on which experimental data is used. The differences are greatest at low energies, and for the plasma-neutral interaction in the outer heliosheath may exceed 50\%. In this paper we assess a number of existing data sets and formulae for proton-hydrogen charge-exchange. We use a global simulation of the heliosphere to quantify the differences between the currently favored cross-section, and a formulation we suggest that more closely matches the majority of available data. We find that in order to make the resulting two heliospheres the same size, the interstellar proton and hydrogen densities need to be adjusted by 10 to 15\%, which provides a way to link the uncertainty in the cross-section to the uncertainty in the parameters of the pristine interstellar plasma.

Interstellar neutral gas atoms penetrate the heliopause and reach 1~au, where they are detected by IBEX. The flow of neutral interstellar helium through the perturbed interstellar plasma in the outer heliosheath (OHS) results in creation of the secondary population of interstellar He atoms, the so-called Warm Breeze, due to charge exchange with perturbed ions. The secondary population brings the imprint of the OHS conditions to the IBEX-Lo instrument. Based on a global simulation of the heliosphere with measurement-based parameters and detailed kinetic simulation of the filtration of He in the OHS, we find the number density of interstellar He$^+$ population at $(8.98\pm 0.12)\times 10^{-3}$~cm$^{-3}$. With this, we obtain the absolute density of interstellar H$^+$ $5.4\times 10^{-2}$~cm$^{-3}$ and electrons $6.3\times 10^{-2}$~cm$^{-3}$, and ionization degrees of H 0.26 and He 0.37. The results agree with estimates of the Very Local Interstellar Matter parameters obtained from fitting the observed spectra of diffuse interstellar EUV and soft X-Ray background.

The effects of turbulence in the very local interstellar medium (VLISM) have been proposed by Giacalone & Jokipii (2015) to be important in determining the structure of the Interstellar Boundary Explorer (IBEX) ribbon via particle trapping by magnetic mirroring. We further explore this effect by simulating the motion of charged particles in a turbulent magnetic field superposed on a large-scale mean field, which we have considered to be either spatially-uniform or a mean field derived from a 3D MHD simulation. We find that the ribbon is not double-peaked, in contrast to Giacalone & Jokipii (2015). However, the magnetic mirror force still plays an important role in trapping particles. Furthermore, the ribbon$'$s thickness is considerably larger if the large-scale mean field is draped around the heliosphere. Voyager 1 observations in the VLISM show a turbulent field component that is stronger than previously thought, which we test in our simulation. We find that the inclusion of turbulent fluctuations at scales ${\gtrsim}$100 au and power consistent with Voyager 1 observations produces a ribbon whose large-scale structure is inconsistent with IBEX observations. However, restricting the fluctuations to ${\sim}$10 au or smaller produces a smoother ribbon structure similar to IBEX observations. Different turbulence realizations produce different small-scale features ${\lesssim}10{\deg}$ in the ribbon, but its large-scale structure is robust if the maximum fluctuation size is ${\sim}$10 au. This suggests that the magnetic field structure at scales ${\lesssim}$10 au is determined by the heliosphere$'$s interaction with the VLISM and cannot entirely be represented by homogeneous interstellar turbulence.

The heliosphere is formed due to interaction between the solar wind (SW) and local interstellar medium (LISM). The shape and position of the heliospheric boundary, the heliopause, in space depend on the parameters of interacting plasma flows. The interplay between the asymmetrizing effect of the interstellar magnetic field and charge exchange between ions and neutral atoms plays an important role in the SW-LISM interaction. By performing three-dimensional, MHD plasma / kinetic neutral atom simulations, we determine the width of the outer heliosheath - the LISM plasma region affected by the presence of the heliosphere - and analyze quantitatively the distributions in front of the heliopause. It is shown that charge exchange modifies the LISM plasma to such extent that the contribution of a shock transition to the total variation of plasma parameters becomes small even if the LISM velocity exceeds the fast magnetosonic speed in the unperturbed medium. By performing adaptive mesh refinement simulations, we show that a distinct boundary layer of decreased plasma density and enhanced magnetic field should be observed on the interstellar side of the heliopause. We show that this behavior is in agreement with the plasma oscillations of increasing frequency observed by the plasma wave instrument onboard Voyager 1. We also demonstrate that Voyager observations in the inner heliosheath between the heliospheric termination shock and the heliopause are consistent with dissipation of the heliospheric magnetic field. The choice of LISM parameters in this analysis is based on the simulations that fit observations of energetic neutral atoms performed by IBEX.

Both physical arguments and simulations of the global heliosphere indicate that the tailward heliopause is flattened considerably in the direction perpendicular to both the incoming flow and the large-scale interstellar magnetic field. Despite this fact, all of the existing global analytical models of the outer heliosheath's magnetic field assume a circular cross section of the heliotail. To eliminate this inconsistency, we introduce a mathematical procedure by which any analytically or numerically given magnetic field can be deformed in such a way that the cross sections along the heliotail axis attain freely prescribed, spatially dependent values for their total area and aspect ratio. The distorting transformation of this method honors both the solenoidality condition and the stationary induction equation with respect to an accompanying flow field, provided that both constraints were already satisfied for the original magnetic and flow fields prior to the transformation. In order to obtain realistic values for the above parameters, we present the first quantitative analysis of the heliotail's overall distortion as seen in state-of-the-art three-dimensional hybrid MHD-kinetic simulations.

The solar wind carves a bubble in the surrounding interstellar medium (ISM), known as the heliosphere. Charged interstellar dust grains (ISDG) encountering the heliosphere may be diverted around the heliopause or penetrate it depending on their charge-to-mass ratio. We present new calculations of trajectories of ISDG in the heliosphere, and the dust density distributions that result. We include up-to-date grain charging calculations using a realistic UV radiation field and full 3-D magnetohydrodynamic fluid + kinetic models for the heliosphere. Models with two different (constant) polarities for the solar wind magnetic field (SWMF) are used, with the grain trajectory calculations done separately for each polarity. Small grains a_gr ~ 0.01 micron are completely excluded from the inner heliosphere. Large grains, a_gr ~ 1.0 micron pass into the inner solar system and are concentrated near the Sun by its gravity. Trajectories of intermediate size grains depend strongly on the SWMF polarity. When the field has magnetic north pointing to ecliptic north, the field de-focuses the grains resulting in low densities in the inner heliosphere, while for the opposite polarity the dust is focused near the Sun. The ISDG density outside the heliosphere inferred from applying the model results to in situ dust measurements is inconsistent with local ISM depletion data for both SWMF polarities, but is bracketed by them. This result points to the need to include the time variation in the SWMF polarity during grain propagation. Our results provide valuable insights for interpretation of the in situ dust observations from Ulysses.

The Interstellar Boundary Explorer (IBEX) spacecraft is providing the first all-sky maps of the energetic neutral atoms (ENAs) produced by charge-exchange between interstellar neutral \HI\ atoms and heliospheric solar wind and pickup ions in the heliosphere boundary regions. The 'edge' of the interstellar cloud presently surrounding the heliosphere extends less than 0.1 pc in the upwind direction, terminating at an unknown distance, indicating that the outer boundary conditions of the heliosphere could change during the lifetime of the IBEX satellite. Using reasonable values for future outer heliosphere boundary conditions, ENA fluxes are predicted for one possible source of ENAs coming from outside of the heliopause. The ENA production simulations use three-dimensional MHD plasma models of the heliosphere that include a kinetic description of neutrals and a Lorentzian distribution for ions. Based on this ENA production model, it is then shown that the sensitivities of the IBEX 1.1 keV skymaps are sufficient to detect the variations in ENA fluxes that are expected to accompany the solar transition into the next upwind cloud. Approximately 20% of the IBEX 1.1 keV pixels appear capable of detecting the predicted model differences at the $ 3 \sigma$ level, with these pixels concentrated in the Ribbon region. Regardless of the detailed ENA production model, the success of the modeled \BdotR\ directions in reproducing the Ribbon locus, together with our results, indicate that the Ribbon phenomenon traces the variations in the heliosphere distortion caused by the relative pressures of the interstellar magnetic and gaseous components.

The distribution of interstellar dust grains (ISDG) observed in the Solar System depends on the nature of the interstellar medium-solar wind interaction. The charge of the grains couples them to the interstellar magnetic field (ISMF) resulting in some fraction of grains being excluded from the heliosphere while grains on the larger end of the size distribution, with gyroradii comparable to the size of the heliosphere, penetrate the termination shock. This results in a skewing the size distribution detected in the Solar System. We present new calculations of grain trajectories and the resultant grain density distribution for small ISDGs propagating through the heliosphere. We make use of detailed heliosphere model results, using three-dimensional (3-D) magnetohydrodynamic/kinetic models designed to match data on the shape of the termination shock and the relative deflection of interstellar neutral H and He flowing into the heliosphere. We find that the necessary inclination of the ISMF relative to the inflow direction results in an asymmetry in the distribution of the larger grains (0.1 micron) that penetrate the heliopause. Smaller grains (0.01 micron) are completely excluded from the Solar System at the heliopause.

Modeling of the global heliosphere seeks to investigate the interaction of the solar wind with the partially ionized local interstellar medium. Models that treat neutral hydrogen self-consistently and in great detail, together with the plasma, but that neglect magnetic fields, constitute a sub-category within global heliospheric models. There are several different modeling strategies used for this sub-category in the literature. Differences and commonalities in the modeling results from different strategies are pointed out. Plasma-only models and fully self-consistent models from four research groups, for which the neutral species is modeled with either one, three, or four fluids, or else kinetically, are run with the same boundary parameters and equations. They are compared to each other with respect to the locations of key heliospheric boundary locations and with respect to the neutral hydrogen content throughout the heliosphere. In many respects, the models' predictions are similar. In particular, the locations of the termination shock agree to within 7% in the nose direction and to within 14% in the downwind direction. The nose locations of the heliopause agree to within 5%. The filtration of neutral hydrogen from the interstellar medium into the inner heliosphere, however, is model dependent, as are other neutral results including the hydrogen wall. These differences are closely linked to the strength of the interstellar bow shock. The comparison also underlines that it is critical to include neutral hydrogen into global heliospheric models.

We investigate heliosheath energetic neutral atom (ENA) fluxes at keV energies, by assuming that the heliosheath proton distribution can be approximated by a kappa-distribution. The choice of the kappa parameter derives from observational data of the solar wind (SW). This has direct applications to the upcoming IBEX mission. We will look at all-sky ENA maps within the IBEX energy range (10 eV to 6 keV), as well as ENA energy spectra in several directions. We find that the use of kappa, as opposed to a Maxwellian, gives rise to greatly increased ENA fluxes above 1 keV, while medium energy fluxes are somewhat reduced. We show how IBEX data can be used to estimate the spectral slope in the heliosheath, and that the use of kappa reduces the differences between ENA maps at different energies. We also investigate the effect introducing a kappa-distribution has on the global interaction between the SW and the local interstellar medium (LISM), and find that there is generally an increase in energy transport from the heliosphere into the LISM, due to the modified profile of ENA's energies. This results in a termination shock that moves out by 4 AU, a heliopause that moves in by 9 AU and a bow shock 25 AU farther out, in the nose direction.

New solar wind data from the Voyager 1 and Voyager 2 spacecraft, together with the SOHO SWAN measurements of the direction that neutral hydrogen enters into the inner heliosheath and neutral helium measurements provided by multiple observations are expected to provide more reliable constraints on the ionization ratio of the local interstellar medium (LISM) and the direction and magnitude of the interstellar magnetic field (ISMF). In this paper we use currently the most sophisticated numerical model of the heliospheric interface, which is based on an MHD treatment of the ion flow and kinetic modeling of neutral particles, to analyze an ISMF-induced asymmetry of the heliosphere in the presence of the interplanetary magnetic field and neutral particles. It is shown that secondary hydrogen atoms modify the LISM properties leading to its shock-free deceleration at the heliopause. We determine the deflection of hydrogen atoms from their original trajectory in the unperturbed LISM and show that it occurs not only in the plane defined by the ISMF and LISM velocity vectors, but also, to a lesser extent, perpendicular to this plane. We also consider the possibility of using 2-3 kHz radio emission data to further constrain the ISMF direction.