In this Letter, we provide constraints on the direction and magnitude of the pristine (i.e., unperturbed by the interaction with the Sun) local interstellar magnetic field. The constraints are based on analysis of the interstellar magnetic field components at the heliopause measured by magnetometer instruments on board Voyager 1 and 2 spacecraft. The analysis was performed with the help of our kinetic-magnetohydrodynamical (MHD) model of the global heliosphere. The model shows that the solar-induced disturbances of the interstellar magnetic field are extended relatively far from the Sun up to 400-500 AU. The field is draped around the heliopause and compressed. By comparison of the model results with Voyager data we found that the model provides results comparable with the data for the interstellar magnetic field of $B_{LISM}$ = 3.7-3.8 $\mu$G in magnitude and directed towards $\approx$125$^\circ$ in longitude, and $\approx$37$^\circ$ in latitude in the heliographic inertial (HGI) coordinate system.

In this paper we consider the distribution of the interstellar dust in the vicinity of the star under an influence of the stellar gravitation and radiation pressure. This study is applicable to the stars with relatively weak stellar wind and strong radiation, when the stellar radiation swept out the interstellar dust much further from the star than the position of the bow shock created by an interaction of the stellar and interstellar plasma flows. In this case the number density of dust for a certain dust grains radius can be calculated analytically based on the classical `cold model'. The dust density distribution for the mixture of dust grains with different radii is calculated. We also calculated intensity maps of the thermal infrared emission at 24~$\mu$m from dust due to heating by the stellar radiation. It is shown that the obtained maps of the infrared emission extremely depends on the model parameters: material of dust grains, dust size distribution assumed in the interstellar medium, the approach used to calculate the dust temperature. The bright distinct arc at the intensity maps is seen for graphite dust grains and almost disappears for silicates. Absolute values of intensity in the case of graphite are several order of magnitudes larger than for silicates due to more intensive heating of graphite. Possible application of the presented theory is proposed as an algorithm for analysis of the observational images of the infrared arc around the star.

In this work, we present for the first time the Lyman $\alpha$ intensities measured by Voyager 1/UVS in 2003-2014 (at 90-130 AU from the Sun). During this period Voyager 1 measured the Lyman $\alpha$ emission in the outer heliosphere at an almost fixed direction close to the upwind (that is towards the interstellar flow). The data show an unexpected behavior in 2003-2009: the ratio of observed intensity to the solar Lyman $\alpha$ flux is almost constant. Numerical modelling of these data is performed in the frame of a state-of-art self-consistent kinetic-MHD model of the heliospheric interface (Izmodenov & Alexashov, 2015). The model results, for various interstellar parameters, predict a monotonic decrease of intensity not seen in the data. We propose two possible scenarios that explain the data qualitatively. The first is the formation of a dense layer of hydrogen atoms near the heliopause. Such a layer would provide an additional backscattered Doppler shifted Lyman $\alpha$ emission, which is not absorbed inside the heliosphere and may be observed by Voyager. About 35 R of intensity from the layer is needed. The second scenario is an external non-heliospheric Lyman $\alpha$ component, which could be galactic or extragalactic. Our parametric study shows that $\sim$25 R of additional emission leads to a good qualitative agreement between the Voyager 1 data and the model results.

High-resolution mid-infrared observations carried out by the Spitzer Space Telescope allowed one to resolve the fine structure of many astrospheres. In particular, they showed that the astrosphere around the B0.7 Ia star kappa Cas (HD 2905) has a clear-cut arc structure with numerous cirrus-like filaments beyond it. Previously, we suggested a physical mechanism for the formation of such filamentary structures. Namely, we showed theoretically that they might represent the non-monotonic spatial distribution of the interstellar dust in astrospheres (viewed as filaments) caused by interaction of the dust grains with the interstellar magnetic field disturbed in the astrosphere due to colliding of the stellar and interstellar winds. In this paper, we invoke this mechanism to explain the structure of the astrosphere around kappa Cas. We performed 3D magnetohydrodynamic modelling of the astrosphere for realistic parameters of the stellar wind and space velocity. The dust dynamics and the density distribution in the astrosphere were calculated in the framework of a kinetic model. It is found that the model results with the classical MRN size distribution of dust in the interstellar medium do not match the observations, and that the observed filamentary structure of the astrosphere can be reproduced only if the dust is composed mainly of big (micron-sized) grains. Comparison of the model results with observations allowed us to estimate parameters (number density and magnetic field strength) of the surrounding interstellar medium.

High-resolution mid-infrared observations of astrospheres show that many of them have filamentary (cirrus-like) structure. Using numerical models of dust dynamics in astrospheres, we suggest that their filamentary structure might be related to specific spatial distribution of the interstellar dust around the stars, caused by a gyrorotation of charged dust grains in the interstellar magnetic filed. Our numerical model describes the dust dynamics in astrospheres under an influence of the Lorentz force and assumption of a constant dust charge. Calculations are performed for the dust grains with different sizes separately. It is shown that non-monotonic spatial dust distribution (viewed as filaments) appears for dust grains with the period of gyromotion comparable with the characteristic time scale of the dust motion in the astrosphere. Numerical modelling demonstrates that number of filaments depends on charge-to-mass ratio of dust.

Opher et al. (2015), Drake et al. (2015) have shown that the heliospheric magnetic field results in formation of two-jet structure of the solar wind flow in the inner heliosheath, i.e. in the subsonic region between the heliospheric termination shock and the heliopause. In this scenario the heliopause has a tube-like topology as compared with a sheet-like topology in the most models of the global heliosphere (e.g. Izmodenov and Alexashov, 2015). In this paper we explore the two-jet scenario for a simplified astrosphere in which 1) the star is at rest with respect to the circumstellar medium, 2) radial magnetic field is neglected as compared with azimuthal component, 3) the stellar wind outflow is assumed to be hypersonic (both the Mach number and the Alfv\'enic Mach number are much greater than unity at the inflow boundary). We have shown that the problem can be formulated in dimensionless form, in which the solution depends only on one dimensionless parameter epsilon that is reciprocal of the Alfv\'enic Mach number at the inflow boundary. This parameter is proportional to stellar magnetic field. We present the numerical solution of the problem for various values of epsilon. Three first integrals of the governing ideal MHD equations are presented, and we make use of them in order to get the plasma distribution in the jets. Simple relations between distances to the termination shock, astropause and the size of the jet are established. These relations allow us to determine the stellar magnetic field from the geometrical pattern of the jet-like astrosphere.

In this paper, we perform numerical modeling of the interstellar hydrogen fluxes measured by IBEX-Lo during orbit 23 (spring 2009) using a state-of-the-art kinetic model of the interstellar neutral hydrogen distribution in the heliosphere. This model takes into account the temporal and heliolatitudinal variations of the solar parameters as well as non-Maxwellian kinetic properties of the hydrogen distribution due to charge exchange in the heliospheric interface. We found that there is a qualitative difference between the IBEX-Lo data and the modeling results obtained with the three-dimensional, time-dependent model. Namely, the model predicts a larger count rate in energy bin~2 (20-41 eV) than in energy bin~1 (11-21 eV), while the data shows the opposite case. We perform study of the model parameter effects on the IBEX-Lo fluxes and the ratio of fluxes in two energy channels. We shown that the most important parameter, which has a major influence on the ratio of the fluxes in the two energy bins, is the solar radiation pressure. The parameter fitting procedure shows that the best agreement between the model result and the data occurs in the case when the ratio of the solar radiation pressure to the solar gravitation, $\mu_0$, is 1.26$^{+0.06}_{-0.076}$, and the total ionization rate of hydrogen at 1 AU is $\beta_{E,0}=3.7^{+0.39}_{-0.35}\times 10^{-7}$~s$^{-1}$. We have found that the value of $\mu_0$ is much larger than $\mu_0=0.89$, which is the value derived from the integrated solar Lyman-alpha flux data for the period of time studied. We discuss possible reasons for the differences.

This paper provides a detailed description of the latest version of our model of the solar wind (SW) interaction with the local interstellar medium (LISM). This model has already been applied to the analysis of Lyman-alpha absorption spectra toward nearby stars and for analyses of Solar and Heliospheric Observatory/SWAN data. Katushkina et al. (this issue) used the model results to analyze IBEX-Lo data. At the same time, the details of this model have not yet been published. This is a three-dimensional (3D) kinetic-magnetohydrodynamical (MHD) model that takes into account SW and interstellar plasmas (including $\alpha$ particles in SW and helium ions in LISM), the solar and interstellar magnetic fields, and the interstellar hydrogen atoms. The latitudinal dependence of SW and the actual flow direction of the interstellar gas with respect to the Sun are also taken into account in the model. It was very essential that our numerical code had been developed in such a way that any numerical diffusion or reconnection across the heliopause had not been allowed in the model. The heliospheric current sheet is a rotational discontinuity in the ideal MHD and can be treated kinematically. In the paper, we focus in particular on the effects of the heliospheric magnetic field and on the heliolatitudinal dependence of SW.

We show that the electron thermal conductivity may strongly affect the heliosheath plasma flow and the global pattern of the solar wind (SW) interaction with the local interstellar medium (LISM). In particular, it leads to strong reduction of the inner heliosheath thickness that makes possible to explain (qualitatively) why Voyager 1 (V1) has crossed the heliopause at unexpectedly small heliocentric distance of 122 AU. To estimate the effect of thermal conductivity we consider a limiting case when thermal conduction is very effective. To do that we assume the plasma flow in the entire heliosphere is nearly isothermal. Due to this effect, the heliospheric distance of the termination shock has increased by about 15 AU in V1 direction compared to the adiabatic case with gamma = 5/3. The heliospheric distance of the heliopause has decreased by about 27 AU. As a result, the thickness of the inner heliosheath in the model has decreased by about 42 AU and become equal to 32 AU.

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.

Recent analysis of the interstellar helium fluxes measured in 2009-2010 at Earth orbit by the Interstellar Boundary Explorer (IBEX) has suggested that the interstellar velocity (both direction and magnitude) is inconsistent with that derived previously from Ulysses/GAS observations made in the period from 1990 to 2002 at 1.5-5.5 AU from the Sun. Both results are model-dependent and models that were used in the analyses are different. In this paper, we perform an analysis of the Uysses/GAS and IBEX-Lo data using our state-of-the-art 3D time-dependent kinetic model of interstellar atoms in the heliosphere. For the first time we analyze Ulysses/GAS data from year 2007, the closest available Ulysses/GAS observations in time to the IBEX observations. We show that the interstellar velocity derived from the Ulysses 2007 data is consistent with previous Ulysses results and does not agree with the velocity derived from IBEX. This conclusion is very robust since, as is shown in the paper, it does not depend on the ionization rates adopted in theoretical models. We conclude that Ulysses data are not consistent with the new LISM velocity vector from IBEX. In contrast, IBEX data, in principle, could be explained with the LISM velocity vector derived from the Ulysses data. This is possible for the models with the interstellar temperature increased from 6300 K to 9000 K. There is a need to perform further study of possible reasons for the broadening of the helium signal core measured by IBEX. This could be an instrumental effect or due to unconsidered physical processes.

After the termination shock (TS) crossing, the Voyager 2 spacecraft has been observing strong variations of the magnetic field and solar wind parameters in the heliosheath. Anomalous cosmic rays, electrons, and galactic cosmic rays present strong intensity fluctuations. Several works suggested that the fluctuations might be attributed to spatial variations within the heliosheath. Additionally, the variability of the solar wind in this region is caused by different temporal events that occur near the Sun and propagate to the outer heliosphere. To understand the spatial and temporal effects in the heliosheath, it is important to study these effects separately. In this work we explore the role of shocks as one type of temporal effects in the dynamics of the heliosheath. Although currently plasma in the heliosheath is dominated by solar minima conditions, with increasing solar cycle shocks associated with transients will play an important role. We used a 3D MHD multi-fluid model of the interaction between the solar wind and the local interstellar medium to study the propagation of a pair of forward-reverse shocks in the supersonic solar wind, interaction with the TS, and propagation to the heliosheath. We found that in the supersonic solar wind the interaction region between the shocks expands, the shocks weaken and decelerate. The fluctuation amplitudes of the plasma parameters vary with heliocentric distance. The interaction of the pair of shocks with the TS creates a variety of new waves and discontinuities in the heliosheath, which produce a highly variable solar wind flow. The collision of the forward shock with the heliopause causes a reflection of fast magnetosonic waves inside the heliosheath.

We present a two-component gasdynamic model of an interstellar cloud embedded in a hot plasma. It is assumed that the cloud consists of atomic hydrogen gas, interstellar plasma is quasineutral. Hydrogen atoms and plasma protons interact through a charge exchange process. Magnetic felds and radiative processes are ignored in the model. The influence of heat conduction within plasma on the interaction between a cloud and plasma is studied. We consider the extreme case and assume that hot plasma electrons instantly heat the plasma in the interaction region and that plasma flow can be described as isothermal. Using the two-component model of the interaction of cold neutral cloud and hot plasma, we estimate the lifetime of interstellar clouds. We focus on the clouds typical for the cluster of local interstellar clouds embedded in the hot Local Bubble and give an estimate of the lifetime of the Local interstellar cloud where the Sun currently travels. The charge transfer between highly charged plasma ions and neutral atoms generates X-ray emission. We assume typical abundance of heavy ions for the Local Bubble plasma and estimate the X-ray emissivity due to charge exchange from the interface between cold neutral cloud and hot plasma. Our results show that charge exchange X-ray emission from the neutral-plasma interfaces can be a non-negligible fraction of the observed X-ray emission.

We present a new 3d self-consistent two-component (plasma and neutral hydrogen) model of the solar wind interaction with the local interstellar medium (LISM). This model (K-MHD) combines the MHD treatment of the solar wind and the ionized LISM component, with a kinetic model of neutral interstellar hydrogen (LISH). The local interstellar magnetic field (BLISM) intensity and orientation are chosen based on an early analysis of the heliosheath flows (Opher et al. 2009). The properties of the plasma and neutrals obtained using the (K-MHD) model are compared to previous multi-fluid (Opher et al. 2009) and kinetic models (Izmodenov et al. 2005). The new treatment of LISH revealed important changes in the heliospheric properties not captures by the multi-fluid model. These include a decrease in the heliocentric distance to the termination shock (TS), a thinner heliosheath and a reduced deflection angle ({\theta}) of the heliosheath flows. The asymmetry of the termination shock, however, seems to be unchanged by the kinetic aspect of the LISH.

We explore imprints of the solar wind interaction with the local interstellar medium on spectral properties of the backscattered solar Ly-alpha radiation. We employ newly developed effective model for the interstellar H atom velocity distribution function in the heliosphere (Katushkina and Izmodenov 2010, 2011). This model takes into account both global effects of H atom perturbations at the heliospheric boundaries and local (i.e. within 10-20 AU from the Sun) effects of the solar ionization, charge exchange, solar gravitation and radiation pressure. Backscattered solar Ly-alpha profiles and their zero, first and second moments were computed at 1 AU for the anti-solar directions of line-of-sight. Then the moments are compared with the moments calculated by using more simplified one-component and two-component hot models. The comparison shows that the Ly-alpha intensities calculated in the frame of new model are somewhat smaller than those calculated by means of the simplified models. Small differences in the first moment (i.e. line-shift) are also observed in the downwind direction. The largest difference between new model and hot models is seen in the second moment of the backscattered Ly-alpha profile. This moment is also called as line-of-sight temperature or line-width of the backscattered spectra. The spectral width depends on the width of the H atom velocity distribution function within few AU from the Sun. We show that the line-widths calculated on the basis of the our model is qualitatively different from those of simplified models.

This paper deals with the modeling of the interstellar hydrogen atoms (H atoms) distribution in the heliosphere. We study influence of the heliospheric interface, that is the region of the interaction between solar wind and local interstellar medium, on the distribution of the hydrogen atoms in vicinity of the Sun. The distribution of H atoms obtained in the frame of the self-consistent kinetic-gasdynamic model of the heliospheric interface is compared with a simplified model which assumes Maxwellian distribution of H atoms at the termination shock and is called often as 'hot' model. This comparison shows that the distribution of H atoms is significantly affected by the heliospheric interface not only at large heliocentric distances, but also in vicinity of the Sun at 1-5 AU. Hence, for analysis of experimental data connected with direct or undirect measurements of the interstellar atoms one necessarily needs to take into account effects of the heliospheric interface. In this paper we propose a new model that is relatively simple but takes into account all major effects of the heliospheric interface. This model can be applied for analysis of backscattered Ly-alpha radiation data obtained on board of different spacecraft.

We present new kinetic-gasdynamic model of the solar wind interaction with the local interstellar medium. The model incorporates several processes suggested by McComas et al. (2009) for the origin of the heliospheric ENA ribbon -- the most prominent feature seen in the all sky maps of heliospheric ENAs discovered by the Interstellar Boundary Explorer (IBEX). The ribbon is a region of enhanced fluxes of ENAs crossing almost the entire sky. Soon after the ribbon's discovery it was realized (McComas et al., 2009) that the enhancement of the fluxes could be in the directions where the radial component of the interstellar magnetic field around the heliopause is close to zero (Schwadron et al., 2009). Our model includes secondary charge exchange of the interstellar H atoms with the interstellar pickup protons outside the heliopause and is a further advancement of the kinetic-gasdynamic model by Malama et al. (2006) where pickup protons were treated as a separate kinetic component. Izmodenov et al. (2009) have shown in the frame of Malama's model that the interstellar pickup protons outside the heliopause maybe a significant source of ENAs at energies above 1 keV. The difference between the current work and that of Izmodenov et al. (2009) is in the assumption of no-scattering for newly created pickup protons outside the heliopause. In this limit the model produces a feature qualitatively similar to the ribbon observed by IBEX.

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.

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.

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

We use newly developed 3D kinetic MHD models of the heliosphere to predict heliospheric H I Lyman-alpha absorption for various lines of sight. These predictions are compared with actual Lyman-alpha spectra from the Hubble Space Telescope, many of which have yielded previous detections of heliospheric absorption. We find that the absorption predicted by the models is weakly affected by both the magnitude and orientation of the assumed ISM magnetic field. Models with B=1.25-2.5 micro-Gauss and an angle of alpha=15-45 deg with respect to the upwind direction of the ISM flow generally provide the best fits to the data, but the sensitivity of the Lyman-alpha absorption to many model input parameters makes it difficult to fully characterize the region of parameter space allowed by the data. We also use the models to assess the degree to which heliospheric asymmetries induced by the ISM field should be apparent in Lyman-alpha absorption. An ISM field that is skewed with respect to the ISM flow vector results in substantial azimuthal asymmetries in both the hydrogen wall and heliosheath, but only the heliosheath asymmetries yield potentially detectable asymmetries in Lyman-alpha absorption; and then only in downwind directions, where comparison with the data is complicated by few actual absorption detections and an insufficient model grid extent.

We assess what information HST observations of stellar Ly-alpha lines can provide on the heliosheath, the region of the heliosphere between the termination shock and heliopause. To search for evidence of heliosheath absorption, we conduct a systematic inspection of stellar Ly-alpha lines reconstructed after correcting for ISM absorption (and heliospheric/astrospheric absorption, if present). Most of the stellar lines are well centered on the stellar radial velocity, as expected, but the three lines of sight with the most downwind orientations relative to the ISM flow (Chi1 Ori, HD 28205, and HD 28568) have significantly blueshifted Ly-alpha lines. Since it is in downwind directions where heliosheath absorption should be strongest, the blueshifts are almost certainly caused by previously undetected heliosheath absorption. We make an initial comparison between the heliosheath absorption and the predictions of a pair of heliospheric models. A model with a complex multi-component treatment of plasma within the heliosphere predicts less absorption than a model with a simple single-fluid treatment, which leads to better agreement with the data. Finally, we find that nonplanetary energetic neutral atom (ENA) fluxes measured by the ASPERA-3 instrument on board Mars Express, which have been interpreted as being from the heliosheath, are probably too high to be consistent with the relative lack of heliosheath absorption seen by HST. This would argue for a local interplanetary source for these ENAs instead of a heliosheath source.

We study the EUV/soft X-ray emission generated by charge transfer between solar wind heavy ions and interstellar H and He neutral atoms in the inner Heliosphere. We present heliospheric maps and spectra for stationary solar wind, depending on solar cycle phase, solar wind anisotropies and composition, line of sight direction and observer position. A time-dependant simulation of the X-ray intensity variations due to temporary solar wind enhancement is compared to XMM Newton recorded data of the Hubble Deep Field North observation (Snowden et al. 2004). Results show that the heliospheric component can explain a large fraction of the line intensity below 1.3 keV, strongly attenuating the need for soft X-ray emission from the Local Interstellar Bubble.

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

We present a new model of the heliospheric interface - the region of the solar wind interaction with the local interstellar medium. This new model performs a multi-component treatment of charged particles in the heliosphere. All charged particles are divided into several co-moving types. The coldest type, with parameters typical of original solar wind protons, is considered in the framework of fluid approximation. The hot pickup proton components created from interstellar H atoms and heliospheric ENAs by charge exchange, electron impact ionization and photoionization are treated kinetically. The charged components are considered self-consistently with interstellar H atoms, which are described kinetically as well. To solve the kinetic equation for H atoms we use the Monte Carlo method with splitting of trajectories, which allows us 1) to reduce statistical uncertainties allowing correct interpretation of observational data, 2) to separate all H atoms in the heliosphere into several populations depending on the place of their birth and on the type of parent protons.

The Sun is moving through a warm ($\sim$6500 K) and partly ionized local interstellar cloud (LIC) with a velocity of $\sim$26 km/s. Recent measurements of the ionization of the LIC (Wolff et al., 1999) suggest that interstellar helium in the vicinity of the Sun is 30-40 % ionized, while interstellar hydrogen is less ionized. Consequently, interstellar helium ions contribute up to 50% of the total dynamic pressure of the ionized interstellar component. Up to now interstellar helium ions have been ignored in existing models of the heliospheric interface. In this paper we present results of a new model of the solar wind interaction with the interstellar medium, which takes into account interstellar helium ions. Using results of this model we find that the heliopause, termination and bow shocks are closer to the Sun when compared to the model results that ignore $He$ ions. The influence of interstellar helium ions is partially compensated by solar wind alpha particles, which are taken into account in our new model as well. Finally, we apply our new model to place constraints on the plausible location of the termination shock.

Physical processes in the tail of the solar wind interaction region with the partially ionized local interstellar medium are investigated in a framework of the self-consistent kinetic-gas dynamic model. It is shown that the charge exchange process of the hydrogen atoms with the plasma protons results in suppression of the gas dynamic instabilities and disappearance the contact discontinuity at sufficiently (~3000 AU) large distances from the Sun. The solar wind plasma temperature decreases and, ultimately, the parameters of the plasma and hydrogen atoms approach to the corresponding parameters of the unperturbed interstellar medium at large heliocentric distances.

The Very Local Interstellar Medium (VLISM) neutral hydrogen and proton densities are still not precisely known even after three decades of deep space research and the existence of the EUV and other diagnostic data obtained by Pioneer 10/11, Voyager 1/2 and other spacecraft. The EUV data interpretation, in particular, has suffered because of inadequate neutral hydrogen-plasma models, difficulty of calculating the multiply scattered Lyman alpha glow and calibration uncertainties. Recently, all these difficulties have been significantly reduced. In the present work we have used the latest state of the art supersonic VLISM neutral hydrogen-plasma and Monte Carlo radiative transfer model, incorporating neutral density, temperature, and velocity variations, actual solar line shape, realistic redistribution function, Doppler and aberration effects. This work presents the methodology of the radiative transfer code and the first results of the comparison of the model predictions with the Pioneer 10 data. Monte Carlo radiative transfer calculations were carried out for five neutral hydrogen- plasma models and compared with Pioneer data. The first results are quite encouraging. We found that the VLISM ionization ratio is between 0.2 and 0.5 and that the VLISM neutral hydrogen density is less than 0.25 cm^{-3}. The present calculation suggests that the Pioneer 10 photometer derived intensities (Rayleighs) need to be increased by a factor of 2. If this model- derived calibration is used then the difference between Pioneer 10 and Voyager 2 intensity values is reduced to about 2.2. The model, neutral hydrogen density = 0.15 cm${-3} and proton density = 0.07 cm^{-3}, is found to best fit the Pioneer 10 data.

This paper reviews the theoretical approaches and existing models of the solar wind interaction with the Local Interstellar Cloud (LIC). Models discussed take into account the multi-component nature of the solar wind and local interstellar medium. Basic results of the modeling and their possible applications to interpretation of space experiments are summarized. Open questions of global modeling of the solar wind/LIC interaction and future perspectives are discussed.

In this paper, we study heliospheric Ly-$\alpha$ absorption toward nearby stars in different lines of sight. We use the Baranov-Malama model of the solar wind interaction with a two-component (charged component and H atoms) interstellar medium. Interstellar atoms are described kinetically in the model. The code allows us to separate the heliospheric absorption into two components, produced by H atoms originating in the hydrogen wall and heliosheath regions, respectively. We study the sensitivity of the heliospheric absorption to the assumed interstellar proton and H atom number densities. These theoretical results are compared with interstellar absorption toward six nearby stars observed by the Hubble Space Telescope.

We use an updated self-consistent kinetic/gasdynamic model of the solar wind/interstellar flow interaction to compute the Lyman-alpha absorption by hydrogen atoms of both interstellar and solar origin generated in the heliospheric interface. We apply this model to the direction of the star Sirius. We show that the neutralized, compressed solar wind from the heliosheath explains the extra absorption at positive redshifts observed in the Lyman-alpha line with the HST/GHRS. This removes the need for the interstellar hot gas previously proposed to explain this extra absorption. We also show that extra absorption in the blue wing can possibly be explained by H atoms formed in an astrosphere around Sirius, providing the stellar wind is at least as massive and fast as the solar wind.

The solar system is moving through the partially ionized local interstellar cloud (LIC). The ionized matter of the LIC interacts with the expanding solar wind forming the heliospheric interface. The neutral component (interstellar atoms) penetrates through the heliospheric interface into the heliosphere, where it is measured directly "in situ" as pick-up ions and neutral atoms (and as anomalous cosmic rays) or indirectly through resonant scattering of solar Ly-alpha. When crossing the heliospheric interface, interstellar atoms interact with the plasma component through charge exchange. This interaction leads to changes of both neutral gas and plasma properties. The heliospheric interface is also the source of radio emissions which have been detected by the Voyager since 1983. In this paper, we have used a kinetic model of the flow of the interstellar atoms with updated values of velocity, temperature, and density of the circumsolar interstellar hydrogen and calculated how all quantities which are directly associated to the observations vary as a function of the interstellar proton number density. These quantities are the degree of filtration, the temperature and the velocity of the interstellar H atoms in the inner heliosphere, the distances to the bow shock (BS), heliopause, and termination shock, and the plasma frequencies in the LIC, at the BS and in the maximum compression region around the heliosphere which constitutes the "barrier" for radio waves formed in the interstellar medium. Comparing the model results with recent pickup ion data, Ly- alpha measurements, and low-frequencies radio emissions, we have searched for a number density of protons in the local interstellar cloud compatible with all observations.