The New Horizons Solar Wind Around Pluto (NH SWAP) instrument has provided the first direct observations of interstellar and He pickup ions (PUIs) at distances between ∼11.26 and 38 au in the solar wind. The observations demonstrate that the distant solar wind beyond the hydrogen ionization cavity is indeed mediated by PUIs. The creation of PUIs modifies the underlying low-frequency turbulence field responsible for their own scattering. The dissipation of these low-frequency fluctuations serves to heat the solar wind plasma, and accounts for the observed non-adiabatic solar wind temperature profile and a possible slow temperature increase beyond ∼30 au. We develop a very general theoretical model that incorporates PUIs, solar wind thermal plasma, the interplanetary magnetic field, and low-frequency turbulence to describe the evolution of the large-scale solar wind, PUIs, and turbulence from 1–84 au, the structure of the perpendicular heliospheric termination shock, and the transmission of turbulence into the inner heliosheath, extending the classical models of Holzer and Isenberg. A detailed comparison of the theoretical model solutions and observations derived from the Voyager 2 and NH SWAP data sets shows excellent agreement between the two for reasonable physical parameters.

The 2D + slab superposition model of solar wind turbulence has its theoretical foundations in nearly incompressible magnetohydrodynamics (NI MHD) in the plasma beta ∼1 or ≪1 regimes. Solar wind turbulence measurements show that turbulence in the inertial range is anisotropic, for which the superposition model offers a plausible explanation. We provide a detailed theoretical analysis of the spectral characteristics of the Elsässer variables in the 2D + NI/slab model. We find that (1) the majority 2D component has a power spectrum in perpendicular wavenumber k ⊥; (2) the strongly imbalanced minority NI/slab turbulence has power spectra and , where k z is aligned with the mean magnetic field; (3) NI/slab turbulence can exhibit a double-power-law spectrum, with the steeper part being G*(k) ∼ k −5/3 and corresponding to strong turbulence and the flatter spectrum satisfying G*(k) ∼ k −3/2 and corresponding to weak turbulence; (4) there is a critical balance regime for NI/slab turbulence that satisfies and ; and (5) the forward and backward Elsässer power spectra can have different spectral forms provided that the triple-correlation times for each are different. We use the spectral analysis to compute the total power spectra in frequency parallel to the solar wind flow for the superposition model, showing that strongly imbalanced turbulence yields an f −5/3 spectrum for all angles between the mean flow and magnetic field, and that double power laws are possible when the nonlinear and Alfvén timescales are both finite.

The origin, structure, and propagation characteristics of a switchback are compelling questions posed by Parker Solar Probe (PSP) observations of velocity spikes and magnetic field reversals. By assuming interchange reconnection between coronal loop and open magnetic field, we show that this results in the generation of upward (into the heliosphere) and downward complex structures propagating at the fast magnetosonic speed (i.e., the Alfvén speed in the low plasma beta corona) that can have an arbitrary radial magnetic field deflection, including “S-shaped.” We derive the evolution equation for the switchback radial magnetic field as it propagates through the inhomogeneous supersonic solar corona. An analytic solution for arbitrary initial conditions is used to investigate the properties of a switchback propagating from launch ∼6 to ∼35 R ⊙ where PSP observed switchbacks during its first encounter. We provide a detailed comparison to an example event, showing that the magnetic field and plasma solutions are in accord with PSP observations. For a simple single switchback, the model predicts either a single or a double-humped structure; the former corresponding to PSP observing either the main body or the flanks of the switchback. The clustering of switchbacks and their sometimes complicated structure may be due to the formation of multiple closely spaced switchbacks created by interchange reconnection with numerous open and loop magnetic field lines over a short period. We show that their evolution yields a complex, aggregated group of switchbacks that includes “sheaths” with large-amplitude radial magnetic field and velocity fluctuations.

The possibility that charged particles are accelerated statistically in a “sea” of small-scale interacting magnetic flux ropes in the supersonic solar wind is gaining credence. In this Letter, we extend the Zank et al. statistical transport theory for a nearly isotopic particle distribution by including an escape term corresponding to particle loss from a finite acceleration region. Steady-state 1D solutions for both the accelerated particle velocity distribution function and differential intensity are derived. We show Ulysses observations of an energetic particle flux enhancement event downstream of a shock near 5 au that is inconsistent with the predictions of classical diffusive shock acceleration (DSA) but may be explained by local acceleration associated with magnetic islands. An automated Grad-Shafranov reconstruction approach is employed to identify small-scale magnetic flux ropes behind the shock. For the first time, the observed energetic particle “time-intensity” profile and spectra are quantitatively compared with theoretical predictions. The results show that stochastic acceleration by interacting magnetic islands accounts successfully for the observed (i) peaking of particle intensities behind the shock instead of at the shock front as standard DSA predicts; (ii) increase in the particle flux amplification factor with increasing particle energy; (ii) increase in distance between the particle intensity peak and the shock front with increasing energy; and (iv) hardening of particle power-law spectra with increasing distance downstream of the shock.

