This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

We present a simplified model of the outer heliosheath to help interpret the observations of interstellar neutral helium by the IBEX-Lo instruments. We assume that the measured particles are composed of the superposition of a primary beam population, with the properties of the local interstellar medium, and a secondary population, created by charge exchange between the primary beam neutrals and the ions that have been deflected as they approach the heliopause. We extract information on the large-scale shape of the heliopause by comparing the helium flux measured at IBEX along four different look directions with simple models of deflected plasma flow around hypothetical obstacles of different aspect ratios to the flow. As a first step in this paper, we model the deflected plasma flow with the analytical solutions for compressible gas flow around a series of oblate ellipsoidal obstacles. Our comparisons between the model results and the observations indicate that the heliopause is very blunt in the vicinity of the heliospheric nose, especially compared to a Rankine half-body or cometary shape. The upstream heliopause seems to be highly elongated in the directions parallel to the interstellar magnetic field, and relatively more compact and symmetric in the directions transverse to that field. The IBEX-Lo helium observations are not consistent with a heliopause elongated in directions parallel to the solar rotation axis.

As the local interstellar plasma flows past our heliosphere, it is slowed and deflected around the magnetic obstacle of the heliopause. The interstellar magnetic field, frozen into this plasma, then becomes draped around the heliopause in a characteristic manner. We derive the analytical solution for this draped magnetic field in the limit of weak field intensity, assuming an ideal potential flow around the heliopause, which we model as a Rankine half-body. We compare the structure of the model magnetic field with observed properties of the IBEX ribbon and with in situ observations at the Voyager 1 spacecraft. We find reasonable qualitative agreement, given the idealizations of the model. This agreement lends support to the secondary ENA model of the IBEX ribbon and to the interpretation that Voyager 1 has crossed the heliopause. We also predict that the magnetic field measured by Voyager 2 after it crosses the heliopause will not be significantly rotated away from the direction of the undisturbed interstellar field.

The narrow ribbon of enhanced energetic neutral atom flux observed by the Interstellar Boundary Explorer (IBEX) spacecraft has prompted numerous ideas to explain its structure and properties. One of these ideas is the "neutral solar wind" scenario, which identifies the source particles as pickup protons in the local interstellar medium originating in solar wind charge-exchange interactions. This scenario has been thought to require unrealistically weak pitch-angle scattering of the pickup protons to explain the narrow structure. Recently, Schwadron & McComas (2013) suggested that this structure could result from a spatial retention of the pickup protons, rather than from a restricted pitch-angle distribution. Here, we present a physically motivated, quantitative mechanism to produce such a spatial configuration. This mechanism is based on the "dominant turbulence" assumption, which can be applied where the production of new pickup protons is slow, and has been used to successfully explain the level of turbulent heating observed in the outer solar wind. This formalism predicts a pickup isotropization process which adds or subtracts energy from the ambient turbulent fluctuations, depending on the initial pitch angle of the pickup protons. We show that a simple model of this process can yield a ribbon structure in qualitative agreement with the observations. The results of this simple model are not yet quantitatively satisfactory, but we suggest several improvements which may reduce the quantitative discrepancy.

We extend previous theories of stochastic ion heating to account for the motion of ions along the magnetic field. We derive an analytic expression for the ion-to-proton perpendicular temperature ratio in the solar wind for any ion species, assuming that stochastic heating is the dominant ion heating mechanism. This expression describes how this temperature ratio depends upon the average velocity of the ions along the magnetic field direction and the ratio of the parallel proton pressure to the magnetic pressure. We compare our model with previously published measurements of alpha particles and protons from the WIND spacecraft. We find that stochastic heating offers a promising explanation for these measurements when the fractional cross helicity and Alfven ratio at the proton-gyroradius scale have values that are broadly consistent with solar-wind measurements.

