In observations of diffuse emissions like, e.g., the Lyman-$\alpha$ heliospheric glow, contributions to the observed signal from point sources (e.g., stars) are considered as a contamination. There are relatively few brightest point sources that are usually properly resolved and can be subtracted or masked. We present results of analysis of the distribution of point sources using UV sky-survey maps from the SOHO/SWAN instrument and spectrophotometry data from the IUE satellite. The estimated distribution suggests that the number of these sources increases with decreasing intensity. Below a certain threshold, these sources cannot be resolved against the diffuse signal from the backscatter glow, that results in a certain physical background from unresolved point sources. Detection, understanding and subtraction of the point-source background has implications for proper characterization of diffuse emissions and accurate comparison with models. Stars are also often used as standard candles for in-flight calibration of satellite UV observations, thus proper understanding of signal contributions from the point sources is important for the calibration process. We present a general approach to quantify the background radiation level from unresolved point sources in UV sky-survey maps. In the proposed method, a distribution of point sources as a function of their intensity is properly integrated to compute the background signal level. These general considerations are applied to estimate the unresolved-point-sources background in the SOHO/SWAN observations that on average amounts to $28.9$ R. We discuss also the background radiation anisotropies and general questions related to modeling the point-source contributions to diffuse UV-emission observations.

The solar wind (SW) and the extreme ultraviolet (EUV) radiation modulate fluxes of interstellar and heliospheric particles inside the heliosphere both in time and in space. Understanding this modulation is necessary to correctly interpret measurements of particles of interstellar origin inside the heliosphere. We present a revision of heliospheric ionization rates and provide the Sun-Heliosphere Observation-based Ionization Rates (SHOIR) model based on the currently available data. We calculate the total ionization rates using revised SW and solar EUV data. We study the in-ecliptic variation of the SW parameters, the latitudinal structure of the SW speed and density, and the reconstruction of the photoionization rates. The revision affects the most the SW out of the ecliptic plane during solar maximum and the estimation of the photoionization rates. The latter, due to a change of the reference data. The revised polar SW is slower and denser during the solar maximum of SC 24. The currently estimated total ionization rates are higher than the previous ones for H, O, and Ne, and lower for He. The changes for the in-ecliptic total ionization rates are less than 10% for H and He, up to 20% for O, and up to 35% for Ne. Additionally, the changes are not constant in time and vary as a function of time and latitude.

We present a modification of a model of solar cycle evolution of the solar Lyman-alpha line profile, along with a sensitivity study of interstellar neutral H hydrogen to uncertainties in radiation pressure level. The line profile model, originally developed by Kowalska-Leszczynska et al. 2018a, is parametrized by the composite solar Lyman-alpha flux, which recently was revised Machol et al. 2019. We present modified parameters of the previously-developed model of solar radiation pressure for neutral hydrogen and deuterium atoms in the heliosphere. The mathematical function used in the model, as well as the fitting procedure, remain unchanged. We show selected effects of the model modification on ISN H properties in the heliosphere and we discuss the sensitivity of these quantities to uncertainties in the calibration of the composite Lyman-alpha series.

The interstellar neutral (ISN) gas enters the heliosphere and is detected at a few au from the Sun, as demonstrated by Ulysses and the nterstellar Boundary Explorer (IBEX) missions. Ulysses observed ISN gas from different vantage points in a polar orbit from 1994 to 2007, while IBEX has been observing in an Earth orbit in a fixed direction relative to the Sun from 2009. McComas et al. 2018 reported about an IMAP-Lo detector on board the Interstellar Mapping and Acceleration Probe (IMAP), with an ability to track the ISN flux in the sky. We present observation geometries for ISN gas for a detector with the capability to adjust the boresight direction along the Earth orbit over a year within a multichoice ISN observation scheme. We study science opportunities from the observations as a function of time during a year and the phase of solar activity. We identify observation geometries and determine the observation seasons separately for various ISN species and populations. We find that using an adjustable viewing direction allows for ISN gas observations in the upwind hemisphere, where the signal is not distorted by gravitational focusing, in addition to the viewing of ISN species throughout the entire year. Moreover, we demonstrate that with appropriately adjusted observation geometries, primary and secondary populations can be fully separated. Additionally, we show that atoms of ISN gas on indirect trajectories are accessible for detection, and we present their impact on the study of the ionization rates for ISN species.

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.

Interaction of the solar wind with interstellar matter involves, among other, charge exchange between interstellar neutral atoms and plasma, which results in the creation of secondary population of interstellar neutral (ISN) atoms. The secondary population of interstellar He was detected by Interstellar Boundary Explorer (IBEX), but interpretation of these measurements was mostly based on an approximation that the primary interstellar neutral population and the secondary population were non-interacting homogeneous Maxwell-Boltzmann functions in the outer heliosheath. We simulate the distribution function in the outer heliosheath and inside the heliopause using method of characteristics with statistical weights obtained from solutions of the production and loss equations for the secondary atoms due to charge-exchange collisions in the outer heliosheath. We show that the two-Maxwellian approximation for the distribution function of neutral He is not a good approximation within the outer heliosheath but a reasonable one inside the termination shock. This is due to a strong selection effect: the He atoms able to penetrate inside the termination shock are a small, peculiar subset of the entire secondary He population. Nevertheless, the two-Maxwellian approximation reproduces the density distribution of ISN He inside the termination shock well and enables a realistic reproduction of the orientation of the plane defined by the Sun's velocity vector through the local interstellar matter and the vector of unperturbed interstellar magnetic field.

