Interplanetary Pick-Up Ion Acceleration A Study of Anisotropic Phase-Space Diffusion

Interplanetary pick-up ions originate from ionizations of neutral interstellar atoms in the heliosphere. Over the past periods it was generally expected that after pick-up by the frozen-in solar wind magnetic fields these ions quickly isotropize in velocity space by strong pitch- angle scattering, they do, however, not assimilate to the ambient solar wind ions. Meanwhile careful investigations of pick-up ion data obtained with the plasma analyzers on AMPTE and ULYSSES could clearly reveal that, especially at periods of flow-aligned fields, noticeably anisotropic distributions must prevail. To better understand the evolutionary tracks of pick-up ions in interplanetary phase-space we carried out an injection study which takes into account all relevant convection and diffusion processes, i.e. describing pitch angle scattering, adiabatic cooling, drifts and energy diffusion. As demonstrated here particles injected at 1 AU establish a distribution function with substantial anisotropies up to distances beyond 6 AU. Only under the action of fairly strong isotropic turbulence levels a trend towards isotropy can be recognized. The bulk velocity of the injected pick-up ions turns out to be remarkably smaller than the solar wind velocity. It also is obvious that pick-ups are strongly spread out from that solar wind plasma parcel into which they were originally implanted. As one consequence it must be concluded that the derivation of interstellar He gas parameters, using He pick-up ion flux data, require appreciable caution. Due to anisotropic spatial diffusion the location of the LISM helium cone axis, i.e. the LISM wind vector, and the LISM helium temperature are hidden in the associated He⁺pick-up ion flux patterns. This revised version was published online in July 2006 with corrections to the Cover Date.     

Physics and Gasdynamics of the Heliospheric Interface

During 30 years, a big theoretical effort to understand the physical processes in the heliospheric interface has followed the pioneer papers by Parker (1961) and Baranov et al. (1971). The heliospheric interface is a shell formed by the solar wind interaction with the ionized component of the circumsolar local interstellar medium (LISM). For fully ionized supersonic interstellar plasma two-shocks (the termination shock and the bow shock) and a contact discontinuity (the heliopause) are formed in the solar wind/LISM interaction. However, LISM consists of at least of three components additional to plasma: H-atoms, galactic cosmic rays and magnetic field. The interstellar atoms that penetrate into the solar wind, are ionized there and form pickup ions. A part of the pickup ions is accelerated to high energies of anomalous cosmic rays (ACRs). ACRs may modify the plasma flow upstream the termination shock and in the heliosheath. In this short review I summarize current understanding of the physical and gasdynamical processes in the heliospheric interface, outline unresolved problems and future perspectives. This revised version was published online in July 2006 with corrections to the Cover Date.     

The “FIP Effect” and the Origins of Solar Energetic Particles and of the Solar Wind

We find that the element abundances in solar energetic particles (SEPs) and in the slow solar wind (SSW), relative to those in the photosphere, show different patterns as a function of the first ionization potential (FIP) of the elements. Generally, the SEP and SSW abundances reflect abundance samples of the solar corona, where low-FIP elements, ionized in the chromosphere, are more efficiently conveyed upward to the corona than high-FIP elements that are initially neutral atoms. Abundances of the elements, especially C, P, and S, show a crossover from low to high FIP at \({\approx}\,10~\mbox{eV}\) in the SEPs but \({\approx}\,14~\mbox{eV}\) for the solar wind. Naively, this seems to suggest cooler plasma from sunspots beneath active regions. More likely, if the ponderomotive force of Alfvén waves preferentially conveys low-FIP ions into the corona, the source plasma that eventually will be shock-accelerated as SEPs originates in magnetic structures where Alfvén waves resonate with the loop length on closed magnetic field lines. This concentrates FIP fractionation near the top of the chromosphere. Meanwhile, the source of the SSW may lie near the base of diverging open-field lines surrounding, but outside of, active regions, where such resonance does not exist, allowing fractionation throughout the chromosphere. We also find that energetic particles accelerated from the solar wind itself by shock waves at corotating interaction regions, generally beyond 1 AU, confirm the FIP pattern of the solar wind.     

