Strong Preferential Ion Heating is Limited to within the Solar Alfven Surface

The decay of the solar wind helium to hydrogen temperature ratio due to Coulomb thermalization can be used to measure how far from the Sun strong preferential ion heating occurs. Previous work has shown that a zone of preferential ion heating, resulting in mass-proportional temperatures, extends about $20-40 R_\odot$ from the Sun on average. Here we look at the motion of the outer boundary of this zone with time and compare it to other physically meaningful distances. We report that the boundary moves in lockstep with the Alfv\'en point over the solar cycle, contracting and expanding with solar activity with a correlation coefficient of better than 0.95 and with an RMS difference of $4.23 R_\odot$. Strong preferential ion heating apparently is predominatly active below the Alfv\'en point. To definitively identify the underlying preferential heating mechanisms, it will be necessary to make in situ measurements of the local plasma conditions below the Alfv\'en surface. We predict Parker Solar Probe (PSP) will be the first spacecraft to directly observe this heating in action, but only a couple of years after launch as activity increases, the zone expands, and PSP's perihelion drops.


INTRODUCTION
Ions in the solar corona and solar wind are too hot, a puzzle since the start of the space age and the first in situ observations of solar wind plasma by spacecraft. Since electrons carry the heat flux and are the ultimate sink of turbulent energy, one would expect electrons to be hotter than ions, and for this difference to grow more extreme with distance from the Sun; however, ions are hotter in the corona and have similar temperatures to electrons at 1 au. Explaining how ions are heated in the corona and solar wind remains a major challenge in the field. Spectroscopic observations allow us to remotely observe the onset of this preferential ion heating and its consequences. A few tenths of a solar radius R above the Sun's photosphere, rising temperatures and falling densities greatly diminish the Coulomb collisions that enforce thermal equilibrium. Unidentified plasma heating mechanisms in this region couple to different ion species with varying efficiency, preferentially heating some ions more rapidly than others. Extreme differences in ion temperatures develop, with some species reaching temperatures beyond 100 MK (Parker 1988;Kohl et al. 1998;Landi & Cranmer 2009). The ratio of heavy ion species temperature to proton temperature T s /T p is observed to reach and even exceed the mass ratio m s /m p . This suggests a kinetic heating process involving interactions with waves or fluctuations with a characteristic velocity, as ions have equal thermal speeds when they have mass proportional temperatures. Remote observations provide some insight into the mechanisms injecting energy at the base of the corona (McIntosh et al. 2011;Grant et al. 2018) but these observations are not sufficient to distinguish between the various mechanisms that have been proposed to lead to preferential ion heating throughout the near-Sun environment, including wave damping, turbulent dissipation, shocks, reconnection, nano-flares, and velocity filtration (see reviews by Ofman (2010), Hansteen & Velli (2012) and Cranmer (2012)). At some distance, strong preferential ion heating ceases, and ion temperature differences in the solar wind begin to decay with increasing time, as infrequent Coulomb collisions begin to thermalize the plasma (Neugebauer 1976;Hernandez et al. 1987;Tracy et al. 2015). We will refer to the heating process active near the Sun that results in mass proportional, and even super-mass proportional temperatures, as strong preferential heating. In solar wind far from the Sun, mass proportional temperatures are only observed when the frequency of Coulomb collisions is low. Helium and heavier ions in the solar wind with high Coulomb collision rates are at most tens of percent hotter than protons, which is either an indication of the temperature measurement error of the instrument or a sign that only much weaker preferential ion heating occurs in interplanetary space Tracy et al. 2016).
Recently we demonstrated a technique for using solar wind observations at 1 au to determine how far from the Sun the strong preferential ion heating occurs (Kasper et al. (2017), referred to as Paper I). We proposed that there is a zone close to the Sun where ion species experience strong preferential heating, and that within this zone of preferential heating ions reach an equilibrium temperature with an unspecified heating mechanism resulting in different steady temperature ratios for different ion species relative to protons. The start of this zone is seen in the spectroscopic observations just a few 0.1R above the photosphere. We further assumed that there is an outer boundary of the zone, at a distance R b from the Sun, beyond which strong preferential heating ends, and ions are either heated equally or at some much weaker preferential rate. Beyond R b , Coulomb relaxation, or the accumulated impact of many small angle Coulomb scattering interactions between ions, dominates over any weak preferential heating, and will slowly act to drive species toward equal temperatures. For an intuitive sense of this process, consider Spitzer (1962) who showed that if two species have a difference in temperature ∆T and exchange thermal energy via Coulomb scattering at frequency ν c , with only one species having a notable change in temperature, then the temperature difference will change with time as Rearranging and defining the temperature excess ≡ T α /T p − 1, where T α and T p are the temperatures of fully ionized helium and hydrogen yields where A c is the Coulomb age, or the number of Coulomb thermalization times that have elapsed from when the plasma crossed R b to when it was observed in space, and o is the steady state excess temperature ratio developed below R b . Both helium and minor ions seen at 1 au exhibit exponential-like decay in temperature excess with A c (Tracy et al. 2015;Kasper et al. 2017). We can locate R b by using the exponential decay of (A c ) as a clock to measure the time it took the solar wind to move from the outer boundary of the zone of preferential heating to the observing spacecraft. In this paper, we report for the first time the temporal dependence of R b , and find that this outer boundary of strong preferential ion heating is well correlated with the Alfvén critical surface. Implications for in situ observations of preferential heating mechanisms are discussed.

