This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Observations at 1 au have confirmed that enhancements in measured energetic particle fluxes are statistically associated with "rough" magnetic fields, i.e., fields having atypically large spatial derivatives or increments, as measured by the Partial Variance of Increments (PVI) method. One way to interpret this observation is as an association of the energetic particles with trapping or channeling within magnetic flux tubes, possibly near their boundaries. However, it remains unclear whether this association is a transport or local effect; i.e., the particles might have been energized at a distant location, perhaps by shocks or reconnection, or they might experience local energization or re-acceleration. The Parker Solar Probe (PSP), even in its first two orbits, offers a unique opportunity to study this statistical correlation closer to the corona. As a first step, we analyze the separate correlation properties of the energetic particles measured by the \isois instruments during the first solar encounter. The distribution of time intervals between a specific type of event, i.e., the waiting time, can indicate the nature of the underlying process. We find that the \isois observations show a power-law distribution of waiting times, indicating a correlated (non-Poisson) distribution. Analysis of low-energy \isois data suggests that the results are consistent with the 1 au studies, although we find hints of some unexpected behavior. A more complete understanding of these statistical distributions will provide valuable insights into the origin and propagation of solar energetic particles, a picture that should become clear with future PSP orbits.

Direct evidence of an inertial-range turbulent energy cascade has been provided by spacecraft observations in heliospheric plasmas. In the solar wind, the average value of the derived heating rate near 1 au is $\sim 10^{3}\, \mathrm{J\,kg^{-1}\,s^{-1}}$, an amount sufficient to account for observed departures from adiabatic expansion. Parker Solar Probe (PSP), even during its first solar encounter, offers the first opportunity to compute, in a similar fashion, a fluid-scale energy decay rate, much closer to the solar corona than any prior in-situ observations. Using the Politano-Pouquet third-order law and the von K\'arm\'an decay law, we estimate the fluid-range energy transfer rate in the inner heliosphere, at heliocentric distance $R$ ranging from $54\,R_{\odot}$ (0.25 au) to $36\,R_{\odot}$ (0.17 au). The energy transfer rate obtained near the first perihelion is about 100 times higher than the average value at 1 au. This dramatic increase in the heating rate is unprecedented in previous solar wind observations, including those from Helios, and the values are close to those obtained in the shocked plasma inside the terrestrial magnetosheath.

The solar wind shows periods of highly Alfv\'enic activity, where velocity fluctuations and magnetic fluctuations are aligned or anti-aligned with each other. It is generally agreed that solar wind plasma velocity and magnetic field fluctuations observed by Parker Solar Probe (PSP) during the first encounter are mostly highly Alfv\'enic. However, quantitative measures of Alfv\'enicity are needed to understand how the characterization of these fluctuations compares with standard measures from prior missions in the inner and outer heliosphere, in fast wind and slow wind, and at high and low latitudes. To investigate this issue, we employ several measures to quantify the extent of Alfv\'enicity -- the Alfv\'en ratio $r_A$, {normalized} cross helicity $\sigma_c$, {normalized} residual energy $\sigma_r$, and the cosine of angle between velocity and magnetic fluctuations $\cos\theta_{vb}$. We show that despite the overall impression that the Alfv\'enicity is large in the solar wind sampled by PSP during the first encounter, during some intervals the cross helicity starts decreasing at very large scales. These length-scales (often $> 1000 d_i$) are well inside inertial range, and therefore, the suppression of cross helicity at these scales cannot be attributed to kinetic physics. This drop at large scales could potentially be explained by large-scale shears present in the inner heliosphere sampled by PSP. In some cases, despite the cross helicity being constant down to the noise floor, the residual energy decreases with scale in the inertial range. These results suggest that it is important to consider all these measures to quantify Alfv\'enicity.

During the Parker Solar Probe's (PSP) first perihelion pass, the spacecraft reached within a heliocentric distance of \(\sim 37~R_\odot\) and observed numerous magnetic and flow structures characterized by sharp gradients. To better understand these intermittent structures in the young solar wind, an important property to examine is their degree of correlation in time and space. To this end, we use the well-tested Partial Variance of Increments (PVI) technique to identify intermittent events in FIELDS and SWEAP observations of magnetic and proton-velocity fields (respectively) during PSP's first solar encounter, when the spacecraft was within 0.25 au from the Sun. We then examine distributions of waiting times between events with varying separation and PVI thresholds. We find power-law distributions for waiting times shorter than a characteristic scale comparable to the correlation time, suggesting a high degree of correlation that may originate in a clustering process. Waiting times longer than this characteristic time are better described by an exponential, suggesting a random memory-less Poisson process at play. These findings are consistent with near-Earth observations of solar wind turbulence. The present study complements the one by Dudok de Wit et al. (2020, present volume), which focuses on waiting times between observed "switchbacks" in the radial magnetic field.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