The Multifluid Character of the `Baranov' Interface

Since about three decades now it is clearly recognized that the interaction of the solar system with the ambient interstellar medium flow mainly is characterized by its hydrodynamic nature invoking structures like the inner shock, the heliopause and the outer shock with plasma sheath regions in between. After the pioneering works by Eugene Parker and Vladimir Baranov the main outlines of this interaction scenario were established, while some discussion on location and geometry of these structures is still going on till now. Fundamentally new aspects of this interaction problem have meanwhile appeared calling for new and more consistent calculations. The revisions of the earlier interaction concept starts with the neutral LISM gas component passing through the solar system. At the occasion of ionizations of this component a medium-energetic plasma component in form of keV-energetic pick-up ions is created. This component changes the distant solar wind properties by mass-, momentum-, and energy-loading, by wave generation and lowering the solar wind Mach numbers. Furthermore pick-up ions serve as a seed population for a high-energetic plasma population with energies between 10 and 100 MeV/nuc called anomalous cosmic rays. This latter component by means of its pressure gradient not only modifies the solar wind flow but also modulates its termination shock. In this paper it is shown how the main features of the enlarged interaction scenario change if the above mentioned multifluid character of the scenario is taken into account. While now we present a `multicolour vision' of the interacting heliosphere, it should never be forgotten that these modern views only were possible due to the fundamental `black-and-white vision' already presented by Baranov in the seventieths. This revised version was published online in July 2006 with corrections to the Cover Date.     

An axisymmetric magnetohydrodynamic model for the interaction of the solar wind with the local interstellar medium

We numerically analyze a magnetohydrodynamic, steady-state model for the interaction of a spherically symmetric solar wind with a three-component local interstellar medium (LISM), which is composed of plasma, hydrogen atoms, and a magnetic field. The magnetic field is assumed to be parallel to the velocity in the LISM. In this case, the model is axisymmetric. We study the effects of magnetic field on the plasma-flow geometry and on the distribution of hydrogen-atom parameters. In particular, we show that the presence of hydrogen atoms does not affect the qualitative change in the shape of the bow shock, the heliopause, and the solar-wind shock with increasing strength of the interstellar magnetic field. The presence of a magnetic field in the LISM can strongly affect the parameters of the energetic hydrogen atoms originated in the solar wind, although its effect on the “hydrogen wall” observed with the GHRS instrument onboard the HST spacecraft (Linsky and Wood 1996) is marginal.     

Signatures of the Interplanetary Helium Cone Reflected by Pick-up Ions

As known for a long time, interstellar wind neutral helium atoms deeply penetrate into the inner heliosphere and, when passing through the solar gravity field, form a strongly pronounced helium density cone in the downwind direction. Helium atoms are photoionized and picked-up by the solar wind magnetic field, but as pick-up ions they are not simply convected outwards with the solar wind in radial directions as assumed in earlier publications. Rather they undergo a complicated diffusion-convection process described here by an appropriate kinetic transport equation taking into account adiabatic cooling and focusing, pitch angle scattering and energy diffusion. In this paper, we solve this equation for He⁺pick-up ions which are injected into the solar wind mainly in the region of the helium cone. We show the resulting He⁺pick-up ion density profile along the orbit of the Earth in many respects differs from the density profile of the neutral helium cone: depending on solar-wind-entrained Alfvénic turbulence levels, the density maximum when looking from the Earth to the Sun is shifted towards the right side of the cone, the ratio of peak-densities to wing-densities varies and a left-to-right asymmetry of the He⁺-density profile is pronounced. Derivation of interstellar helium parameters from these He⁺-structures, such as the local interstellar medium (LISM) wind direction, LISM velocity and LISM temperature, are very much impeded. In addition, the pitch-angle spectrum of He⁺pick-up ions systematically becomes more anisotropic when passing from the left to the right wing of the cone structure. All effects mentioned are more strongly pronounced in high velocity solar wind compared to the low velocity solar wind.     