We study the distribution of the interstellar neutral (ISN) gas density and the pick-up ion (PUI) density of hydrogen, helium, neon, and oxygen in the heliosphere for heliocentric distances from inside 1 au up to the solar wind termination shock (TS), both in and out of the ecliptic plane. We discuss similarities and differences in the large-scale structures of the ISN gas and PUIs formed in the heliosphere between various species. We discuss the distribution of ISN gas and PUI densities for two extreme phases of the solar activity cycle, it is the solar minimum and the solar maximum. We identify the location of the ISN gas density cavity of various species. We study the relative abundance ratios of Ne/O, H/He, Ne/He, and O/He for ISN gas and PUIs densities and their variation with location in the heliosphere. We also discuss the modulation of relative abundance ratios of ISN gas and PUIs along the TS. We conclude that the preferable locations for detection of He$^+$ and Ne$^+$ PUIs are in the downwind hemisphere within 1~au, whereas for H$^+$ and O$^+$ PUIs the preferable locations for detection are for distances from Jupiter to Pluto orbits.

Interstellar neutral hydrogen (ISN H) gas penetrates freely the heliopause. Inside the inner heliosheath, the charge-exchange interaction of this gas with the shocked solar wind and pickup ions creates energetic neutral atoms (ENAs). ISN H is strongly depleted inside the termination shock but a fraction reaches the Earth's orbit. In these regions of the heliosphere, ISN H is the source population for interstellar pickup ions and for the heliospheric backscatter glow. The Globally Distributed Flux (GDF) of ENAs created in the inner heliosheath has been sampled directly by Interstellar Boundary Explorer (IBEX). Based on these measurements, we calculate the density of the GDF ENA population at the Earth's orbit. We find that this number density is between $10^{-4}$ and $10^{-3}$ cm$^{-3}$, i.e., comparable in magnitude to the number density of ISN H in the downwind portion of the Earth's orbit. Half of this atom population have energies less than $\sim 80$ eV. This GDF population of neutral hydrogen is likely to provide a significant contribution to the intensity of heliospheric glow in the downwind hemisphere, may be the source of the inner source of hydrogen pickup ions, and may be responsible for the excess of production of pickup ions found in the analysis of magnetic wave events induced by the proton pickup process in the downwind region at 1 au from the Sun.

Following the derivation of a more accurate model of the evolution of the solar Lyman-$\alpha$ line with the changing solar activity by Kowalska-leszczynska et al. 2018 (IKL18) than the formerly used model by Tarnopolski et al. 2009 (ST09), we investigate potential consequences that adoption of the resulting refined model of radiation pressure has for the model distribution of interstellar neutral (ISN) H in the inner heliosphere and on the interpretation of selected observations. We simulated the ISN H densities using the two alternative radiation pressure models and identical models of all other factors affecting the ISN H distribution. We found that during most of the solar cycle, the IKL18 model predicts larger densities of ISN H and PUIs than ST09 in the inner heliosphere, especially in the downwind hemisphere. However, the density of ISN H at the termination shock estimated by Bzowski et al. 2008 obtained using ST09 does not need revision, and the detection of ISN D by IBEX is supported. However, we point out the existence of a considerable absorption of a portion of the solar Lyman-$\alpha$ spectral flux inside the heliosphere. Therefore, the model of radiation pressure for ISN H is still likely to need revision, and hence the available models of ISN H are not self-consistent.

Recent studies of interstellar neutral (ISN) hydrogen observed by the Interstellar Boundary Explorer (IBEX) suggested that the present understanding of the radiation pressure acting on hydrogen atoms in the heliosphere should be revised. There is a significant discrepancy between theoretical predictions of the ISN H signal using the currently used model of the solar Lyman-alpha profile by Tarnopolski et al. 2009 (TB09) and the signal due to ISN H observed by IBEX-Lo. We developed a new model of evolution of the solar Lyman-alpha profile that takes into account all available observations of the full-disk solar Lyman-alpha profiles from SUMER/SOHO, provided by Lemaire et al. 2015 (L15), covering practically the entire 23rd solar cycle. The model has three components that reproduce different features of the profile. The main shape of the emission line that is produced in the chromosphere is modeled by the kappa function; the central reversal due to absorption in the transition region is modeled by the Gauss function; the spectral background is represented by the linear function. The coefficients of all those components are linear functions of the line-integrated full-disk Lyman-alpha irradiance, which is the only free parameter of the model. The new model features potentially important differences in comparison with the model by TB09, which was based on a limited set of observations. This change in the understanding of radiation pressure, especially during low solar activity, may significantly affect the interstellar H and D distributions in the inner heliosphere and their derivative populations.

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.

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.

