Source author record

W. H. Matthaeus

W. H. Matthaeus appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2025arXiv

Factors Controlling the Statistics of Magnetic Reconnection in MHD Turbulence

We study the statistics of dynamical quantities associated with magnetic reconnection events embedded in a sea of strong background magnetohydrodynamic (MHD) turbulence using direct numerical simulations. We focus on the relationship of the reconnection properties to the statistics of global turbulent fields. For the first time, we show that the distribution in turbulence of reconnection rates (determined by upstream fields) is strongly correlated with the magnitude of the global turbulent magnetic field at the correlation scale. The average reconnection rates, and associated dissipation rates, during turbulence are thus much larger than predicted by using turbulent magnetic field fluctuation amplitudes at the dissipation or kinetic scales. Magnetic reconnection may therefore be playing a major role in energy dissipation in astrophysical and heliospheric turbulence.

preprint2022arXiv

Energetic Particle Perpendicular Diffusion: Simulations and Theory in Noisy Reduced Magnetohydrodynamic Turbulence

The transport of energetic charged particles (e.g., cosmic rays) in turbulent magnetic fields is usually characterized in terms of the diffusion parallel and perpendicular to a large-scale (or mean) magnetic field. The nonlinear guiding center theory (NLGC) has been a prominent perpendicular diffusion theory. A recent version of this theory, based on random ballistic spreading of magnetic field lines and a backtracking correction (RBD/BC), has shown good agreement with test particle simulations for a two-component magnetic turbulence model. The aim of the present study is to test the generality of the improved theory by applying it to the noisy reduced magnetohydrodynamic (NRMHD) turbulence model, determining perpendicular diffusion coefficients that are compared with those from the field line random walk (FLRW) and unified nonlinear (UNLT) theories and our test particle simulations. The synthetic NRMHD turbulence model creates special conditions for energetic particle transport, with no magnetic fluctuations at higher parallel wavenumbers so there is no resonant parallel scattering if the particle Larmor radius $R_{\rm L}$ is even slightly smaller than the minimum resonant scale. This leads to non-monotonic variation in the parallel mean free path $λ_\parallel$ with $R_{\rm L}$. Among the theories considered, only RBD/BC matches simulations within a factor of two over the range of parameters considered. This accuracy is obtained even though the theory depends on $λ_\parallel$ and has no explicit dependence on $R_{\rm L}$. In addition, the UNLT theory often provides accurate results and even the FLRW limit provides a very simple and reasonable approximation in many cases.

preprint2022arXiv

Sub-Alfvenic Solar Wind observed by PSP: Characterization of Turbulence, Anisotropy, Intermittency, and Switchback

In the lower solar coronal regions where the magnetic field is dominant, the Alfven speed is much higher than the wind speed. In contrast, the near-Earth solar wind is strongly super-Alfvenic, i.e., the wind speed greatly exceeds the Alfven speed. The transition between these regimes is classically described as the "Alfven point" but may in fact occur in a distributed Alfven critical region. NASA's Parker Solar Probe (PSP) mission has entered this region, as it follows a series of orbits that gradually approach more closely to the sun. During its 8th and 9th solar encounters, at a distance of 16 solar radii from the Sun, PSP sampled four extended periods in which the solar wind speed was measured to be smaller than the local Alfven speed. These are the first in-situ detections of sub-Alfvenic solar wind in the inner heliosphere by PSP. Here we explore properties of these samples of sub-Alfvenic solar wind, which may provide important previews of the physical processes operating at lower altitude. Specifically, we characterize the turbulence, anisotropy, intermittency, and directional switchback properties of these sub-Alfvenic winds and contrast these with the neighboring super-Alfvenic periods.

preprint2021arXiv

Magnetic reconnection as an energy cascade process

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 cascade 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 itself is fundamentally an energy cascade process.

preprint2021arXiv

PSP/IS$\odot$IS Observation of a Solar Energetic Particle Event Associated With a Streamer Blowout Coronal Mass Ejection During Encounter 6

In this paper we examine a low-energy SEP event observed by IS$\odot$IS's Energetic Particle Instrument-Low (EPI-Lo) inside 0.18 AU on September 30, 2020. This small SEP event has a very interesting time profile and ion composition. Our results show that the maximum energy and peak in intensity is observed mainly along the open radial magnetic field. The event shows velocity dispersion, and strong particle anisotropies are observed throughout the event showing that more particles are streaming outward from the Sun. We do not see a shock in the in-situ plasma or magnetic field data throughout the event. Heavy ions, such as O and Fe were detected in addition to protons and 4He, but without significant enhancements in 3He or energetic electrons. Our analysis shows that this event is associated with a slow streamer-blowout coronal mass ejection (SBO-CME) and the signatures of this small CME event are consistent with those typical of larger CME events. The time-intensity profile of this event shows that PSP encountered the western flank of the SBO-CME. The anisotropic and dispersive nature of this event in a shockless local plasma give indications that these particles are most likely accelerated remotely near the Sun by a weak shock or compression wave ahead of the SBO-CME. This event may represent direct observations of the source of low-energy SEP seed particle population.

