von Karman correlation similarity in solar wind magnetohydrodynamic turbulence

A major development underlying hydrodynamic turbulence theory is the similarity decay hypothesis due to von Karman and Howarth, here extended empirically to plasma turbulence in the solar wind. In similarity decay the second-order correlation experiences a continuous transformation based on a universal functional form and a rescaling of energy and characteristic length. Solar wind turbulence follows many principles adapted from classical fluid turbulence, but previously this similarity property has not been examined explicitly. Here, we analyze an ensemble of Elsässer autocorrelation functions computed from Advanced Composition Explorer data at 1 astronomical unit (AU), and demonstrate explicitly that the two-point correlation functions undergo a collapse to a similarity form of the type anticipated from von Karman's hypothesis applied to weakly compressive magnetohydrodynamic turbulence. This provides a firm empirical basis for employing the similarity decay hypothesis to the Elsässer correlations that represent the incompressive turbulence cascade. This approach is of substantial utility in space turbulence data analysis, and for adopting von Karman-type heating rates in global and subgrid-scale dynamical modeling.

Energy transfer in reconnection and turbulence

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. These results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.

Probing the energetic particle environment near the Sun

NASA's Parker Solar Probe mission1 recently plunged through the inner heliosphere of the Sun to its perihelia, about 24 million kilometres from the Sun. Previous studies farther from the Sun (performed mostly at a distance of 1 astronomical unit) indicate that solar energetic particles are accelerated from a few kiloelectronvolts up to near-relativistic energies via at least two processes: 'impulsive' events, which are usually associated with magnetic reconnection in solar flares and are typically enriched in electrons, helium-3 and heavier ions2, and 'gradual' events3,4, which are typically associated with large coronal-mass-ejection-driven shocks and compressions moving through the corona and inner solar wind and are the dominant source of protons with energies between 1 and 10 megaelectronvolts. However, some events show aspects of both processes and the electron-proton ratio is not bimodally distributed, as would be expected if there were only two possible processes5. These processes have been very difficult to resolve from prior observations, owing to the various transport effects that affect the energetic particle population en route to more distant spacecraft6. Here we report observations of the near-Sun energetic particle radiation environment over the first two orbits of the probe. We find a variety of energetic particle events accelerated both locally and remotely including by corotating interaction regions, impulsive events driven by acceleration near the Sun, and an event related to a coronal mass ejection. We provide direct observations of the energetic particle radiation environment in the region just above the corona of the Sun and directly explore the physics of particle acceleration and transport.

Kinetic Range Spectral Features of Cross Helicity Using the Magnetospheric Multiscale Spacecraft

We study spectral features of ion velocity and magnetic field correlations in the magnetosheath and in the solar wind using data from the Magnetospheric Multiscale (MMS) spacecraft. High-resolution MMS observations enable the study of the transition of these correlations between their magnetofluid character at larger scales into the subproton kinetic range, previously unstudied in spacecraft data. Cross-helicity, angular alignment, and energy partitioning is examined over a suitable range of scales, employing measurements based on the Taylor frozen-in approximation as well as direct two-spacecraft correlation measurements. The results demonstrate signatures of alignment at large scales. As kinetic scales are approached, the alignment between v and b is destroyed by demagnetization of protons.

Magnetospheric Multiscale Observation of Plasma Velocity-Space Cascade: Hermite Representation and Theory

Plasma turbulence is investigated using unprecedented high-resolution ion velocity distribution measurements by the Magnetospheric Multiscale mission (MMS) in the Earth's magnetosheath. This novel observation of a highly structured particle distribution suggests a cascadelike process in velocity space. Complex velocity space structure is investigated using a three-dimensional Hermite transform, revealing, for the first time in observational data, a power-law distribution of moments. In analogy to hydrodynamics, a Kolmogorov approach leads directly to a range of predictions for this phase-space transport. The scaling theory is found to be in agreement with observations. The combined use of state-of-the-art MMS data sets, novel implementation of a Hermite transform method, and scaling theory of the velocity cascade opens new pathways to the understanding of plasma turbulence and the crucial velocity space features that lead to dissipation in plasmas.

Turbulence generation during the head-on collision of Alfvénic wave packets