Time-series photometry from space in the ultraviolet can be presently done with only a few platforms, none of which is able to provide wide-field long-term high-cadence photometry. We present a concept of UVSat, a twin space telescope which will be capable to perform this kind of photometry, filling an observational niche. The satellite will host two telescopes, one for observations in the ultraviolet, the other for observations in the optical band. We also briefly show what science can be done with UVSat.

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.

We simulated the signal due to neutral He atoms, observed by Interstellar Boundary Explorer (IBEX), assuming that charge exchange collisions between neutral He atoms and He+ ions operate everywhere between the heliopause and a distant source region in the local interstellar cloud, where the neutral and charged components are in thermal equilibrium. We simulated several test cases of the plasma flow within the outer heliosheath and investigated the signal generation for plasma flows both in the absence and in the presence of the interstellar magnetic field. We found that a signal in the portion of IBEX data identified as due to the Warm Breeze does not arise when a homogeneous plasma flow in front of the heliopause is assumed, but it appears immediately when any reasonable disturbance in its flow due to the presence of the heliosphere is assumed. We obtained a good qualitative agreement between the data selected for comparison and the simulations for a model flow with the velocity vector of the unperturbed gas and the direction and intensity of magnetic field adopted from recent determinations. We conclude that direct-sampling observations of neutral He atoms at 1~AU from the Sun are a sensitive tool for investigating the flow of interstellar matter in the outer heliosheath, that the Warm Breeze is indeed the secondary population of interstellar helium, as it was hypothesized earlier, and that the WB signal is consistent with the heliosphere distorted from axial symmetry by the interstellar magnetic field.

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.

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.

We simulated the modulation of the interstellar neutral (ISN) He, Ne, and O density and pick-up ion (PUI) production rate and count rate along the Earth's orbit over the solar cycle from 2002 to 2013 to verify if solar cycle-related effects may modify the inferred ecliptic longitude of the ISN inflow direction. We adopted the classical PUI model with isotropic distribution function and adiabatic cooling, modified by time- and heliolatitude-dependent ionization rates and non-zero injection speed of PUIs. We found that the ionization losses have a noticeable effect on the derivation of the ISN inflow longitude based on the Gaussian fit to the crescent and cone peak locations. We conclude that the non-zero radial velocity of the ISN flow and the energy range of the PUI distribution function that is accumulated are of importance for a precise reproduction of the PUI count rate along the Earth orbit. However, the temporal and latitudinal variations of the ionization in the heliosphere, and particularly their variation on the solar cycle time-scale, may significantly modify the shape of PUI cone and crescent and also their peak positions from year to year and thus bias by a few degrees the derived longitude of the ISN gas inflow direction.

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.

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 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.

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.

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.

The journey of the Sun through space carries the solar system through a dynamic interstellar environment that is presently characterized by Mach 1 motion between the heliosphere and the surrounding interstellar medium (ISM). The interaction between the heliosphere and ISM is an evolving process due to the variable solar wind and to interstellar turbulence. Frisch et al. presented a meta-analysis of the historical data on the interstellar wind flowing through the heliosphere and concluded that temporal changes in the ecliptic longitude of the wind were statistically indicated by the data available in the refereed literature at the time of that writing. Lallement and Bertaux disagree with this result, and suggested, for instance, that a key instrumental response function of IBEX-Lo was incorrect and that the STEREO pickup ion data are unsuitable for diagnosing the flow of interstellar neutrals through the heliosphere. Here we show that temporal variations in the interstellar wind through the heliosphere are consistent with our knowledge of local ISM. The statistical analysis of the historical helium wind data is revisited, and a recent correction of a typographical error in the literature is incorporated into the new fits. With this correction, and including no newer IBEX results, these combined data still indicate that a change in the longitude of the interstellar neutral wind over the past forty years is statistically likely, but that a constant flow longitude is now also statistically possible. It is shown that the IBEX instrumental response function is known, and that the STEREO pickup ion data have been correctly utilized in this analysis.

The NIS He gas has been observed at a few AU to the Sun almost from the beginning of the space age. To model its flow an estimate of the loss rates due to ionization by solar extreme-ultraviolet (EUV) flux is needed. The EUV irradiance has been measured directly from mid 1990-ties, but with high temporal and spectral resolution only from 2002. Beforehand only EUV proxies are available. A new method of reconstruction of the Carrington rotation averaged photoionization rates for neutral interstellar helium (NIS He) in the ecliptic at 1 AU to the Sun before 2002 is presented. We investigate the relation between the solar rotation averaged time series of the ionization rates for NIS He at 1 AU derived from TIMED measurements of EUV irradiance and the solar 10.7 cm flux (F10.7) only. We perform a weighted iterative fit of a nonlinear model to data split into sectors. The obtained formula allows to reconstruct the solar rotation averages of photoionization rates for He between ~1947 and 2002 with an uncertainty ranging from less than 10% during solar minimum up to 20% for solar maximum.

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.