preprint2021arXiv

Suprathermal Ion Energy spectra and Anisotropies near the Heliospheric Current Sheet crossing observed by the Parker Solar Probe during Encounter 7

We present observations of >10-100 keV/nucleon suprathermal (ST) H, He, O, and Fe ions associated with crossings of the heliospheric current sheet (HCS) at radial distances <0.1 au from the Sun. Our key findings are: 1) very few heavy ions are detected during the 1st full crossing, the heavy ion intensities are reduced during the 2nd partial crossing and peak just after the 2nd crossing; 2) ion arrival times exhibit no velocity dispersion; 3) He pitch-angle distributions track the magnetic field polarity reversal and show up to ~10:1 anti-sunward, field-aligned flows and beams closer to the HCS that become nearly isotropic further from the HCS; 4) the He spectrum steepens either side of the HCS and the He, O, and Fe spectra exhibit power-laws of the form ~E^4-6; and 5) maximum energies EX increase with the ion's charge-to-mass (Q/M) ratio as EX/EH proportional to [(QX/MX)]^alpha where alpha~0.65-0.76, assuming that the average Q-states are similar to those measured in gradual and impulsive solar energetic particle events at 1 au. The absence of velocity dispersion in combination with strong field-aligned anisotropies closer to the HCS appears to rule out solar flares and near-sun coronal mass ejection-driven shocks. These new observations present challenges not only for mechanisms that employ direct parallel electric fields and organize maximum energies according to E/Q, but also for local diffusive and magnetic reconnection-driven acceleration models. Re-evaluation of our current understanding of the production and transport of energetic ions is necessary to understand this near-solar, current-sheet-associated population of ST ions.

preprint2020arXiv

Energetic Particle Increases Associated with Stream Interaction Regions

The Parker Solar Probe was launched on 2018 August 12 and completed its second orbit on 2019 June 19 with perihelion of 35.7 solar radii. During this time, the Energetic particle Instrument-Hi (EPI-Hi, one of the two energetic particle instruments comprising the Integrated Science Investigation of the Sun, ISOIS) measured seven proton intensity increases associated with stream interaction regions (SIRs), two of which appear to be occurring in the same region corotating with the Sun. The events are relatively weak, with observed proton spectra extending to only a few MeV and lasting for a few days. The proton spectra are best characterized by power laws with indices ranging from -4.3 to -6.5, generally softer than events associated with SIRs observed at 1 au and beyond. Helium spectra were also obtained with similar indices, allowing He/H abundance ratios to be calculated for each event. We find values of 0.016-0.031, which are consistent with ratios obtained previously for corotating interaction region events with fast solar wind < 600 km s-1. Using the observed solar wind data combined with solar wind simulations, we study the solar wind structures associated with these events and identify additional spacecraft near 1 au appropriately positioned to observe the same structures after some corotation. Examination of the energetic particle observations from these spacecraft yields two events that may correspond to the energetic particle increases seen by EPI-Hi earlier.

preprint2020arXiv

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.

preprint2020arXiv

Reconnection from a Turbulence Perspective

The spectral properties associated with laminar, anti-parallel reconnection are examined using a 2.5D kinetic particle in cell (PIC) simulation. Both the reconnection rate and the energy spectrum exhibit three distinct phases: an initiation phase where the reconnection rate grows, a quasi-steady phase, and a declining phase where both the reconnection rate and the energy spectrum decrease. During the steady phase, the energy spectrum exhibits approximately a double power-law behavior, with a slope near -5/3 at wavenumbers smaller than the inverse ion inertial length, and a slope steeper than -8/3 for larger wavenumbers up to the inverse electron inertial length. This behavior is consistent with a Kolmogorov energy cascade and implies that laminar reconnection may fundamentally be an energy cascade process. Consistent with this idea is that the reconnection rate exhibits a rough correlation with the energy spectrum at wave numbers near the inverse ion inertial length. The 2D spectrum is strongly anisotropic with most energy associated with the wave vector direction normal to the current sheet. Reconnection acts to isotropize the energy spectrum, reducing the Shebalin angle from an initial value of 70 degrees to about 48 degrees (nearly isotropic) by the end of the simulation.