The description of the Moffatt and Parker problem recently revisited by O. Pezzi et al. [Astrophys. J. 834, 166 (2017)1538-435710.3847/1538-4357/834/2/166] is here extended by analyzing the features of the turbulence produced by the interaction of two colliding Alfvénic wave packets in a kinetic plasma. Although the approach based on the presence of linear modes features is still helpful in characterizing some low-energy fluctuations, other signatures, which go beyond the pure linear modes analysis, are recovered, such as the significant weakening of clear dispersion relations and the production of zero frequency fluctuations.

Energy transfer channels and turbulence cascade in Vlasov-Maxwell turbulence

Analysis of the Vlasov-Maxwell equations from the perspective of turbulence cascade clarifies the role of electromagnetic work, and reveals the importance of the pressure-strain relation in generating internal energy. Particle-in-cell simulation demonstrates the relative importance of the several energy exchange terms, indicating that the traceless pressure-strain interaction "Pi-D" is of particular importance for both electrons and protons. The Pi-D interaction and the second tensor invariants of the strain are highly localized in similar spatial regions, indicating that energy transfer occurs preferentially in coherent structures. The collisionless turbulence cascade may be fruitfully explored by study of these energy transfer channels, in addition to examining transfer across spatial scales.

Explosive Particle Dispersion in Plasma Turbulence

Particle dynamics are investigated in plasma turbulence, using self-consistent kinetic simulations, in two dimensions. In the steady state, the trajectories of single protons and proton pairs are studied, at different values of plasma β (ratio between kinetic and magnetic pressure). For single-particle displacements, results are consistent with fluids and magnetic field line dynamics, where particles undergo normal diffusion for very long times, with higher β's being more diffusive. In an intermediate time range, with separations lying in the inertial range, particles experience an explosive dispersion in time, consistent with the Richardson prediction. These results, obtained for the first time with a self-consistent kinetic model, are relevant for astrophysical and laboratory plasmas, where turbulence is crucial for heating, mixing, and acceleration processes.

Ensemble Space-Time Correlation of Plasma Turbulence in the Solar Wind

Single point measurement turbulence cannot distinguish variations in space and time. We employ an ensemble of one- and two-point measurements in the solar wind to estimate the space-time correlation function in the comoving plasma frame. The method is illustrated using near Earth spacecraft observations, employing ACE, Geotail, IMP-8, and Wind data sets. New results include an evaluation of both correlation time and correlation length from a single method, and a new assessment of the accuracy of the familiar frozen-in flow approximation. This novel view of the space-time structure of turbulence may prove essential in exploratory space missions such as Solar Probe Plus and Solar Orbiter for which the frozen-in flow hypothesis may not be a useful approximation.

Intermittency, nonlinear dynamics and dissipation in the solar wind and astrophysical plasmas

An overview is given of important properties of spatial and temporal intermittency, including evidence of its appearance in fluids, magnetofluids and plasmas, and its implications for understanding of heliospheric plasmas. Spatial intermittency is generally associated with formation of sharp gradients and coherent structures. The basic physics of structure generation is ideal, but when dissipation is present it is usually concentrated in regions of strong gradients. This essential feature of spatial intermittency in fluids has been shown recently to carry over to the realm of kinetic plasma, where the dissipation function is not known from first principles. Spatial structures produced in intermittent plasma influence dissipation, heating, and transport and acceleration of charged particles. Temporal intermittency can give rise to very long time correlations or a delayed approach to steady-state conditions, and has been associated with inverse cascade or quasi-inverse cascade systems, with possible implications for heliospheric prediction.

Anisotropy in solar wind plasma turbulence

A review of spectral anisotropy and variance anisotropy for solar wind fluctuations is given, with the discussion covering inertial range and dissipation range scales. For the inertial range, theory, simulations and observations are more or less in accord, in that fluctuation energy is found to be primarily in modes with quasi-perpendicular wavevectors (relative to a suitably defined mean magnetic field), and also that most of the fluctuation energy is in the vector components transverse to the mean field. Energy transfer in the parallel direction and the energy levels in the parallel components are both relatively weak. In the dissipation range, observations indicate that variance anisotropy tends to decrease towards isotropic levels as the electron gyroradius is approached; spectral anisotropy results are mixed. Evidence for and against wave interpretations and turbulence interpretations of these features will be discussed. We also present new simulation results concerning evolution of variance anisotropy for different classes of initial conditions, each with typical background solar wind parameters.