METHODOLOGY
This work uses measurements from the Solar Wind Experiment (Ogilvie et al. 1995) and Magnetic Field Investigation (Lepping et al. 1995) instruments on the NASA Wind spacecraft. The same data selection criterion were used as were presented in Paper I with the added exclusion of data collected prior to October 27, 1997. Before this date, a different observation mode yielded larger measurement uncertainties for T α .
Our model for A c and R b uses radial power law exponents δ and σ to capture how overall temperature and speed vary with distance from the Sun, T (r) = T 0 (r/R) −δ and U (r) = U 0 (r/R) −σ . The solar wind proton density and magnetic field amplitude scale as and where the angular frequency of the Sun's rotation in the equatorial plane is ω = 2.7 × 10 −6 rad s −1 . Quantities with a subscript 0 are values measured at R = 1au. Our full model for (A c ) is more complicated than Eqn 2 because it allows both species temperatures to change, with the relative rates of heating and cooling determined by the relative mass density F ≡ m α n α /m p n p . An increase in F causes a faster equalization of T p and T α . We also account for the variation of ν c due to these relative temperatures changing with distance; see Eqn. 11 in Paper I for the full expression.
In Paper I, we used a range of δ that bracketed published values seen by Helios (Marsch et al. 1982;Hellinger et al. 2011) and assumed σ was zero. A best fit for observed as a function of A c , with R b and o as free parameters was found. The analysis in Paper I and in this work is limited to solar wind intervals where the speed is between 300 and 500 km/s to ensure good data coverage of both high and low collisional age plasma. Using all solar wind speeds does not qualitatively affect the results presented here. Over the entire Wind mission the model fits the observations with a Pearson's chi-squared test of χ 2 /dof of less than 2, and can predict the mean for a given A c with an RMS error of less than 10%. If δ is specified, the best fit uncertainty in R b is much less than one R . For every 0.1 increase in δ, R b drops 8.8 R closer to the Sun. For the range of δ reported in the literature R b could be between 20 − 40R . Figure 1. Observed (Ac) for three year-and-a-half intervals with increasing solar activity, colored symbols. Solid lines are predicted (Ac) for each interval, factoring in the increasing observed mass density ratio F but holding R b = 25 . At solar minimum, 25R is an overestimate of R b , with the observed decaying faster than the model. As activity increased R b is underestimated, suggesting it is moving outwards from the Sun.