The Parker Solar Probe (PSP) primary mission extends seven years and consists of 24 orbits of the Sun with descending perihelia culminating in a closest approach of ($\sim 9.8~R_\odot$). In the course of these orbits PSP will pass through widely varying conditions, including anticipated large variations of turbulence properties such as energy density, correlation scales and cross helicities. Here we employ global magnetohydrodynamics simulations with self-consistent turbulence transport and heating \citep{usmanov2018} to preview likely conditions that will be encountered by PSP, by assuming suitable boundary conditions at the coronal base. The code evolves large-scale parameters -- such as velocity, magnetic field, and temperature -- as well as turbulent energy density, cross helicity, and correlation scale. These computed quantities provide the basis for evaluating additional useful parameters that are derivable from the primary model outputs. Here we illustrate one such possibility in which computed turbulence and large-scale parameters are used to evaluate the accuracy of the Taylor "frozen-in" hypothesis along the PSP trajectory. Apart from the immediate purpose of anticipating turbulence conditions that PSP will encounter, as experience is gained in comparisons of observations with simulated data, this approach will be increasingly useful for planning and interpretation of subsequent observations.

The solar corona and young solar wind may be characterized by critical surfaces -- the sonic, Alfv\'en, and first plasma-$\beta$ unity surfaces -- that demarcate regions where the solar wind flow undergoes certain crucial transformations. Global numerical simulations and remote sensing observations offer a natural mode for the study of these surfaces at large scales, thus providing valuable context for the high-resolution in-situ measurements expected from the soon-to-be-launched Parker Solar Probe (PSP). The present study utilizes global three-dimensional magnetohydrodynamic simulations of the solar wind to characterize the critical surfaces and investigate the flow in propinquitous regions. Effects of solar activity are incorporated by varying source magnetic dipole tilts and employing magnetogram-based boundary conditions. A magnetohydrodynamic turbulence model is self-consistently coupled to the bulk flow equations, enabling investigation of turbulence properties of the flow in the vicinity of critical regions. The simulation results are compared with a variety of remote sensing observations. A simulated PSP trajectory is used to provide contextual predictions for the spacecraft in terms of the computed critical surfaces. Broad agreement is seen in the interpretation of the present results in comparison with existing remote sensing results, both from heliospheric imaging and from radio scintillation studies. The trajectory analyses show that the period of time that PSP is likely to spend inside the $\beta=1$, sonic and Alfv\'en surfaces depends sensitively on the degree of solar activity and the tilt of the solar dipole and location of the heliospheric current sheet.

MMS observations recently confirmed that crescent-shaped electron velocity distributions in the plane perpendicular to the magnetic field occur in the electron diffusion region near reconnection sites at Earth's magnetopause. In this paper, we re-examine the origin of the crescent-shaped distributions in the light of our new finding that ions and electrons are drifting in opposite directions when displayed in magnetopause boundary-normal coordinates. Therefore, ExB drifts cannot cause the crescent shapes. We performed a high-resolution multi-scale simulation capturing sub-electron skin depth scales. The results suggest that the crescent-shaped distributions are caused by meandering orbits without necessarily requiring any additional processes found at the magnetopause such as the highly asymmetric magnetopause ambipolar electric field. We use an adiabatic Hamiltonian model of particle motion to confirm that conservation of canonical momentum in the presence of magnetic field gradients causes the formation of crescent shapes without invoking asymmetries or the presence of an ExB drift. An important consequence of this finding is that we expect crescent-shaped distributions also to be observed in the magnetotail, a prediction that MMS will soon be able to test.

Continuous plasma coherent emission is maintained by repetitive Langmuir collapse driven by the nonlinear evolution of a strong electron two-stream instability. The Langmuir waves are modulated by solitary waves in the linear stage, and by electrostatic whistler waves in the nonlinear stage. Modulational instability leads to Langmuir collapse and electron heating that fills in cavitons. The high pressure is released via excitation of a short wavelength ion acoustic mode that is damped by electrons and that re-excites small-scale Langmuir waves---this process closes a feedback loop that maintains the continuous coherent emission.