Intermittent Dissipation and Heating in 3D Kinetic Plasma Turbulence

High resolution, fully kinetic, three dimensional (3D) simulation of collisionless plasma turbulence shows the development of turbulence characterized by sheetlike current density structures spanning a range of scales. The nonlinear evolution is initialized with a long wavelength isotropic spectrum of fluctuations having polarizations transverse to an imposed mean magnetic field. We present evidence that these current sheet structures are sites for heating and dissipation, and that stronger currents signify higher dissipation rates. The analyses focus on quantities such as J·E, electron, and proton temperatures, and conditional averages of these quantities based on local electric current density. Evidently, kinetic scale plasma, like magnetohydrodynamics, becomes intermittent due to current sheet formation, leading to the expectation that heating and dissipation in astrophysical and space plasmas may be highly nonuniform. Comparison with previous results from 2D kinetic simulations, as well as high frequency solar wind observational data, are discussed.

Magnetic field reversals and long-time memory in conducting flows

Employing a simple ideal magnetohydrodynamic model in spherical geometry, we show that the presence of either rotation or finite magnetic helicity is sufficient to induce dynamical reversals of the magnetic dipole moment. The statistical character of the model is similar to that of terrestrial magnetic field reversals, with the similarity being stronger when rotation is present. The connection between long-time correlations, 1/f noise, and statistics of reversals is supported, consistent with earlier suggestions.

Proton kinetic effects and turbulent energy cascade rate in the solar wind

The first observed connection between kinetic instabilities driven by proton temperature anisotropy and estimated energy cascade rates in the turbulent solar wind is reported using measurements from the Wind spacecraft at 1 AU. We find enhanced cascade rates are concentrated along the boundaries of the (β∥, T⊥/T∥) plane, which includes regions theoretically unstable to the mirror and firehose instabilities. A strong correlation is observed between the estimated cascade rate and kinetic effects such as temperature anisotropy and plasma heating, resulting in protons 5-6 times hotter and 70%-90% more anisotropic than under typical isotropic plasma conditions. These results offer new insights into kinetic processes in a turbulent regime.

Von Kármán energy decay and heating of protons and electrons in a kinetic turbulent plasma

Decay in time of undriven weakly collisional kinetic plasma turbulence in systems large compared to the ion kinetic scales is investigated using fully electromagnetic particle-in-cell simulations initiated with transverse flow and magnetic disturbances, constant density, and a strong guide field. The observed energy decay is consistent with the von Kármán hypothesis of similarity decay, in a formulation adapted to magnetohydrodyamics. Kinetic dissipation occurs at small scales, but the overall rate is apparently controlled by large scale dynamics. At small turbulence amplitudes the electrons are preferentially heated. At larger amplitudes proton heating is the dominant effect. In the solar wind and corona the protons are typically hotter, suggesting that these natural systems are in the large amplitude turbulence regime.

Inhomogeneous kinetic effects related to intermittent magnetic discontinuities

A connection between kinetic processes and two-dimensional intermittent plasma turbulence is observed using direct numerical simulations of a hybrid Vlasov-Maxwell model, in which the Vlasov equation is solved for protons, while the electrons are described as a massless fluid. During the development of turbulence, the proton distribution functions depart from the typical configuration of local thermodynamic equilibrium, displaying statistically significant non-Maxwellian features. In particular, temperature anisotropy and distortions are concentrated near coherent structures, generated as the result of the turbulent cascade, such as current sheets, which are nonuniformly distributed in space. Here, the partial variance of increments (PVI) method has been employed to identify high magnetic stress regions within a two-dimensional turbulent pattern. A quantitative association between non-Maxwellian features and coherent structures is established.

Intermittent dissipation at kinetic scales in collisionless plasma turbulence

High resolution kinetic simulations of collisionless plasma driven by shear show the development of turbulence characterized by dynamic coherent sheetlike current density structures spanning a range of scales down to electron scales. We present evidence that these structures are sites for heating and dissipation, and that stronger current structures signify higher dissipation rates. Evidently, kinetic scale plasma, like magnetohydrodynamics, becomes intermittent due to current sheet formation, leading to the expectation that heating and dissipation in astrophysical and space plasmas may be highly nonuniform and patchy.

