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.

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.

Observations of energetic neutral atoms (ENAs) allow for remote studies of the plasma condition in the heliosphere and the neighboring local interstellar medium. The first IBEX results revealed an arc-like enhancement of the ENA intensity in the sky, known as the ribbon. The ribbon was not expected from the heliospheric models prior to the IBEX launch. One of the proposed explanations of the ribbon is the mechanism of the secondary ENA emission. The ribbon reveals energy-dependent structure in the relative intensity along its circumference and in the position. Namely, the ribbon geometric center varies systematically by about 10$^\circ$ in the energy range 0.7-4.3 keV. Here, we show by analytic modeling that this effect is a consequence of the helio-latitudinal structure of the solar wind reflected in the secondary ENAs. Along with a recently measured distance to the ribbon source just beyond the heliopause, our findings support the connection of the ribbon with the local interstellar magnetic field by the mechanism of the secondary ENA emission. However, the magnitude of the center shift in the highest IBEX energy channel is much larger in the observations than expected from the modeling. This may be due to another, not currently recognized, process of ENA generation.

Direct-sampling observations of interstellar neutral (ISN) He by Interstellar Boundary Explorer (IBEX) provide valuable insight into the physical state of and processes operating in the interstellar medium ahead of the heliosphere. The ISN He atom signals are observed at the four lowest ESA steps of the IBEX-Lo sensor. The observed signal is a mixture of the primary and secondary components of ISN He and H. Previously, only data from one of the ESA steps have been used. Here, we extended the analysis to data collected in the three lowest ESA steps with the strongest ISN He signal, for the observation seasons 2009-2015. The instrument sensitivity is modeled as a linear function of the atom impact speed onto the sensor's conversion surface separately for each ESA step of the instrument. We found that the sensitivity increases from lower to higher ESA steps, but within each of the ESA steps it is a decreasing function of the atom impact speed. This result may be influenced by the hydrogen contribution, which was not included in the adopted model, but seems to exist in the signal. We conclude that the currently accepted temperature of ISN He and velocity of the Sun through the interstellar medium do not need a revision, and we sketch a plan of further data analysis aiming at investigating ISN H and a better understanding of the population of ISN He originating in the outer heliosheath.

Observations of energetic neutral atoms (ENAs) are a fruitful tool for remote diagnosis of the plasma in the heliosphere and its vicinity. So far, instruments detecting ENAs from the heliosphere were configured for observations of hydrogen atoms. Here, we estimate emissions of ENAs of the heavy chemical elements helium, oxygen, nitrogen, and neon. A large portion of the heliospheric ENAs is created in the inner heliosheath from neutralized interstellar pick-up ions (PUIs). We modeled this process and calculated full-sky intensities of ENAs for energies 0.2-130 keV/nuc. We found that the largest fluxes among considered species are expected for helium, smaller for oxygen and nitrogen, and smallest for neon. The obtained intensities are 50-10^6 times smaller than the hydrogen ENA intensities observed by IBEX. The detection of heavy ENAs will be possible if a future ENA detector is equipped with the capability to measure the masses of observed atoms. Because of different reaction cross-sections among the different species, observations of heavy ENAs can allow for a better understanding of global structure of the heliosphere as well as the transport and energization of PUIs in the heliosphere.

Observations of energetic neutral atoms (ENAs) allow for remote sensing of plasma properties in distant regions of the heliosphere. So far, most of the observations have concerned only hydrogen atoms. In this paper, we present perspectives for observations of helium energetic neutral atoms (He ENAs). We calculated the expected intensities of He ENAs created by the neutralization of helium ions in the inner heliosheath and through the secondary ENA mechanism in the outer heliosheath. We found that the dominant source region for He ENAs is the inner heliosheath. The obtained magnitudes of intensity spectra suggest that He ENAs can be observed with future ENA detectors, as those planned on Interstellar Mapping and Acceleration Probe. Observing He ENAs is most likely for energies from a few to a few tens of keV/nuc. Estimates of the expected count rates show that the ratio of helium to hydrogen atoms registered in the detectors can be as low as 1:10^4. Consequently, the detectors need to be equipped with an appropriate mass spectrometer capability, allowing for recognition of chemical elements. Due to the long mean free paths of helium ions in the inner heliosheath, He ENAs are produced also in the distant heliospheric tail. This implies that observations of He ENAs can resolve its structure, which seems challenging from observations of hydrogen ENAs since energetic protons are neutralized before they progress deeper in the heliospheric tail.