Analysis of IBEX measurements of neutral interstellar He flux brought the inflow velocity vector different from the results of earlier analysis of observations from GAS/Ulysses. Recapitulation of results on the helium inflow direction from the past ~40 years suggested that the inflow direction may be changing with time. We reanalyze the old Ulysses data and reprocess them to increase the accuracy of the instrument pointing to investigate if the GAS observations support the hypothesis that the interstellar helium inflow direction is changing. We employ a similar analysis method as in the analysis of the IBEX data. We seek a parameter set that minimizes reduced chi-squared, using the Warsaw Test Particle Model for the interstellar He flux at Ulysses with a state of the art model of neutral He ionization in the heliosphere, and precisely reproducing the observation conditions. We also propose a supplementary method of constraining the parameters based on cross-correlations of parameters obtained from analysis of carefully selected subsets of data. We find that the ecliptic longitude and speed of interstellar He are in a very good agreement with the values reported in the original GAS analysis. We find, however, that the temperature is markedly higher. The 3-seasons optimum parameter set is lambda = 255.3, beta = 6, v = 26.0 km/s, T = 7500 K. We find no evidence that it is varying with time, but the uncertainty range is larger than originally reported. The originally-derived parameters of interstellar He from GAS are in good agreement with presently derived, except for the temperature, which seems to be appreciably higher, in good agreement with interstellar absorption line results. While the results of the present analysis are in marginal agreement with the earlier reported results from IBEX, the most likely values from the two analyses differ for reasons that are still not understood.

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.

In the heliosphere, especially in the inner heliosheath, mass-, momentum-, and energy loading induced by the ionization of neutral interstellar species plays an important, but for some species, especially Helium, an underestimated role. We discuss the implementation of charge exchange and electron impact processes for interstellar neutral Hydrogen and Helium and their implications for further modeling. Especially, we emphasize the importance of electron impact and a more sophisticated numerical treatment of the charge exchange reactions. Moreover, we discuss the non-resonant charge exchange effects. The rate coefficients are discussed and the influence of the cross-sections in the (M)HD equations for different reactions are revised as well as their representation in the collision integrals. Electron impact is in some regions of the heliosphere, particularly in the heliotail, more effective than charge exchange, and the ionization of neutral interstellar Helium contributes about 40% to the mass- and momentum loading in the inner heliosheath. The charge exchange cross-sections need to be modeled with higher accuracy, especially in view of the latest developments in their description. The ionization of Helium and electron impact ionization of Hydrogen needs to be taken into account for the modeling of the heliosheath and, in general, astrosheaths. Moreover, the charge exchange cross-sections need to be handled in a more sophisticated way, either by developing better analytic approximations or by solving the collision integrals numerically.

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.

The abundance of deuterium in the interstellar gas in front of the Sun gives insight into the processes of filtration of neutral interstellar species through the heliospheric interface and potentially into the chemical evolution of the Galactic gas. We investigate the possibility of detection of neutral interstellar deuterium at 1 AU from the Sun by direct sampling by the Interstellar Boundary Explorer (IBEX). We simulate the flux of neutral interstellar D at IBEX for the actual measurement conditions. We assess the number of interstellar D atom counts expected during the first three years of IBEX operation. We also simulate observations expected during an epoch of high solar activity. In addition, we calculate the expected counts of D atoms from the thin terrestrial water layer, sputtered from the IBEX-Lo conversion surface by neutral interstellar He atoms. Most D counts registered by IBEX-Lo are expected to originate from the water layer, exceeding the interstellar signal by 2 orders of magnitude. However, the sputtering should stop once the Earth leaves the portion of orbit traversed by interstellar He atoms. We identify seasons during the year when mostly the genuine interstellar D atoms are expected in the signal. During the first 3 years of IBEX operations about 2 detectable interstellar D atoms are expected. This number is comparable with the expected number of sputtered D atoms registered during the same time intervals. The most favorable conditions for the detection occur during low solar activity, in an interval including March and April each year. The detection chances could be improved by extending the instrument duty cycle, e.g., by making observations in the special deuterium mode of IBEX-Lo.