preprint2019arXiv

CME -Associated Energetic Ions at 0.23 AU -- Consideration of the Auroral Pressure Cooker Mechanism Operating in the Low Corona as a Possible Energization Process

We draw a comparison between a solar energetic particle event associated with the release of a slow coronal mass ejection close to the sun, and the energetic particle population produced in high current density field-aligned current structures associated with auroral phenomena in planetary magnetospheres. We suggest that this process is common in CME development and lift-off in the corona, and may account for the electron populations that generate Type III radio bursts, as well as for the prompt energetic ion and electron populations typically observed in interplanetary space.

preprint2019arXiv

Energetic Particle Observations from Parker Solar Probe using Combined Energy Spectra from the IS$\odot$IS Instrument Suite

The Integrated Science Investigations of the Sun (IS$\odot$IS) instrument suite includes two Energetic Particle instruments: EPI-Hi, designed to measure ions from ~1-200 MeV/nuc, and EPI-Lo, designed to measure ions from ~20 keV/nuc to ~15 MeV/nuc. We present an analysis of eight energetic proton events observed across the energy range of both instruments during PSP's first two orbits in order to examine their combined energy spectra. Background corrections are applied to help resolve spectral breaks between the two instruments and are shown to be effective. In doing so we demonstrate that, even in the early stages of calibration, IS$\odot$IS is capable of producing reliable spectral observations across broad energy ranges. In addition to making groundbreaking measurements very near the Sun, IS$\odot$IS also characterizes energetic particle populations over a range of heliocentric distances inside 1 au. During the first two orbits, IS$\odot$IS observed energetic particle events from a single corotating interaction region (CIR) at three different distances from the Sun. The events are separated by two Carrington rotations and just 0.11 au in distance, however the relationship shown between proton intensities and proximity of the spacecraft to the source region shows evidence of the importance of transport effects on observations of energetic particles from CIRs. Future IS$\odot$IS observations of similar events over larger distances will help disentangle the effects of CIR-related acceleration and transport. We apply similar spectral analyses to the remaining five events, including four that are likely related to stream interaction regions (SIRs) and one solar energetic particle (SEP) event.

preprint2019arXiv

Enhanced Energy Transfer Rate in Solar Wind Turbulence Observed near the Sun from Parker Solar Probe

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ármán 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.

preprint2019arXiv

Observations of Energetic-Particle Population Enhancements along Intermittent Structures near the Sun from Parker Solar Probe

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.

preprint2019arXiv

Observations of heating along intermittent structures in the inner heliosphere from PSP data

The solar wind proton temperature at 1-au has been found to be correlated with small-scale intermittent magnetic structures, i.e., regions with enhanced temperature are associated with coherent structures such as current sheets. Using Parker Solar Probe data from the first encounter, we study this association using measurements of radial proton temperature, employing the Partial Variance of Increments (PVI) technique to identify intermittent magnetic structures. We observe that the probability density functions of high-PVI events have higher median temperatures than those with lower PVI, The regions in space where PVI peaks were also locations that had enhanced temperatures when compared with similar regions suggesting a heating mechanism in the young solar wind that is associated with intermittency developed by a nonlinear turbulent cascade.n the immediate vicinity.

preprint2019arXiv

Observations of the 2019 April 4 Solar Energetic Particle Event at the Parker Solar Probe

A solar energetic particle event was detected by the Integrated Science Investigation of the Sun (ISOIS) instrument suite on Parker Solar Probe (PSP) on 2019 April 4 when the spacecraft was inside of 0.17 au and less than 1 day before its second perihelion, providing an opportunity to study solar particle acceleration and transport unprecedentedly close to the source. The event was very small, with peak 1 MeV proton intensities of ~0.3 particles (cm^2 sr s MeV)^-1, and was undetectable above background levels at energies above 10 MeV or in particle detectors at 1 au. It was strongly anisotropic, with intensities flowing outward from the Sun up to 30 times greater than those flowing inward persisting throughout the event. Temporal association between particle increases and small brightness surges in the extreme-ultraviolet observed by the Solar TErrestrial RElations Observatory, which were also accompanied by type III radio emission seen by the Electromagnetic Fields Investigation on PSP, indicates that the source of this event was an active region nearly 80 degrees east of the nominal PSP magnetic footpoint. This suggests that the field lines expanded over a wide longitudinal range between the active region in the photosphere and the corona.

preprint2019arXiv

Seed Population Pre-Conditioning and Acceleration Observed by Parker Solar Probe