Maps of Energetic Neutral Atom (ENA) fluxes obtained from Interstellar Boundary Explorer (IBEX) observations revealed a bright structure extending over the sky, subsequently dubbed the IBEX ribbon. The ribbon had not been expected from the existing models and theories prior to IBEX, and a number of mechanisms have since been proposed to explain the observations. In these mechanisms, the observed ENAs emerge from source plasmas located at different distances from the Sun. Since each part of the sky is observed by IBEX twice during the year from opposite sides of the Sun, the apparent position of the ribbon as observed in the sky is shifted due to parallax. To determine the ribbon parallax, we found the precise location of the maximum signal of the ribbon observed in each orbital arc. The obtained apparent positions were subsequently corrected for the Compton-Getting effect, gravitational deflection, and radiation pressure. Finally, we selected a part of the ribbon where its position is similar between the IBEX energy passbands. We compared the apparent positions obtained from the viewing locations on the opposite sides of the Sun, and found that they are shifted by a parallax angle of $0.41^\circ\pm0.15^\circ$, which corresponds to a distance of $140^{+84}_{-38}$ AU. This finding supports models of the ribbon with the source located just outside the heliopause.

With the velocity vector and temperature of the pristine interstellar neutral (ISN) He recently obtained with high precision from a coordinated analysis summarized by McComas et al.2015b, we analyzed the IBEX observations of neutral He left out from this analysis. These observations were collected during the ISN observation seasons 2010---2014 and cover the region in the Earth's orbit where the Warm Breeze persists. We used the same simulation model and a very similar parameter fitting method to that used for the analysis of ISN He. We approximated the parent population of the Warm Breeze in front of the heliosphere with a homogeneous Maxwell-Boltzmann distribution function and found a temperature of $\sim 9\,500$ K, an inflow speed of 11.3 km s$^{-1}$, and an inflow longitude and latitude in the J2000 ecliptic coordinates $251.6^\circ$, $12.0^\circ$. The abundance of the Warm Breeze relative to the interstellar neutral He is 5.7\% and the Mach number is 1.97. The newly found inflow direction of the Warm Breeze, the inflow directions of ISN H and ISN He, and the direction to the center of IBEX Ribbon are almost perfectly co-planar, and this plane coincides within relatively narrow statistical uncertainties with the plane fitted only to the inflow directions of ISN He, ISN H, and the Warm Breeze. This co-planarity lends support to the hypothesis that the Warm Breeze is the secondary population of ISN He and that the center of the Ribbon coincides with the direction of the local interstellar magnetic field. The common plane for the direction of inflow of ISN gas, ISN H, the Warm Breeze, and the local interstellar magnetic field %includes the Sun and is given by the normal direction: ecliptic longitude $349.7^\circ \pm 0.6^\circ$ and latitude $35.7^\circ \pm 0.6^\circ$ in the J2000 coordinates, with the correlation coefficient of 0.85.

