Interstellar Hydrogen Atoms Research Articles

Effects of Heliolatitudinal Anisotropy of Solar Far-ultraviolet/Extreme-ultraviolet Emissions on Lyα Helioglow

We present a study of the influence of solar UV anisotropy on the heliospheric backscatter helioglow generated by resonant scattering of solar Lyα photons on interstellar hydrogen atoms around the Sun. Simulations based on the WawHelioGlow model suggest that the response of the helioglow pole-to-ecliptic ratio to the anisotropy is linear, but 15% of the anisotropy (polar darkening) generates 30%–40% change in the ratio in the solar minimum and 15%–20% in the solar maximum. We attribute this difference to an interplay between the solar UV anisotropy and the latitudinal structure of the solar wind in solar minima. The solar UV anisotropy also increases the helioglow intensity from the downwind direction by 5%–10%, due to the influence of the anisotropy on the ionization losses and trajectories of atoms passing by the Sun in polar regions. Consequently, midlatitude regions (in the heliographic and ecliptic coordinates) are least affected by the UV anisotropy. By comparison of the simulation results with observations of the Solar and Heliospheric Observatory/SWAN satellite instrument, we derive the day-by-day time evolution of the solar Lyα anisotropy for the north and south poles over two solar cycles from 1996 to 2022. The inferred anisotropy is ∼5%–10% in solar minima and ∼15%–25% in solar maxima, the northern anisotropy being stronger than the southern one. Our study suggests that in solar minima a highly structured solar wind is associated with relatively small solar UV anisotropy, while in solar maxima the solar wind is more isotropic but a substantial solar UV anisotropy appears.

Effects of charge exchange on plasma flow in the heliosheath and astrosheathes

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.

Absorption of the Lyα Radiation in the Heliosphere

Abstract Absorption of the Lyα radiation on interstellar neutral hydrogen (ISN H) atoms in the heliosphere is a potentially important effect to account for in precise gas distribution simulations. In this paper, we develop a method to estimate the magnitude of absorption of solar Lyα radiation inside the solar wind termination shock and to include absorption effects in the Warsaw Test Particle Model (WTPM) by an appropriate modification of radiation pressure. We perform calculations of absorption effects on a 3D grid in the heliosphere and present a set of parameters to model absorption effects for the mean solar activity conditions. We show that absorption can change by up to 3%, depending on the solar activity level. Using a modified version of WTPM, we calculate the expected signal from IBEX-Lo and show that absorption may modify the simulated flux by up to 8%.

New Horizons Detection of the Local Galactic Lyman-α Background

Abstract Since 2007 the Alice spectrograph on the New Horizons (NH) spacecraft has been used to periodically observe the Lyman-α (Lyα) emissions of the interplanetary medium (IPM), which mostly result from resonant scattering of solar Lyα emissions by interstellar hydrogen atoms passing through the solar system. Three observations of IPM Lyα along a single great circle were made during the NH cruise to Pluto, and these have been supplemented by observations along six great circles (spread over the sky at 30° intervals), acquired one month before and one day after the NH flyby of Pluto, and on a further five occasions since then, out to just over 47 au from the Sun. These data indicate a distant Lyα background of 43 ± 3 Rayleigh brightness (equivalent to 56 ± 4 nW m−2 sr−1), which is present in all directions (i.e., not only in the upstream direction, as previously reported). This result is found independently by: (1) the falloff with distance from the Sun of the IPM Lyα brightness observed by NH–Alice in several directions on the sky, and (2) the residual between the observed brightness and a model brightness accounting for the resonantly scattered solar Lyα component alone. The repeated observations show that this distant Lyα background is constant and uniform over the sky, and represents the local Galactic Lyα background. The observations show no strong correlation with the cloud structure of the local IPM. The observed brightness constrains the absorption coefficient of interstellar dust at Lyα to 0.2 ± 0.01 kpc−1.

Imprints of the secondary interstellar hydrogen atoms at 1 au

ABSTRACT In this paper, we search for the possible imprints of the secondary interstellar hydrogen atoms created at the heliospheric boundary in the full-sky maps of the hydrogen fluxes at the Earth orbit. By using our three-dimensional time-dependent kinetic model, the maps of the hydrogen fluxes are calculated for different phases of the solar cycle and different energy ranges. It is shown that the flux maps obtained during the solar minimum conditions for the energy range 1–20 eV have specific features such as blobs and tails, which are pronounced due to a signal of the secondary component of the interstellar atoms. We investigate how these features depend on parameters of the secondary atoms far away from the Sun and found that the geometry and shape of the tails depend on the averaged velocity and kinetic temperatures of the secondary population. The results of the paper provide a strategy where and when we need to look in order to detect the secondary component of the interstellar hydrogen separately from the primary component at 1 au. This can be important for future space missions devoted to the exploration of the heliospheric boundary.