A series of solar energetic particle (SEP) events were observed at Parker Solar Probe (PSP) by the Integrated Science Investigation of the Sun (\ISOIS) during the period from April 18, 2019 through April 24, 2019. The PSP spacecraft was located near 0.48 au from the Sun on Parker spiral field lines that projected out to 1 au within $\sim 25^\circ$ of near Earth spacecraft. These SEP events, though small compared to historically large SEP events, were amongst the largest observed thus far in the PSP mission and provide critical information about the space environment inside 1 au during SEP events. During this period the Sun released multiple coronal mass ejections (CMEs). One of these CMEs observed was initiated on April 20, 2019 at 01:25 UTC, and the interplanetary CME (ICME) propagated out and passed over the PSP spacecraft. Observations by the Electromagnetic Fields Investigation (FIELDS) show that the magnetic field structure was mostly radial throughout the passage of the compression region and the plasma that followed, indicating that PSP did not directly observe a flux rope internal to the ICME, consistent with the location of PSP on the ICME flank. Analysis using relativistic electrons observed near Earth by the Electron, Proton and Alpha Monitor (EPAM) on the Advanced Composition Explorer (ACE) demonstrates the presence of electron seed populations (40--300 keV) during the events observed. The energy spectrum of the \ISOIS~ observed proton seed population below 1 MeV is close to the limit of possible stationary state plasma distributions out of equilibrium. \ISOIS~ observations reveal the \revise{enhancement} of seed populations during the passage of the ICME, which \revise{likely indicates a key part} of the pre-acceleration process that occurs close to the Sun.

preprint2019arXiv

Transition from ion-coupled to electron-only reconnection: Basic physics and implications for plasma turbulence

Using kinetic particle-in-cell (PIC) simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion Larmor radius, the ions do not respond to the reconnection dynamics leading to ''electron-only'' reconnection with very large quasi-steady reconnection rates. The transition to more traditional ''ion-coupled'' reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasi-steady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfvén speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths.

preprint2016arXiv

Explosive Particle Dispersion in Plasma Turbulence