The formation of the observed core-halo feature in the solar wind electron velocity distribution function is a long-time puzzle. In this letter based on the current knowledge of nanoflares we show that the nanoflare-accelerated electron beams are likely to trigger a strong electron two-stream instability that generates kinetic Alfv\'en wave and whistler wave turbulence, as we demonstrated in a previous paper. We further show that the core-halo feature produced during the origin of kinetic turbulence is likely to originate in the inner corona and can be preserved as the solar wind escapes to space along open field lines. We formulate a set of equations to describe the heating processes observed in the simulation and show that the core-halo temperature ratio of the solar wind is insensitive to the initial conditions in the corona and is related to the core-halo density ratio of the solar wind and to the quasi-saturation property of the two-stream instability at the time when the exponential decay ends. This relation can be extended to the more general core-halo-strahl feature in the solar wind. The temperature ratio between the core and hot components is nearly independent of the heliospheric distance to the sun. We show that the core-halo relative drift previously reported is a relic of the fully saturated two stream instability. Our theoretical results are consistent with the observations while new tests for this model are provided.

Based on global conservation principles, magnetohydrodynamic (MHD) relaxation theory predicts the existence of several equilibria, such as the Taylor state or global dynamic alignment. These states are generally viewed as very long-time and large-scale equilibria, which emerge only after the termination of the turbulent cascade. As suggested by hydrodynamics and by recent MHD numerical simulations, relaxation processes can occur during the turbulent cascade that will manifest themselves as local patches of equilibrium-like configurations. Using multi-spacecraft analysis techniques in conjunction with Cluster data, we compute the current density and flow vorticity and for the first time demonstrate that these localized relaxation events are observed in the solar wind. Such events have important consequences for the statistics of plasma turbulence.

The observed sub-proton scale turbulence spectrum in the solar wind raises the question of how that turbulence originates. Observations of keV energetic electrons during solar quite-time suggest them as possible source of free energy to drive the turbulence. Using particle-in-cell simulations, we explore how free energy in energetic electrons, released by an electron two-stream instability drives Weibel-like electromagnetic waves that excite wave-wave interactions. Consequently, both kinetic Alfv\'enic and whistler waves are excited that evolve through inverse and forward magnetic energy cascades.

Electron scale solar wind turbulence has attracted great interest in recent years. Clear evidences have been given from the Cluster data that turbulence is not fully dissipated near the proton scale but continues cascading down to the electron scales. However, the scaling of the energy spectra as well as the nature of the plasma modes involved at those small scales are still not fully determined. Here we survey 10 years of the Cluster search-coil magnetometer (SCM) waveforms measured in the solar wind and perform a statistical study of the magnetic energy spectra in the frequency range [$1, 180$]Hz. We show that a large fraction of the spectra exhibit clear breakpoints near the electon gyroscale $\rho_e$, followed by steeper power-law like spectra. We show that the scaling below the electron breakpoint cannot be determined unambiguously due to instrumental limitations that will be discussed in detail. We compare our results to recent ones reported in other studies and discuss their implication on the physical mechanisms and the theoretical modeling of energy dissipation in the SW.

The nonlinear development of the strong Buneman instability and the associated fast electron heating in thin current layers with $\Omega_e/\omega_{pe} <1$ are explored. Phase mixing of the electrons in wave potential troughs and a rapid increase in temperature are observed during the saturation of the instability. We show that the motion of trapped electrons can be described using a Hamiltonian formalism in the adiabatic approximation. The process of separatrix crossing as electrons are trapped and de-trapped is irreversible and guarantees that the resulting electron energy gain is a true heating process.

The nature of solar wind (SW) turbulence below the proton gyroscale is a topic that is being investigated extensively nowadays. Although recent observations gave evidence of the dominance of Kinetic Alfv\'en Waves (KAW) at sub-ion scales with $\omega<{\omega_{ci}}$, other studies suggest that the KAW mode cannot carry the turbulence cascade down to electron scales and that the whistler mode (i.e., $\omega>\omega_{ci}$) is more relevant. Here, we propose to study key properties of the short wavelength plasma modes under realistic SW conditions, typically $\beta_i\gtrsim \beta_e\sim 1$ and for high oblique angles of propagation $80^\circ\leq \Theta_{\bf kB}<90^\circ$ as observed from the Cluster data. The linear properties of the plasma modes under these conditions are poorly known, which contrasts with the well-documented cold plasma limit and/or moderate oblique angles of propagation ($\Theta_{\bf kB} <80^\circ$). Based on linear solutions of the Vlasov kinetic theory, we discuss the relevance of each plasma mode (fast, Bernstein, KAW, whistler) in carrying the energy cascade down to electron scales. We show, in particular, that the shear Alfv\'en mode extends at scales $k\rho_i\gtrsim1$ following either a whistler mode ($\omega>\omega_{ci}$) or a KAW mode (with $\omega<\omega_{ci}$) depending on the anisotropy $k_\parallel/ k_\perp$. This contrasts with the well-accepted idea that the whistler branch develops as a continuation at high frequencies of the fast magnetosonic mode. We show, furthermore, that the whistler branch is more damped than the KAW one, which makes the latter a more relevant candidate to carry the energy cascade down to electron scales. We discuss how these new findings may facilitate resolution of the controversy concerning the nature of the small scale turbulence, and we discuss the implications for present and future spacecraft wave measurements in the SW.