NASA's Interstellar Boundary Explorer (IBEX) mission has recently opened a new window on the interstellar medium (ISM) by imaging neutral atoms. One "bright" feature in the sky is the interstellar wind flowing into the solar system. Composed of remnants of stellar explosions as well as primordial gas and plasma, the ISM is by no means uniform. The interaction of the local ISM with the solar wind shapes our heliospheric environment with hydrogen being the dominant component of the very local ISM. In this paper, we report on direct sampling of the neutral hydrogen of the local ISM over four years of IBEX observations. The hydrogen wind observed at 1 AU has decreased and nearly disappeared as the solar activity has increased over the last four years; the signal at 1 AU has dropped off in 2012 by a factor of ~8 to near background levels. The longitudinal offset has also increased with time presumably due to greater radiation pressure deflecting the interstellar wind. We present longitudinal and latitudinal arrival direction measurements of the bulk flow as measured over four years beginning at near solar minimum conditions. The H distribution we observe at 1 AU is expected to be different from that outside the heliopause due to ionization, photon pressure, gravity, and filtration by interactions with heliospheric plasma populations. These observations provide an important benchmark for modeling of the global heliospheric interaction. Based on these observations we suggest a further course of scientific action to observe neutral hydrogen over a full solar cycle with IBEX.

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.

The 3D structure of solar wind and its evolution in time is needed for heliospheric modeling and interpretation of energetic neutral atoms observations. We present a model to retrieve the solar wind structure in heliolatitude and time using all available and complementary data sources. We determine the heliolatitude structure of solar wind speed on a yearly time grid over the past 1.5 solar cycles based on remote-sensing observations of interplanetary scintillations, in situ out-of-ecliptic measurements from Ulysses, and in situ in-ecliptic measurements from the OMNI-2 database. Since the in situ information on the solar wind density structure out of ecliptic is not available apart from the Ulysses data, we derive correlation formulae between solar wind speed and density and use the information on the solar wind speed from interplanetary scintillation observations to retrieve the 3D structure of solar wind density. With the variations of solar wind density and speed in time and heliolatitude available we calculate variations in solar wind flux, dynamic pressure and charge exchange rate in the approximation of stationary H atoms.

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.

The goal of the Fully Online Datacenter of Ultraviolet Emissions (FONDUE) Working Team of the International Space Science Institute in Bern, Switzerland, was to establish a common calibration of various UV and EUV heliospheric observations, both spectroscopic and photometric. Realization of this goal required an up-to-date model of spatial distribution of neutral interstellar hydrogen in the heliosphere, and to that end, a credible model of the radiation pressure and ionization processes was needed. This chapter describes the solar factors shaping the distribution of neutral interstellar H in the heliosphere. Presented are the solar Lyman-alpha flux and the solar Lyman-alpha resonant radiation pressure force acting on neutral H atoms in the heliosphere, solar EUV radiation and the photoionization of heliospheric hydrogen, and their evolution in time and the still hypothetical variation with heliolatitude. Further, solar wind and its evolution with solar activity is presented in the context of the charge exchange ionization of heliospheric hydrogen, and in the context of dynamic pressure variations. Also the electron ionization and its variation with time, heliolatitude, and solar distance is presented. After a review of all of those topics, we present an interim model of solar wind and the other solar factors based on up-to-date in situ and remote sensing observations of solar wind. Results of this effort will further be utilised to improve on the model of solar wind evolution, which will be an invaluable asset in all heliospheric measurements, including, among others, the observations of Energetic Neutral Atoms by the Interstellar Boundary Explorer (IBEX).

The brightest and most surprising feature in the first all-sky maps of Energetic Neutral Atoms (ENA) emissions (0.2-6 keV) produced by the Interstellar Boundary Explorer (IBEX) is an almost circular ribbon of a ~140{\deg} opening angle, centered at (l,b) = (33{\deg}, 55{\deg}), covering the part of the celestial sphere with the lowest column densities of the Local Interstellar Cloud (LIC). We propose a novel interpretation of the IBEX results based on the idea of ENA produced by charge-exchange between the neutral H atoms at the nearby edge of the LIC and the hot protons of the Local Bubble (LB). These ENAs can reach the Sun's vicinity because of very low column density of the intervening LIC material. We show that a plane-parallel or slightly curved interface layer of contact between the LIC H atoms (n_H = 0.2 cm^-3, T = 6000-7000 K) and the LB protons (n_p = 0.005 cm^-3, T ~ 10^6 K), together with indirect contribution coming from multiply-scattered ENAs from the LB, may be able to explain both the shape of the ribbon and the observed intensities provided that the edge is < (500-2000) AU away, the LIC proton density is (correspondingly) < (0.04-0.01) cm^-3, and the LB contains ~1% of non-thermal protons over the IBEX energy range. If this model is correct, then IBEX, for the first time, has imaged in ENAs a celestial object from beyond the confines of the heliosphere and can directly diagnose the plasma conditions in the LB.

We discuss a consolidation of determinations of the density of neutral interstellar H at the nose of the termination shock carried out with the use of various data sets, techniques, and modeling approaches. In particular, we focus on the determination of this density based on observations of H pickup ions on Ulysses during its aphelion passage through the ecliptic plane. We discuss in greater detail a novel method of determination of the density from these measurements and review the results from its application to actual data. The H density at TS derived from this analysis is equal to 0.087 \pm 0.022 cm-3, and when all relevant determinations are taken into account, the consolidated density is obtained at 0.09 \pm 0.022 cm-3. The density of H in CHISM based on literature values of filtration factor is then calculated at 0.16 \pm 0.04 cm-3.

