Anomalous Cosmic Rays and Heliospheric Energetic Particles

We present a review of Anomalous Cosmic Rays (ACRs), including the history of their discovery and recent insights into their acceleration and transport in the heliosphere. We focus on a few selected topics including a discussion of mechanisms of their acceleration, escape from the heliosphere, their effects on the dynamics of the heliosheath, transport in the inner heliosphere, and their solar cycle dependence. A discussion concerning their name is also presented towards the end of the review. We note that much is known about ACRs and perhaps the term Anomalous Cosmic Ray is not particularly descriptive to a non specialist. We suggest that the more-general term: “Heliospheric Energetic Particles”, which is more descriptive, for which ACRs and other energetic particle species of heliospheric origin are subsets, might be more appropriate.

The cosmic-ray spectrum at 1 AU. This figure is from Giacalone et al. (2012), and originally from Jokipii (1990)

Introduction
In the early 1970's, analysis of the spectrum of galactic cosmic rays (GCR) revealed an unusual, or "anomalous", enhancement at lower energies. In particular, the GCR intensity did not decrease with decreasing energy, as was expected based on our understanding that low-energy GCRs were unable to reach 1 AU due to their interaction with the solar wind. The anomalous enhancements were most notable in the spectra of oxygen and helium. A solar origin to these enhancements was not favored (Garcia-Munoz et al. 1973;Hovestadt et al. 1973;McDonald et al. 1974). These have since become known as anomalous cosmic rays (ACRs). Figure 1 shows the full cosmic ray spectrum, including both GCRs and ACRs. GCRs are those indicated with the black squares, for which the spectrum has an approximately power law dependence on energy above a few GeV, and turns over, or is "modulated" at energies below about 1 GeV. This modulation is caused by the interaction of GCRs, most of which are produced in supernovae remnants very far form the Sun, with our heliosphere. Low energy ( GeV) GCRs have considerable difficulty reaching the inner heliosphere due to their interaction with the solar wind. ACRs are visible as an additional enhancement in the spectrum of oxygen (cross symbols), as indicated in the figure with the arrow, at about 100 MeV, which is a bit below the modulated part of the GCR oxygen spectrum. ACRs are enhancements in the spectrum of helium, nitrogen, oxygen, neon, protons, but not in carbon (Klecker 1995). The intensity of ACRs were observed to increase with distance from the Sun, suggesting their source is in the distant heliosphere (McDonald et al. 1974).
It was recognized soon after their discovery, by Fisk et al. (1974), that ACRs are accelerated interstellar pickup ions. These are interstellar atoms which have become ionized through either charge exchange with the solar wind, or photo-ionization by solar UV. Interstellar atoms enter the solar system because the Sun moves relative to the local interstellar medium (e.g. Frisch 1996). Upon ionization, a pickup ion is subject to forces which cause it to drift away from the Sun towards the outer heliosphere. Moreover, a freshly ionized pickup ion has an energy of about ∼ 1 keV relative to a frame of reference moving with the solar wind and must be significantly accelerated to achieve the energy associated with ACRs. It was first noted by Pesses et al. (1981) that this energization could be achieved by the interaction with the solar wind termination shock. The acceleration mechanism, in this case, is known as diffusive shock acceleration discovered in the late 1970's (Axford et al. 1977;Krymsky 1977;Bell 1978;Blandford and Ostriker 1978). This theory is based on solutions to the well-known cosmic-ray transport equation (Parker 1965) for the distribution function of particles interacting with a shock-like plasma discontinuity. Jokipii (1986) solved this equation giving the first quantitative calculation for the distribution of ACRs in the heliosphere, including acceleration at the termination shock and transport in the heliosphere.

