Pickup Protons Research Articles

Transport Equation of Kappa Distributions in the Heliosphere

Abstract In this paper, we develop the transport equation of kappa, the fundamental thermodynamic parameter that labels kappa distributions of particle velocities. Using the recently developed concept of entropy defect, we are able to formulate the transport equation of kappa as a function of a general, positive or negative, rate of entropy change. Then, we derive the particular case of exchanging plasma ions with low-dimensionality, newly born pickup protons, which interact and decrease the entropy of the flow of otherwise kappa-distributed plasma protons. Finally, we apply the transport equation of kappa to the solar wind plasma protons, which leads to the radial profile of kappa values, as well as the evolution of the kappa distributions through the heliosphere. The results show that the solar wind kappa decreases with increasing heliocentric distance, corresponding to plasmas residing in stationary states far from classical thermal equilibrium. Moreover, in the outer heliosphere and the heliosheath, kappa reaches its lowest values and is spread across the far-equilibrium region of 1.5 < κ < 2.5, which coincides with independent observations provided by NASA’s Interstellar Boundary Explorer mission.

Adiabatic energy change in the inner heliosheath: how does it affect the distribution of pickup protons and energetic neutral atom fluxes?

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.

Upstream proton cyclotron waves at Mars during the passage of solar wind stream interaction regions

Proton cyclotron waves (PCWs) upstream from Mars are generated by ion–ion instabilities due to the interaction between the solar wind and the pickup protons that originate from the extended hydrogen (H) exosphere of Mars, indicating a loss of H to interplanetary space. Their occurrences and properties associated with nominal solar wind conditions have been widely studied, but relatively little is known about PCW activity under extreme solar wind conditions. In this work we characterize the upstream PCWs during the passage of solar wind stream interaction regions (SIRs) utilizing in situ observations by the Mars Atmosphere and Volatile EvolutioN spacecraft from October 2014 to August 2021. The 46 analyzed SIR events are a subset of the events that occurred in the perihelion season of Martian years 32–35. We find that the PCW occurrence rate is increased by a factor of about 1.8 during the SIR phase relative to the pre-SIR and post-SIR phases. Furthermore, the PCW activity tends to occur more frequently in the leading portion of the SIR. The PCWs detected during the SIR phase have more pronounced wave characteristics, namely larger wave amplitudes, smaller propagation angles with respect to the background magnetic field, and closer-to-circular polarization.

MAVEN Proton Observations Near the Martian Moon Phobos: Does Phobos Backscatter Solar Wind Protons?

AbstractESA's Mars Express (MEX) may have observed twice solar wind protons backscattered by Phobos. However, these detections remain uncertain and call for an independent confirmation. Here, we analyze the proton measurements collected by the Suprathermal And Thermal Ion Composition experiment onboard NASA's Mars Atmosphere and Volatile EvolutioN mission (MAVEN) during the 15 closest encounters of MAVEN with Phobos (closest approach between 80 and 370 km). We use a model of magnetic connectivity to the Martian bow shock and an exospheric pickup ion model to identify the origin of the protons observed close to the Martian moon. We find no evidence of protons backscattered by Phobos. Instead, foreshock‐shock foot protons and Martian exospheric pickup protons are often observed near Phobos. These results put into new context not only the past MEX detections, but also the future observation attempts of the upcoming JAXA Martian Moons Exploration mission.

In Situ Observation of Solar-flare-induced Proton Cyclotron Waves Upstream from Mars

Proton cyclotron waves (PCWs) upstream from Mars are usually interpreted as waves generated by ion/ion instabilities due to the interaction between the solar wind plasma and the pickup protons, originating from the extended hydrogen (H) exosphere of Mars. Their generation mainly depends on the solar wind properties and the relative density of the newborn protons with respect to the background solar wind. Under stable solar wind conditions, a higher solar irradiance leads to both increased exospheric H density and ionization rate of H atoms, and therefore a higher relative density, which tends to increase the linear wave growth rate. Here we show that the solar irradiance is likely to contribute significantly to PCW generation. Specifically, we present observations from the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft indicating that, around the peak of the X8.2 flare on 2017 September 10, the increased solar irradiance gave rise to higher pickup H+ fluxes, which in turn excited PCWs. This result has implications for inferring the loss of hydrogen to space in early Martian history with more intense and frequent X-class flares as well as their contributions to the total loss.

