This White Paper outlines the importance of addressing the fundamental science theme <<How are charged particles energized in space plasmas>> through a future ESA mission. The White Paper presents five compelling science questions related to particle energization by shocks, reconnection,waves and turbulence, jets and their combinations. Answering these questions requires resolving scale coupling, nonlinearity and nonstationarity, which cannot be done with existing multi-point observations. In situ measurements from a multi-point, multi-scale L-class plasma observatory consisting of at least 7 spacecraft covering fluid, ion and electron scales are needed. The plasma observatory will enable a paradigm shift in our comprehension of particle energization and space plasma physics in general, with very important impact on solar and astrophysical plasmas. It will be the next logical step following Cluster, THEMIS and MMS for the very large and active European space plasmas community. Being one of the cornerstone missions of the future ESA Voyage 2035-2050 science program, it would further strengthen the European scientific and technical leadership in this important field.

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.

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.

We present a comprehensive study of energetic neutral atoms (ENAs) of 10 eV to 2.5 keV from the downwind hemisphere of the heliosphere. These ENAs are believed to originate mostly from pickup protons and solar wind protons in the inner heliosheath. This study includes all low-energy observations made with the Interstellar Boundary Explorer over the first 8 years. Since the protons around 0.1 keV dominate the plasma pressure in the inner heliosheath in downwind direction, these ENA observations offer the unique opportunity to constrain the plasma properties and dimensions of the heliosheath where no in-situ observations are available. We first derive energy spectra of ENA intensities averaged over time for 49 macropixels covering the entire downwind hemisphere. The results confirm previous studies regarding integral intensities and the roll-over around 0.1 keV energy. With the expanded dataset we now find that ENA intensities at 0.2 and 0.1 keV seem to anti-correlate with solar activity. We then derive the product of total plasma pressure and emission thickness of protons in the heliosheath to estimate lower limits on the thickness of the inner heliosheath. The temporally averaged ENA intensities support a rather spherical shape of the termination shock and a heliosheath thickness between 150 and 210 au for most regions of the downwind hemisphere. Around the nominal downwind direction of 76{\deg} ecliptic longitude, the heliosheath is at least 280 au thick. There, the neutral hydrogen density seems to be depleted compared to upwind directions by roughly a factor of 2.

We present a simplified model of the outer heliosheath to help interpret the observations of interstellar neutral helium by the IBEX-Lo instruments. We assume that the measured particles are composed of the superposition of a primary beam population, with the properties of the local interstellar medium, and a secondary population, created by charge exchange between the primary beam neutrals and the ions that have been deflected as they approach the heliopause. We extract information on the large-scale shape of the heliopause by comparing the helium flux measured at IBEX along four different look directions with simple models of deflected plasma flow around hypothetical obstacles of different aspect ratios to the flow. As a first step in this paper, we model the deflected plasma flow with the analytical solutions for compressible gas flow around a series of oblate ellipsoidal obstacles. Our comparisons between the model results and the observations indicate that the heliopause is very blunt in the vicinity of the heliospheric nose, especially compared to a Rankine half-body or cometary shape. The upstream heliopause seems to be highly elongated in the directions parallel to the interstellar magnetic field, and relatively more compact and symmetric in the directions transverse to that field. The IBEX-Lo helium observations are not consistent with a heliopause elongated in directions parallel to the solar rotation axis.

The Interstellar Boundary Explorer (IBEX) has now operated in space for 7 years and returned nearly continuous observations that have led to scientific discoveries and reshaped our entire understanding of the outer heliosphere and its interaction with the local interstellar medium. Here we extend prior work, adding the 2014-2015 data for the first time, and examine, validate, initially analyze, and provide a complete 7-year set of Energetic Neutral Atom (ENA) observations from ~0.1 to 6 keV. The data, maps, and documentation provided here represent the 10th major release of IBEX data and include improvements to various prior corrections to provide the citable reference for the current version of IBEX data. We are now able to study time variations in the outer heliosphere and interstellar interaction over more than half a solar cycle. We find that the Ribbon has evolved differently than the globally distributed flux (GDF), with a leveling off and partial recovery of ENAs from the GDF, owing to solar wind output flattening and recovery. The Ribbon has now also lost its latitudinal ordering, which reflects the breakdown of solar minimum solar wind conditions and exhibits a greater time delay than for the surrounding GDF. Together, the IBEX observations strongly support a secondary ENA source for the Ribbon, and we suggest that this be adopted as the nominal explanation of the Ribbon going forward.

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.

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.

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.

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.