We derive an equation whose solutions describe self-consistent states of marginal stability for a proton-electron plasma interacting with parallel-propagating ion cyclotron waves. Ion cyclotron waves propagating through this marginally stable plasma will neither grow nor damp. The dispersion relation of these waves, {\omega} (k), smoothly rises from the usual MHD behavior at small |k| to reach {\omega} = {\Omega}p as k \rightarrow \pm\infty. The proton distribution function has constant phase-space density along the characteristic resonant surfaces defined by this dispersion relation. Our equation contains a free function describing the variation of the proton phase-space density across these surfaces. Taking this free function to be a simple "box function", we obtain specific solutions of the marginally stable state for a range of proton parallel betas. The phase speeds of these waves are larger than those given by the cold plasma dispersion relation, and the characteristic surfaces are more sharply peaked in the v\bot direction. The threshold anisotropy for generation of ion cyclotron waves is also larger than that given by estimates which assume bi-Maxwellian proton distributions.

Whether the phenomenology governing MHD turbulence is Kolmogorov or Iroshnikov-Kraichnan (IK) remains an open question, theoretically as well as observationally. The ion heating profile observed in the solar wind provides a quantitative, if indirect, observational constraint on the relevant phenomenology. Recently, a solar wind heating model based on Kolmogorov spectral scaling has produced reasonably good agreement with observations, provided the effect of turbulence generation due to pickup ions is included in the model. Without including the pickup ion contributions, the Kolmogorov scaling predicts a proton temperature profile that decays too rapidly beyond a radial distance of 15 AU. In the present study, we alter the heating model by applying an energy cascade rate based on IK scaling, and show that the model yields higher proton temperatures, within the range of observations, with or without the inclusion of the effect due to pickup ions. Furthermore, the turbulence correlation length based on IK scaling seems to follow the trend of observations better.

In recent years, a phenomenological solar wind heating model based on a turbulent energy cascade prescribed by the Kolmogorov theory has produced reasonably good agreement with observations on proton temperatures out to distances around 70 AU, provided the effect of turbulence generation due to pickup ions is included in the model. In a recent study [Ng et al., J. Geophys. Res., 115, A02101 (2010)], we have incorporated in the heating model the energy cascade rate based on Iroshnikov-Kraichnan (IK) scaling. We showed that the IK cascade rate can also produce good agreement with observations, with or without the inclusion of pickup ions. This effect was confirmed both by integrating the model using average boundary conditions at 1 AU, and by applying a method [Smith et al., Astrophys. J., 638, 508 (2006)] that uses directly observed values as boundary conditions. The effects due to pickup ions is found to be less important for the IK spectrum, which is shallower than the Kolmogorov spectrum. In this paper, we will present calculations of the pickup ions effect in more details, and discuss the physical reason why a shallower spectrum generates less waves and turbulence.

In February 2011 we proposed a NASA Explorer Mission of Opportunity program to develop and operate a large-aperture ultraviolet coronagraph spectrometer called the Coronal Physics Investigator (CPI) as an attached International Space Station (ISS) payload. The primary goal of this program is to identify and characterize the physical processes that heat and accelerate the primary and secondary components of the fast and slow solar wind. In addition, CPI can make key measurements needed to understand CMEs. UVCS/SOHO allowed us to identify what additional measurements need to be made to answer the fundamental questions about how solar wind streams are produced, and CPI's next-generation capabilities were designed specifically to make those measurements. Compared to previous instruments, CPI provides unprecedented sensitivity, a wavelength range extending from 25.7 to 126 nm, higher temporal resolution, and the capability to measure line profiles of He II, N V, Ne VII, Ne VIII, Si VIII, S IX, Ar VIII, Ca IX, and Fe X, never before seen in coronal holes above 1.3 solar radii. CPI will constrain the properties and effects of coronal MHD waves by (1) observing many ions over a large range of charge and mass, (2) providing simultaneous measurements of proton and electron temperatures to probe turbulent dissipation mechanisms, and (3) measuring amplitudes of low-frequency compressive fluctuations. CPI is an internally occulted ultraviolet coronagraph that provides the required high sensitivity without the need for a deployable boom, and with all technically mature hardware including an ICCD detector. A highly experienced Explorer and ISS contractor, L-3 Com Integrated Optical Systems and Com Systems East, will provide the tracking and pointing system as well as the instrument, and the integration to the ISS.