Cometary ions detected by the Cassini spacecraft 6.5 au downstream of Comet 153P/Ikeya–Zhang

During March–April 2002, while between the orbits of Jupiter and Saturn, the Cassini spacecraft detected a significant enhancement in pickup proton flux. The most likely explanation for this enhancement was the addition of protons to the solar wind by the ionization of neutral hydrogen in the corona of comet 153P/Ikeya–Zhang. This comet passed relatively close to the Sun–Cassini line during that period, allowing pickup ions to be carried to Cassini by the solar wind. This pickup proton flux could have been further modulated by the passage of the interplanetary counterparts of coronal mass ejections past the comet and spacecraft. The radial distance of 6.5 Astronomical Units (au) traveled by the pickup protons, and the implied total tail length of >7.5 au make this cometary ion tail the longest yet measured.

On the Energization of Pickup Ions Downstream of the Heliospheric Termination Shock by Comparing 0.52–55 keV Observed Energetic Neutral Atom Spectra to Ones Inferred from Proton Hybrid Simulations

We present an unprecedented comparison of ∼0.52–55 keV energetic neutral atom (ENA) heliosheath measurements, remotely sensed by the Interstellar Boundary Explorer (IBEX) mission and the Ion and Neutral Camera (INCA) on the Cassini mission, with modeled ENAs inferred from interstellar pickup protons that have been accelerated at the termination shock, using hybrid simulations, to assess the pickup ion energetics within the heliosheath. This is the first study to use hybrid simulations that are able to accurately model the acceleration of ions to tens of keV energies, which is essential in order to model ENA fluxes in the heliosheath, covering the full energy range observed by IBEX and CASSINI/INCA. The observed ENA intensities are an average value over the time period from 2009 to the end of 2012, along the Voyager 2 (V2) trajectory. The hybrid simulations upstream of the termination shock, where V2 crossed, are constrained by observations. We report an energy-dependent discrepancy between observed and simulated ENA fluxes, with the observed ENA fluxes being persistently higher than the simulated ones. Our analysis reveals that the termination shock may not accelerate pickup ions to sufficient energies to account for the observed ENA fluxes. We, thus, suggest that the further acceleration of these pickup ions is most likely occurring within the heliosheath, via additional physical processes like turbulence or magnetic reconnection. However, the redistribution of energy inside the heliosheath remains an open question.

Energetic pickup proton population downstream of the termination shock as revealed by IBEX-Hi data

ABSTRACT Pickup protons originate as a result of the ionization of hydrogen atoms in the supersonic solar wind, forming the suprathermal component of protons in the heliosphere. While they are being picked by the heliospheric magnetic field and convected into the heliosheath, the pickup protons may suffer stochastic acceleration from the solar wind turbulence in the region from the Sun up to the heliospheric termination shock, where they can also experience shock-drift acceleration or reflection from the cross-shock potential. These processes create a high-energy tail in the pickup ion energy distribution. The properties of this energetic pickup proton population are still not well defined, in spite of the fact that they are vital for models that simulate energetic neutral atom fluxes. We consider two scenarios for the pickup proton velocity distribution downstream of the heliospheric termination shock (a filled shell with an energetic power-law tail, and bi-Maxwellian). Based on a numerical kinetic model and observations of the energetic neutral atom fluxes from the inner heliosheath by the IBEX-Hi instrument, we characterize the pickup proton distribution and provide estimations of the properties of the energetic pickup proton population downstream of the termination shock.

Hybrid Simulations of Interstellar Pickup Protons Accelerated at the Solar-wind Termination Shock at Multiple Locations