Interstellar matter and the location of the shock front

The initially supersonic flow of the solar wind passes through a magnetic shock front where its velocity is supposed to be reduced to subsonic values. The location of this shock front is primarily determined by the energy density of the external interstellar magnetic field and the momentum density of the solar wind plasma. Interstellar hydrogen penetrating into the heliosphere undergoes charge exchange processes with the solar wind protons and ionization processes by the solar EUV radiation. This results in an extraction of momentum from the solar wind plasma. Changes of the geometry and the location of the shock front due to this interaction are studied in detail and it is shown that the distance of the magnetic shock front from the Sun decreases from 200 to 80 AU for an increase of the interstellar hydrogen density from 0.1 to 1.0 cm⁻³. The geometry of the shock front is essentially spherical with a pronounced embayment in the direction opposite to the approach of interstellar matter which depends very much on the temperature of the interstellar gas. Due to the energy loss by the interaction with neutral matter the solar wind plasma reduces its velocity with increasing distance from the Sun. This modifies Parker's solution of a constant solar wind velocity.     

Low energy neutral atoms imaging of the Moon

Imaging of low-energy neutral atoms (LENAs) in the vicinity of the Moon can provide wide knowledge of the Moon from the viewpoint of plasma physics and planetary physics. At the surface of the Moon, neutral atoms are mainly generated by photon-stimulated desorption, micrometeorite vaporization and sputtering by solar wind protons. LENAs, the energetic neutral atoms with energy range of 10-500 eV, are mainly created by sputtering of solar wind particles. We have made quantitative estimates of sputtered LENAs from the Moon surface. The results indicate that LENAs can be detected by a realistic instrument and that the measurement will provide the global element maps of sputtered particles, which substantially reflect the surface composition, and the magnetic anomalies. We have also found that LENAs around dark regions, such as the permanent shadow inside craters in the pole region, can be imaged. This is because the solar wind ions can penetrate shaded regions due to their finite gyro-radius and the pressure gradient between the solar wind and the wake region. LENAs also extend our knowledge about the magnetic anomalies and associated mini-magnetosphere systems, which are the smallest magnetospheres as far as one knows. It is thought that no LENAs are generated from mini-magnetosphere regions because no solar wind may penetrate inside them. Imaging such void areas of LENAs will provide another map of lunar magnetic anomalies.     

Anisotropy of > 35keV ions in corotating particle events at 1 AU

The anisotropy of 35–1000 keV ions in two corotating particle events associated with high-speed solar wind streams at 1 AU is examined in terms of the diffusion-convection propagation model using data from the Energetic Proton Anisotropy Spectrometer on ISEE-3. The calculated diffusive anisotropy in the solar wind frame is found to be sunward and closely field-aligned, with a nearly energy-independent magnitude of ～ 40%. For one stream, using the Voyager 2 data of Decker et al. (1981), a positive gradient of ～ 100%/AU is found for ? 50 keV ions between 1 and 4 AU. The observations do not appear to support the scatter-free propagation model and indicate that ions with energies as low as a few tens of keV may be in diffusive equilibrium with the solar wind in this class of events.     

Time variations of the Heliospheric Termination Shock Radius

In this work we study the temporal variation of the heliospheric termination shock radius using in-ecliptic combined measurements from different spacecraft at 1 AU near the Earth. The results show that the radius of the heliospheric termination shock varies in time with a period of 11 years. During some 11-year time periods the shock radius anti-correlates with the solar cycle activity, specifically with the sunspot number. The average radius is approximately 115 AU with minimum value 80 AU and maximum value 150 AU. These values are the upper limits since we do not take into account the charge exchange effects of solar wind with the interstellar neutrals. We also compare the results with those from other spacecraft (Helios 1 and Voyager 2). We find that Helios 1 measurements give almost the same result as the one obtained from measurements at 1 AU while Voyager 2 measurements give a heliospheric termination shock radius approximately 15 AU lower.     