Energy Dependent Evolution of Energy Spectrum and Intensity Across the Termination Shock
Prior to the Voyager crossings of the termination shock, ACRs were observed to have large, positive radial gradients, indicating the source was more distant (McDonald et al. 1974). It was also noted that the low-energy end of ACR (H and He) spectra revealed a large increase with radial distance compared to the higher energy end (Stone et al. 1999). Assuming that the particles diffuse from the outer heliosphere, this observation suggests the lower-energy particles have a smaller diffusion coefficient than the higher energy ones. It was also expected that the energy spectrum would become a power law, at least up to the typical ACR energy of several tens to about 200 MeV (total energy), at the time of the crossing of the termination shock because it was thought that the shock was the source of the particles (e.g. Pesses et al. 1981). This was the prediction of the models (e.g. Jokipii 1986;Potgieter and Moraal 1988;Steenberg and Moraal 1996). However, neither Voyager observed the expected power-law energy spectrum upon crossing the termination shock, and instead observed that in the 3-30 MeV energy range, the spectrum did not increase at the shock, remaining "modulated" (Stone et al. 2005. Later, when each spacecraft was deeper into the heliosheath, the spectrum eventually filled in with a nearly power law dependence on energy with a spectral index of about −1.6 . Clearly, the particles above a few MeV were not accelerated locally where each spacecraft crossed the shock. Instead, their fluxes continued to rise as the spacecraft moved into the heliosheath. This suggests they are accelerated by some other mechanism than shock acceleration, or are accelerated elsewhere on the shock. This is discussed further below. In contrast, at energies below a few MeV, both Voyagers observed significant enhancements of the particle flux at the time of the shock crossing (Decker et al. 2005(Decker et al. , 2008. This was found to be consistent with the acceleration of pickup ions locally at the termination shock (Giacalone and Decker 2010; see also Florinski et al. 2009).

Characteristics of Termination Shock Particles
As noted in the previous section, the term "TSP" refers to particles whose flux was observed to peak locally at the crossing of the termination shock (Decker et al. 2005(Decker et al. , 2008, and were also observed to increase in intensity for a few years before the shock crossings (see Fig. 13). Unusual TSP anisotropies were observed by V1 in 2001-2003, prior to the shock crossing. These anisotropies were directed away from the Sun and were nearly aligned with the Parker spiral magnetic field. At the time, this was thought to be unusual because it suggested a source inside the location of the spacecraft, and not beyond it. In fact, it was suggested that V1 had crossed the termination shock (Krimigis et al. 2003), which, at the time, was assumed to be roughly spherically shaped. However, it was soon realized that the observations could be understood if the surface of the shock was blunt shaped, while the Parker spiral magnetic field lines were very nearly perpendicular to the radial direction . A key prediction of this model was that when V2 crossed the shock, the anisotropies would be observed to be moving in the opposite direction (sunward) since the V2 crossing would be on the other side of the sub-solar point of the blunt-shaped shock compared to V1. This was indeed observed (Decker et al. 2005(Decker et al. , 2008, confirming the picture put forth by . The current paradigm is that the termination shock is blunt shaped. This is also seen in large-scale fluid simulations of the heliosphere (e.g. Zank 1999).
TSPs were also observed to be highly variable with intermittent increases in intensity. It was suggested that these variations, and the significant field-aligned anisotropies, could be understood as the spacecraft crossing through magnetic field lines that were either filled or devoid of energetic particles depending on the magnetic connectivity to the source (assumed to be the termination shock) (Giacalone and Jokipii 2006). In addition, it was also observed that some of the intermittent intensity enhancements occurred at times of specific magnetic polarity reversals associated with crossings of the heliospheric current sheet (HCS) (Richardson et al. 2006). This was interpreted by Giacalone and Burgess (2010) as the result of the interaction and acceleration of pickup ions at the point of intersection of the termination shock and HCS of a particular polarity such that the accelerated particles drift upstream, away from the shock, along the HCS.
In the heliosheath, particles below an MeV or so, were observed to be remarkably uniform in intensity as seen by both Voyagers 1 and 2, despite the large difference in their trajectories leading to a > 100 AU separation. At the heliopause crossing, their intensity dropped to background levels, indicating their rapid escape from the heliosphere at this point . This is discussed further in Sect. 4. Regarding the uniform intensity in the heliosheath, this suggests that their acceleration at the shock is also uniform along the shock surface. This was supported by numerical simulation of pickup ion acceleration at the termination shock at multiple locations (Giacalone et al. 2021).

ACR Charge States
Analysis of observations of ACRs revealed a rather low charge state. Klecker et al. (1980) reported an upper bound for ACR oxygen O of 4, and later,  found that the observations were consistent with a charge state of 1. ACRs with higher charge states have been reported as well, as discussed below, but it is now generally understood that most ACRs below 100-200 MeV total energy are singly charged. Adams and Leising (1991) noted that a singly charged ACR moving in the interstellar medium will lose additional electrons if the source is beyond about 0.2 pc. This can also be thought of as a constraint on the age to be less than about four years. This places a rather stringent constraint on the time scale associated with the acceleration of pickup ions to the observed ACR energies. Jokipii (1992) noted that this time scale could be achieved most naturally for acceleration at the nearly perpendicular termination shock, and other mechanisms were too slow.
Multiply charged ACRs have also been observed (Mewaldt et al. 1996). This work suggests an even more-stringent constraint of the acceleration time scale of about one year. Large-scale global modeling of ACRs, including acceleration at the termination shock, transport in the heliosphere, and further stripping of electrons via interaction with the background plasma using known formulas for the cross-sections, by Jokipii (1996) found that higherenergy ACRs, > 300 MeV, change from being singly charged to multiply charged. This modeling, which agreed well with the observations, was expanded upon by Barghouty et al. (2000) using updated hydrogen-impact ionization rates giving an even better agreement with the observations. An important key takeaway from the historical observations of ACR charge states is that the acceleration time scale must be of the order of a year. In a separate study, analysis of Pioneer-10 observations showed that there was a shift in the correlation between ACR-fluxes with Lyman-alpha data (a proxy for solar variability) that was attributed to the time it takes to accelerate pickup ions to the observed ACR energies (Scherer et al. 1997). The time scale was found to be less than a year, which is consistent with the charge-state observations. Given that the initial energy of the particles is of the order of a keV, the energy of a freshly ionized pickup ion, there must be a nearly 5 order-of-magnitude-increase in the particle energy on a time scale of a year.

Recovery of ACR Intensity Following a Global Merged Interaction Region: Pinpointing the ACR Source Location
In the early 1990's, a large global merged interaction region (GMIR), a large magnetic/plasma disturbance associated with intense solar activity, moved through the heliosphere causing a depletion in the intensities of GCRs and ACRs (McDonald et al. 2000;Webber et al. 2001). The reduction of the GCR intensity is the well-known Forbush-decrease phenomenon. McDonald et al. (2000) noted that the recovery time of the flux of ACRs and GCRs were different, with the ACRs recovering considerably more rapidly than the GCRs. This suggests that the ACR source is nearer to the spacecraft than is the GCR source. In fact, these authors used the recovery time to estimate the distance to the source, coming up with a location of the termination shock between 80-96 AU. This is consistent with the distances at which each Voyager later crossed the termination shock. Moreover, time-dependent simulations of ACRs and GCRs during the passage of a GMIR through the outer heliosphere by , see also le Roux and Fichtner (1999), obtained results that were in agreement with the observations. These calculations assumed an ACR source at the termination shock.

General Considerations: Diffusive Shock Acceleration
Particle acceleration at shocks has been studied for many years, with a major breakthrough in understanding occurring in the late 1970's with the development of diffusive shock acceleration (DSA) theory (Axford et al. 1977;Bell 1978;Blandford and Ostriker 1978;Krymsky 1977). Perri et al, this journal, also addresses the current state of knowledge of this topic and this topic was also recently reviewed by Liu and Jokipii (2021). The theory is based on the fact that charged particles are accelerated by any compression of the plasma, and shocks provide a very strong plasma compression owing to their the change in plasma density over a scale that is of the order of the thermal ion inertial length. In DSA theory, the rate of energy gain depends on how effectively the particles are trapped near the shock, through their diffusion in turbulent magnetic fields, which are either pre-existing (independent of the shock), or associated with plasma instabilities related to the shock. The fundamental equation of DSA theory is known as the Parker equation (Parker 1965), which is also used in cosmic-ray transport theory. The equation gives the phase-space distribution function, f , as a function of position vector x i , momentum magnitude, p and time, t . It is given by: where κ ij is the diffusion tensor which contains components associated with spatial diffusion both along and across the magnetic field (symmetric components), as well as large-scale guiding center drifts (the anti-symmetric components), U i = U is the plasma velocity, and Q is a source. The term that gives acceleration, or deceleration, is the second from the last, and one notes immediately that it is proportional to the divergence of the flow speed, which from the continuity equation in fluid dynamics, is related to the plasma compression. There are two key assumptions in the derivation of this equation, one is that the particle distribution is assumed to be isotropic in pitch angle, and the other is that the electric field, E is that of ideal MHD, given by E = −(1/c)U × B, where B is the magnetic field vector, and c is the speed of light. Equation (1) contains the four main transport effects: advection, spatial diffusion, drifts, and energy change. DSA is essentially a solution to Equation (1), for the case in which the flow velocity vector changes instantaneously at a shock. This is discussed in Giacalone et al. (2012). In steady state, it is found that the downstream distribution function has a power-law dependence on the momentum with a spectral index that depends only on the ratio of the upstream to downstream plasma density. Moreover, the distribution upstream decreases with distance from the shock, exponentially if the diffusion coefficient does not depend on distance. It is constant downstream of the shock, in the steady-state case. We note from Fig. 13 that from about the year 2000 until the crossing of the termination shock by each Voyager, the intensity of 140-220 keV ions increased until the shock crossing and then became constant downstream, as expected from the theory. Of course, as is now well known, the peak was not at the shock at higher energies (ACRs), which is discussed further below.
For the case of the termination shock, the magnetic field is essentially normal to the unit normal to the shock everywhere over the shock's surface owing to its quasi-spherical shape and because the Parker spiral is highly wound in the outer heliosphere. This is also true at each of the heliographic poles if one considers the behavior of the transverse and radial components of the turbulent magnetic field, as suggested by Jokipii and Kota (1989). They noted that while the radial component of the interplanetary magnetic field behaves as 1/r 2 , the transverse component falls off only as 1/r. Thus, even if the transverse field is small compared to the radial component near the Sun, at the termination shock, nearly 100 AU from the Sun, the transverse field dominates. Directly above the poles, this depends somewhat on the choice of parameters for the turbulent fluctuations, but the affected portion of the termination shock is quite small, less than 0.01% of the surface. Thus, the field is almost entirely in the transverse direction at the termination shock, even at the heliographic poles. In this case, particles gain energy primary by drifting along the shock due to the gradient in the magnetic field, in the direction of the motional electric field. In the absence of particle diffusion, this process is known as shock-drift acceleration (Armstrong et al. 1985). This process is actually a subset of DSA, since drift along the shock, and energy gain through drift along the electric field, is contained within DSA theory.
For the time-dependent case, and also for the case in which there are losses due to geometrical considerations (such a finite-sized shock), the power-law part of the distribution is only at the lowest energies. The characteristic energy where the spectrum deviates from a power-law depends on a number of factors, such as the age of the shock, the size and geometry of the shock, and the rate of acceleration, among other things. For the case of the termination shock, ignoring solar-cycle related effects, the problem is essentially time independent, but the shock is of finite size. In fact, it was emphasized by Jokipii and Giacalone (1998) that the maximum energy is limited by the finite drift distance along the shock. It was noted that the electric potential difference from pole-to-equator is about 240 MeV/charge, which, if fully realized for protons, would give a maximum energy near to the observed maximum ACR energy. The departure from the power law at the shock, and the maximum energy, can also be understood as the energy at which the rate of acceleration at the shock becomes comparable to the rate of energy loss in the expanding solar wind (note that the last term in Equation (1) gives energy loss for the solar wind). Potgieter and Moraal (1988) emphasized the importance of shock curvature, noting that the power-law result was obtained only for a planar shock and the spectral rollover occurs at the energy for which the relevant diffusive length scale is of the same order as the radius of curvature of the shock. The topic of the spectral rollover was also addressed by Steenberg and Moraal (1999).

Acceleration at the Blunt-Shaped Termination Shock
It is now widely accepted that the termination shock is not spherically shaped, even on average, as was assumed in the early models of cosmic ray acceleration and transport, but, instead, has a blunt shape. This shape arises because of the flow of the interstellar plasma, providing a dynamic pressure that is highest in the region where the interstellar flow velocity is opposite to that of the solar wind. This is the so-called "nose" region, and at this point, the termination shock is closer to the Sun than it is at other locations. The blunt-shaped termination shock was first suggested by analysis of TSP anisotropies ).
Although not emphasized earlier, it can also be seen in large-scale MHD simulations of the heliosphere (e.g. Zank 1999).
The blunt-shaped termination shock has a significant and important effect on the acceleration of particles. While the shock normal varies with longitude, the Parker spiral magnetic field is essentially normal to the radial direction. Thus, as was noted by McComas and Schwadron (2006), this leads to the situation in which the intersection point of any given magnetic line of force with the shock moves along the shock from its first crossing, at the nose, towards the flanks of the heliosphere. Since it takes time to accelerate particles, only the lowest-energy particles, accelerated the most rapidly, are expected to peak at this point. It takes up to a year to accelerate to the highest energies, during which time the intersection point will have moved far towards the flanks of the heliosphere. Thus, one expects that at the flanks is where the ACRs will have the highest intensity.
The conceptual, intuitive picture (shown in the left panel of Fig. 2, discussed below) suggested by McComas and Schwadron (2006) has been quantitatively confirmed by Kóta and Jokipii (2008) and . Moreover, it provides a natural explanation for the observations. For instance, at low energies, the TSPs were observed to peak at the shock and subsequently had a relatively constant intensity. In contrast, the highest-energy particles, the ACRs, continued to rise into the heliosheath which is also consistent with the picture above. This is because when the peak fluxes are in the flanks, the particles subsequently are transported within the heliosheath. Along the trajectory of the Voyagers, the peak will not be at the shock, but farther into the heliosheath, as observed (cf. McComas et al. 2019).
Moreover, analysis of ACR anisotropies by , showing a diffusive streaming of ACRs coming from the flanks of the heliosphere support this picture.
It should also be noted that similar effects may arise for the case in which the shock interacts with large-scale plasma turbulence. In this case, the shock surface may become rippled, which can effect the particle acceleration at the termination shock (e.g. Li and Zank 2006). If the scale of the ripples is larger than the characteristic scale associated with the particle diffusion, effects similar to those discussed above for the blunt-shock case arise (e.g. Kóta 2010; Guo et al. 2010). Similarly, even at a planar shock, the turbulent plasma upstream of the shock may lead to meandering magnetic lines of force which also produce similar effects as considered by Kóta and Jokipii (2008). It is likely that these effects also play an important role in the acceleration of heliospheric energetic particles, such as TSPs and ACRs, at the termination shock.

Acceleration at Low Energies: The Injection Problem
We know that ACRs and TSPs originate as interstellar pickup ions whose energy upon their creation is of the order of ∼ 1 keV. It is well accepted that DSA theory is applicable for particles with energies much larger than this, but it is also applicable at lower energies, such as those associated with pickup ions, provided there is sufficient scattering to maintain an isotropic distribution and, also, to remain near the shock. It is now generally known that shocks accelerate low energy particles, even a fraction of those that are part of the thermal distribution, to high energies (Ellison and Eichler 1984;Scholer 1990;Giacalone et al. 1992;Giacalone 2005b). The question arises of how this can happen given that their random velocity is similar to the bulk speed of the plasma which advects them downstream. Since shocks convert energy entering the shock -in the form of dynamic pressure -into thermal energy, the (bulk) plasma is subsonic downstream of the shock. Thus, a fraction of the downstream plasma can return to the shock and enter back into the upstream flow. Thus, even shock-heated thermal plasma is capable of participating in the shock-acceleration process, and at the termination shock, the dominant species is shock-heated interstellar pickup ions. Accelerated particles likely play a critical role in the energy dissipation at shocks, and the efficiency of the acceleration depends on how much of the incident energy at the shock can be proportioned to the high-energy tail.
For the case of interest here, i.e. pickup ion acceleration at the termination shock, one must also consider that the average magnetic field is nearly perpendicular to the shock normal, even for the blunt-shaped shock. It has been previously thought that such shocks have a difficulty to accelerate low-energy ions (e.g. Ellison et al. 1995). However, it has been shown that such shocks are indeed capable of accelerating low (and even thermal) energy particles provided there is a sufficient level of pre-existing magnetic fluctuations present in the plasma that passes through the shock (e.g. Giacalone 2005b). Also, as can clearly be seen in Fig. 13, the peak intensity of 140-220 keV ions, admittedly higher energy than freshly ionized pickup ions, but still much lower than ACR energies, is at the shock. Recently, Lario et al. (2019) analysed observations of energetic particles at interplanetary shocks, including observations which include the thermal, suprathermal, and high-energy tail of the distribution, and found that very low-energy ions are accelerated directly at the shock, even including at least one event in which the shock-normal angle was close to ninety degrees. Neergaard- Parker and Zank (2012) and Neergaard-Parker et al. (2014) quantified this further by assuming a form for source distribution of particles in diffusive shock acceleration, and fit the shock-accelerated spectrum to that observed for a number of interplanetary shocks. They found that in many cases, acceleration of solar-wind, at near thermal energies, accounted for the intensity of the observed high-energy tail. Thus, observations of interplanetary shocks, as well as that at the termination shock, generally reveal that there is no problem accelerating low energy ions at quasi-perpendicular shocks. Thus, it is reasonable to say that it is a misconception to suggest otherwise.
Hybrid simulations of the interaction of pickup ions with the termination shock, showing efficient acceleration, were performed initially by Giacalone and Decker (2010), and more recently by Giacalone et al. (2021). The simulations treat all the ions (solar wind and pickup) kinetically and effectively solve for the kinetic shock physics on ion scales. It was found that the pickup ions are thermalized at the shock, and a high-energy tail forms from the shockheated distribution. This high-energy tail has a nearly power-law dependence on energy up to about 50-100 keV, where the simulated spectrum falls off sharply. This rapid fall off of the spectrum is caused by the limitations imposed by the very small simulation domain and finite simulation time. By comparing the results with V2/LECP observations, and by using simulation parameters consistent with the observed plasma and field parameters, it was found that at about 50 keV, the flux of accelerated pickup ions agreed well with the observed TSP intensities downstream of the shock. This suggests strongly that TSPs are, in fact, accelerated interstellar pickup ions, representing a key step in the process which ultimately leads to ACRs. More-recently, Zirnstein et al. (2021) performed a similar study using a test-particle approach of particle acceleration at the termination shock in the presence of magnetic turbulence, consistent with that observed by Voyager 2. The results were compared to IBEX observations of energetic neutral atoms. It was found that a high-energy, power-law tail formed from the pickup ion distribution due to shock acceleration, with a spectral index consistent with that inferred from IBEX observations.
Recently, Giacalone et al. (2021) extended the study of Giacalone and Decker (2010) to study the acceleration of pickup ions at other locations along the termination shock in addition to that where V2 crossed. The right panel of Fig. 2 is taken from this paper. The left panel shows a representation of blunt-shaped termination shock (thick black line) and representative Parker-spiral magnetic field lines (gray), and a color-code representation of the ACR intensity with red being the highest and yellow the lowest. Three locations along the shock are identified in the left panel with purple, green, and blue circle symbols. This figure is an adaptation of the figure shown in McComas and Schwadron (2006). The locations identified in this figure relate to the hybrid simulated spectra, in the same colors, in the right panel. In the right panel, solid lines are hybrid-simulation pickup ion differential intensity spectra as a function of energy, and dashed lines are those for the solar wind. Figure 2 also shows V2/LECP data just after the crossing of the termination shock (black square symbols). The simulated spectra at about 50-100 keV agree well with the observations. Above this energy, the simulated spectra fall off dramatically, which is due to limits imposed by the use of the small simulation domain and finite simulation time, and are not comparable to the observations at these energies. As noted above, the fact that the simulated spectra match the observations at 50-100 keV indicates that TSPs are accelerated pickup ions. This figure also indicates that there is little variation of the 50-100 keV intensity at widely different locations along the termination shock. Given the TSPs were observed to be rather uniform throughout the heliosheath, prior to their going away completely at the heliopause, and that they are uniform along the shock as well, suggests that TSPs are likely quite uniform throughout the heliosheath. For more details, see Giacalone et al. (2021).
With regards to the physics of particle acceleration at low energis, the primary concern is how particles are trapped near the shock. Thermal particles, or those just slightly more energetic than the bulk plasma, are largely dominated by convection. However, as the plasma is heated across the shock, some of these particles are capable of returning upstream. For  Giacalone et al. (2021). The black squares in the right figure are Voyager 2 LECP observations taken just after the crossing of the termination shock (see also Giacalone and Decker 2010, for more details) shocks that move normal to the magnetic field, the particles must move across the mean magnetic field in order to stay near the shock. This has led to the suggestion that nearlyperpendicular shocks have difficulty in accelerating low-energy particles. However, magnetic field-line meandering, caused by large-scale turbulence in the solar wind, is one means by which low-energy ions, even those near the thermal population, can remain near the shock and be efficiently accelerated (cf. Giacalone 2005a,b). Is has also been suggested that the required process for the injection of particles into DSA can also be related to the reflection of ions at the shock (Chalov and Fahr 1996;Dworsky and Fahr 2000;le Roux et al. 2000).  and  investigated ion reflections from the shock due to pitch-angle scattering driven by Alfvènic turbulence and demonstrated that suprathermal pick-up ions are not necessary, but that the injection already works for the energetic part of the regular thermal solar wind ions. This process works especially effectively for parallel MHD shock configurations where about 18% of the incoming solar wind ions are reflected back to the upstream side and do enter the ACR acceleration process. This pitch-angle-induced ion reflection could in principle continue working even further downstream of the shock since Alfvénic turbulence amplitudes grow downstream (Zank et al., this journal).

Acceleration Mechanisms in the Heliosheath
The topic of the acceleration of cosmic rays in the inner heliosheath received attention after the Voyager spacecraft crossed the termination shock. As noted above, they found that only the intensities of low-energy particles peaked at the shock, while those of the ACRs increased further in the downstream region and that their spectra exhibited multiple powerlaw spectral slopes (Decker et al. 2005;. At the time, this was not well understood since classical DSA models were not readily applicable (le Roux and Fichtner 1997).
Explanations for this unexpected behaviour include refined DSA models (as discussed above) (McComas and Schwadron 2006;le Roux and Webb 2009;Giacalone and Decker 2010;Senanayake et al. 2015), second-order Fermi processes such as momentum diffusion (Moraal et al. 2006;Ferreira et al. 2007;Fisk and Gloeckler 2009;Kallenbach et al. 2009;Strauss et al. 2010a), and magnetic reconnection (Lazarian and Opher 2009;Drake et al. 2010;Zhao et al. 2019). A discussion of these mechanisms in the context of ACRs was included in the comprehensive review by Giacalone et al. (2012). While, so far, there is still some debate regarding the acceleration mechanism of ACRs, there is agreement with regard to the seed population, which is the pick-up ion population.
A further acceleration of the seed population of ACRs in the inner heliosheath is also of interest in the context of energetic neutral atoms (ENAs). The main ENA component observed by the IBEX spacecraft is produced via charge exchange of solar wind and pickup protons with interstellar neutral atoms in this region (e.g., Siewert et al. 2013). ENAs 1 keV-6 keV need a pre-acceleration of protons, which is usually assumed to be momentum diffusion (Fahr et al. 2016;Zirnstein et al. 2018). For such modeling, kappa distributions have been considered, which are discussed in the chapter by Perri et al. in this book. For ENAs with higher energies, as observed with SoHO/HSTOF, charge exchange with ACRs has been considered (Czechowski et al. 2001) but without taking into account their acceleration in the inner heliosheath. This modeling is, therefore, left as a task for the future. Strauss et al. (2010b) developed a numerical model for the acceleration and propagation of anomalous Oxygen in the heliosheath. Apart from diffusive shock acceleration at the termination shock (TS), their model also included momentum diffusion and adiabatic energy changes in the heliosheath (see also Langner et al. 2006;Ferreira et al. 2007). The TS compression ratio and particle injection efficiency were obtained from appropriate hydrodynamic models (Scherer et al. 2006). Their model satisfactorily reproduced the modulated spectrum of anomalous Oxygen observed at the TS and the unfolding into the heliosheath. They concluded that a combination of momentum diffusion and adiabatic energy increase, under reasonable assumptions of the solar wind speed in the heliosheath, produce a feasible continuous acceleration process to explain the main features of the full anomalous Oxygen spectrum in the heliosheath. Strauss et al. (2010a) emphasized that the observed anomalous Oxygen spectrum near the TS reveals a power law form above the energy roll-over point, instead of an exponential spectral cut-off at energies > 20 MeV/nuc (Webber et al. 2007). Using the numerical model mentioned above, they showed that this deviation from an expected exponential form could be explained taking into account the acceleration of multiply-charged anomalous Oxygen up to 70 MeV/nuc above which Galactic Oxygen dominates. They illustrated how the roll-over energy shifts to higher energies when the higher charge states for anomalous Oxygen were incorporated into their simulations, while the low energy portion of these spectra (1-10 MeV/nuc) remains relatively unaltered. This causes in the process an increase in high energy anomalous Oxygen intensities throughout the heliosphere, also at the Earth. This was displayed in the modulated spectra at the Earth in comparison with relevant observations near the Earth from several space missions (Marsden et al. 1999;Hill et al. 2002;Cummings et al. 2009). The topic of multiply-charged anomalous cosmic rays was discussed previously in this chapter (see Jokipii 1996;Mewaldt et al. 1996;. Strauss and Potgieter (2010) modeled ACR Oxygen spectra and associated spatial gradients applicable to the inner heliosphere. They specifically simulated latitudinal and radial gradients during both heliospheric magnetic field polarity cycles and compared their simulations with observations close to and at the Earth; see their Table 1 for a list of observed and radial gradients. Their global simulation of spatial gradients for anomalous Oxygen is depicted in Fig. 3. This model was a refinement of the numerical model of Strauss et al. (2010b), applicable to the main features of this type of charged particle population which was not addressed earlier.

Transport Within the Heliosheath and Escape from the Heliosphere
The V1 and V2 heliopause (HP) encounters revealed that the magnetic boundary of the solar system was not impervious to HEPs (heliospheric energetic particles, as discussed in Sect. 8). Viewed on a global scale, the HP appeared to be essentially a tangential discontinuity characterized by an approximately twofold increase in the magnetic field strength without a significant change in direction. The situation was much more complex on short (< AU) scales. V1 encountered two narrow regions within the last 0.3 AU before the heliopause where the magnetic field was stronger, and the ACR intensity was depressed by a factor of 3-5 compared with the surrounding heliosheath (Krimigis et al. 2013). The intensities dropped again at the HP itself. The decay time of these particles was rather different at the V1 and V2 crossings. In the case of the V1 crossing, the ACRs were seen to deplete over a time period of several days, but at V2 ACRs persisted for a month or more after the crossing. The extent of this boundary layer was found to be weakly dependent on particle's velocity, but strongly rigidity-dependent ). In addition, ions streaming along the field disappeared first, while those gyrating near the 90 • pitch angle persisted the longest. As a result, within the dropouts and the LISM boundary layer, ACRs developed a "pancake" distribution with a large second-order anisotropy.
The dropouts were interpreted by Krimigis et al. (2013) as interstellar magnetic flux tubes penetrating the HP, into the heliosheath, as a consequence of the magnetic interchange instability. The stability criteria were derived theoretically by Florinski (2015). The instability is driven by a force exerted on the magnetic field by ions moving along curved trajectories corresponding to the convexity of the HP, and requires a negative pressure gradient across the boundary. Because interchange is inhibited by a degree of magnetic rotation across the HP, it was suggested that a flux tube would undergo a kink and thus become aligned with the field on the opposite side. Florinski et al. (2013) simulated the anisotropies of 5 MeV protons near the HP by treating it as a shear layer with a large decrease in the amplitudes of magnetic fluctuations from the heliospheric to the interstellar side. They assumed that magnetic shear suppressed field line meandering across the HP, so that only particles with large pitch angle were able to cross into the VLISM (very local interstellar medium) by virtue of their large gyro-radius, whereas particles streaming parallel to the local field remained trapped inside the heliosheath. The model successfully explained the observed ∼ 5 MeV ion spatial profiles and pitch-angle distributions reported by Krimigis et al. (2013). These ideas were further developed by  who investigated the rigidity dependence of the ACR depletion layer. They proposed that gradient magnetic-field drift was responsible for spatially separating particles according to their rigidity because the gradient drift velocity is proportional to the gyroradius. This mechanism requires an increase in the field strength toward the south, which is expected based on the magnetic field draping pattern from large-scale MHD simulations (e.g., Izmodenov and Alexashov 2020). Figure 4 compares simulated intensities of 5-MeV protons at two different pitch angles. Note, in particular, a more gradual decrease in particles gyrating near the 90 • pitch angle compared with those streaming parallel to the field lines. An alternative to this model was developed by Strauss and Fichtner (2014), who suggested that perpendicular diffusion was large at 90 • and decreased to zero as the pitch-angle cosine, μ → ±1. Subsequently, a theory was developed in support of this scenario. The theory predicted that perpendicular diffusion in the presence of compressive longitudinal fluctuations would indeed have the required dependence on pitch angle (Strauss et al. 2016).
Measuring the anisotropy of ACRs associated with the gradient in the density of their guiding centers, that was very prominent within the dropouts, allowed  to derive an independent estimate of the radial component of the plasma velocity in the heliosheath, on the heliospheric side of the HP. This observation was possible owing to the relative orientation of the multiple Voyager LET telescopes. The value for the radial speed obtained was consistent with zero or slightly positive, but under 5 km/s. This result is consistent with the notion that the HP is a tangential discontinuity, impervious to the low-energy solar-wind ions. Nonetheless, a sizable fraction of the higher energy ions were evidently able to overcome the magnetic shear layer and escape into the VLISM. Note that Gloeckler and Fisk (2016) concluded that the sum of the pressures of the solar wind, pickup ions, ACRs, plus the magnetic pressure on the heliospheric side of the HP exceeds the combined pressures of the interstellar plasma and magnetic field, which appears to contradict the notion of a pressurebalanced boundary. This difficulty is naturally resolved by adding the escaped heliospheric ions to the pressure balance on the interstellar side. Guo et al. (2018Guo et al. ( , 2019 investigated the effects of this pressure loss on the size of the heliosphere. In their model, the energetic particles were treated as a massless fluid whose energy density E c satisfied the equation (e.g., Zank et al. 1993) where p c is the energetic particle pressure, u is the plasma velocity, κ is the energy-averaged diffusion coefficient, and α is the rate of injection from the low-energy pickup ion population that is not distinguished from the thermal plasma. It was found that the ACRs play a significant role in the structure of the outer heliosphere in some situations, being largely dependent on the momentum and energy transfer from the solar wind plasma to the ACRs, the diffusive shock acceleration at the termination shock, and the subsequent loss of ACRs across the heliopause and their rapid escape into the interstellar medium. Under favorable conditions, characterized by a large fraction of energy conversion α and a high enough diffusion coefficient κ in the solar wind, the ACRs were found to reduce the width of the heliosheath by up to ∼ 18 AU, as seen in Fig. 5. Consequently, these results indicate that the effect of ACRs is a potential key factor for the global structure of the outer heliosphere in a numerical model that could partially explain the thickness of the heliosheath based on the timing of the HP and termination shock encounters of the two Voyager probes.

Solar Cycle Dependence
Observations of both ACRs and GCRs near Earth span multiple solar cycles. By 1 AU, both ACR and GCR intensities have been heavily modulated because of their transport in the solar   Mewaldt et al. 1993 for data references). Figure updated from Leske et al. (2013) wind, by an amount which varies greatly with the phase, and to some extent, the magnetic polarity of the solar cycle. As seen in Fig. 6, during periods of low sunspot number, solar minimum, and when the solar magnetic field points outward in the northern hemisphere (A > 0), the intensities of positively-charged ACR and GCR ions tend to plateau for several years, creating a "flat-topped" profile. For the alternate magnetic polarity (still during solar minimum), when the field points outward in the southern hemisphere (A < 0), the highest intensities persist for a much briefer period. This behavior is understood to arise from the fact that during A < 0 solar minima, particles drift into the inner heliosphere along the heliospheric current sheet (HCS), while under A > 0 conditions particles enter from the polar regions of the heliosphere and drift outward along the HCS (Jokipii and Thomas 1981;Kota and Jokipii 1983). The result is that in addition to the well-known 11-year cosmic-ray cycle pattern, there is also a 22-year cosmic ray cycle which results from the swapping of the solar magnetic polarity with consecutive solar minima. During the 2008-2009 A < 0 solar minimum period, cosmic-ray modulation was significantly less and GCR intensities were correspondingly higher than were seen in earlier cycles (Mewaldt et al. 2010); however, the ACR intensities remained below their 1996 values (Leske et al. 2013). The solar wind density, dynamic pressure and magnetic field strength were all quite weak Smith and Balogh 2008). Yet, the HCS, which Fig. 7 Energy spectra from the SIS (at lower energies) and CRIS (at higher energies) instruments on ACE for elements N, O, and Ne with strong ACR components (left panels) and C, Si, and Fe that are essentially pure GCRs at 1 AU (right panels). For each element, spectra are shown in two periods, one near the 1997 A > 0 solar minimum (lighter dashed curves) and one at the 2018-2020 A > 0 solar minimum (darker solid curves), with ratios of the spectra at the two time periods shown in the bottom panels. All GCRs are seen to be significantly enhanced by a similar amount in the latest solar minimum relative to the earlier minimum throughout this energy range, while ACRs are somewhat depleted. For a similar comparison of spectra in the 2009 A < 0 minimum with those in the 1997 A > 0 period, see Leske et al. (2013) might be expected to be quite flat during the weakest part of the solar cycle, remained tilted. It only became "flattened" in 2009. This lead to speculation that since the ACR source intensity in the outer heliosphere is expected to be stronger at low heliographic latitudes during A < 0 minima periods (Jokipii 1986;Florinski et al. 2004), while the GCRs are more uniformly distributed at the boundary of the heliosphere, that perhaps ACR intensities at 1 AU in near ecliptic plane (at Earth) were more sensitive to the higher HCS tilt angle than GCRs and thus not as greatly enhanced (Leske et al. 2013). Effects of the HCS tilt angle would be expected to be less important during the present A > 0 cycle, however.
Although GCR intensities in the last two solar minima have been higher than previously observed, the same is not the case for the ACR intensities. Note from Fig. 6 that the ACR intensities are seen to return to similar peak intensities at successive solar maxima. Differences between the ACR and GCR behavior for different ion species and energies in the last two A > 0 solar minima are seen in more detail in the spectra shown in Fig. 7. The intensities of all GCR elements from C to Fe were enhanced by ∼ 30% in 2020 over their values in 1997 at energies from at least 50 to over 200 MeV/nucleon. For elements that do not have a strong ACR component, such as those in the right panels of Fig. 7, this enhancement extends all the way down to 10 MeV/nucleon. In contrast, the ACR elements N, O, and Ne all show no change or even a slight depletion in their 2020 peak intensities compared with their 1997 values, as seen in the lower left panel of Fig. 7. This is clearly not simply an energy-dependent effect on the spectra since it is not present in GCRs (right panel at the same energies), nor can it be attributed to instrumental effects since the ACRs and GCRs were measured with the same instrument.
Both ACRs and GCRs arrive at 1 AU from the outer heliosphere and are subject to similar transport, and subsequent modulation, since they move in the same field and plasma. They generally have the same waxing and waning with the solar cycle, as indicated by the similarities in their intensity variations as the solar cycle progresses (e.g. Fig. 6). However, the ACRs are "locally" accelerated, as described in the previous sections of this review, while the GCRs are not, and come from far away from the heliosphere. The GCR behavior shows that there must be less modulation in the current solar minimum period compared to the last A > 0 minimum, but the ACR intensities seem rather unaffected. This is not presently well understood. On the one hand, the GCR intensity varied considerably across the HP (see Rankin, this journal), suggesting that perhaps much of the GCR modulation occurs beyond the source of the ACRs. On the other hand, the production of ACRs depends on their acceleration rate, which depends on the solar magnetic field, shock strength, and other factors, and these likely play a critical role. This is discussed further below. Moreover, the low-energy source of the ACRs, interstellar pickup ions, may have a lower intensity due to reduction in the ionization rate of inflowing interstellar neutral atoms because of the reduced solar output in the recent weak solar maxima (Sokół et al. 2019). At present, there is no widely accepted explanation for the observed differences in the ACR and GCR intensity variations (at Earth) with the solar cycle.
In addition to modifying the ACR intensity, changes in the ACR source might alter their composition and/or energy spectra. For example, the production of pick-up ions depends on both charge exchange collisions with the solar wind and photoionization by solar UV radiation. To ionize neutral Ne, photoionization is about 100 times more significant than charge exchange, while for oxygen the two processes contribute nearly equally (Sokół et al. 2019). Thus, if the solar wind density and speed, which governs charge exchange, varies by a different amount than the solar UV emission that determines the photoionization rate, the elemental composition of the pick-up ions and hence ACRs would be expected to change. Alternatively, if the ACR acceleration process is mass-dependent, changes in the plasma and fields near where ACRs are accelerated might affect the resulting composition of the ACRs. Similarly, reduced acceleration efficiency might preferentially deplete the higher energies, resulting in a softer spectrum. Changes in diffusion might alter the production of higher ionic charge states via stripping of the accelerated particles (Jokipii 1996), again changing the measured elemental spectra since the ions with higher charge states are observed at higher energies (Mewaldt et al. 1996;Klecker et al. 1997;Selesnick et al. 1997). These hypothetical composition and spectral changes would be much more subtle and harder to detect at 1 AU than the differences in the overall ACR intensity, but perhaps limits to these effects could be obtained from the existing data to help constrain models.
A reasonable explanation for the observed discrepancy in ACR-GCR solar-cycle variations is the potentially less effective acceleration of ACRs at the TS, as noted above. ACRs reaching Earth must originate from higher-energy ACRs at the TS because they lose energy during their transport to Earth; thus the observed ACR flux at Earth is sensitive to the flux at the high-energy part of the source spectrum at the TS. Moraal and Stoker (2010) were the first to point out that while the likely larger diffusion coefficient responsible for milder (or weaker) modulation will, at the same time, also lead to a slower rate of acceleration and a lower cut-off energy of the ACR source spectrum. This is because the rate of acceleration is inversely proportional to the diffusion coefficient κ (e.g. Drury 1983). Note also that it can be shown that for the simplest spherical configuration, the powerlaw exponent associated with acceleration of particles at the spherical TS, takes the form: The first term is the result from the usual planar-shock case, while the second term arises because of the spherical geometry, and inclusion of adiabatic cooling. The net result is that the spectrum steepens if either κ increases or the shock radius, R sh decreases. The later may be related to the effects of varying solar-wind dynamic pressure.
Finally, the weaker heliospheric magnetic field results in faster particle drift along the shock front, which again leads to less effective acceleration. ACRs drift along the shock face as they are accelerated, hence the maximum energy of accelerated ACRs should be in the range of the electric potential difference between the heliospheric equator and pole, as we discussed previously (e.g. Jokipii 1996). This electric potential difference is determined by the 26 day rotational velocity of the Sun multiplied by the total open flux of the solar magnetic field, hence a weaker field should lower the maximum energy of ACRs.
At present we cannot single out the precise cause of the difference in solar-cycle behavior of the ACRs and GCRs. The mechanisms discussed above can, and likely do, act together, the exploration of their respective contribution needs further research.

Radial and Latitudinal ACR Gradients in the Heliosphere
An key prediction from cosmic-ray transport models is that radial and latitudinal gradients of ACRs depend strongly on the solar magnetic cycle. The simulated radial gradients using the model described in the previous section at 2, 5, 10 and 20 AU, and a polar angle of 45 degrees, for anomalous Oxygen for A > 0 and A < 0 magnetic cycles are shown in Fig. 8. For comparison, observations of non-local radial gradient are shown as indicated. The value of 45%/AU for the radial gradient of anomalous Oxygen close and inside 1 AU were found by Marquardt et al. (2018) to be qualitatively consistent with the predictions from the model.
Observational evidence for positive and negative latitudinal gradients came from the two Pioneer and Voyager missions in the outer heliosphere (McKibben et al. 1979;Cummings et al. 1987;Christon et al. 1986). McKibben et al. (1979) and Cummings et al. Fig. 9 Overview of latitudinal gradients observed during the first and second Ulysses fast latitude scan (1987) reported values for the latitudinal gradients ranging from 2.1 to 3.1%/ • and −2.2 to −1.6%/ • for ACR Helium at about 15 MeV/nuc. during an A > 0-and A < 0-solar magnetic cycles, respectively. For ACR oxygen Cummings et al. (1987) found even larger values for the latitudinal gradient (−3.7 to −2.9%/ • ). However their analysis was restricted to heliographic latitudes of less than 30 • .
The Ulysses spacecraft mission provided instruments (Simpson et al. 1992;Simpson et al. 1995;Lanzerotti and Maclennan 1995;Keppler et al. 1992) with which to measure cosmic ray latitudinal gradients in the heliographic polar regions. Trattner et al. (1995bTrattner et al. ( ,a, 1996, Lanzerotti and Maclennan (1995) and Heber et al. (1999) report latitudinal gradients varying between 0.39%/ • to a 2.12%/ • in an A > 0-solar magnetic epoch. Figure 9, which comes from Heber (2001), displays the latitudinal gradient as a function of rigidity for both ACRs and GCRs during Ulysses' first latitudinal scan. It was also found that the ACR oxygen gradient during the A < 0 magnetic cycle was 5 times smaller than during the A > 0-solar magnetic cycle (Cummings et al. 2009). These authors concluded that ACRs do not propagate into the inner heliosphere by drifting along the HCS. In contrast, latitudinal gradients in the outer heliosphere during the A < 0-solar magnetic epoch are generally considerably larger, indicating that the transport of ACRs is dominated in the inner heliosphere by diffusion, rather than drift along the HCS. Simpson et al. (1996) also reported a north-south asymmetry in the cosmic rays observations from Ulysses after it had observed both polar regions of the Sun. They noted that while the latitudinal gradients in each hemisphere were quite similar in each hemisphere, and the modulation was approximately symmetric, the plane of symmetry was offset southward, by about 10 degrees from the heliographic equator. To our knowledge, this has not been fully explained.

Observations of Anomalous Cosmic Rays Inside of 1 AU
The study of cosmic rays in the near-Sun environment is a territory that is just beginning to be explored. The first measurements from ∼ 0.3 to 1.0 AU were taken by the Helios spacecraft over 45 years ago, with no follow-up until Parker Solar Probe. Re-visiting the Helios data, Marquardt et al. (2018) investigated the spatial distribution of ACR oxygen during the 1974-1977 solar minimum and found a significant radial gradient: 48 ± 12% AU −1 for energies of 9 to 29 MeV/nuc, roughly 3 times larger than measured by Pioneer 10 from 1 to 10 AU around the same time period. For instance, Webber et al. (1981) found values of the SIS and CRIS instruments on ACE (light and dark green triangles, respectively) were used to correct for modulation due to long-term variations in solar output. Data were taken during the ∼ 2018.7 to ∼ 2019.9 and ∼ 2018.7 to ∼ 2021.2 time-frames, respectively. GCRs at 1 AU were simulated using the HelMod online calculator (www.helmod.org; see text for further details). Figure adapted from Rankin et al. (2021) (Helium) and Rankin et al. (in press) (Oxygen) 15 ± 3% AU −1 from 1 to 10 AU for 9.5 to 24 MeV/nuc. These findings were not well reproduced by models. For example the model of Strauss and Potgieter (2010) for ACR oxygen in the heliosphere, including close to 1 AU, predicted the largest radial gradient is about 2 AU beyond the orbit of Earth, suggesting a smaller gradient inside 1 AU.
The first studies of cosmic rays on Parker Solar Probe were conducted by Rankin et al. (2021) and Rankin et al. (2022) who analyzed ACR helium and oxygen into 0.17 and 0.096 AU, respectively, using measurements made by the EPI-Hi sensor of the ISOIS instrument suite (McComas et al. 2016). The resulting event-subtracted spectra are shown in Fig. 10. Like Marquardt et al. (2018), they found that the gradients inside 1 AU were much larger than those at 1 AU and beyond. For example, Rankin et al. (2021) reported a helium radial intensity gradient of 25 ± 5% AU −1 in the ∼ 4 to ∼ 45 MeV/nuc energy range, vs. ∼ 10% AU −1 in the inner heliosphere under similar solar cycle conditions (see their Table 1). They also determined the gradient's magnitude for specific discrete energy intervals, as shown in Fig. 11. Results were acquired by performing log fits in linear space, assuming a differential radial gradient of the form: Distinguishing amongst the different populations (i.e., SEPs, ACRs, and GCRs) is challenging. To differentiate between SEPs and cosmic rays, most studies utilize data taken during times of low, or quiet solar activity, ideally during sunspot minimum when cosmic ray fluxes are high. Distinguishing ACRs from GCRs often proves more challenging as many cosmic ray instruments do not distinguish amongst charge states. When their spectra are well defined, as evidenced in Oxygen (Fig. 10b), the two populations can be readily separated by energy. However, for species such as Helium (Fig. 10a), the distinction is not as clear. These results reflect the combination of ACRs and GCRs over their stated energy ranges. Rankin et al. (2021) additionally evaluated the behavior of the gradients of ACRs after GCR subtraction. Using the HELMOD model at 1 AU (Bobik et al. 2012;Boschini et al. 2020, version 4.0.1, 2021www.helmod.org) to approximate the near-Sun behavior of GCRs, they Fig. 11 Radial gradients of PSP/ISOIS helium as a function of energy. Data were taken from ∼ 2018.7 to ∼ 2019.9 using the LET1 and HET B telescopes on the EPI-Hi sensor. After subtracting solar energetic particle events and de-trending from long-term changes in solar modulation (using SOHO/EPHIN as a 1-AU baseline; see Fig. 10a), fits were applied to 3 forms of data: (i) averaged over Carrington longitudes in the spacecraft frame (Carrington Ave.), (ii) averaged over daily time scales (Daily Ave.), and (iii) averaged in 0.01 AU radial increments (0.01-AU bins). The results for (ii) and (iii) are mostly consistent with each other, while a possibly larger gradient is evident in (i), for which variations caused by solar wind streams have been minimized. The fluxes used in this analysis reflect the combined contribution from all sources at a given energy (e.g., ACRs dominate but the GCRs have not yet been subtracted). Figure from Rankin et al. (2021) reported a dramatic increase in the magnitude of the ACR gradient: from 25.2 ± 4.1% AU −1 to 34.3 ± 5.6% AU −1 in the 4.0 to 32 MeV/nuc energy range, and 26.3 ± 5.8% AU −1 to 44.7 ± 10.2% AU −1 in the 13 to 45 MeV/nuc range, after GCR background subtraction.
The ACR Oxygen radial gradient reported recently by Rankin et al. (in press) was even larger in magnitude: 49.4 ± 8.0% AU −1 (6.7 to 27 MeV nuc −1 ; GCR's subtracted and averaged in 0.025-AU radial bins), and showed surprising agreement with the comparableenergy value obtained by Helios ∼ 45 years ago (Marquardt et al. 2018), indicative of a constant radial gradient from 0.1 to 1 AU. For a summary of ACR Oxygen radial gradients measured over the past 5 solar cycles, the reader is referred to Fig. 4 from Rankin et al. (2022). Figure 12 presents a comparison of the observed ACR Oxygen gradients during multiple solar minima to results produced by an updated version of the model by Strauss and Potgieter (2010). Beyond ∼ 2 AU, the model and observations show strong agreement. However, inside ∼ 2 AU, they start to diverge. This was noted above. While the model predicts a significant decrease toward the Sun, following the typical r −2 Parker Spiral behavior, the observations appear more consistent with ∼ r −1 (as seen in blue dotted line in the left panel of Fig. 12). This and the remarkable consistency between Parker Solar Probe and HELIOS data led Rankin et al. (2022) to suggest a potential explanation: transverse magnetic fluctuations, whose intensity decays as ∼ r −1 , could potentially overwhelm the radial component near the Sun, which decays according to ∼ r −2 , and alter the structure of the magnetic field at scales from ∼ 0.1 to several AU, in such a way as to significantly influence the transport of cosmic rays (e.g. Fisk and Schwadron 1995;Jokipii 2001).
Overall, the new ACR data inside 1 AU presents an intriguing challenge for theorists and modelers investigating particle transport near the Sun. We suggest the interested reader also see Englebrecht et al. (this journal), for more details and recent advances in ACR, GCR, and SEP transport modeling.

Fig. 12
Comparison of simulated and observed radial gradients of ACR Oxygen for qA > 0 (left) and qA < 0 (right) polarity cycles. Model results from Strauss and Potgieter (2010) are included, demonstrating clear challenges for reproducing measured gradients inside ∼ 2 AU. Interestingly, results near the Sun reported by Rankin et al. (2021) and Marquardt et al. (2018) appear more consistent with a r −1 behavior rather than that of the standard Parker Spiral assumed here. Figure from Rankin et al. (2022) 8 Is the Term "Anomalous" Still Appropriate?
Given that we have learned much about ACRs, it is reasonable to address the question of whether the term "Anomalous Cosmic Ray" is useful. On the one hand, ACRs has become a well known acronym in the heliophysics community, and anomalous cosmic rays are described in many articles and textbooks. Thus, it might cause considerable confusion to introduce a new term. However, to those not familiar with the physics of the outer heliosphere, the term "anomalous" is itself confusing, and not particularly descriptive on its own. By its strict definition, it implies that they are not understood or different from expectations. On the other hand, it is widely accepted that ACRs are interstellar pickup ions that are accelerated in the outer heliosphere. So, they are clearly not anomalous in this context. Moreover, the term "cosmic ray" could cause a bit of confusion to the non-expert as well. Most definitions of this include that the particle move at near the speed of light. However, ACRs generally move only a fraction of the speed of light (∼ 0.4c for a 300 MeV proton). Moreover, ACRs are not even the most-energetic species produced in the solar system since the Sun occasionally creates particles with energies of a few GeV, which is more than that of ACRs. The term "solar cosmic ray" was in use some decades ago, but is not as commonly used today.
It should be noted that despite all we have learned about ACRs over the past (nearly) 50 years, there remain important puzzles, particularly concerning the mechanism of acceleration, that are useful to discuss with regards to a consideration of their designation. Prior to the Voyager crossings of the termination shock in 2004 (Voyager 1) and 2007 (Voyager 2), pickup-ion acceleration at the solar-wind termination shock was the paradigm for ACRs. Observations by interplanetary spacecraft prior to this agreed well with the predictions of the theory. The paradigm was challenged when the ACR spectrum did not "unfold", as expected, when these spacecraft crossed the termination shock (Stone et al. 2005. However, it is important to note that charged particles with energies < 1 MeV did indeed reach  Krimigis and Decker (2015), but only panel B made it to the final publication. It is reproduced here with permission from R. Decker (private communication) a maximum near the time of the termination shock crossing (Decker et al. 2005(Decker et al. , 2008. In fact, they were also observed to be enhanced well before the shock crossing, as can be seen in Fig. 2. Panel (B) of this figure shows the flux of 140-220 keV ions observed by the Voyager LECP experiment for the entire mission since the launch of each spacecraft. Voyager 1 (V1) is in black, and Voyager 2 (V2) in red. The top panel (A) shows the radial location and latitude of each spacecraft, panel (C) shows the flux of > 270 MeV ions, and the bottom panel (D) shows the average sunspot number. The decline in the intensity of 140-220 keV ions after launch is very noteworthy. The dominant source of these particles, during this period, are solar-related energetic particles and their intensity decreases with distance from the Sun. Starting in about 2002, the flux of particles in this energy range began to rise again, reaching a peak at the crossing of the termination shock in late 2004 by Voyager 1 and mid 2007 by Voyager 2, as can be seen in this figure. These particles are known as Termination Shock Particles (TSPs), and are the are likely progenitors of ACRs, as are interstellar pickup ions.
A potentially useful analogy is the acronym "SEP", for solar energetic particles. SEP is a general name given to energetic particles related to processes occurring at the Sun, including those associated with solar flares, coronal mass ejections (CME), corotating (or stream) interaction regions, etc. In the case of CME-related SEP events, which are among the most intense and reach the highest energies, it is well understood that the shock driven by the CME accelerates the particles. However, only at lower energies is it common to see an intensity peak at the time of arrival of the shock. At energies above a few MeV, the peak is usually well before the arrival of the shock. Scientists call the ones that peak at the shock ESPs (Energetic Storm Particles). These particles are locally accelerated. We do not have a commonly used name for the high-energy ones that do not peak at the shock when it crosses 1 AU. Yet, they are all called SEPs. This is very analogous to the case of acceleration at the termination shock, as described above. The high-energy part, what we call ACRs, did not peak at the shock, but the low-energy part, what we call TSPs, did. Thus, it is reasonable to refer to them all as heliospheric energetic particles, or HEPs. We note that this suggestion also appeared in the article by Wimmer-Schweingruber and Bochsler (2001).
The acronym HEP is to be understood as a general term which includes anomalous cosmic rays and termination shock particles, and possibly others. More specifically, they are energetic particles accelerated in the heliosphere, away from the Sun, and are understood to be distinctly different from SEPs. Interstellar pickup ions could reasonably be included with this general term because they represent a suprathermal population, distinctly different from the solar wind thermal population. Neutral species might also be considered in the general definition since the term "particle" does not restrict it to ionized species. In fact, the Sun produces both energetic neutral atoms and energetic neutrons, which are sometimes generally referred to as SEPs. By extension, one might also consider energetic neutral atoms produced by the processes in the heliosphere to be also included in the general HEP term.

Summary
We have presented a review of selected topics in the observation, theory, and modeling of anomalous cosmic rays. We reviewed a number of topical subjects, focusing on recent advances in our understanding of the origin and distribution of these particles in the heliosphere, and even their escape. Certainly we were unable to cover all aspects of this topic; however, we have a very extensive list of references at the end of this chapter, which provides the interested reader considerably more sources of information.
We also raised the question of whether the term "anomalous" is still appropriate in describing these particles. Although it is challenging to introduce a new term once it is well established, as ACRs are in this case, we noted that we have learned much about these particles to the point that this term is no longer descriptive or useful, particularly to the broader scientific community. We introduced the term "heliospheric energetic particles" (HEP), of which ACRs are a subset. HEPs represent any suprathermal particle species which originate in the heliosphere, away from the Sun. The latter qualifier is to distinguish them from solar energetic particles, or SEPs. HEPs might reasonably include interstellar pickup ions, in addition to termination shock particles, and ACRs, but exclude GCRs which originate far from the heliosphere. ACRs are the high-energy component of HEPs, while termination shock particles, and interstellar pickup ions are lower-energy components. We suggest that this more-general term is more appropriate, and descriptive to a broader audience, given what we have learned about energetic particles in the heliosphere from nearly five decades of research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.