We determine the magnetic helicity, along with the magnetic energy, at high latitudes using data from the Ulysses mission. The data set spans the time period from 1993 to 1996. The basic assumption of the analysis is that the solar wind is homogeneous. Because the solar wind speed is high, we follow the approach first pioneered by Matthaeus et al. (1982, Phys. Rev. Lett. 48, 1256) by which, under the assumption of spatial homogeneity, one can use Fourier transforms of the magnetic field time series to construct one-dimensional spectra of the magnetic energy and magnetic helicity under the assumption that the Taylor frozen-in-flow hypothesis is valid. That is a well-satisfied assumption for the data used in this study. The magnetic helicity derives from the skew-symmetric terms of the three-dimensional magnetic correlation tensor, while the symmetric terms of the tensor are used to determine the magnetic energy spectrum. Our results show a sign change of magnetic helicity at wavenumber k~2 AU^{-1} (or frequency nu~2 uHz) at distances below 2.8 AU and at k~30 AU^{-1} (or nu~25 uHz) at larger distances. At small scales the magnetic helicity is positive at northern heliographic latitudes and negative at southern latitudes. The positive magnetic helicity at small scales is argued to be the result of turbulent diffusion reversing the sign relative to what is seen at small scales at the solar surface. Furthermore, the magnetic helicity declines toward solar minimum in 1996. The magnetic helicity flux integrated separately over one hemisphere amounts to about 10^{45} Mx^2/cycle at large scales and to a 3 times lower value at smaller scales.

Semi-analytic models are a powerful tool for studying the formation of galaxies. However, these models inevitably involve a significant number of poorly constrained parameters that must be adjusted to provide an acceptable match to the observed universe. In this paper, we set out to quantify the degree to which observational data-sets can constrain the model parameters. By revealing degeneracies in the parameter space we can hope to better understand the key physical processes probed by the data. We use novel mathematical techniques to explore the parameter space of the GALFORM semi-analytic model. We base our investigation on the Bower et al. 2006 version of GALFORM, adopting the same methodology of selecting model parameters based on an acceptable match to the local bJ and K luminosity functions. The model contains 16 parameters that are poorly constrained, and we investigate this parameter space using the Model Emulator technique, constructing a Bayesian approximation to the GALFORM model that can be rapidly evaluated at any point in parameter space. By combining successive waves of emulation, we show that only 0.26% of the initial volume is of interest for further exploration. However, within this region we show that the Bower et al. 2006 model is only one choice from an extended sub-space of model parameters that can provide equally acceptable fits. We explore the geometry of this region and begin to explore the physical connections between parameters that are exposed by this analysis. We also consider the impact of adding additional observational data to further constrain the parameter space.

(Abridged abstract) A theory of incompressible MHD turbulence recently developed by Boldyrev predicts the existence of a scale dependent angle of alignment between velocity and magnetic field fluctuations that is proportional to the lengthscale of the fluctuations to the power 1/4. In this study, plasma and magnetic field data from the Wind spacecraft are used to investigate the angle between velocity and magnetic field fluctuations in the solar wind as a function of the timescale of the fluctuations and to look for the power law scaling predicted by Boldyrev.

We present spectra of the extended 12CO and 13CO J=1-0 emission along the major axes of 17 nearby galaxies. Spatial variations in the ratio of CO and 13CO integrated intensities, R, are found in nearly every galaxy observed. There is an overall variation in R of 20-40% from the inner 2 kpc to the disk. Roughly a third of the survey galaxies have such gradients in R detected above the 2-sigma confidence level. Though some galaxies show a lower central value of R, on average R inside 2 kpc is 10-30% higher than R outside of 2 kpc. An increase in the CO/13CO intensity ratio from disk to nucleus may imply that the conversion factor between CO intensity and H2 column density, X, is lower in galactic nuclei. Also variations in physical conditions, most notably the gas kinetic temperature, affect both R and X. Abundance variations probably do not cause the gradient in R, though we do not rule out a decrease in effective cloud column densities in galactic nuclei possibly caused by destructive starburst superwinds. A modest rise in temperature (less than a factor of 2 or 3) from outside a 2 kpc radius towards the nucleus can easily account for the observed gradient. These results support previous work implying that X is lower in the center of the Milky Way and probably most galactic nuclei. Therefore calculating H2 masses using the standard Galactic X-factor, especially within the central few kpc of galaxies, overestimates the true mass by factors of a few. The standard X-factor still appears to be appropriate for galactic disks.