RESULTS
Temporal variation in (A c ) and the outer boundary R b can be seen in a relatively raw view of the observations. Fig. 1 compares the observed (A c ) (symbols) with the expected decay (lines) if R b is taken to be constant but the mass density ratio F is updated to account for its observed solar cycle dependence. Over these three 1.5year intervals F grows from 0.06 to 0.14 (Kasper et al. 2012), and with more helium, the two species can reach thermal equilibrium faster. Most of the observations fall below the curve in solar minimum (light blue diamonds), and are generally all above the curve as activity increases (red squares). This can be explained as an underestimate of A c in solar minimum because R b is closer to the Sun than we assumed, and an overestimate of A c in times of high activity because R b has moved closer to the observer than assumed. The elevated at large A c in 2008 appears to be to due a higher uncertainty in the temperature of helium caused by lower helium densities which does not impact this analysis.
We next calculate a best fit value for R b as a function of year over the entire Wind mission, and compare it to a proxy for solar activity and to several critical surfaces surrounding the Sun, where the bulk solar wind speed transitions from below to above some characteristic wave speed. It is also known that critical surfaces around the Sun also have temporal dependencies associated with the solar cycle (Katsikas et al. 2010). We focus on two of these surfaces, related to two fundamental waves speeds in a magnetized plasma, the Alfvén speed v A and the sound speed v s . For a radial solar wind profile U (r) the radial location of the Alfvén and sound critical surfaces, R A and R s , are where v A (r) = U (r) and v s (r) = U (r). Beyond these critical surfaces, an Alfvén or sound wave respectively cannot travel back to the Sun. The Alfvén speed is calculated as Using an Alfvén speed with only the proton mass density or the total ion mass density does not significantly alter the results presented below. For the sound critical surface R s , we use Eqn 11 from Katsikas et al. (2010) evaluated at equatorial latitudes. For each measurement at 1au, assuming a particular radial scaling for the solar wind temperature (δ) and velocity (σ), we determine R A and R s . The distribution of R A as a function of time is shown in Fig. 2 as a column-normalized two-dimensional histogram; R A typically has values around 25R with significant expansion and contraction between solar minimum and solar maximum. The outer boundary of the zone of preferential heating R b is calculated using the same scheme as in Paper I, except instead of segregating by solar wind speed, the data is sub-divided into one-and-a-half year intervals. The mean value of R b as a function of time from 1998 to 2017 is plotted in Fig. 3 for values for radial scalings of proton temperature and speed consistent with previous measurements of the solar wind, (δ = 0.814 and σ = −0.05, c.f. Hellinger et al. (2011Hellinger et al. ( , 2013). We observe significant variation in R b , ranging from ∼ 10R at solar minimum to ∼ 35R at solar maximum. This variation occurs for both the relatively strong cycle 23 and the weaker cycle 24. The one-sigma error in R b is quite narrow, on the order of 0.5R . To compare with R b , we also calculate average values and RMS variations of R A and R s , plotted in red and blue in Fig. 3. The sound critical surface R s is much closer to the Sun than R b and does not have the same temporal variation with solar activity. The Alfvén critical surface's temporal variation is well correlated with R b , with a RMS difference between the two distances of less than 4.23R and a Spearman rank correlation of 0.956 with significance of 3 × 10 −6 . The correlation is significantly better than with sunspot number, a typical indicator of solar activity (0.842), or the sound surface (−0.367). Another physically meaningful distance is where the ratio of thermal and magnetic pressures, β(r) = 8πn p (r)k B T p (r)B(r) −2 , crosses some critical value; we denote this surface as R β . Choosing different critical β values of less than unity, we find that R β and R b are well correlated, with a Spearman rank correlation coefficient of 0.820, though not as well correlated as R b and R A . Using standard minimization techniques, we determine that a critical β value of 0.023 has the smallest RMS distance between R b and R β . We note, however, that this may be a transitive effect; the mean value of β at R A is approximately 0.03, with only minimal temporal variation. Therefore, as R b is correlated with R A and β(R A ) ≈ 0.03, we expect the distance R(β ≈ 0.03) will be correlated with R b .
We repeat the process of calculating R b , R A , R s , and R β for a range of radial power law exponents for proton temperature and speed, δ ∈ [0.75, 0.95] and σ ∈ [−0.1, 0] compatible with radial trends extracted from Helios. The Spearman rank correlation coefficients between R b and the three critical surfaces have little variation due to the power law exponents (not shown); for all values of δ and σ considered, R A is by far the best correlated surface. Further, we find that there is no meaningful global minimum value for δ and σ in the RMS difference between R b and R A , ∆R A,b nor any global maximum for their correlation coefficient. Rather, there is a family of solutions for which ∆R A,b is minimized. Specifically, there is a unique value of δ for any σ resulting in a minimum ∆R A,b (σ, δ) of 5R , as illustrated in Fig. 4. The δ leading to the minimum ∆R A,b obeys the equation δ = 0.813 − 1.037σ. These preferred solutions follow closely, but are not identical to the δ = 2/3−4/3σ scaling suggested by the radial dependence of the collisional age integral discussed in Paper I. The minor disagreement between these trends is likely caused by uncertainty introduced by our model for radial variation in densities and Alfvén speeds. Venzmer & Bothmer (2018) have produced limits on values for σ and δ based upon recent analysis of Helios radial trends which bound the ∆R A,b minimum extracted from the outer boundary analysis. Using the mean value from that study, σ = −0.05, we predict that δ = 0.85.

DISCUSSION
These results suggest that the outer boundary of the zone of strong preferential ion heating R b is the Alfvén critical surface R A , and that the zone and the Alfvén surface expand in lockstep as solar activity changes. This leads to the question of why any preferential heating mechanism would be affected by a transition across this surface. In the expanding solar wind, some fraction of outward propagating wave-power is reflected back toward the Sun due to large-scale gradients in background quantities. Below R A , these waves can travel all the way back to the Sun and interact with outward propagating waves, leading to wave-reflection driven turbulence (Matthaeus et al. 1999;Perez & Chandran 2013). Above R A , backward propagating waves in the plasma frame are advected forward in the Sun's reference frame. The abundance of counter-propagating waves below R A can dramatically enhance local preferential ion heating, either due to wave-particle interactions  or alterations to the background turbulence (Velli et al. 1989;Matthaeus et al. 1999;Dmitruk et al. 2001;Verdini & Velli 2007;Cranmer et al. 2007;Chandran & Hollweg 2009;Verdini et al. 2012). No theoretical predictions of dissipation have suggested a sharp change in the preferential heating exactly at R A . Like crossing the event horizon of a black hole, there is no sudden change experienced when crossing R A , even as the plasma becomes causally disconnected from the Sun. While it is possible that there is a discrete termination of preferential heating mechanisms, it is more likely that the preferential heating gradually shuts off after the plasma passes R A . Given that we find ∆R A,b to always be greater than ≈ 5R , this difference may serve as a estimate for the thickness over which the preferential heating ceases.
Could there be some form of discontinuity or transition at the Alfvén point? Weber & Davis (1967) proposed that the corona co-rotates with the Sun out to the Alfvén point, with a sudden drop in rotational speed as the solar wind Alfvén mach number exceeded unity. Perhaps the diminished role of magnetic tension and the rotation of the magnetic field either alters the local turbulence or enhances reflection back towards the Sun. Remote observations where the solar wind transitions character from striated to flocculated (DeForest et al. 2016) report distances of 44 − 88R , slightly beyond R A , indicating that this transition region is not coterminous with R b . The Alfvén surface has been identified in numerical MHD simulations (Chhiber et al. 2018) as the region where large-scale magnetohydrodynamic turbulence first manifests, potentially changing the mechanisms heating the plasma.
We emphasize that we are only predicting that strong preferential ion heating shuts off beyond R A , not that all ion heating is terminated outside of this zone. Weak preferential heating of minor ions may also continue outside of this zone, leading to no more than a ten percent difference between proton and alpha temperature, as reported for highly collisional solar wind at 1au Tracy et al. 2016); this temperature difference between species is also consistent with the temperature measurement error of the particle instruments (Kasper et al. 2006).
Parker Solar Probe (PSP), launched in August 2018, is the first spacecraft to enter the near-Sun environment, with an initial perihelion of 35R in 2018 and final perihelia of 9.86R starting in 2024 (Fox et al. 2016). The first scientific objective of PSP is to "trace the flow of energy that heats and accelerates the solar corona and solar wind". By closing to within 10R of the Sun's surface, PSP will have a high probability of observing non-thermal heating in action with its electromagnetic field ) and ion and electron plasma  instruments. The explicit assumption has been that we are much more likely to observe this heating in action below the Alfvén point because we are in the magnetic atmosphere of the Sun, casually connected to the Sun, or simply because it may be easier to map the plasma to its sources that close to the Sun.
With the results reported here we can make a specific prediction for how PSP may observe and reveal preferential ion heating in action, and the underlying physics, for the first time. We have found that not only is there a zone of preferential ion heating surrounding the Sun that extends tens of R from the Sun, but the outer boundary R b of this zone expands with solar activity, closely tracking the location of the Alfvén surface and likely intimately connected to changes in the nature of the plasma and waves across this surface. In Fig. 5, we advance our calculation of R b using Wind observations by 11 and 22 years to project forward data from the last two solar cycles into the PSP mission timeframe. Against this projection we plot the minimum distance of approach to the Sun by PSP, which steps closer to the Sun via six Venus gravity assists. We find that at launch in solar minimum, PSP's perihelion is too high, and R b is too low, for PSP to enter the zone. However, in late 2020 as PSP's perihelion lowers, the preferential heating zone and the Alfvén point will extend outwards and cross the trajectory of the spacecraft. This prediction varies slightly when using R b from the relatively strong cycle 23 or the weaker cycle 24. See also Chhiber et al. (2019) for a discussion of predictions from global MHD simulations for the Alfvén critical surface.
As R b rarely exceeds 35R , it is likely that no previous spacecraft has sampled this region of preferential heating. As PSP approaches the zone, the derived distance to R b calculated using the method described in this work should decrease. By comparing local measurements of the Alfvén speed and the solar wind bulk speed, we should be able to determine when PSP crosses the Alfvén point. At that point we predict PSP will be able to detect if has reached its expected asymptotic value as well as signatures local heating processes in this region. By doing so, PSP should be able to fulfill its first scientific objective of characterizing how the solar corona and solar wind are heated. While PSP starts with a perihelion above R b , we suggest the spacecraft and boundary will cross in late-2020, permitting the first direct observations of preferential ion heating.