With the plethora of detailed results from heliospheric missions and at the advent of the first mission dedicated IBEX, we have entered the era of precision heliospheric studies. Interpretation of these data require precision modeling, with second-order effects quantitatively taken into account. We study the influence of the non-flat shape of the solar Ly-alpha line on the distribution of neutral interstellar H in the inner heliosphere. Based on available data, we (i) construct a model of evolution for the solar Ly-alpha line profile with solar activity, (ii) modify an existing test-particle code used to calculate the distribution of neutral interstellar H in the inner heliosphere so that it takes the dependence of radiation pressure on radial velocity into account, and (iii) compare the results of the old and new version. Discrepancies between the classical and Doppler models appear between ~5 and ~3 AU and increase towards the Sun from a few percent to a factor of 1.5 at 1 AU. The classical model overestimates the density everywhere except for a ~60-degr cone around the downwind direction, where a density deficit appears. The magnitude of the discrepancies appreciably depends on the phase of the solar cycle, but only weakly on the parameters of the gas at the termination shock. For in situ measurements of neutral atoms performed at ~1 AU, the Doppler correction will need to be taken into account, because the modifications include both the magnitude and direction of the local flux by a few km/s and degrees, respectively, which, when unaccounted for, would introduce an error of a few km/s and degrees in determination of the magnitude and direction of the bulk velocity vector at the termination shock.

By reevaluating a 13-month stretch of Ulysses SWICS H pickup ion measurements near 5 AU close to the ecliptic right after the previous solar minimum, this paper presents a determination of the neutral interstellar H density at the solar wind termination shock and implications for the density and ionization degree of hydrogen in the LIC. The density of neutral interstellar hydrogen at the termination shock was determined from the local pickup ion production rate as obtained close to the cut-off in the distribution function at aphelion of Ulysses. As shown in an analytical treatment for the upwind axis and through kinetic modeling of the pickup ion production rate at the observer location, with variations in the ionization rate, radiation pressure, and the modeling of the particle behavior, this analysis turns out to be very robust against uncertainties in these parameters and the modeling. Analysis using current heliospheric parameters yields the H density at the termination shock equal to $0.087\pm0.022$ cm$^{-3}$, including observational and modeling uncertainties.

Aims: Investigate/Study de-charging of solar wind C, N, O, Mg, Si and S ions and assess fluxes of resulting ENA in the heliosphere. Methods: The model treats the heavy ions as test particles convected by (and in a particular case also diffusing through) a hydrodynamically calculated background plasma flow, from 1 AU to the termination shock (TS), to heliosheath (HS) and finally to heliospheric tail (HT). The ions undergo radiative and dielectronic recombinations, charge exchanges, photo- and electron impact ionizations with plasma particles, interstellar neutral atoms (calculated on a Monte-Carlo model) and solar photons. Results: Highly-charged heavy coronal ions flowing with the solar wind undergo successive de-ionizations, mainly in the heliosheath, leading to charge-states much lower than in the supersonic solar wind. If Coulomb scattering is the main ion energy loss mechanism, the end product of these deionizations are fluxes of ENA of ~1 keV/nucleon originating in the upwind heliosheath that for C, Mg, Si and S may constitute sources of pickup ions (PUI) significantly exceeding the interstellar supply. Conclusions: Discussed processes result in (i) distinct difference of ion charge q in the supersonic solar wind (approximately q >= +Z/2, Z = atomic number) compared to that in the HS (approximately 0 <= q <= +Z/2)), (ii) probable concentration of singly ionized atoms (q = +1) in the heliosheath towards the heliopause (HP) and in the HT, (iii) possible significant production of ENA in the HS offering natural explanation for production of PUI, and -- after acceleration at the TS -- anomalous cosmic rays (ACR) of species (like C, Mg, Si, S) unable to enter the heliospheric cavity from outside because of their total ionization in the local interstellar medium.

Context: With the forthcoming launch of a NASA SMEX mission IBEX devoted to imaging of heliospheric interface by in-situ detection of Energetic Neutral Atoms (ENA) an important issue becomes recognizing of transport of these atoms from the termination shock of the solar wind to Earth orbit. Aims: Investigate modifications of energy and of survival probability of the H ENA detectable by IBEX (0.01 -- 6 keV) between the termination shock and Earth orbit taking into account the influence of the variable and anisotropic solar wind and solar EUV radiation. Methods: Energy change of the atoms is calculated by numerical simulations of orbits of the H ENA atoms from ~100 AU from the Sun down to Earth orbit, taking into account solar gravity and Lyman-$\alpha$ radiation pressure, which is variable in time and depends on radial velocity of the atom. To calculate survival probabilities of the atoms against onization, a detailed 3D and time-dependent model of H ENA ionization based on observations of the solar wind and EUV ionizing radiation is constructed, and wth the use of this model probabilities of survival of the atoms are calculated by numerical integration along the previously calculated orbits. Results: Owing to the radiation pressure, H ENA reach the Earth orbit practically without energy and direction change except the atoms with energy lower than 0.1 keV during high solar activity. For a given energy at Earth orbit one expects fluctuations of survival probability from ~20% at 0.01 keV down to just a few percent at 6 keV and a modulation of survival probability as a function of the location at Earth orbit, ecliptic latitude of the arrival direction, and the phase of solar cycle with an amplitude of a few dozen percent for 0.1 keV atoms at solar minimum to a few percent for 6 keV atoms at solar maximum.