Magnetic moment nonconservation in magnetohydrodynamic turbulence models

The fundamental assumptions of the adiabatic theory do not apply in the presence of sharp field gradients or in the presence of well-developed magnetohydrodynamic turbulence. For this reason, in such conditions the magnetic moment μ is no longer expected to be constant. This can influence particle acceleration and have considerable implications in many astrophysical problems. Starting with the resonant interaction between ions and a single parallel propagating electromagnetic wave, we derive expressions for the magnetic moment trapping width Δμ (defined as the half peak-to-peak difference in the particle magnetic moments) and the bounce frequency ω(b). We perform test-particle simulations to investigate magnetic moment behavior when resonance overlapping occurs and during the interaction of a ring-beam particle distribution with a broadband slab spectrum. We find that the changes of magnetic moment and changes of pitch angle are related when the level of magnetic fluctuations is low, δB/B(0) = (10(-3),10(-2)), where B(0) is the constant and uniform background magnetic field. Stochasticity arises for intermediate fluctuation values and its effect on pitch angle is the isotropization of the distribution function f(α). This is a transient regime during which magnetic moment distribution f(μ) exhibits a characteristic one-sided long tail and starts to be influenced by the onset of spatial parallel diffusion, i.e., the variance <(Δz)(2)> grows linearly in time as in normal diffusion. With strong fluctuations f(α) becomes completely isotropic, spatial diffusion sets in, and the f(μ) behavior is closely related to the sampling of the varying magnetic field associated with that spatial diffusion.

Kinetic signatures and intermittent turbulence in the solar wind plasma

A connection between kinetic processes and intermittent turbulence is observed in the solar wind plasma using measurements from the Wind spacecraft at 1 A.U. In particular, kinetic effects such as temperature anisotropy and plasma heating are concentrated near coherent structures, such as current sheets, which are nonuniformly distributed in space. Furthermore, these coherent structures are preferentially found in plasma unstable to the mirror and firehose instabilities. The inhomogeneous heating in these regions, which is present in both the magnetic field parallel and perpendicular temperature components, results in protons at least 3-4 times hotter than under typical stable plasma conditions. These results offer a new understanding of kinetic processes in a turbulent regime, where linear Vlasov theory is not sufficient to explain the inhomogeneous plasma dynamics operating near non-Gaussian structures.

Intermittency and local heating in the solar wind

Evidence for nonuniform heating in the solar wind plasma near current sheets dynamically generated by magnetohydrodynamic (MHD) turbulence is obtained using measurements from the ACE spacecraft. These coherent structures only constitute 19% of the data, but contribute 50% of the total plasma internal energy. Intermittent heating manifests as elevations in proton temperature near current sheets, resulting in regional heating and temperature enhancements extending over several hours. The number density of non-Gaussian structures is found to be proportional to the mean proton temperature and solar wind speed. These results suggest magnetofluid turbulence drives intermittent dissipation through a hierarchy of coherent structures, which collectively could be a significant source of coronal and solar wind heating.

Anisotropic third-moment estimates of the energy cascade in solar wind turbulence using multispacecraft data

The first direct determination of the inertial range energy cascade rate, using an anisotropic form of Yaglom's law for magnetohydrodynamic turbulence, is obtained in the solar wind with multispacecraft measurements. The two-point mixed third-order structure functions of Elsässer fluctuations are integrated over a sphere in magnetic field-aligned coordinates, and the result is consistent with a linear scaling. Therefore, volume integrated heating and cascade rates are obtained that, unlike previous studies, make only limited assumptions about the underlying spectral geometry of solar wind turbulence. These results confirm the turbulent nature of magnetic and velocity field fluctuations in the low frequency limit, and could supply the energy necessary to account for the nonadiabatic heating of the solar wind.

Emergence of very long time fluctuations and 1/f noise in ideal flows