Slowdown and Heating of Interstellar Neutral Helium by Elastic Collisions Beyond the Heliopause.

Direct sampling of interstellar neutral (ISN) atoms close to the Sun enables studies of the very local interstellar medium (VLISM) around the heliosphere. The primary population of ISN helium atoms has, until now, been assumed to reflect the pristine VLISM conditions at the heliopause. Consequently, the atoms observed at 1 au by the Interstellar Boundary Explorer (IBEX) were used to determine the VLISM temperature and velocity relative to the Sun, without accounting for elastic collisions with other species outside the heliopause. Here, we evaluate the effect of these collisions on the primary ISN helium population. We follow trajectories of helium atoms and track their collisions with slowed plasma and interstellar hydrogen atoms ahead of the heliopause. Atoms typically collide a few times in the outer heliosheath, and only ~1.5% of the atoms are not scattered at all. We use calculated differential cross sections to randomly choose scattering angles in these collisions. We estimate that the resulting primary ISN helium atoms at the heliopause are slowed down by ~0.45 km s-1 and heated by ~1100 K compared to the pristine VLISM. The resulting velocity distribution is asymmetric and shows an extended tail in the antisunward direction. Accounting for this change in the parameters derived from IBEX observations gives the Sun's relative speed of 25.85 km s-1 and temperature of 6400 K in the pristine VLISM. Finally, this paper serves as a source of the differential cross sections for elastic collisions with helium atoms.

The multiphase and magnetized neutral hydrogen seen by LOFAR

Faraday tomography of polarimetric observations at low frequency in the radio is a unique tool for studying the structure of the magneto-ionic diffuse interstellar medium (ISM) based on Faraday depth. LOFAR data below 200 MHz have revealed a plethora of features in polarization, whose origin remains unknown. Previous studies have highlighted the remarkable association of such features with tracers of the magnetized-neutral ISM, such as interstellar dust and atomic hydrogen (HI). However, the physical conditions responsible for the correlation between magneto-ionic and neutral media have not been clarified yet. In this Letter we further investigate the correlation between LOFAR data and the HI spectroscopic observations at 21 cm from the Effelsberg-Bonn HI Survey (EBHIS). We focus on the multiphase properties of the HI gas. We present the first statistical study on the morphological correlation between LOFAR tomographic data and the cold (CNM), lukewarm (LNM), and warm (WNM) neutral medium HI phases. We use the Regularized Optimization for Hyper-Spectral Analysis approach to decompose the HI phases based on a Gaussian decomposition of the HI spectra. We study four fields of view – Fields 3C196, A, B, and C – and find, in at least the first two, a significant correlation between the LOFAR and EBHIS data using the histograms of oriented gradients (HOG) feature. The absence of a correlation in Fields B and C is caused by a low signal-to-noise ratio in polarization. The observed HOG correlation in Fields 3C196 and A is associated with all HI phases and it is surprisingly dominant in the CNM and LNM phases. We discuss possible mechanisms that would explain the correlation between CNM, LNM, and WNM with polarized emission at Faraday depths up to 10 rad m−2. Our results show how the complex structure of the ionic medium seen by the LOFAR tomographic data is tightly related to phase transition in the diffuse and magnetized neutral ISM traced by HI spectroscopic data.

Density of Neutral Hydrogen in the Sun's Interstellar Neighborhood

Abstract Interstellar neutral atoms, unlike charged particles, freely penetrate the heliosphere, allowing us to sample the physical state of the interstellar matter directly. Most interstellar hydrogen atoms are ionized before reaching the inner heliosphere and become energetic protons picked up by the solar wind and transported away from the Sun. Consequently, observations of interstellar hydrogen atoms by missions operating within a few astronomical units from the Sun are subject to significant systematic uncertainties. We analyze observations from the Solar Wind Around Pluto instrument on New Horizons, the first experiment to provide extensive measurements of the picked-up protons far from the Sun. Analyzing the density of these protons, we find an interstellar neutral hydrogen density at the termination shock of 0.127 ± 0.015 cm−3, i.e., ∼40% higher than previously thought. We show that the Voyager observations of the slowdown of the solar wind further support this value. This result resolves a problem of why energetic neutral atom fluxes, created from pickup ions by charge exchange with hydrogen atoms, are roughly twice that expected from numerical models. Our result also implies higher charge exchange rates at the heliospheric boundaries and, consequently, a less asymmetric shape of the heliosphere. Based on a previous study of the atom filtration in the heliospheric boundaries, we estimate the neutral hydrogen density in the unperturbed local interstellar medium of 0.195 ± 0.033 cm−3. This value agrees with astrophysical observations of the interstellar clouds in the Sun proximity.