We study the influence of the non-flat shape of the solar Lyman-alpha line on the distribution of neutral interstellar hydrogen in the inner heliosphere and assess importance of this effect for interpretation of heliospheric in situ measurements. Based on available data, construct a model of evolution of the solar Lyman-alpha line profile with solar activity. Modify an existing test-particle code calculating distribution of neutral interstellar hydrogen in the inner heliosphere to take into account the dependence of radiation pressure on radial velocity. Discrepancies between the classical and Doppler models appear at ~ 5 AU and increase towards the Sun from a few percent to a factor of 1.5 at 1 AU. The classical model overestimates density everywhere except a ~ 60 deg cone around the downwind direction, where a density deficit appears. The magnitude of discrepancies depends appreciably on the phase of solar cycle, but only weakly on the parameters of the gas at the termination shock. For in situ measurements of neutral atoms performed at ~ 1 AU, as those planned for IBEX, the Doppler correction will need to be taken into account, because the modifications include both the magnitude and direction of the local flux by a few km/s and degree, which, when unaccounted for, would bring an error of a few degrees and a few km/s in determination of the bulk velocity vector at the termination shock. The Doppler correction is appreciable for in situ observations of neutral H populations and their derivatives performed a few AU from the Sun.

Context: Reconnaissance of feasibility of detection of interstellar neutral D by the forthcoming NASA SMEX mission IBEX. Aims: To investigate by numerical simulations the absolute density and flux at Earth orbit of neutral interstellar D and to check its detectability by IBEX. Methods: The simulations were performed with the use of the Warsaw 3D time-dependent test-particle model of neutral interstellar gas in the inner heliosphere, specially adapted to the case of D, and of the state-of-the-art models of the ionization field and radiation pressure. The modeling returned density, bulk velocity, and flux of interstellar D along at the Earth locations during the solar cycle. Particular attention was paid to the time interval corresponding to the planned operations of IBEX. Results: Simulations predict a large enhancement of D abundance at Earth orbit with respect to the abundance at the termination shock. The energy of the D atoms at IBEX will be within the energetic sensitivity band of its Lo instrument except a short time interval between September and November each year. Because of the specific observing geometry of IBEX, there will be one opportunity each year to search for I/S D, when Earth is near ecliptic longitude 136 deg, i.e. in February. Assuming the TS abundance of D identical as in the Local Cloud equal to 1.56E-5, and the density of H at TS 0.11/cm2/s one obtains the expected relative flux about 0.015/cm2/s, which corresponds to the local absolute flux about 0.007/cm2/s. The dependence of the expected flux on the phase of solar cycle is relatively weak. The flux scales proportionally to the density of deuterium at the termination shock and depends only weakly on the bulk velocity and temperature of the gas in this region.

We show that multiply ionized coronal C, N, O, Mg, Si, S ions carried by the solar wind and neutralized by consecutive electron captures from neutral interstellar atoms constitute an important new source of neutral atoms in the inner heliosheath, with energies up to ~ 1 keV/n. In the model we developed, the heavy ions are treated as test particles carried by hydrodynamic plasma flow (with a Monte-Carlo description of interstellar neutrals) and undergoing all relevant atomic processes determining the evolution of all charge-states of considered species (radiative and dielectronic recombination, charge exchange, photo-, and electron impact ionization). The total strength of the source is from ~10^6 g/s for S to ~10^8 g/s for O, deposited as neutrals below the heliopause. These atoms should provide, as they drift to supersonic wind region, important sources of PUIs and eventually ACRs, especially for species that are excluded from entering the heliosphere because of their ionization in the LISM. The expected corresponding ENA fluxes at 1 AU are in the range 10^-4 - 10^0 at./(cm^2 s sr), depending on the species and direction (Table 2).

Dynamics of H, D, and heavy Energetic Neutral Atoms (ENA) between the termination shock and 1 AU is discussed in the context of the forthcoming NASA SMEX mission IBEX. In particular, effects of the velocity-dependent radiation pressure on atomic trajectories are considered and ionization losses between TS and 1 AU are studied. It is shown, among others, that most of the dynamical effects and ionization losses are induced within a few AU from the Sun, which translates to the time domain into $\sim 1 - 3$ solar rotations before detection. This loosens considerably time requirements for tracking the ionization and radiation pressure history to just prior 3 months. ENA seem excellent tracers of the processes within the heliospheric interface, with the transport effects between the termination shock and detector relatively mild and easy to account for.