This paper shows the connection between three previously observed but seemingly unrelated phenomena in hydrodynamic (HD) and magnetohydrodynamic (MHD) turbulent flows, involving the emergence of fluctuations occurring on very long time scales: the low-frequency 1/f noise in the power frequency spectrum, the delayed ergodicity of complex valued amplitude fluctuations in wave number space, and the spontaneous flippings or reversals of large-scale fields. Direct numerical simulations of ideal MHD and HD are employed in three space dimensions, at low resolution, for long periods of time, and with high accuracy to study several cases: different geometries, presence of rotation and/or a uniform magnetic field, and different values of the associated conserved global quantities. It is conjectured that the origin of all these long-time phenomena is rooted in the interaction of the longest wavelength fluctuations available to the system, with fluctuations at much smaller scales. The strength of this nonlocal interaction is controlled either by the existence of conserved global quantities with a back-transfer in Fourier space or by the presence of a slow manifold in the dynamics.

Large-scale behavior and statistical equilibria in rotating flows

We examine long-time properties of the ideal dynamics of three-dimensional flows, in the presence or not of an imposed solid-body rotation and with or without helicity (velocity-vorticity correlation). In all cases, the results agree with the isotropic predictions stemming from statistical mechanics. No accumulation of excitation occurs in the large scales, although, in the dissipative rotating case, anisotropy and accumulation, in the form of an inverse cascade of energy, are known to occur. We attribute this latter discrepancy to the linearity of the term responsible for the emergence of inertial waves. At intermediate times, inertial energy spectra emerge that differ somewhat from classical wave-turbulence expectations and with a trace of large-scale excitation that goes away for long times. These results are discussed in the context of partial two dimensionalization of the flow undergoing strong rotation as advocated by several authors.

Similarity decay of enstrophy in an electron fluid

A similarity decay law is proposed for enstrophy of a one-signed-vorticity fluid in a circular free-slip domain. It excludes the metastable equilibrium enstrophy which cannot drive turbulence, and approaches Batchelor's t(-2) law for strong turbulence. Measurements of the decay of a turbulent electron fluid agree well with the predictions of the decay law for a variety of initial conditions.

Waiting-time distributions of magnetic discontinuities: clustering or Poisson process?

Using solar wind data from the Advanced Composition Explorer spacecraft, with the support of Hall magnetohydrodynamic simulations, the waiting-time distributions of magnetic discontinuities have been analyzed. A possible phenomenon of clusterization of these discontinuities is studied in detail. We perform a local Poisson's analysis in order to establish if these intermittent events are randomly distributed or not. Possible implications about the nature of solar wind discontinuities are discussed.

Hydrodynamic relaxation of an electron plasma to a near-maximum entropy state

Dynamical relaxation of a pure electron plasma in a Malmberg-Penning trap is studied, comparing experiments, numerical simulations and statistical theories of weakly dissipative two-dimensional (2D) turbulence. Simulations confirm that the dynamics are approximated well by a 2D hydrodynamic model. Statistical analysis favors a theoretical picture of relaxation to a near-maximum entropy state with constrained energy, circulation, and angular momentum. This provides evidence that 2D electron fluid relaxation in a turbulent regime is governed by principles of maximum entropy.

Magnetic reconnection in two-dimensional magnetohydrodynamic turbulence

Systematic analysis of numerical simulations of two-dimensional magnetohydrodynamic turbulence reveals the presence of a large number of X-type neutral points where magnetic reconnection occurs. We examine the statistical properties of this ensemble of reconnection events that are spontaneously generated by turbulence. The associated reconnection rates are distributed over a wide range of values and scales with the geometry of the diffusion region. Locally, these events can be described through a variant of the Sweet-Parker model, in which the parameters are externally controlled by turbulence. This new perspective on reconnection is relevant in space and astrophysical contexts, where plasma is generally in a fully turbulent regime.

Anisotropic magnetohydrodynamic spectral transfer in the diffusion approximation

A theoretical model of spectral transfer for anisotropic magnetohydrodynamic (MHD) turbulence is introduced, approximating energy transport in wave vector (k) space as a nonlinear diffusion process, extending previous isotropic k -space diffusion theories for hydrodynamics and MHD. This formal closure at the spectral equation level may be useful in space and astrophysical applications.

Ergodicity of ideal Galerkin three-dimensional magnetohydrodynamics and Hall magnetohydrodynamics models