Direct sampling of neutral interstellar (NIS) atoms by the Interstellar Boundary Explorer (IBEX) can potentially provide a complementary method for studying element abundances in the Local Interstellar Cloud and processes in the heliosphere interface.}{We set the stage for abundance-aimed in-depth analysis of measurements of NIS He, Ne, and O by IBEX and determine systematic differences between abundances derived from various calculation methods and their uncertainties.}{Using a model of ionization rates of the NIS species in the heliosphere, based on independent measurements of the solar wind and solar EUV radiation, we develop a time-dependent method of calculating the survival probabilities of NIS atoms from the termination shock (TS) of the solar wind to IBEX. With them, we calculate densities of these species along the Earth's orbit and simulate the fluxes of NIS species as observed by IBEX. We study pairwise ratios of survival probabilities, densities and fluxes of NIS species at IBEX to calculate correction factors for inferring the abundances at TS.}{The analytic method to calculate the survival probabilities gives acceptable results only for He and Ne during low solar activity. For the remaining portions of the solar cycle, and at all times for O, a fully time dependent model should be used. Electron impact ionization is surprisingly important for NIS O. Interpreting the IBEX observations using the time dependent model yields the LIC Ne/O abundance of $0.16\pm40%$. The uncertainty is mostly due to uncertainties in the ionization rates and in the NIS gas flow vector.}{The Ne/He, O/He and Ne/O ratios for survival probabilities, local densities, and fluxes scaled to TS systematically differ and thus an analysis based only on survival probabilities or densities is not recommended, except the Ne/O abundance for observations at low solar activity.

Post-launch boresight of the IBEX-Lo instrument onboard the Interstellar Boundary Explorer (IBEX) is determined based on IBEX-Lo Star Sensor observations. Accurate information on the boresight of the neutral gas camera is essential for precise determination of interstellar gas flow parameters. Utilizing spin-phase information from the spacecraft attitude control system (ACS), positions of stars observed by the Star Sensor during two years of IBEX measurements were analyzed and compared with positions obtained from a star catalog. No statistically significant differences were observed beyond those expected from the pre-launch uncertainty in the Star Sensor mounting. Based on the star observations and their positions in the spacecraft reference system, pointing of the IBEX satellite spin axis was determined and compared with the pointing obtained from the ACS. Again, no statistically significant deviations were observed. We conclude that no systematic correction for boresight geometry is needed in the analysis of IBEX-Lo observations to determine neutral interstellar gas flow properties. A stack-up of uncertainties in attitude knowledge shows that the instantaneous IBEX-Lo pointing is determined to within $\sim 0.1\degr$ in both spin angle and elevation using either the Star Sensor or the ACS. Further, the Star Sensor can be used to independently determine the spacecraft spin axis. Thus, Star Sensor data can be used reliably to correct the spin phase when the Star Tracker (used by the ACS) is disabled by bright objects in its field-of-view. The Star Sensor can also determine the spin axis during most orbits and thus provides redundancy for the Star Tracker.

Neutral Interstellar Helium (NISHe) is almost unaffected at the heliospheric interface with the interstellar medium and freely enters the solar system. It provides some of the best information on the characteristics of the interstellar gas in the Local Interstellar Cloud. The Interstellar Boundary Explorer (IBEX) is the second mission to directly detect NISHe. We present a comparison between recent IBEX NISHe observations and simulations carried out using a well-tested quantitative simulation code. Simulation and observation results compare well for times when measured fluxes are dominated by NISHe (and contributions from other species are small). Differences between simulations and observations indicate a previously undetected secondary population of neutral helium, likely produced by interaction of interstellar helium with plasma in the outer heliosheath. Interstellar neutral parameters are statistically different from previous in situ results obtained mostly from the GAS/Ulysses experiment, but they do agree with the local interstellar flow vector obtained from studies of interstellar absorption: the newly-established flow direction is ecliptic longitude 79.2 deg, latitude -5.1 deg, the velocity is \sim 22.8 km/s, and the temperature is 6200 K. These new results imply a markedly lower absolute velocity of the gas and thus significantly lower dynamic pressure on the boundaries of the heliosphere and different orientation of the Hydrogen Deflection Plane compared to prior results from Ulysses. A different orientation of this plane also suggests a new geometry of the interstellar magnetic field and the lower dynamic pressure calls for a compensation by other components of the pressure balance, most likely a higher density of interstellar plasma and strength of interstellar magnetic field.

Hydrogen gas is the dominant component of the local interstellar medium. However, due to ionization and interaction with the heliosphere, direct sampling of neutral hydrogen in the inner heliosphere is more difficult than sampling the local interstellar neutral helium, which penetrates deep into the heliosphere. In this paper we report on the first detailed analysis of the direct sampling of neutral hydrogen from the local interstellar medium. We confirm that the arrival direction of hydrogen is offset from that of the local Helium component. We further report the discovery of a variation of the penetrating Hydrogen over the first two years of IBEX observations. Observations are consistent with hydrogen experiencing an effective ratio of outward solar radiation pressure to inward gravitational force greater than unity ({\mu}>1); the temporal change observed in the local interstellar hydrogen flux can be explained with solar variability.

We present a recent analysis with 1996 SOHO CELIAS CTOF data, which reveals correlations of He+ pickup ion fluxes and spectra with the magnetic field strength and solar wind density. The motivation is to better understand the ubiquitous large variations in both pickup ion fluxes and their velocity distributions found in interstellar pickup ion datasets. We concentrate on time periods of that can be associated with compression regions in the solar wind. Along with enhancements of the overall pickup ion fluxes, adiabatic heating and acceleration of the pickup ions are also observed in these regions. Transport processes that lead to the observed compressions and related heating or acceleration are discussed. A shift in velocity space associated with traveling interplanetary compression regions is observed, and a simple model presented to explain this phenomenon based on the conserved magnetic adiabatic moment.