Angular Scattering in Charge Exchange: Issues and Implications for Secondary Interstellar Hydrogen

Abstract Interstellar neutral atoms provide a remote diagnostic of plasma in the outer heliosheath and the very local interstellar medium via charge exchange collisions that convert ions into atoms and vice versa. So far, most studies of interstellar atoms assumed that daughter hydrogen atoms directly inherit the kinetic properties of parent protons. This assumption neglects angular scattering of the interacting particles. However, for low relative velocities, as expected for charge exchanges in the outer heliosheath, this scattering is significant. In this study, we present how the parameters of daughter populations depend on the relative velocity and temperatures of parent populations. For this purpose, we numerically compute collision terms with and without this scattering. We find that the secondary population of interstellar hydrogen atoms, for the parent populations with the relative bulk velocity of 20 km s−1 and equal temperatures of 7500 K, has ∼2 km s−1 higher bulk velocity if the scattering is taken into account. Additionally, temperatures are higher by ∼2400 K and ∼1200 K in parallel and perpendicular directions to the relative motion of parent populations, respectively. Moreover, a significant departure of secondary atoms from the Maxwell–Boltzmann distribution is expected for high relative velocities of parent populations. This process affects the distribution and density of interstellar atoms in the heliosphere and production of pickup ions. Thus, we show that angular scattering in charge exchange collisions is important to include in analyses of interstellar neutral atoms and pickup ions observed at 1 au and in the outer heliosphere.

My C.G.S.I.S.A.H. Theory of Dark Matter

Theory and observations concerning the cosmic reionization epoch are briefly discussed in the context of recent observations attributed to dark matter. A case is made that cold ground state interstellar atomic hydrogen of average density of about one atom per cubic centimeter (1.67 × 10-21 kg·m-3 or 1.67 × 10-24 g·cm-3) appears to be the most likely candidate to explain these observations.

Synthetic Observations of 21 cm H i Line Profiles from Inhomogeneous Turbulent Interstellar H i Gas with Magnetic Fields

Abstract We carried out synthetic observations of interstellar atomic hydrogen at 21 cm wavelength by utilizing the magnetohydrodynamic numerical simulations of the inhomogeneous turbulent interstellar medium. The cold neutral medium (CNM) shows a significantly clumpy distribution with a small volume filling factor of 3.5%, whereas the warm neutral medium (WNM) has a distinctly different and smooth distribution with a large filling factor of 96.5%. In projection on the sky, the CNM exhibits a highly filamentary distribution with a subparsec width, whereas the WNM shows a smooth, extended distribution. In the H i optical depth, the CNM is dominant and the contribution of the WNM is negligibly small. The CNM has an area covering factor of 30% in projection, while the WNM has a covering factor of 70%. This means that the emission–absorption measurements toward radio continuum compact sources tend to sample the WNM with a probability of 70%, yielding a smaller H i optical depth and a smaller H i column density than those of the bulk H i gas. The emission–absorption measurements, which are significantly affected by the small-scale large fluctuations of the CNM properties, are not suitable for characterizing the bulk H i gas. Larger-beam emission measurements that are able to fully sample the H i gas will provide a better tool for that purpose, if a reliable proxy for hydrogen column density, possibly dust optical depth and gamma rays, is available. The present results provide a step toward precise measurements of the interstellar hydrogen with ∼10% accuracy. This will be crucial in interstellar physics, including identification of the proton–proton interaction in gamma-ray supernova remnants.

The role of charge-exchange cross-section for pickup protons and neutrals in the inner heliosheath