We explore the problem of the ergodicity of magnetohydrodynamics and Hall magnetohydrodynamics in three-dimensional, ideal Galerkin systems that are truncated to a finite number of Fourier modes. We show how single Fourier modes follow the Gibbs ensemble prediction, and how the ergodicity of the phase space is restored for long-time Galerkin solutions. Running time averages and two-time correlation functions show, at long times, a convergence towards zero of time averaged single Fourier modes. This suggests a delayed approach to, rather than a breaking of, ergodicity. Finally, we present some preliminary ideas concerning the origin of the associated time scales.

Depression of nonlinearity in decaying isotropic MHD turbulence

Spectral method simulations show that undriven magnetohydrodynamic turbulence spontaneously generates coherent spatial correlations of several types, associated with local Beltrami fields, directional alignment of velocity and magnetic fields, and antialignment of magnetic and fluid acceleration components. These correlations suppress nonlinearity to levels lower than what is obtained from Gaussian fields, and occur in spatial patches. We suggest that this rapid relaxation leads to non-Gaussian statistics and spatial intermittency.

Rapid alignment of velocity and magnetic field in magnetohydrodynamic turbulence

We show that local directional alignment of the velocity and magnetic field fluctuations occurs rapidly in magnetohydrodynamics for a variety of parameters and is seen both in direct numerical simulations and in solar wind data. The phenomenon is due to an alignment between magnetic field and gradients of either pressure or kinetic energy, and is similar to alignment of velocity and vorticity in Navier-Stokes turbulence. This rapid and robust relaxation process leads to a local weakening of nonlinear terms.

Low-frequency 1/f fluctuations in hydrodynamic and magnetohydrodynamic turbulence

We investigate the occurrence of 1/f spectra of low-frequency fluctuations in numerical simulations of three-dimensional hydrodynamic and magnetohydrodynamic turbulence driven by a random forcing with a controlled correlation time. A range of one decade of 1/f spectrum is observed when a strong background magnetic field is present. The frequency spectra of individual Fourier modes is also analyzed and it is observed that the 1/f range is present in the largest available wavelength mode for the magnetohydrodynamic simulations with and without a background magnetic field and it is not observed (or is less clear) for the hydrodynamic case. The presence of 1/f spectra of low-frequency fluctuations is also analyzed for two-dimensional magnetohydrodynamic and hydrodynamic turbulence simulations and it is observed in both cases. The origin of these long period fluctuations is discussed.

Spatial correlation of solar-wind turbulence from two-point measurements

Interplanetary turbulence, the best studied case of low frequency plasma turbulence, is the only directly quantified instance of astrophysical turbulence. Here, magnetic field correlation analysis, using for the first time only proper two-point, single time measurements, provides a key step in unraveling the space-time structure of interplanetary turbulence. Simultaneous magnetic field data from the Wind, ACE, and Cluster spacecraft are analyzed to determine the correlation (outer) scale, and the Taylor microscale near Earth's orbit.

Spectral distribution of the cross helicity in the solar wind

There are a variety of theoretical and observational indications that fluctuation energy in astrophysical and space plasma turbulence is distributed anisotropically in space relative to the magnetic field direction. The cross helicity, represented by correlations between velocity and magnetic field fluctuations, enters a magnetohydrodynamic description on equal footing with the energy, but its anisotropy has not been examined in the same degree of detail. Here we employ Advanced Coronal Explorer data to examine the rotational symmetry of the cross helicity. We find that the normalized cross helicity is associated more or less equally with all angular components of the fluctuations. This favors turbulence models that allow for cross communication between parallel and perpendicular wave numbers, suggesting that "wavelike" and "turbulencelike" fluctuations are strongly coupled.

One-point statistics of the induced electric field in quasinormal magnetofluid turbulence

We study one-point statistical properties of the induced turbulent electric field for a magnetohydrodynamic (MHD) plasma under the quasinormal approximation. Assuming exact Gaussianity for both the velocity field and the magnetic field, and different degrees of correlations between their Cartesian components, we derive the probability distribution function (PDF) for the Cartesian components of the electric field e(i). We show that the PDF reduces in some canonical cases to an exponential function of the form exp(-/e(i)/). To study deviations from these results in the more realistic case in which the velocity and magnetic fields are not exactly normal but quasinormal instead, we perform three-dimensional numerical simulations of the MHD equations at moderate Reynolds numbers. For turbulent relaxation from an initial condition, we find that the analytical results give a very good first-order approximation to the computed PDF.