Particle dynamics are investigated in plasma turbulence, using self-consistent kinetic simulations, in two dimensions. In steady state, the trajectories of single protons and proton-pairs are studied, at different values of plasma "beta" (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 "beta" 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.

preprint2016arXiv

Fading Coronal Structure and the Onset of Turbulence in the Young Solar Wind

Above the top of the solar corona, the young slow solar wind transitions from low-beta, magnetically structured flow dominated by radial structures, to high-beta, less structured flow dominated by hydrodynamics. This transition, long inferred via theory, is readily apparent in the sky region close to 10 degrees from the Sun, in processed, background-subtracted solar wind images. We present image sequences collected by the STEREO/HI1 instrument in 2008 Dec, covering apparent distances from approximately 4 to 24 degrees from the center of the Sun and spanning this transition in large-scale morphology of the wind. We describe the observation and novel techniques to extract evolving image structure from the images, and we use those data and techniques to present and quantify the clear textural shift in the apparent structure of the corona and solar wind in this altitude range. We demonstrate that the change in apparent texture is due both to anomalous fading of the radial striae that characterize the corona, and to anomalous relative brightening of locally dense puffs of solar wind that we term "flocculae;" and show that these phenomena are inconsistent with smooth radial flow, but consistent with onset of hydrodynamic or MHD instabilities leading to a turbulent cascade in the young solar wind.

preprint2016arXiv

On the compressibility effect in test particle acceleration by magnetohydrodynamic turbulence

The effect of compressibility in charged particle energization by magnetohydrodynamic (MHD) fields is studiedin the context of test particle simulations. This problem is relevant to the solar wind and the solar corona due to the compressible nature of the flow in those astrophysical scenarios. We consider turbulent electromagnetic fields obtained from direct numerical simulations of the MHD equations with a strong background magnetic field. In order to explore the flow compressibilty effect over the particle dynamics we performed different numerical experiments: an incompressible case, and two weak compressible cases with Mach number M = 0.1 and M = 0.25. We analyze the behavior of protons and electrons in those turbulent fields, which are well known to form aligned current sheets in the direction of the guide magnetic field. What we call protons and electrons are test particles with scales comparable to (for protons) and much smaller than (for electrons) the dissipative scale of MHD turbulence, maintaining the correct mass ratio me /mi. For these test particles we show that compressibility enhances the efficiency of proton acceleration, and that the energization is caused by perpendicular electric fields generated between currents sheets. On the other hand, electrons remain magnetized and display an almost adiabatic motion, with no effect of compressibility observed. Another set of numerical experiments takes into account two fluid modifications, namely electric field due to Hall effect and electron pressure gradient. We show that the electron pressure has an important contribution to electron acceleration allowing highly parallel energization. In contrast, no significant effect of these additional terms is observed for the protons.

preprint2016arXiv

On the spatio-temporal behavior of magnetohydrodynamic turbulence in a magnetized plasma

Using direct numerical simulations of three-dimensional magnetohydrodynamic (MHD) turbulence the spatio-temporal behavior of magnetic field fluctuations is analyzed. Cases with relatively small, medium and large values of a mean background magnetic field are considered. The (wavenumber) scale dependent time correlation function is directly computed for different simulations, varying the mean magnetic field value. From this correlation function the time decorrelation is computed and compared with different theoretical times, namely, the local non-linear time, the random sweeping time, and the Alfvénic time, the latter being a wave effect. It is observed that time decorrelations are dominated by sweeping effects, and only at large values of the mean magnetic field and for wave vectors mainly aligned with this field time decorrelations are controlled by Alfvénic effects.

preprint2015arXiv

Kinetic Alfvén waves generation by large-scale phase-mixing

One view of the solar-wind turbulence is that the observed highly anisotropic fluctuations at spatial scales near the proton inertial length $d_p$ may be considered as Kinetic Alfvén waves (KAWs). In the present paper, we show how phase-mixing of large-scale parallel propagating Alfvén waves is an efficient mechanism for the production of KAWs at wavelengths close to $d_p$ and at large propagation angle with respect to the magnetic field. Magnetohydrodynamic (MHD), Hall-Magnetohydrodynamic (HMHD), and hybrid Vlasov-Maxwell (HVM) simulations modeling the propagation of Alfvén waves in inhomogeneous plasmas are performed. In linear regime, the role of dispersive effects is singled out by comparing MHD and HMHD results. Fluctuations produced by phase-mixing are identified as KAWs through a comparison of polarization of magnetic fluctuations and wave group velocity with analytical linear predictions. In the nonlinear regime, comparison of HMHD and HVM simulations allows to point out the role of kinetic effects in shaping the proton distribution function. We observe generation of temperature anisotropy with respect to the local magnetic field and production of field-aligned beams. The regions where the proton distribution function highly departs from thermal equilibrium are located inside the shear layers, where the KAWs are excited, this suggesting that the distortions of the proton distribution are driven by a resonant interaction of protons with KAW fluctuations. Our results are relevant in configurations where magnetic field inhomogeneities are present, as, for example, in the solar corona where the presence of Alfvén waves has been ascertained.

preprint2014arXiv

Characterization of the Turbulent Magnetic Integral Length in the Solar Wind: From 0.3 to 5 Astronomical Units

The solar wind is a structured and complex system, in which the fields vary strongly over a wide range of spatial and temporal scales. As an example, the turbulent activity in the wind affects the evolution in the heliosphere of the integral turbulent scale or correlation length [λ], usually associated with the breakpoint in the turbulent-energy spectrum that separates the inertial range from the injection range. This large variability of the fields demands a statistical description of the solar wind. In this work, we study the probability distribution function (PDF) of the magnetic autocorrelation lengths observed in the solar wind at different distances from the Sun. We use observations from Helios, ACE, and Ulysses spacecraft. We distinguish between the usual solar wind and one of its transient components (Interplanetary Coronal Mass Ejections, ICMEs), and study also solar wind samples with low and high proton beta [β_p ]. We find that in the last 3 regimes the PDF of λ is a log-normal function, consistent with the multiplicative and non-linear processes that take place in the solar wind, the initial λ (before the Alfvénic point) being larger in ICMEs.

preprint2014arXiv

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.

preprint2014arXiv

Magnetic Reconnection and Intermittent Turbulence in the Solar Wind

A statistical relationship between magnetic reconnection, current sheets and intermittent turbulence in the solar wind is reported for the first time using in-situ measurements from the Wind spacecraft at 1 AU. We identify intermittency as non-Gaussian fluctuations in increments of the magnetic field vector, $\mathbf{B}$, that are spatially and temporally non-uniform. The reconnection events and current sheets are found to be concentrated in intervals of intermittent turbulence, identified using the partial variance of increments method: within the most non-Gaussian 1% of fluctuations in $\mathbf{B}$, we find 87%-92% of reconnection exhausts and $\sim$9% of current sheets. Also, the likelihood that an identified current sheet will also correspond to a reconnection exhaust increases dramatically as the least intermittent fluctuations are removed from the dataset. Hence, the turbulent solar wind contains a hierarchy of intermittent magnetic field structures that are increasingly linked to current sheets, which in turn are progressively more likely to correspond to sites of magnetic reconnection. These results could have far reaching implications for laboratory and astrophysical plasmas where turbulence and magnetic reconnection are ubiquitous.

preprint2014arXiv

Nonlinear and Linear Timescales near Kinetic Scales in Solar Wind Turbulence

The application of linear kinetic treatments to plasma waves, damping, and instability requires favorable inequalities between the associated linear timescales and timescales for nonlinear (e.g., turbulence) evolution. In the solar wind these two types of timescales may be directly compared using standard Kolmogorov-style analysis and observational data. The estimated local nonlinear magnetohydrodynamic cascade times, evaluated as relevant kinetic scales are approached, remain slower than the cyclotron period, but comparable to, or faster than, the typical timescales of instabilities, anisotropic waves, and wave damping. The variation with length scale of the turbulence timescales is supported by observations and simulations. On this basis the use of linear theory - which assumes constant parameters to calculate the associated kinetic rates - may be questioned. It is suggested that the product of proton gyrofrequency and nonlinear time at the ion gyroscales provides a simple measure of turbulence influence on proton kinetic behavior.

preprint2014arXiv

Quantification of the strength of inertial waves in a rotating turbulent flow

We quantify the strength of the waves and their impact on the energy cascade in rotating turbulence by studying the wave number and frequency energy spectrum, and the time correlation functions of individual Fourier modes in numerical simulations in three dimensions in periodic boxes. From the spectrum, we find that a significant fraction of the energy is concentrated in modes with wave frequency $ω\approx 0$, even when the external forcing injects no energy directly into these modes. However, for modes for which the period of the inertial waves $τ_ω$ is faster than the turnover time $τ_\textrm{NL}$, a significant fraction of the remaining energy is concentrated in the modes that satisfy the dispersion relation of the waves. No evidence of accumulation of energy in the modes with $τ_ω= τ_\textrm{NL}$ is observed, unlike what critical balance arguments predict. From the time correlation functions, we find that for modes with $τ_ω< τ_\textrm{sw}$ (with $τ_\textrm{sw}$ the sweeping time) the dominant decorrelation time is the wave period, and that these modes also show a slower modulation on the timescale $τ_\textrm{NL}$ as assumed in wave turbulence theories. The rest of the modes are decorrelated with the sweeping time, including the very energetic modes modes with $ω\approx 0$.

preprint2014arXiv

Test-particle acceleration in a hierarchical three-dimensional turbulence model

The acceleration of charged particles is relevant to the solar corona over a broad range of scales and energies. High-energy particles are usually detected in concomitance with large energy release events like solar eruptions and flares, nevertheless acceleration can occur at smaller scales, characterized by dynamical activity near current sheets. To gain insight into the complex scenario of coronal charged particle acceleration, we investigate the properties of acceleration with a test-particle approach using three-dimensional magnetohydrodynamic (MHD) models. These are obtained from direct solutions of the reduced MHD equations, well suited for a plasma embedded in a strong axial magnetic field, relevant to the inner heliosphere. A multi-box, multi-scale technique is used to solve the equations of motion for protons. This method allows us to resolve an extended range of scales present in the system, namely from the ion inertial scale of the order of a meter up to macroscopic scales of the order of $10\,$km ($1/100$th of the outer scale of the system). This new technique is useful to identify the mechanisms that, acting at different scales, are responsible for acceleration to high energies of a small fraction of the particles in the coronal plasma. We report results that describe acceleration at different stages over a broad range of time, length and energy scales.

preprint2013arXiv

Proton Kinetic Effects in Vlasov and Solar Wind Turbulence

Kinetic plasma processes have been investigated in the framework of solar wind turbulence, employing Hybrid Vlasov-Maxwell (HVM) simulations. The dependency of proton temperature anisotropy T_{\perp}/T_{\parallel} on the parallel plasma beta β_{\parallel}, commonly observed in spacecraft data, has been recovered using an ensemble of HVM simulations. By varying plasma parameters, such as plasma beta and fluctuation level, the simulations explore distinct regions of the parameter space given by T_{\perp}/T_{\parallel} and β_{\parallel}, similar to solar wind sub-datasets. Moreover, both simulation and solar wind data suggest that temperature anisotropy is not only associated with magnetic intermittent events, but also with gradient-type structures in the flow and in the density. This connection between non-Maxwellian kinetic effects and various types of intermittency may be a key point for understanding the complex nature of plasma turbulence.

preprint2013arXiv

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 (MHD). Kinetic dissipation occurs at small scales, but the overall rate is apparently controlled by large scale dynamics. At small turbulence amplitude 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 large amplitude turbulence regime.

preprint2012arXiv

Interchange reconnection in a turbulent Corona

Magnetic reconnection at the interface between coronal holes and loops, so-called interchange reconnection, can release the hotter, denser plasma from magnetically confined regions into the heliosphere, contributing to the formation of the highly variable slow solar wind. The interchange process is often thought to develop at the apex of streamers or pseudo-streamers, near Y and X-type neutral points, but slow streams with loop composition have been recently observed along fanlike open field lines adjacent to closed regions, far from the apex. However, coronal heating models, with magnetic field lines shuffled by convective motions, show that reconnection can occur continuously in unipolar magnetic field regions with no neutral points: photospheric motions induce a magnetohydrodynamic turbulent cascade in the coronal field that creates the necessary small scales, where a sheared magnetic field component orthogonal to the strong axial field is created locally and can reconnect. We propose that a similar mechanism operates near and around boundaries between open and closed regions inducing a continual stochastic rearrangement of connectivity. We examine a reduced magnetohydrodynamic model of a simplified interface region between open and closed corona threaded by a strong unipolar magnetic field. This boundary is not stationary, becomes fractal, and field lines change connectivity continuously, becoming alternatively open and closed. This model suggests that slow wind may originate everywhere along loop-coronal hole boundary regions, and can account naturally and simply for outflows at and adjacent to such boundaries and for the observed diffusion of slow wind around the heliospheric current sheet.

preprint2012arXiv

Intermittent Dissipation and Local Heating in the Solar Wind

Evidence for inhomogeneous heating in the interplanetary 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.

preprint2012arXiv

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 AU. In particular, kinetic effects such as temperature anisotropy and plasma heating are concentrated near coherent structures, such as current sheets, which are non-uniformly 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.

preprint2012arXiv

Local Anisotropy, Higher Order Statistics, and Turbulence Spectra

Correlation anisotropy emerges dynamically in magnetohydrodynamics (MHD), producing stronger gradients across the large-scale mean magnetic field than along it. This occurs both globally and locally, and has significant implications in space and astrophysical plasmas, including particle scattering and transport, and theories of turbulence. Properties of local correlation anisotropy are further documented here by showing through numerical experiments that the effect is intensified in more localized estimates of the mean field. The mathematical formulation of this property shows that local anisotropy mixes second-order with higher order correlations. Sensitivity of local statistical estimates to higher order correlations can be understood in connection with the stochastic coordinate system inherent in such formulations. We demonstrate this in specific cases, illustrate the connection to higher order statistics by showing the sensitivity of local anisotropy to phase randomization, and thus establish that the local structure function is not a measure of the energy spectrum. Evidently the local enhancement of correlation anisotropy is of substantial fundamental interest, and this phenomenon must be understood in terms of higher order correlations, fourth-order and above.

preprint2012arXiv

Magnetic moment non-conservation in magnetohydrodynamic turbulence models

The fundamental assumptions of the adiabatic theory do not apply in presence of sharp field gradients as well as in 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 moment) and the bounce frequency $ω_b$. We perform test-particle simulations to investigate magnetic moment behavior when resonances overlapping occurs and during the interaction of a ring-beam particle distribution with a broad-band slab spectrum. We find that magnetic moment dynamics is strictly related to pitch angle $α$ for a low level of magnetic fluctuation, $δ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(α)$ isotropizes completely, spatial diffusion sets in and $f(μ)$ behavior is closely related to the sampling of the varying magnetic field associated with that spatial diffusion.