The process of deceleration of the solar wind downstream of the termination shock is studied on the basis of a one-dimensional multi-component model. It is assumed that the solar wind consists of thermal protons, electrons and interstellar pickup protons. The protons interact with interstellar hydrogen atoms by charge-exchange. Two cases are considered. In the first one, the charge-exchange cross-section for thermal protons and hydrogen atoms is the same as for pickup protons and atoms. Under this condition, there is a strong dependence of the solar wind velocity on the downstream temperature of pickup protons. When the proton temperature is close to 10 keV, the change in the velocity with the distance from the termination shock is similar to that measured on the Voyager 1 spacecraft: linear velocity decrease is accompanied by an extended transition region with near-zero velocity. However, with a more careful approach to the choice of the charge-exchange cross-section, the situation changes dramatically. The strong dependence of the solar wind speed on the pickup proton temperature disappears and the transition region in the heliosheath disappears as well, at least at reasonable distances from the TS.

Prescriptions on antiproton cross section data for precise theoretical antiproton flux predictions

After the breakthrough from the satellite-borne PAMELA detector, the flux of cosmic-ray (CR) antiprotons has been provided with unprecedented accuracy by AMS-02 on the International Space Station. Its data spans an energy range from below 1 GeV up to 400 GeV and most of the data points contain errors below the amazing level of 5%. The bulk of the antiproton flux is expected to be produced by the scatterings of CR protons and helium off interstellar hydrogen and helium atoms at rest. The modeling of these interactions, which requires the relevant production cross sections, induces an uncertainty in the determination of the antiproton source term that can even exceed the uncertainties in the CR $\bar{p}$ data itself. The aim of the present analysis is to determine the uncertainty required for $p+p\rightarrow \bar{p} + X$ cross section measurements such that the induced uncertainties on the $\bar{p}$ flux are at the same level. Our results are discussed both in the center-of-mass reference frame, suitable for collider experiments, and in the laboratory frame, as occurring in the Galaxy. We find that cross section data should be collected with accuracy better that few percent with proton beams from 10 GeV to 6 TeV and a pseudorapidity $\eta$ ranging from 2 to almost 8 or, alternatively, with $p_T$ from 0.04 to 2 GeV and $x_R$ from 0.02 to 0.7. Similar considerations hold for the $p$He production channel. The present collection of data is far from these requirements. Nevertheless, they could, in principle, be reached by fixed target experiments with beam energies in the reach of CERN accelerators.

MODELING THE SOLAR WIND AT THE ULYSSES, VOYAGER, AND NEW HORIZONS SPACECRAFT

ABSTRACT The outer heliosphere is a dynamic region shaped largely by the interaction between the solar wind and the interstellar medium. While interplanetary magnetic field and plasma observations by the Voyager spacecraft have significantly improved our understanding of this vast region, modeling the outer heliosphere still remains a challenge. We simulate the three-dimensional, time-dependent solar wind flow from 1 to 80 astronomical units (au), where the solar wind is assumed to be supersonic, using a two-fluid model in which protons and interstellar neutral hydrogen atoms are treated as separate fluids. We use 1 day averages of the solar wind parameters from the OMNI data set as inner boundary conditions to reproduce time-dependent effects in a simplified manner which involves interpolation in both space and time. Our model generally agrees with Ulysses data in the inner heliosphere and Voyager data in the outer heliosphere. Ultimately, we present the model solar wind parameters extracted along the trajectory of the New Horizons spacecraft. We compare our results with in situ plasma data taken between 11 and 33 au and at the closest approach to Pluto on 2015 July 14.

Ionization of Interstellar Hydrogen Beyond the Termination Shock

Models of solar wind interaction with the surrounding interstellar medium usually disregard ionization of interstellar hydrogen atoms beyond the solar wind termination shock. If and when included, the effects of ionization in the heliospheric interface region are often obscured by complexities of the interaction. This work assesses the importance of interstellar hydrogen ionization in the heliosheath. Photoionization could be accounted for in a straightforward way. In contrast, electron impact ionization is largely unknown because of poorly understood energy transfer to electrons at the termination shock and beyond. We first estimate the effect of photoionization and then use it as a yardstick to assess the role of electron impact ionization. The physical estimates show that ionization of interstellar hydrogen may lead to significant mass loading in the inner heliosheath which would slow down plasma flowing toward the heliotail and deplete populations of nonthermal protons, with the corresponding effect on heliospheric fluxes of energetic neutral atoms.

On the evolution of the κ distribution of protons in the inner heliosheath