Abstract We estimate the intensity of interstellar pickup protons accelerated to ∼50 keV at various locations along the solar-wind termination shock, using two-dimensional hybrid simulations. Parameters for the solar wind, interstellar pickup ions (PUIs), and magnetic field just upstream of the termination shock at one flank of the heliosphere, and at the location in the downwind (or tail-ward) direction are based on a solar-wind/pickup-ion/turbulence model. The parameters upstream of the shock where Voyager 2 crossed are based on observations. The simulation is limited in size, and therefore cannot accurately model the distribution to energies much beyond ∼50 keV. This is sufficient to study the origin of the high-energy tail of the distribution, which is the low-energy portion of the anomalous cosmic-ray spectrum. We also extrapolate our results to other locations along the termination shock, such as the other flank, and the poles of the heliosphere. We find that the intensity of ∼10–50 keV accelerated pickup protons is remarkably similar at all three locations we simulated, suggesting that particles in this energy range are relatively uniformly distributed along the termination shock, and are likely quite uniform throughout the entire heliosheath. In addition, we find significant differences in the distribution in the 0.5–1 keV energy range for energetic neutral atoms coming from the tail region of the heliosphere compared to that at the nose or flank look directions. This is because the peak in the PUI distribution is at a higher energy there.

Probing the thermodynamic conditions of the heliosheath plasma by shock wave propagation

Context. The pressure equilibrium between the inner heliosheath and the outer heliosheath (referred to as the local interstellar medium) is an eminent theoretical and practical problem; theoretical, because the relevant pressure carriers have to be identified, and practical, because data must be gathered in order to confirm such a pressure equilibrium. The problem is closely connected with the stability of the heliopause, that is, of the tangential discontinuity between these two counterflowing media, and is of utmost importance for understanding the stability of the whole circumsolar plasma structure. Aims. In this paper we analyze the thermodynamic conditions of the multi-fluid plasma between the solar wind termination shock and the heliopause determining the total heliosheath pressure. We look into this problem from a theoretical standpoint and revisit theoretical descriptions of the solar wind plasma after its passage over the solar wind termination shock, thereafter forming the subsonic heliosheath region. Methods. Hereby we take into account the 3D magnetohydrodynamics shock conditions and the resulting 3D temperature structure of the downstream plasma flow. We use a kind of seismological procedure to probe the heliosheath plasma by inquiring into the propagation conditions of traveling shock wave perturbations in this predetermined 3D heliosheath plasma structure. We discuss the fact that the front geometry of such a traveling shock wave most probably does not remain spherical, if it was to begin with, due to asymmetric shock propagation conditions. In contrast, the wave front is likely to become strongly deformed into an upwind bulge. Results. Concerning the plasma pressure, in addition to solar wind and pick-up proton pressures, we have to take into account the solar wind electron pressure which as a surprise turns out to be of comparable magnitude. As a consequence, the characteristic propagation speed of the traveling shock wave in the weakly magnetized heliosheath plasma is given as a mixed speed expressed by the sound speeds of the protons and the electrons. We describe local low-energy proton density signatures that can be found in Voyager-2 proton data as a consequence of traveling shock wave passages and show that the total local plasma pressure can be directly derived from them.

Turbulent Heating in Solar Wind Thermodynamics

Abstract This paper considers the concept of wave-particle thermodynamic equilibrium in order to improve our understanding of the role of turbulent heating in the solar wind proton plasma. The thermodynamic equilibrium in plasmas requires the energy of a plasmon—the quantum of plasma fundamental oscillation—to be balanced by the proton-magnetized plasma energy, that is, the magnetic field and proton kinetic/thermal energy. This equilibrium has already been confirmed in several prior analyses, but also in this paper, by analyzing (i) multi-spacecraft data sets along the radial profile of the inner heliosphere, and (ii) representative data sets of a variety of 27 different space and astrophysical plasmas. Recently, it was shown that the slow mode of the near-Earth solar wind plasma is characterized by a missing energy source that is necessary for keeping the energy balance in the plasmon–proton-magnetized plasma. Here we show strong evidence that this missing energy is the turbulent energy heating the solar wind. In particular, we derive and compare the radial and velocity profiles of this missing energy and the turbulent energy in the inner heliosphere, also considering other minor contributions, such as the temperature of pickup protons. The connection of the missing plasmon–proton energy with the turbulent energy provides a new method for estimating and cross-examining the turbulent energy in space and astrophysical plasmas, while it confirms the universality of the involved new Planck-type constant that implies a large-scale quantization.