Corotating Shock Waves and the Solar-wind Source of Energetic Ion bundances: Power Laws in $$/Q$Q$

We find that element abundances in energetic ions accelerated by shock waves formed at corotating interaction regions (CIRs) mirror the abundances of the solar wind modified by a decreasing power-law dependence on the mass-to-charge ratio \(A\)/\(Q\) of the ions. This behavior is similar in character to the well-known power-law dependence on \(A\)/\(Q\) of abundances in large gradual solar energetic particles (SEP). The CIR ions reflect the pattern of \(A\)/\(Q\), with \(Q\) values of the source plasma temperature or freezing-in temperature of 1.0?–?1.2 MK typical of the fast solar wind in this case. Thus the relative ion abundances in CIRs are of the form \((A\mbox{/}Q)^{a}\), where \(a\) is nearly always negative and evidently decreases with distance from the shocks, which usually begin beyond 1 AU. For one unusual historic CIR event where \(a \approx 0\), the reverse shock wave of the CIR seems to occur at 1 AU, and these abundances of the energetic ions become a direct proxy for the abundances of the fast solar wind.     

ACE Observations of the Bastille Day 2000 Interplanetary Disturbances

We present ACE observations for the six-day period encompassing the Bastille Day 2000 solar activity. A high level of transient activity at 1 AU, including ICME-driven shocks, magnetic clouds, shock-accelerated energetic particle populations, and solar energetic ions and electrons, are described. We present thermal ion composition signatures for ICMEs and magnetic clouds from which we derive electron temperatures at the source of the disturbances and we describe additional enhancements in some ion species that are clearly related to the transient source. We describe shock acceleration of 0.3–2.0 MeV nucl⁻¹ protons and minor ions and the relative inability of some of the shocks to accelerate significant energetic ion populations near 1 AU. We report the characteristics of < 20 MeV nucl⁻¹ solar energetic ions and < 0.32 MeV electrons and attempt to relate the release of energetic electrons to particular source regions.     

The Injection, Acceleration, and Dynamical Influence of Interstellar Pickup Ions at the Solar Wind Termination Shock

The interaction of interstellar pickup ions with the solar wind termination shock is reviewed and assessed. The pickup ions mass and momentum load the wind and increase its pressure, effects which decrease the strength of the shock and its distance from the Sun. The pickup hydrogen may contribute substantially to the "reflected" ion population, which should provide most of the dissipation at the supercritical quasi-perpendicular shock. A fraction of the pickup ions impinging on the shock is "injected" into the process of diffusive shock acceleration to form the anomalous cosmic ray component. An injection mechanism which accounts for the apparent absence of solar wind ions in the anomalous component is "shock surfing", in which pickup ions which approach the shock slowly may be trapped between the upstream Lorentz force and the shock potential and accelerated in the motional electric field beyond the energy threshold for diffusive shock acceleration. However, the simplest interpretation of shock surfing would favor less massive pickup ion species, in contradiction with Voyager observations of anomalous component composition. A possible extension of the shock surfing mechanism is considered, as well as other injection mechanisms. Finally, the pressure of the anomalous component may modify the structure of the termination shock, which in turn may influence injection rates. This revised version was published online in July 2006 with corrections to the Cover Date.     

The solar wind

The current status of our understanding of the nature and origin of the solar wind is briefly reviewed, with emphasis being placed on the need for wave-particle interactions to account for the main energy source as well as details of the particle distribution functions. There has been considerable progress in the theoretical treatment of various aspects of the physics of the solar wind but a complete understanding is not yet in sight. Arguments concerning the ultimate fate of the solar wind are reviewed, in particular those concerning the distance to the shock wave which marks the termination of supersonic flow. This is of particular significance in view of recent observations suggesting that the termination might occur at about 50 AU from the Sun.     

The solar wind interaction with Venus through the eyes of the Pioneer Venus Orbiter