Understanding the physical processes responsible for accelerating the solar wind requires detailed measurements of the collisionless plasma in the extended solar corona. Some key clues about these processes have come from instruments that combine the power of an ultraviolet (UV) spectrometer with an occulted telescope. This combination enables measurements of ion emission lines far from the bright solar disk, where most of the solar wind acceleration occurs. Although the UVCS instrument on SOHO made several key discoveries, many questions remain unanswered because its capabilities were limited. This white paper summarizes these past achievements and also describes what can be accomplished with next-generation instrumentation of this kind.

Resonant interactions between ions and Alfv\'en/ion-cyclotron (A/IC) waves may play an important role in the heating and acceleration of the fast solar wind. Although such interactions have been studied extensively for "parallel" waves, whose wave vectors ${\bf k}$ are aligned with the background magnetic field ${\bf B}_0$, much less is known about interactions between ions and oblique A/IC waves, for which the angle $\theta$ between ${\bf k}$ and ${\bf B}_0$ is nonzero. In this paper, we present new numerical results on resonant cyclotron interactions between protons and oblique A/IC waves in collisionless low-beta plasmas such as the solar corona. We find that if some mechanism generates oblique high-frequency A/IC waves, then these waves initially modify the proton distribution function in such a way that it becomes unstable to parallel waves. Parallel waves are then amplified to the point that they dominate the wave energy at the large parallel wave numbers at which the waves resonate with the particles. Pitch-angle scattering by these waves then causes the plasma to evolve towards a state in which the proton distribution is constant along a particular set of nested "scattering surfaces" in velocity space, whose shapes have been calculated previously. As the distribution function approaches this state, the imaginary part of the frequency of parallel A/IC waves drops continuously towards zero, but oblique waves continue to undergo cyclotron damping while simultaneously causing protons to diffuse across these kinetic shells to higher energies. We conclude that oblique A/IC waves can be more effective at heating protons than parallel A/IC waves, because for oblique waves the plasma does not relax towards a state in which proton damping of oblique A/IC waves ceases.

The Balmer line profiles of nonradiative supernova remnant shocks provide the means to measure the post-shock proton velocity distribution. While most analyses assume a Maxwellian velocity distribution, this is unlikely to be correct. In particular, neutral atoms that pass through the shock and become ionized downstream form a nonthermal distribution similar to that of pickup ions in the solar wind. We predict the H alpha line profiles from the combination of pickup protons and the ordinary shocked protons, and we consider the extent to which this distribution could affect the shock parameters derived from H alpha profiles. The Maxwellian assumption could lead to an underestimate of shock speed by up to about 15%. The isotropization of the pickup ion population generates wave energy, and we find that for the most favorable parameters this energy could significantly heat the thermal particles. Sufficiently accurate profiles could constrain the strength and direction of the magnetic field in the shocked plasma, and we discuss the distortions from a Gaussian profile to be expected in Tycho's supernova remnant.

We present a recent analysis with 1996 SOHO CELIAS CTOF data, which reveals correlations of He+ pickup ion fluxes and spectra with the magnetic field strength and solar wind density. The motivation is to better understand the ubiquitous large variations in both pickup ion fluxes and their velocity distributions found in interstellar pickup ion datasets. We concentrate on time periods of that can be associated with compression regions in the solar wind. Along with enhancements of the overall pickup ion fluxes, adiabatic heating and acceleration of the pickup ions are also observed in these regions. Transport processes that lead to the observed compressions and related heating or acceleration are discussed. A shift in velocity space associated with traveling interplanetary compression regions is observed, and a simple model presented to explain this phenomenon based on the conserved magnetic adiabatic moment.