Connection of Turbulence with Polytropic Index in the Solar Wind Proton Plasma

This paper improves our understanding of the interplay of the proton plasma turbulent heating sources of the expanding solar wind in the heliosphere. Evidence is shown of the connections between the polytropic index, the rate of the heat absorbed by the solar wind, and the rate of change of the turbulent energy, which heats the solar wind in the inner and outer heliosphere. In particular, we: (i) show the theoretical connection of the rate of a heat source, such as the turbulent energy, with the polytropic index and the thermodynamic process; (ii) calculate the effect of the pick-up protons in the total proton temperature and the relationship connecting the rate of heating with the polytropic index; (iii) derive the radial profiles of the solar wind heating in the outer and inner heliosphere; and (iv) use the radial profile of the turbulent energy in the solar wind proton plasma in the heliosphere, in order to show its connection with the radial profiles of the polytropic index and the heating of the solar wind.

The behavior of electrons at the heliospheric shock transition

The jump of the temperature of thermal protons and the compression at the solar wind termination shock, measured by the Voyager 2 spacecraft during its crossing of the shock, can be explained if one assumes that the downstream temperature of electrons essentially higher than that of the thermal protons. The upstream temperature of the electrons has little effect on this result. For instance, a tenfold increase in the upstream electron temperature relative to the thermal proton temperature results in the temperature jump change by only 10% and the compression change by 1%. Within the framework of the considered model, changing the relative content of pickup protons over a wide range does not affect the temperature jump and compression at the termination shock.

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.

The Downwind Hemisphere of the Heliosphere: Eight Years of IBEX-Lo Observations

Abstract We present a comprehensive study of energetic neutral atoms (ENAs) of 10 eV to 2.5 keV from the downwind hemisphere of the heliosphere. These ENAs are believed to originate mostly from pickup protons and solar-wind protons in the inner heliosheath. This study includes all low-energy observations made with the Interstellar Boundary Explorer over the first eight years. Because the protons around 0.1 keV dominate the plasma pressure within the inner heliosheath in downwind direction, these ENA observations offer the unique opportunity to constrain the plasma properties and dimensions of the heliosheath where no in situ observations are available. We first derive energy spectra of ENA intensities averaged over time for 49 macropixels covering the entire downwind hemisphere. The results confirm previous studies regarding integral intensities and the roll-over around 0.1 keV energy. With the expanded data set, we now find that ENA intensities at 0.2 and 0.1 keV seem to anti-correlate with solar activity. We then derive the product of total plasma pressure and emission thickness of protons in the heliosheath to estimate lower limits on the thickness of the inner heliosheath. The temporally averaged ENA intensities support a rather spherical shape of the termination shock and a heliosheath thickness between 150 and 210 au for most regions of the downwind hemisphere. Around the nominal downwind direction of 76° ecliptic longitude, the heliosheath is at least 280 au thick. There, the neutral hydrogen density seems to be depleted compared to upwind directions by roughly a factor of 2.

A Test of the Interstellar Boundary EXplorer Ribbon Formation in the Outer Heliosheath