The modeling of the heliosphere requires continuous three-dimensional solar wind data. The in-situ out-of-ecliptic measurements are very rare, so that other methods of solar wind detection are needed. We use the remote-sensing data of the solar wind speed from observations of interplanetary scintillation (IPS) to reconstruct spatial and temporal structures of the solar wind proton speed from 1985 to 2013. We developed a method of filling the data gaps in the IPS observations to obtain continuous and homogeneous solar wind speed records. We also present a method to retrieve the solar wind density from the solar wind speed, utilizing the invariance of the solar wind dynamic pressure and energy flux with latitude. To construct the synoptic maps of the solar wind speed we use the decomposition into spherical harmonics of each of the Carrington rotation map. To fill the gaps in time we apply the singular spectrum analysis to the time series of the coefficients of spherical harmonics. We obtained helio-latitudinal profiles of the solar wind proton speed and density over almost three recent solar cycles. The accuracy in the reconstruction is, due to computational limitations, about 20%. The proposed methods allow us to improve the spatial and temporal resolution of the model of the solar wind parameters presented in our previous paper (Sok\'o{\l} et al. 2013) and give a better insight into the time variations of the solar wind structure. Additionally, the solar wind density is reconstructed more accurately and it fits better to the in-situ measurements from Ulysses.

This paper is one of three companion papers presenting the results of our in-depth analysis of the interstellar neutral helium (ISN He) observations carried out using the IBEX-Lo during the first six Interstellar Boundary Explorer (IBEX) observation seasons. We derive corrections for losses due to the limited throughput of the interface buffer and determine the IBEX spin-axis pointing. We develop an uncertainty system for the data, taking into account the resulting correlations between the data points. This system includes uncertainties due to Poisson statistics, background, spin-axis determination, systematic deviation of the boresight from the prescribed position, correction for the interface buffer losses, and the expected Warm Breeze (WB) signal. Subsequently, we analyze the data from 2009 to examine the role of various components of the uncertainty system. We show that the ISN He flow parameters are in good agreement with the values obtained by the original analysis. We identify the WB as the principal contributor to the global $\chi^2$ values in previous analyses. Other uncertainties have a much milder role and their contributions are comparable to each other. The application of this uncertainty system reduced the minimum $\chi^2$ value 4-fold. The obtained $\chi^2$ value, still exceeding the expected value, suggests that either the uncertainty system may still be incomplete or the adopted physical model lacks a potentially important element, which is likely an imperfect determination of the WB parameters. The derived corrections and uncertainty system are used in the accompanying paper by Bzowski et al. in an analysis of the data from six seasons.

We analyzed observations of interstellar neutral helium (ISN~He) obtained from the Interstellar Boundary Explorer (IBEX) satellite during its first six years of operation. We used a refined version of the ISN~He simulation model, presented in the companion paper by Sokol_et al. 2015, and a sophisticated data correlation and uncertainty system and parameter fitting method, described in the companion paper by Swaczyna et al 2015. We analyzed the entire data set together and the yearly subsets, and found the temperature and velocity vector of ISN~He in front of the heliosphere. As seen in the previous studies, the allowable parameters are highly correlated and form a four-dimensional tube in the parameter space. The inflow longitudes obtained from the yearly data subsets show a spread of ~6 degree, with the other parameters varying accordingly along the parameter tube, and the minimum chi-square value is larger than expected. We found, however, that the Mach number of the ISN~He flow shows very little scatter and is thus very tightly constrained. It is in excellent agreement with the original analysis of ISN~He observations from IBEX and recent reanalyses of observations from Ulysses. We identify a possible inaccuracy in the Warm Breeze parameters as the likely cause of the scatter in the ISN~He parameters obtained from the yearly subsets, and we suppose that another component may exist in the signal, or a process that is not accounted for in the current physical model of ISN~He in front of the heliosphere. From our analysis, the inflow velocity vector, temperature, and Mach number of the flow are equal to lambda_ISNHe = 255.8 +/- 0.5 degree, beta_ISNHe = 5.16 +/- 0.10 degree, T_ISNHe = 7440 +/- 260 K, v_ISNHe = 25.8 +/- 0.4$ km/s, and M_ISNHe = 5.079 +/- 0.028, with uncertainties strongly correlated along the parameter tube.

We have developed a refined and optimized version of the Warsaw Test Particle Model of interstellar neutral gas in the heliosphere, specially tailored for analysis of IBEX-Lo observations. The former version of the model was used in the analysis of neutral He observed by IBEX that resulted in an unexpected conclusion that the interstellar neutral He flow vector was different than previously thought and that a new population of neutral He, dubbed the Warm Breeze, exists in the heliosphere. It was also used in the reanalysis of Ulysses observations that confirmed the original findings on the flow vector, but suggested a significantly higher temperature. The present version model has two strains targeted for different applications, based on an identical paradigm, but differing in the implementation and in the treatment of ionization losses. We present the model in detail and discuss numerous effects related to the measurement process that potentially modify the resulting flux of ISN He observed by IBEX, and identify those of them that should not be omitted in the simulations to avoid biasing the results. This paper is part of a coordinated series of papers presenting the current state of analysis of IBEX-Lo observations of ISN He. Details of the analysis method are presented by Swaczyna et al. 2015, and results of the analysis are presented by Bzowski et al. 2015.

Neutral interstellar helium has been observed by the Interstellar Boundary Explorer (IBEX) since 2009 with a signal-to-noise ratio well above 1000. Because of the geometry of the observations, the signal observed from January to March each year is the easiest to identify. However, as we show via simulations, the portion of the signal in the range of intensities from 10^{-3} to 10^{-2} of the peak value, previously mostly left out from the analysis, may bring important information about the details of the distribution function of interstellar He gas in front of the heliosphere. In particular, these observations may inform us about possible departures of the parent interstellar He population from equilibrium. We compare the expected distribution of the signal for the canonical assumption of a single Maxwell-Boltzmann population with the distributions for a superposition of the Maxwell-Boltzmann primary population and the recently discovered Warm Breeze, and for a single primary population given by a kappa function. We identify the regions on the sky where the differences between those cases are expected to be the most visible against the background. We discuss the diagnostic potential of the fall peak of the interstellar signal, reduced by a factor of 50 due to the Compton-Getting effect but still above the detection limit of IBEX. We point out the strong energy dependence of the fall signal and suggest that searching for this signal in the data could bring an independent assessment of the low-energy measurement threshold of the IBEX-Lo sensor.

We investigate the signals from neutral He atoms observed from Earth orbit in 2010 by IBEX. The full He signal observed during the 2010 observation season can be explained as a superposition of pristine neutral interstellar He gas and an additional population of neutral He that we call the Warm Breeze. The Warm Breeze is approximately two-fold slower and 2.5 times warmer than the primary interstellar He population, and its density in front of the heliosphere is ~7% that of the neutral interstellar helium. The inflow direction of the Warm Breeze differs by ~19deg from the inflow direction of interstellar gas. The Warm Breeze seems a long-term feature of the heliospheric environment. It has not been detected earlier because it is strongly ionized inside the heliosphere, which brings it below the threshold of detection via pickup ion and heliospheric backscatter glow observations, as well as by the direct sampling of GAS/Ulysses. Possible sources for the Warm Breeze include (1) the secondary population of interstellar helium, created via charge exchange and perhaps elastic scattering of neutral interstellar He atoms on interstellar He+ ions in the outer heliosheath, or (2) a gust of interstellar He originating from a hypothetic wave train in the Local Interstellar Cloud. A secondary population is expected from models, but the characteristics of the Warm Breeze do not fully conform to modeling results. If, nevertheless, this is the explanation, IBEX-Lo observations of the Warm Breeze provide key insights into the physical state of plasma in the outer heliosheath. If the second hypothesis is true, the source is likely to be located within a few thousand of AU from the Sun, which is the propagation range of possible gusts of interstellar neutral helium with the Warm Breeze characteristics against dissipation via elastic scattering in the Local Cloud.

Context. Since 1996, during periods of low solar activity, the HSTOF instrument onboard the SOHO satellite has been measuring weak fluxes of He atoms of 28-58 keV/n (helium energetic neutral atoms - He ENA). The probable source region is the inner heliosheath. Aims. We aim to understand the emission mechanism of He ENA based on knowledge of heliosheath spatial extent and plasma content resulting from Voyager 1 & 2 measurements in the period after termination shock crossings. Methods. He ENA are generated by charge-exchange neutralization of energetic helium ions on interstellar neutral H and He. Energy spectra of helium ions in the heliosheath are calculated by following the evolution of their velocity distribution functions when carried by and undergoing binary interactions with plasma constituents of a background flow whose particle populations are modeled to approximately render post-termination-shock Voyager data. Results. The observed HSTOF He ENA form a higher energy part of general heliospheric He ENA fluxes and can be explained by the proposed mechanism to within 2{\sigma} error. The main factor determining the level of emission (and its uncertainty) is the energy spectrum of He^+ pickup ions in post-termination shock plasmas.

Full sky maps of energetic neutral atoms (ENA) obtained with the Interstellar Boundary Explorer revealed a bright, arc-like Ribbon. We compare possible, though as yet undetected, He ENA emission in two models of the Ribbon origin. The models were originally developed for hydrogen ENA. In the first one, ENA are produced outside the heliopause from the ionized neutral solar wind in the direction where the local interstellar magnetic field is perpendicular to the line-of-sight. The second model proposes production at the contact layer between the Local Interstellar Cloud (LIC) and the Local Bubble (LB). The models are redesigned to helium using relevant interactions between atoms and ions. Resulting intensities are compared with possible emission of helium ENA from the heliosheath. In the first model, the expected intensity is ~ 0.014 (cm^2 s sr keV)^-1, i.e., of the order of the He emission from the heliosheath, whereas in the second, the LIC/LB contact layer model, the intensity is ~ (2-7) (cm^2 s sr keV)^-1, i.e., a few hundred times larger. If the IBEX Ribbon needs a source population of He ENA leaving the heliosphere, it should not be visible in He ENA fluxes mainly because of the insufficient supply of the parent He ENA originating from the neutralized solar wind {\alpha}-particles. We conclude that full-sky measurements of He ENA could give promising prospect for probing the Local Interstellar Medium at the distance of a few thousand AU and create possibility of distinction between above mentioned models of Ribbon origin.

Aims. A model of heliosheath density and energy spectra of alpha-particles and He+ ions carried by the solar wind is developed. Neutralization of heliosheath He+ ions, mainly by charge exchange (CX) with neutral interstellar H and He atoms, gives rise to ~0.2 - ~100 keV fluxes of energetic neutral He atoms (He ENA). Such fluxes, if observed, would give information about plasmas in the heliosheath and heliospheric tail. Methods. Helium ions crossing the termination shock (TS) constitute suprathermal (test) particles convected by (locally also diffusing through) hydrodynamically calculated background plasma flows (three versions of flows are employed). The He ions proceed from the TS towards heliopause (HP) and finally to the heliospheric tail (HT). Calculations of the evolution of alpha- and He+ particle densities and energy spectra include binary interactions with background plasma and interstellar atoms, adiabatic heating (cooling) resulting from flow compression (rarefaction), and Coulomb scattering on background plasma. Results. Neutralization of suprathermal He ions leads to the emergence of He ENA fluxes with energy spectra modified by the Compton-Getting effect at emission and ENA loss during flight to the Sun. Energy-integrated He ENA intensities are in the range ~0.05 - ~50 cm^-2 s^-1 sr^-1 depending on spectra at the TS (assumed kappa-distributions), background plasma model, and look direction. The tail/apex intensity ratio varies between ~1.8 and ~800 depending on model assumptions. Energy spectra are broad with maxima in the ~0.2 - ~3 keV range depending on the look direction and model. Conclusions. Expected heliosheath He ENA fluxes may be measurable based on the capabilities of the IBEX spacecraft. Data could offer insight into the heliosheath structure and improve understanding of the post-TS solar wind plasmas. HT direction and extent could be assessed.