preprint2011arXiv

Aging of anisotropy of solar wind magnetic fluctuations in the inner heliosphere

We analyze the evolution of the interplanetary magnetic field spatial structure by examining the inner heliospheric autocorrelation function, using Helios 1 and Helios 2 "in situ" observations. We focus on the evolution of the integral length scale (λ) anisotropy associated with the turbulent magnetic fluctuations, with respect to the aging of fluid parcels traveling away from the Sun, and according to whether the measured λis principally parallel (λ_parallel) or perpendicular (λ_perp) to the direction of a suitably defined local ensemble average magnetic field B0. We analyze a set of 1065 24-hour long intervals (covering full missions). For each interval, we compute the magnetic autocorrelation function, using classical single-spacecraft techniques, and estimate λwith help of two different proxies for both Helios datasets. We find that close to the Sun, λ_parallel < λ_perp. This supports a slab-like spectral model, where the population of fluctuations having wavevector k parallel to B0 is much larger than the one with k-vector perpendicular. A population favoring perpendicular k-vectors would be considered quasi-two dimensional (2D). Moving towards 1 AU, we find a progressive isotropization of λand a trend to reach an inverted abundance, consistent with the well-known result at 1 AU that λ_parallel > λ_perp, usually interpreted as a dominant quasi-2D picture over the slab picture. Thus, our results are consistent with driving modes having wavevectors parallel to B0 near Sun, and a progressive dynamical spectral transfer of energy to modes with perpendicular wavevectors as the solar wind parcels age while moving from the Sun to 1 AU.