The Parker Solar Probe (PSP) entered a region of sub-Alfvénic solar wind during encounter 8, and we present the first detailed analysis of low-frequency turbulence properties in this novel region. The magnetic field and flow velocity vectors were highly aligned during this interval. By constructing spectrograms of the normalized magnetic helicity, cross-helicity, and residual energy, we find that PSP observed primarily Alfvénic fluctuations, a consequence of the highly field-aligned flow that renders quasi-2D fluctuations unobservable to PSP. We extend Taylor’s hypothesis to sub- and super-Alfvénic flows. Spectra for the fluctuating forward and backward Elsässer variables ( z ±, respectively) are presented, showing that z + modes dominate z − by an order of magnitude or more, and the z + spectrum is a power law in frequency (parallel wavenumber) f −3/2 ( k ∥ − 3 / 2 ) compared to the convex z − spectrum with f −3/2 ( k ∥ − 3 / 2 ) at low frequencies, flattening around a transition frequency (at which the nonlinear and Alfvén timescales are balanced) to f −1.25 at higher frequencies. The observed spectra are well fitted using a spectral theory for nearly incompressible magnetohydrodynamics assuming a wavenumber anisotropy k ⊥ ∼ k ∥ 3 / 4 , that the z + fluctuations experience primarily nonlinear interactions, and that the minority z − fluctuations experience both nonlinear and Alfvénic interactions with z + fluctuations. The density spectrum is a power law that resembles neither the z ± spectra nor the compressible magnetic field spectrum, suggesting that these are advected entropic rather than magnetosonic modes and not due to the parametric decay instability. Spectra in the neighboring modestly super-Alfvénic intervals are similar.

A primary goal of the Parker Solar Probe (PSP) Mission is to answer the outstanding question of how the solar corona plasma is heated to the high temperatures needed for the acceleration of the solar wind. Various heating mechanisms have been suggested, but one that is gaining increasing credence is associated with the dissipation of low frequency magnetohyrodynamic (MHD) turbulence. However, the MHD turbulence models come in several flavors: one in which outwardly propagating Alfvén waves experience reflection from the large-scale flow and density gradients associated with the solar corona, and the resulting counterpropagating Alfvén waves couple nonlinearly to produce quasi-2D turbulence that dissipates and heats the corona, thereby driving the solar wind. The second approach eschews a dominant outward flux of Alfvén waves but argues instead that quasi-2D turbulence dominates the lower coronal plasma and is generated in the constantly upwelling magnetic carpet, experiencing dissipation as it is advected through the corona and into the solar wind, yielding temperatures in the corona that exceed a million degrees. We review the two turbulence models, describe the modeling that has been done, and relate PSP observations to the basic predictions of both models. Although PSP measurements are made in the super-Alfvénic solar wind, the observations are close to the coronal region, thus providing a glimpse into the likely properties of coronal turbulence. Observations of low-frequency MHD turbulence by PSP in the super-Alfvénic solar wind allow us to place constraints on models of the turbulently heated solar corona that drive the supersonic solar wind.

We propose a turbulence-driven solar wind model for a fast solar wind flow in an open coronal hole where the solar wind flow and the magnetic field are highly aligned. We compare the numerical results of our model with Parker Solar Probe measurements of the fast solar wind flow and find good agreement between them. We find that (1) the majority quasi-2D turbulence is mainly responsible for coronal heating, raising the temperature to about K within a few solar radii, which leads in turn to the acceleration of the solar wind; (2) the heating rate due to quasi-2D turbulence near the coronal base is larger than that due to nearly incompressible/slab turbulence; (3) the quasi-2D energy in forward-propagating modes decreases with increasing distance, while the nearly incompressible/slab energy in forward-propagating modes increases, reaching a peak value at ∼11.7 before decreasing with increasing heliocentric distance; (4) the correlation length increases with increasing distance from the coronal base; and (5) the variance of the density fluctuations decreases as a function of heliocentric distance.

Recent studies of unusual or atypical energetic particle flux events (AEPEs) observed at 1 au show that another mechanism, different from diffusive shock acceleration, can energize particles locally in the solar wind. The mechanism proposed by Zank et al. is based on the stochastic energization of charged particles in regions filled with numerous small-scale magnetic islands (SMIs) dynamically contracting or merging and experiencing multiple magnetic reconnection in the super-Alfvénic solar wind flow. A first- and second-order Fermi mechanism results from compression-induced changes in the shape of SMIs and their developing dynamics. Charged particles can also be accelerated by the formation of antireconnection electric fields. Observations show that both processes often coexist in the solar wind. The occurrence of SMIs depends on the presence of strong current sheets like the heliospheric current sheet (HCS), and related AEPEs are found to occur within magnetic cavities formed by stream–stream, stream–HCS, or HCS–shock interactions that are filled with SMIs. Previous case studies comparing observations with theoretical predictions were qualitative. Here we present quantitative theoretical predictions of AEPEs based on several events, including a detailed analysis of the corresponding observations. The study illustrates the necessity of accounting for local processes of particle acceleration in the solar wind.