ABSTRACT 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 Ly 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 particles in SW and helium ions in LISM), the solar and interstellar magnetic fields, and 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 was developed in such a way that any numerical diffusion or reconnection across the heliopause were not 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.
ABSTRACT
Shock boundary layers are regions bounded by a shock wave on the one side and tangential discontinuity on the other side. These boundary layers are commonly observed in astrophysics. For ...example, they exist in the regions of the interaction of the stellar winds with the surrounding interstellar medium. Additionally, the shock layers (SLs) are often penetrated by the flows of interstellar atoms, as, for instance, in the astrospheres of the stars embedded by the partially ionized interstellar clouds. This paper presents a simple toy model of an SL that aims to qualitatively describe the influence of charge exchange with interstellar hydrogen atoms on the plasma flow in astrospheric SLs. To clearly explore this effect, magnetic fields are neglected, and the geometry is kept as simple as possible. The model explains why the cooling of plasma due to charge exchange in the inner heliosheath leads to an increase in plasma density in front of the heliopause. It also demonstrates that the source of momentum causes changes in the pressure profile within the SL. The paper also discusses the decrease in plasma density near the astropause in the outer SL in the case of layer heating, which is particularly relevant in light of the Voyager measurements in the heliospheric SL.
ABSTRACT
Voyager 1 and 2 crossed the heliopause at ∼122 au in 2012 and ∼119 au in 2018, respectively. It was quite a surprise because the thickness of the inner heliosheath obtained at that time by ...the existing time models of the global heliosphere was significantly larger (by 20–40 au). Until now, the problem of the heliosheath thickness has not been fully resolved. Earlier, in the frame of an oversimplified toy model of nearly isothermal solar wind plasma, it has been shown that the effect of electron thermal conduction may significantly reduce the thickness of the inner heliosheath. In this paper, we present the first results of our 3D kinetic-magnetohydrodynamic (MHD) model of the global heliosphere, where the effect of thermal electron conduction has been considered rigorously. The thermal conduction acts mainly along the magnetic field lines. Classical and saturated thermal fluxes are employed when appropriate. It is shown that the effects of thermal conduction are significant. The thickness of the inner heliospheric is reduced. It is a desired effect since it helps to reconcile the thickness obtained in the model with Voyager data. The other effects are the strong depletion of the heliosheath plasma temperature towards the heliopause and the increase of the plasma temperature in the supersonic solar wind upstream of the termination shock.
ABSTRACT
The main goal of this paper is to explore why observations of many astrospheres (or circumstellar bubbles) show quite stable and smooth structures of astropauses – the tangential ...discontinuities separating the stellar and interstellar winds – while both theory and numerical simulations suggest that tangential discontinuities are unstable due to well known Kelvin–Helmholtz (K-H) instability. It was recognized before that magnetic fields may stabilize the astropauses. In this paper, we explore another mechanism to reduce the K-H instability of the astropauses. This mechanism is a periodic change of the stellar wind dynamic pressure. Fluctuations of the stellar wind parameters are quite expected. For example, the Sun has an 11-yr cycle of global activity although there are also shorter periods of the solar wind fluctuations. We performed the parametric numerical study and demonstrate that the development of the K-H instability depends on the dimensionless parameter χ, which is the ratio of the stellar wind terminal speed and interstellar flow speed. The larger the parameter χ, the larger the fluctuations caused by the K-H instability. It has been shown that the K-H instability is convective which agrees with the previous linear analysis. The stabilization of the astropause by the periodic fluctuations in the stellar wind lead is demonstrated. It is shown that for the solar wind the most effective stabilization occurs when the period of stellar parameter change is about 1–4 yr. For the 11-yr solar cycle, the stabilization effect is weaker.
ABSTRACT
The interstellar boundary explorer (IBEX) has been measuring fluxes of the energetic neutral atoms (ENAs) using the IBEX-Hi (0.3–6 keV) instrument since 2008. We have developed a numerical ...time-depended code to calculate globally distributed flux (GDF) of hydrogen ENAs employing both (1) 3D kinetic-MHD model of the global heliosphere and (2) reconstruction of atom trajectories from 1 au, where they are observed by IBEX, to the point of their origin in the inner heliosheath (IHS). The key factor in the simulation is a detailed kinetic consideration of the pickup ions (PUIs), the supra-thermal component of protons in the heliosphere, which is ‘parental’ to the ENAs and originates in the region of the supersonic solar wind being picked by the heliospheric magnetic field. As a result of our study, we have concluded that (1) the developed model is able to reproduce the geometry of the multilobe structure seen in the IBEX-Hi GDF maps, (2) the GDF is extremely sensitive to the form of the velocity distribution function of PUIs in the IHS, and the accounting for the existence of an additional energetic population of PUIs is essential to explain the data, (3) despite a relatively good agreement, there are some quantitative differences between the model calculations and IBEX-Hi data. Possible reasons for these differences are discussed.
We investigate the signals from neutral helium atoms observed in situ from Earth orbit in 2010 by the Interstellar Boundary Explorer (IBEX). The full helium 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 helium that we call the Warm Breeze. We discuss possible sources for theWarm Breeze, including (1) the secondary population of interstellar helium, created via charge exchange and perhaps elastic scattering of neutral interstellar He atoms on interstellar Hesup + ions in the outer heliosheath, or (2) a gust of interstellar He originating from a hypothetic wave train in the Local Interstellar Cloud. If the second hypothesis is true, the source is likely to be located within a few thousand 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. Whatever the nature of the Warm Breeze, its discovery exposes a critical new feature of our heliospheric environment.
Neutral gas of the local interstellar medium flows through the inner solar system while being deflected by solar gravity and depleted by ionization. The dominating feature in the energetic neutral ...atom Interstellar Boundary Explorer (IBEX) all-sky maps at low energies is the hydrogen, helium, and oxygen interstellar gas flow. The He and O flow peaked around 8 February 2009 in accordance with gravitational deflection, whereas H dominated after 26 March 2009, consistent with approximate balance of gravitational attraction by solar radiation pressure. The flow distributions arrive from a few degrees above the ecliptic plane and show the same temperature for He and O. An asymmetric O distribution in ecliptic latitude points to a secondary component from the outer heliosheath.
ABSTRACT
The hydrogen atoms penetrate the heliosphere from the local interstellar medium, and while being ionized, they form the population of pickup protons. The distribution of pickup protons is ...modified by the adiabatic heating (cooling) induced by the solar wind plasma compression (expansion). In this study, we emphasize the importance of the adiabatic energy change in the inner heliosheath that is usually either neglected or considered improperly. The effect of this process on the energy and spatial distributions of pickup protons and energetic neutral atoms (ENAs), which originate in the charge exchange of pickup protons, has been investigated and quantified using a kinetic model. The model employs the global distributions of plasma and hydrogen atoms in the heliosphere from the simulations of a kinetic-magnetohydrodynamic model of solar wind interaction with the local interstellar medium. The findings indicate that the adiabatic energy change is responsible for the broadening of the pickup proton velocity distribution and the significant enhancement of ENA fluxes (up to ∼5 and ∼20 times in the upwind and downwind directions at energies ∼1–2 keV for an observer at 1 au). It sheds light on the role of adiabatic energy change in explaining the discrepancies between the ENA flux observations and the results of numerical simulations.
Abstract
The heliosphere is the bubble formed by the solar wind as it interacts with the interstellar medium (ISM). The collimation of the heliosheath (HS) flows by the solar magnetic field in the ...heliotail into distinct north and south columns (jets) is seen in recent global simulations of the heliosphere. However, there is disagreement between the models about how far downtail the two-lobe feature persists and whether the ambient ISM penetrates into the region between the two lobes. Magnetohydrodynamic simulations show that these heliospheric jets become unstable as they move down the heliotail and drive large-scale turbulence. However, the mechanism that produces this turbulence had not been identified. Here we show that the driver of the turbulence is the Rayleigh–Taylor (RT) instability produced by the interaction of neutral H atoms streaming from the ISM with the ionized matter in the HS. The drag between the neutral and ionized matter acts as an effective gravity, which causes an RT instability to develop along the axis of the HS magnetic field. A density gradient exists perpendicular to this axis due to the confinement of the solar wind by the solar magnetic field. The characteristic timescale of the instability depends on the neutral H density in the ISM and for typical values the growth rate is ∼3 years. The instability destroys the coherence of the heliospheric jets and magnetic reconnection ensues, allowing ISM material to penetrate the heliospheric tail. Signatures of this instability should be observable in Energetic Neutral Atom maps from future missions such as the Interstellar Mapping and Acceleration Probe (IMAP). The turbulence driven by the instability is macroscopic and potentially has important implications for particle acceleration.
ABSTRACT
In this paper we consider the distribution of interstellar dust in the vicinity of a star under the influence of stellar gravitation and radiation pressure. This study is applicable to stars ...with relatively weak stellar wind and strong radiation, when the stellar radiation has 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 particular dust-grain radius can be calculated analytically based on the classical ‘cold model’. The dust density distribution for a mixture of dust grains with different radii is calculated. We also calculate intensity maps of thermal infrared emission at 24 $\mu$m from dust due to heating by stellar radiation. It is shown that the obtained maps of infrared emission strongly depend on the model parameters: the material of dust grains, the dust-size distribution assumed in the interstellar medium and the approach used to calculate the dust temperature. A bright distinct arc in 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 orders of magnitude larger than for silicates due to more intensive heating of graphite. A possible application of the presented theory is proposed as an algorithm for analysis of the observational images of the infrared arc around the star.