preprint2011arXiv

Morphology, dynamics and plasma parameters of plumes and inter-plume regions in solar coronal holes

Coronal plumes, which extend from solar coronal holes (CH) into the high corona and - possibly - into the solar wind (SW), can now continuously be studied with modern telescopes and spectrometers on spacecraft, in addition to investigations from the ground, in particular, during total eclipses. Despite the large amount of data available on these prominent features and related phenomena, many questions remained unanswered as to their generation and relative contributions to the high-speed streams emanating from CHs. An understanding of the processes of plume formation and evolution requires a better knowledge of the physical conditions at the base of CHs, in plumes and in the surrounding inter-plume regions (IPR). More specifically, information is needed on the magnetic field configuration, the electron densities and temperatures, effective ion temperatures, non-thermal motions, plume cross-sections relative to the size of a CH, the plasma bulk speeds, as well as any plume signatures in the SW. In spring 2007, the authors proposed a study on "Structure and dynamics of coronal plumes and inter-plume regions in solar coronal holes" to the International Space Science Institute (ISSI) in Bern to clarify some of these aspects by considering relevant observations and the extensive literature. This review summarizes the results and conclusions of the study. Stereoscopic observations allowed us to include three-dimensional reconstructions of plumes. Multi-instrument investigations carried out during several campaigns led to progress in some areas, such as plasma densities, temperatures, plume structure and the relation to other solar phenomena, but not all questions could be answered concerning the details of plume generation process(es) and interaction with the SW.

