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.

We report the results of our study of the high-velocity (\approx160 km/s) runaway O star CPD-64 2731 and its associated horseshoe-shaped nebula discovered with the Wide-field Infrared Survey Explorer. Spectroscopic observations with the Southern African Large Telescope and spectral analysis indicate that CPD-64 2731 is a fast-rotating main-sequence O5.5 star with enhanced surface nitrogen abundance. We derive a projected rotational velocity of \approx300 km/s which is extremely high for this spectral type. Its kinematic age of \approx6 Myr, assuming it was born near the Galactic plane, exceeds its age derived from single star models by a factor of two. These properties suggest that CPD-64 2731 is a rejuvenated and spun-up binary product. The geometry of the nebula and the almost central location of the star within it argue against a pure bow shock interpretation for the nebula. Instead, we suggest that the binary interaction happened recently, thereby creating the nebula, with a cavity blown by the current fast stellar wind. This inference is supported by our results of 2D numerical hydrodynamic modelling.

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.

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.

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.

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.

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.

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

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.

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.