Most of the sungrazing comets observed by LASCO at SOHO belong to the Kreutz group of comets and follow trajectories that are tightly clumped in space. Statistical analysis of 9 years of SOHO observations suggests that the true apparition rate of these comets is as high as one every other day. Practically all these comets break up before perihelion passage. Their material is dissociated and ionized, and subsequently transported away from the Sun as pickup ions in the solar wind. Their mean mass flux is about $3.1 10^4 g/s. Since the breakup occurs between 40 and 4 solar radii and the ionization is almost immediate, the expected location of these ions in the phase space is close to the location of the inner source of pickup ions. Assuming radial propagation, the cometary pickups should be observable at Earth between August and January, with peak probability at the end of September. At Ulysses, they should be observable approximately between -25 and 40 degrees ecliptic latitude during the fast latitude scans, the first of which occurred in 1995 and the second in 2001. The population of cometary pickup ions should be augmented by about 40% by solar wind protons as a result of charge exchange with the cometary neutral hydrogen and oxygen atoms and subsequent reionization of the newly-created Energetic Neutral Atoms, streaming with respect to the solar wind. In total, the average flux of the pickup ions related to the sungrazing comets at 1 AU should be about 1.6 10^5 g/s/sr within the detection area (and null outside it). This value is comparable to the flux of the inner source-pickup ions.

We present arguments that at least part of the inner source of pickup ions in the solar wind might be the material released by sungrazing comets. Based on a statistical analysis of sungrazing comets detected over almost eight years of LASCO operation (1996 -- September 09 2004) an overwhelming majority of the observed sungrazers belong to the Kreutz group of comets, follows tightly clumped orbits and break up at $\sim 40 - 4$ solar radii in a well defined region of space. The material released from these comets could be (after ionization) an important portion of the inner source of pickup ions (PUIs), as the local mass flux of the inner source and cometary PUIs seem comparable. We indicate time intervals during the year when the cometary PUIs could be observed from a spacecraft on the Earth's orbit (from the end of July until the end of the year) and show three time intervals when they should be observable by Ulysses (from its launch time until the end of 1990, from the end of November 1994 until mid-May 1995 and from February 2001 until the end of July, 2001). We argue that the PUIs from the inner source should include both singly and doubly charged ions and that this cometary hypothesis alleviates some difficulties (in particular, the issue of hydrogen deficit) in the interpretation of the inner source as solar wind neutralized on dust grains close to the Sun.

Various hydrodynamic models of the heliospheric interface have been presented meanwhile, numerically simulating the interaction of the solar wind plasma bubble with the counterstreaming partially ionized interstellar medium. In these model approaches the resulting interface flows are found by the use of hydrodynamic simulation codes trying to consistently describe the dynamic and thermodynamic coupling of the different interacting fluids of protons, H-atoms and pick-up ions. Within such approaches, the fluids are generally expected to be correctly described by the three lowest velocity moments, i.e., by shifted Maxwellians. We shall show that in these approaches the charge-exchange-induced momentum coupling is treated in an unsatisfactory representation valid only at supersonic differential flow speeds. Though this flaw can be removed by an improved coupling term, we shall further demonstrate that the assumption of shifted Maxwellians in some regions of the interface is insufficiently well fulfilled both for H-atoms and protons. Using a Boltzmann-kinetic description of the proton- and H-atom- distribution functions coupled by charge exchange processes we emphasize the fact that non-negligible deviations from shifted Maxwellians are generated in the interface. This has to be taken into account when interpreting inner heliospheric measurements in terms of interstellar parameters.

It is well known that the Solar System is presently moving through a partially ionized local interstellar medium. This gives rise to a counter-flow situation requiring a consistent description of behaviour of the two fluids -- ions and neutral atoms -- which are dynamically coupled by mutual charge exchange processes. Solutions to this problem have been offered in the literature, all relying on the assumption that the proton fluid, even under evidently nonequilibrium conditions, can be expected to stay in a highly-relaxated distribution function given by mono-Maxwellians shifted by the local proton bulk velocity. Here we check the validity of this assumption, calculating on the basis of a Boltzmann-kinetic approach the actually occurring deviations. As we show, especially for low degrees of ionization, $\xi \le 0.3$, both the H-atoms and protons involved do generate in the heliospheric interface clearly pronounced deviations from shifted Maxwellians with asymmetrically shaped distribution functions giving rise to non-convective transport processes and heat conduction flows. Also in the inner heliosheath region and in the heliotail deviations of the proton distribution from the hydrodynamic one must be expected. This sheds new light on the correctness of current calculations of H-atom distribution functions prevailing in the inner heliosphere and also of the Lyman-$\alpha $ absorption features in stellar spectra due to the presence of the hydrogen wall atoms. Deviations from LTE-functions would be even more pronounced in magnetic interfaces, which via CGL-effects cause temperature anisotropies to arise.