preprint2011arXiv

On the emergence of very long time fluctuations and 1/f noise in ideal flows

This study 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 wavenumber 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 non-local 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.

preprint2010arXiv

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, even though 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.

preprint2010arXiv

Ultraviolet Coronagraph Spectroscopy: A Key Capability for Understanding the Physics of Solar Wind Acceleration

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.

preprint2009arXiv

A turbulence-driven model for heating and acceleration of the fast wind in coronal holes

A model is presented for generation of fast solar wind in coronal holes, relying on heating that is dominated by turbulent dissipation of MHD fluctuations transported upwards in the solar atmosphere. Scale-separated transport equations include large-scale fields, transverse Alfvenic fluctuations, and a small compressive dissipation due to parallel shears near the transition region. The model accounts for proton temperature, density, wind speed, and fluctuation amplitude as observed in remote sensing and in situ satellite data.

preprint2009arXiv

Electron and proton heating by solar wind turbulence

Previous formulations of heating and transport associated with strong magnetohydrodynamic (MHD) turbulence are generalized to incorporate separate internal energy equations for electrons and protons. Electron heat conduction is included. Energy is supplied by turbulent heating that affects both electrons and protons, and is exchanged between them via collisions. Comparison to available Ulysses data shows that a reasonable accounting for the data is provided when (i) the energy exchange timescale is very long and (ii) the deposition of heat due to turbulence is divided, with 60% going to proton heating and 40% into electron heating. Heat conduction, determined here by an empirical fit, plays a major role in describing the electron data.

preprint2009arXiv

The third-order law for increments in magnetohydrodynamic turbulence with constant shear

We extend the theory for third-order structure functions in homogeneous incompressible magnetohydrodynamic (MHD) turbulence to the case in which a constant velocity shear is present. A generalization is found of the usual relation [Politano and Pouquet, Phys. Rev. E, 57 21 (1998)] between third-order structure functions and the dissipation rate in steady inertial range turbulence, in which the shear plays a crucial role. In particular, the presence of shear leads to a third-order law which is not simply proportional to the relative separation. Possible implications for laboratory and space plasmas are discussed.

preprint2008arXiv

Kinetic dissipation and anisotropic heating in a turbulent collisionless plasma

The kinetic evolution of the Orszag-Tang vortex is studied using collisionless hybrid simulations. In the magnetohydrodynamic regime this vortex leads rapidly to broadband turbulence. Significant differences from MHD arise at small scales, where the fluid scale energy dissipates into heat almost exclusively through the magnetic field because the protons are decoupled from the magnetic field. Although cyclotron resonance is absent, the protons heat preferentially in the plane perpendicular to the mean field, as in the corona and solar wind. Effective transport coefficients are calculated.