AbstractThe evolution of the solar wind proton distribution function along the plasma flow downstream of the heliospheric termination shock is studied. Starting from a kinetic phase space transport equation valid in the bulk frame of the plasma flow that takes into account convective changes, cooling, velocity diffusion, and charge exchange‐induced injection and loss, the associated moment equation for the evolution of the pressure of the total proton population consisting of thermal solar wind and suprathermal pickup protons is derived. Assuming that the local joint proton distribution can always be represented by a so‐called κ distribution, an ordinary differential equation for the variation of the parameter κ along a given streamline is obtained. This way, the proton velocity distribution can be computed for the whole inner heliosheath. It is demonstrated that the accelerating effect of velocity diffusion is to be expected to overcompensate the loss effect on the proton distribution at higher velocities due to charge exchange with cold interstellar hydrogen atoms. While this corroborates a value of κ < 1.65 in the inner (upwind) heliosheath, at the same time it reveals that the assumption of a constant κ, which was commonly made in earlier studies, should be abandoned.

Mass loading of bow shock pulsar wind nebulae

We investigate the dynamics of bow shock nebulae created by pulsars moving supersonically through a partially ionized interstellar medium. A fraction of interstellar neutral hydrogen atoms penetrating into the tail region of a pulsar wind will undergo photo-ionization due to the UV light emitted by the nebula, with the resulting mass loading dramatically changing the flow dynamics of the light leptonic pulsar wind. Using a quasi 1-D hydrodynamic model of both non-relativistic and relativistic flow, and focusing on scales much larger than the stand-off distance, we find that if a relatively small density of neutral hydrogen, as low as $10^{-4}\,\text{cm}^{-3}$, penetrate inside the pulsar wind, this is sufficient to strongly affect the tail flow. Mass loading leads to the fast expansion of the pulsar wind tail, making the tail flow intrinsically non-stationary. The shapes predicted for the bow shock nebulae compare well with observations, both in H$\alpha$ and X-rays.

THREE-DIMENSIONAL KINETIC-MHD MODEL OF THE GLOBAL HELIOSPHERE WITH THE HELIOPAUSE-SURFACE FITTING

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.

Lyman-α observations of astrospheres

Abstract. Charge-exchange reactions between outflowing stellar wind protons and interstellar neutral hydrogen atoms entering a stellar astrosphere produce a region of piled-up-decelerated neutral hydrogen called the hydrogen wall. Absorption by this gas, which is observed in stellar Lyman-α emission lines, provides the only viable technique at this time for measuring the mass-loss rates of F–M dwarf stars. We describe this technique, present an alternative way for understanding the relation of mass-loss rate with X-ray emission, and identify several critical issues.

THREE-FLUID, THREE-DIMENSIONAL MAGNETOHYDRODYNAMIC SOLAR WIND MODEL WITH EDDY VISCOSITY AND TURBULENT RESISTIVITY

We have developed a three-fluid, three-dimensional magnetohydrodynamic solar wind model that incorporates turbulence transport, eddy viscosity, turbulent resistivity, and turbulent heating. The solar wind plasma is described as a system of co-moving solar wind protons, electrons, and interstellar pickup protons, with separate energy equations for each species. Numerical steady-state solutions of Reynolds-averaged solar wind equations coupled with turbulence transport equations for turbulence energy, cross helicity, and correlation length are obtained by the time relaxation method in the corotating with the Sun frame of reference in the region from 0.3 to 100 AU (but still inside the termination shock). The model equations include the effects of electron heat conduction, Coulomb collisions, photoionization of interstellar hydrogen atoms and their charge exchange with the solar wind protons, turbulence energy generation by pickup protons, and turbulent heating of solar wind protons and electrons. The turbulence transport model is based on the Reynolds decomposition and turbulence phenomenologies that describe the conversion of fluctuation energy into heat due to a turbulent cascade. In addition to using separate energy equations for the solar wind protons and electrons, a significant improvement over our previous work is that the turbulence model now uses an eddy viscosity approximation for themore » Reynolds stress tensor and the mean turbulent electric field. The approximation allows the turbulence model to account for driving of turbulence by large-scale velocity gradients. Using either a dipole approximation for the solar magnetic field or synoptic solar magnetograms from the Wilcox Solar Observatory for assigning boundary conditions at the coronal base, we apply the model to study the global structure of the solar wind and its three-dimensional properties, including embedded turbulence, heating, and acceleration throughout the heliosphere. The model results are compared with plasma and magnetic field observations on WIND, Ulysses, and Voyager 2 spacecraft.« less