Venus has no internal magnetic dynamo and thus its ionosphere and hot oxygen exosphere dominate the interaction with the solar wind. The solar wind at 0.72 AU has a dynamic pressure that ranges from 4.5 nPa (at solar max) to 6.6 nPa (at solar min), and its flow past the planet produces a shock of typical magnetosonic Mach number 5 at the subsolar point. At solar maximum the pressure in the ionospheric plasma is sufficient to hold off the solar wind at an altitude of 400 km above the surface at the subsolar point, and 1000 km above the terminators. The deflection of the solar wind occurs through the formation of a magnetic barrier on the inner edge of the magnetosheath, or shocked solar wind. Under typical solar wind conditions the time scale for diffusion of the magnetic field into the ionosphere is so long that the ionosphere remains field free and the barrier deflects almost all the incoming solar wind. Any neutral atoms of the hot oxygen exosphere that reach the altitude of the magnetosheath are accelerated by the electric field of the flowing magnetized plasma and swept along cycloidal paths in the antisolar direction. This pickup process, while important for the loss of the Venus atmosphere, plays a minor role in the deceleration and deflection of the solar wind. Like at magnetized planets, the Venus shock and magnetosheath generate hot electrons and ions that flow back along magnetic field lines into the solar wind to form a foreshock. A magnetic tail is created by the magnetic flux that is slowed in the interaction and becomes mass-loaded with thermal ions.The structure of the ionosphere is very much dependent on solar activity and the dynamic pressure of the solar wind. At solar maximum under typical solar wind conditions, the ionosphere is unmagnetized except for the presence of thin magnetic flux ropes. The ionospheric plasma flows freely to the nightside forming a well-developed night ionosphere. When the solar wind pressure dominates over the ionospheric pressure the ionosphere becomes completely magnetized, the flow to the nightside diminishes, and the night ionosphere weakens. Even at solar maximum the night ionosphere has a very irregular density structure. The electromagnetic environment of Venus has not been well surveyed. At ELF and VLF frequencies there is noise generated in the foreshock and shock. At low altitude in the night ionosphere noise, presumably generated by lightning, can be detected. This paper reviews the plasma environment at Venus and the physics of the solar wind interaction on the threshold of a new series of Venus exploration missions.     

Interaction of interstellar pick-up ions with the solar wind

The interaction of interstellar pick-up ions with the solar wind is studied by comparing a model for the velocity distribution function of pick-up ions with actual measurements of He⁺ ions in the solar wind. The model includes the effects of pitch-ang'e diffusion due to interplanetary Alfvén waves, adiabatic deceleration in the expanding solar wind and the radial variation of the source function. It is demonstrated that the scattering mean free path is in the range 0.1 AU and that energy diffusion can be neglected as compared with adiabatic deceleration. The effects of adiabatic focusing, of the radial variation of the neutral density and of a variation of the solar wind velocity with distance from the Sun are investigated. With the correct choice of these parameters we can model the measured energy spectra of the pick-up ions reasonably well. It is shown that the measured differential energy density of the pick-up ions does not vary with the solar wind velocity and the direction of the interplanetary magnetic field for a given local neutral gas density and ionization rate. Therefore, the comparison of the model distributions with the measurements leads to a quantitative determination of the local interstellar gas density.Paper dedicated to Professor Hannes Alfvén on the occasion of his 80th birthday, 30 May 1988.     

Fluctuations in the heliospheric hydrogen distribution induced by generalized and time-dependent interstellar boundary conditions

It is well known that the neutral component of the local interstellar medium can effectively pass through the plasma interface ahead of the solar system and can penetrate deeply into the inner heliosphere. Here we present a newly-developed theoretical approach to describe the distribution function of LISM neutral hydrogen in the heliosphere, also taking into account time-dependent solar and interstellar boundary conditions. For this purpose we start from a Boltzmann-Vlasov equation, Fourier-transformed with respect to space and time coordinates, in connection with correspondingly transformed solar radiation forces and ionization rates, and then arrive at semi-analytic solutions for the transformed hydrogen velocity distribution function. As interstellar boundary conditions we allow for very general, non-Maxwellian and time-dependent distribution functions to account for the case that some LISM turbulence patterns or nonlinear wave-like shock structures pass over the solar system. We consider this theoretical approach to be an ideal instrument for the synoptic interpretation of huge data samples on interplanetary Ly- resonance glow intensities registered from different celestial directions over extended periods of time. In addition we feel that the theoretical approach presented here, when applied to interplanetary resonance glow data, may permit the detection of genuine fluctuations in the local interstellar medium.     