Abstract NASA’s Interstellar Boundary EXplorer (IBEX) mission is imaging energetic neutral atoms (ENAs) propagating to Earth from the outer heliosphere and local interstellar medium (LISM). A dominant feature in all ENA maps is a ribbon of enhanced fluxes that was not predicted before IBEX. While more than a dozen models of the ribbon formation have been proposed, consensus has gathered around the so-called secondary ENA model. Two classes of secondary ENA models have been proposed; the first class assumes weak scattering of the energetic pickup protons in the LISM, and the second class assumes strong but spatially localized scattering. Here we present a numerical test of the “weak scattering” version of the secondary ENA model using our gyro-averaged kinetic model for the evolution of the phase-space distribution of protons in the outer heliosheath. As input for our test, we use distributions of the primary ENAs from our MHD-plasma/kinetic-neutral model of the heliosphere-LISM interaction. The magnetic field spectrum for the large-scale interstellar turbulence and an upper limit for the amplitude of small-scale local turbulence (SSLT) generated by protons are taken from observations by Voyager 1 in the LISM. The hybrid simulations of energetic protons are also used to set the bounding wavenumbers for the spectrum of SSLT. Our test supports the “weak scattering” version. This makes an additional solid step on the way to understanding the origin and formation of the IBEX ribbon and thus to improving our understanding of the interaction between the heliosphere and the LISM.

Listing of 502 Times When the Ulysses Magnetic Fields Instrument Observed Waves Due to Newborn Interstellar Pickup Protons

Abstract In two earlier publications we analyzed 502 intervals of magnetic waves excited by newborn interstellar pickup protons that were observed by the Ulysses spacecraft. Due to the considerable effort required in identifying these events, we provide a list of the times for the 502 wave event intervals previously identified. In the process, we provide a brief description of how the waves were found and what their properties are. We also remind the reader of the conditions that permit the waves to reach observable levels and explain why the waves are not seen more often.

PROTON HEATING BY PICK-UP ION DRIVEN CYCLOTRON WAVES IN THE OUTER HELIOSPHERE: HYBRID EXPANDING BOX SIMULATIONS

ABSTRACT Using a one-dimensional hybrid expanding box model, we investigate properties of the solar wind in the outer heliosphere. We assume a proton–electron plasma with a strictly transverse ambient magnetic field and, aside from the expansion, we take into account the influence of a continuous injection of cold pick-up protons through the charge-exchange process between the solar wind protons and hydrogen of interstellar origin. The injected cold pick-up protons form a ring distribution function, which rapidly becomes unstable, and generate Alfvén cyclotron waves. The Alfvén cyclotron waves scatter pick-up protons to a spherical shell distribution function that thickens over that time owing to the expansion-driven cooling. The Alfvén cyclotron waves heat solar wind protons in the perpendicular direction (with respect to the ambient magnetic field) through cyclotron resonance. At later times, the Alfvén cyclotron waves become parametrically unstable and the generated ion-acoustic waves heat protons in the parallel direction through Landau resonance. The resulting heating of the solar wind protons is efficient on the expansion timescale.

The Sun's Dynamic Influence on the Outer Heliosphere, the Heliosheath, and the Local Interstellar Medium

The Sun has been observed for many years to be a dynamic influence in the heliosphere, and as the Voyager missions have continued long after achieving their original goals of observing the major planets they have provided the first in situ observations of the effects of solar activity in the heliosheath (HS), and the nearest portions of the local Interstellar Medium (LISM). Comparing these observations with models provides key insights. We employ two three-dimensional (3D) time-dependent models that simulate the propagation of shocks, other specific features, and the background solar wind throughout the heliosphere, starting with the solar background and solar event boundary conditions near the Sun at 2.5 Rs. The Hybrid Heliospheric Modeling System with Pickup Protons (HHMS-PI) is a 3D time- dependent Magnetohydrodynamic (MHD) simulation. HAFSS (HAF Solar Surface) is a 3D time-dependent kinematic simulation. Comparing our models with the observations indicates that solar effects are seen in the heliosphere, the HS, and the LISM in in-situ spacecraft measurements of plasma, magnetic field, energetic particles, cosmic rays, and plasma waves. There is quantitative agreement (at ACE, Ulysses, VI, V2) with data (e.g., solar wind, IMF, Ulysses SWICS pickup protons (PUPs)). Propagating shocks are slowed due to PUPs. The 3D locations of solar events and of various spacecraft are key to understanding the 3D propagation and timing of shocks, other specific features, and gradients throughout the heliosphere, HS, and LISM.

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.