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.

We report the discovery of filamentary structures stretched behind the bow-shock-producing high-mass X-ray binary Vela X-1 using the SuperCOSMOS H-alpha Survey and present the results of optical spectroscopy of the bow shock carried out with the Southern African Large Telescope (SALT). The geometry of the detected structures suggests that Vela X-1 has encountered a wedge-like layer of enhanced density on its way and that the shocked material of the layer partially outlines a wake downstream of Vela X-1. To substantiate this suggestion, we carried out 3D magnetohydrodynamic simulations of interaction between Vela X-1 and the layer for three limiting cases. Namely, we run simulations in which (i) the stellar wind and the interstellar medium (ISM) were treated as pure hydrodynamic flows, (ii) a homogeneous magnetic field was added to the ISM, while the stellar wind was assumed to be unmagnetized, and (iii) the stellar wind was assumed to possess a helical magnetic field, while there was no magnetic field in the ISM. We found that although the first two simulations can provide a rough agreement with the observations, only the third one allowed us to reproduce not only the wake behind Vela X-1, but also the general geometry of the bow shock ahead of it.

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.

Aims. A recent debate on the decade-long stability of the interstellar He flow vector, and in particular the flow longitude, has prompted us to check for any variability in the interstellar H flow vector as observed by the SWAN instrument on board SOHO. Methods. We used a simple model-independent method to determine the interstellar H flow longitude, based on the parallax effects induced on the Lyman-{\alpha} intensity measured by SWAN following the satellite motion around the Sun. Results. Our results show that the interstellar H flow vector longitude does not vary significantly from an average value of 252.9$^{\circ}$ $\pm$ 1.4$^{\circ}$ throughout the 20-year span of the SWAN dataset, further strengthening the arguments for the stability of the interstellar gas flow.

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.

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.

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.