Influence of the magnetic field on the density distribution of solar wind protons and cometary ions in the shock layer ahead of cometary ionospheres

The “mass loading” of the solar wind by cometary ions produced by the photoionization of neutral molecules outflowing from the cometary nucleus plays a major role in the interaction of the solar wind with cometary atmospheres. In particular, this process leads to a decrease in the solar wind velocity with a transition from supersonic velocities to subsonic ones through the bow shock. The so-called single-fluid approximation, in which the interacting plasma flows are considered as a single fluid, is commonly used in modeling such an interaction. However, it is occasionally necessary to know the distribution of parameters for the components of the interacting plasma flows. For example, when the flow of the cometary dust component in the interplanetary magnetic field is considered, the dust particle charge, which depends significantly on the composition of the surrounding plasma, needs to be known. In this paper, within the framework of a three-dimensional magnetohydrodynamic model of the solar wind flow around cometary ionospheres, we have managed to separately obtain the density distributions of solar wind protons and cometary ions between the bow shock and the cometary ionopause (in the shock layer). The influence of the interplanetary magnetic field on the position of the point of intersection between the densities with the formation of a region near the ionopause where the proton density is essentially negligible compared to the density of cometary ions is investigated. Such a region was experimentally detected by the Vega-2 spacecraft when investigating Comet Halley in March 1986. The results of the model considered below are compared with some experimental data obtained by the Giotto spacecraft under the conditions of flow around Comets Halley and Grigg–Skjellerup in 1986 and 1992, respectively. Unfortunately, our results of calculations on Comet Churyumov–Gerasimenko are only predictive in character, because the trajectory of the Rosetta spacecraft, which manoeuvred near its surface for several months, is complex.     

Observations of energetic ions near the Venus ionopause

Ions (primarily O⁺) with spacecraft rest frame energies >40 eV have been observed by the Pioneer Venus Neutral Mass Spectrometer. The signature occurs in about 13% of the 700 orbits examined, primarily near the ionopause and at all solar zenith angles. The energetic ions coincide in location with superthermal ions observed by the Ion Mass Spectrometer and more rarely occur in some of the plasma clouds observed by the Electron Temperature Probe. These observations in conjunction with measurements by the Plasma Wave Instrument near the ionopause suggest that the ions are accelerated out of ionospheric plasma by the shocked solar wind through plasma waveparticle interactions.     

Spatial Distribution of Solar Energetic Particles in the Inner Heliosphere

We study the spatial distribution of solar energetic particles (SEPs) throughout the inner heliosphere during six large SEP events from the period 1977 through 1979, as deduced from observations on the Helios 1 and 2, IMP 7 and 8, ISEE 3, and Voyager 1 and 2 spacecraft. Evidence of intensity maxima associated with the expanding shock wave is commonly seen along its central and western flanks, although the region of peak acceleration or “nose” of the shock is sometimes highly localized in longitude. In one event (1 January 1978) a sharp peak in 20?–?30 MeV proton intensities is seen more strongly by Voyager at ～?2 AU than it is by spacecraft at nearby longitudes at ～?1 AU. Large spatial regions, or “reservoirs,” often exist behind the shocks with spatially uniform SEP intensities and invariant spectra that decrease adiabatically with time as their containment volume expands. Reservoirs are seen to sweep past 0.3 AU and can extend out many AU. Boundaries of the reservoirs can vary with time and with particle velocity, rather than rigidity. In one case, a second shock wave from the Sun reaccelerates protons that retain the same hard spectrum as protons in the reservoir from the preceding SEP event. Thus reservoirs can provide not only seed particles but also a “seed spectrum” with a spectral shape that is unchanged by a weaker second shock.