Source author record

Daniel Verscharen

Daniel Verscharen 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

37works
7topics
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

37 published item(s)

preprint2022arXiv

Coherence of ion cyclotron resonance for damping ion cyclotron waves in space plasmas

Ion cyclotron resonance is one of the fundamental energy conversion processes through field-particle interaction in collisionless plasmas. However, the key evidence for ion cyclotron resonance (i.e., the coherence between electromagnetic fields and the ion phase space density) and the resulting damping of ion cyclotron waves (ICWs) has not yet been directly observed. Investigating the high-quality measurements of space plasmas by the Magnetospheric Multiscale (MMS) satellites, we find that both the wave electromagnetic field vectors and the bulk velocity of the disturbed ion velocity distribution rotate around the background magnetic field. Moreover, we find that the absolute gyro-phase angle difference between the center of the fluctuations in the ion velocity distribution functions and the wave electric field vectors falls in the range of (0, 90) degrees, consistent with the ongoing energy conversion from wave-fields to particles. By invoking plasma kinetic theory, we demonstrate that the field-particle correlation for the damping ion cyclotron waves in our theoretical model matches well with our observations. Furthermore, the wave electric field vectors ($δ\mathbf{E'}_{\mathrm {wave,\perp}}$), the ion current density ($δ\mathbf{J}_\mathrm {i,\perp}$) and the energy transfer rate ($δ\mathbf{J}_\mathrm {i,\perp}\cdot δ\mathbf{E'}_{\mathrm {wave,\perp}}$) exhibit quasi-periodic oscillations, and the integrated work done by the electromagnetic field on the ions are positive, indicates that ions are mainly energized by the perpendicular component of the electric field via cyclotron resonance. Therefore, our combined analysis of MMS observations and kinetic theory provides direct, thorough, and comprehensive evidence for ICW damping in space plasmas.

preprint2022arXiv

Direct evidence of magnetic reconnection onset via the tearing instability

Magnetic reconnection is a sporadic process responsible for energy release in space and laboratory plasmas. It is believed that the tearing mode instability may be responsible for the onset of reconnection in the magnetotail. However, due to its elusive nature, there is an absence of in-situ observations of the tearing instability prior to magnetic reconnection in our nearest natural plasma laboratory. Using neural network outlier detection methods in conjunction with Cluster spacecraft data, we find unique electron pitch angle distributions that are consistent with simulation predictions of the tearing instability and the subsequent evolution of plasma electrons and reconnection. We confirm that the events identified via our neural network outlier method are well above the tearing stability threshold based on the criterion detailed in this paper. We find signatures of magnetic reconnection minutes after the majority of tearing observations. Our analysis of the tearing instability provides new insights into the fundamental understanding of the mechanism responsible for reconnection, a process that is ubiquitous in different astrophysical plasma regimes across the universe and in laboratory experiments on Earth.

preprint2022arXiv

Electron-driven instabilities in the solar wind

The electrons are an essential particle species in the solar wind. They often exhibit non-equilibrium features in their velocity distribution function. These include temperature anisotropies, tails (kurtosis), and reflectional asymmetries (skewness), which contribute a significant heat flux to the solar wind. If these non-equilibrium features are sufficiently strong, they drive kinetic micro-instabilities. We develop a semi-graphical framework based on the equations of quasi-linear theory to describe electron-driven instabilities in the solar wind. We apply our framework to resonant instabilities driven by temperature anisotropies. These include the electron whistler anisotropy instability and the propagating electron firehose instability. We then describe resonant instabilities driven by reflectional asymmetries in the electron distribution function. These include the electron/ion-acoustic, kinetic Alfvén heat-flux, Langmuir, electron-beam, electron/ion-cyclotron, electron/electron-acoustic, whistler heat-flux, oblique fast-magnetosonic/whistler, lower-hybrid fan, and electron-deficit whistler instability. We briefly comment on non-resonant instabilities driven by electron temperature anisotropies such as the mirror-mode and the non-propagating firehose instability. We conclude our review with a list of open research topics in the field of electron-driven instabilities in the solar wind.

preprint2022arXiv

Radial evolution of thermal and suprathermal electron populations in the slow solar wind from 0.13 to 0.5 au : Parker Solar Probe Observations

We develop and apply a bespoke fitting routine to a large volume of solar wind electron distribution data measured by Parker Solar Probe (PSP) over its first five orbits, covering radial distances from 0.13 to 0.5 au. We characterise the radial evolution of the electron core, halo and strahl populations in the slow solar wind during these orbits. The fractional densities of these three electron populations provide evidence for the growth of the combined suprathermal halo and strahl populations from 0.13 to 0.17 au. Moreover, the growth in the halo population is not matched by a decrease of the strahl population at these distances, as has been reported for previous observations at distances greater than 0.3 au. We also find that the halo is negligible at small heliocentric distances. The fractional strahl density remains relatively constant ~1 % below 0.2 au, suggesting that the rise in the relative halo density is not solely due to the transfer of strahl electrons into the halo.

preprint2022arXiv

Test Problems for Potential Field Source Surface Extrapolations of Solar and Stellar Magnetic Fields

The potential field source surface (PFSS) equations are commonly used to model the coronal magnetic field of the Sun and other stars. As with any computational model, solving equations using a numerical scheme introduces errors due to discretisation. We present a set of tests for quantifying these errors by taking advantage of analytic solutions to the PFSS equations when the input field is proportional to a single spherical harmonic. From the spherical harmonic solutions we derive analytic equations for magnetic field lines traced through the three dimensional magnetic field solution. We propose these as a set of standard analytic solutions that all PFSS solvers should be tested against to quantify their inherent errors. We apply these tests to the pfsspy software package, showing that it reproduces spherical harmonic solutions well with a slight overestimation of the unsigned open magnetic flux. It is also successful at reproducing analytic field line equations, with errors in field line footpoints typically much less than one degree.

preprint2022arXiv

The Kinetic Expansion of Solar-Wind Electrons: Transport Theory and Predictions for the very Inner Heliosphere

We propose a transport theory for the kinetic evolution of solar-wind electrons in the heliosphere. We derive a gyro-averaged kinetic transport equation that accounts for the spherical expansion of the solar wind and the geometry of the Parker-spiral magnetic field. To solve our three-dimensional kinetic equation, we develop a mathematical approach that combines the Crank--Nicolson scheme in velocity space and a finite-difference Euler scheme in configuration space. We initialize our model with isotropic electron distribution functions and calculate the kinetic expansion at heliocentric distances from 5 to 20 solar radii. In our kinetic model, the electrons evolve mainly through the combination of the ballistic particle streaming, the magnetic mirror force, and the electric field. By applying fits to our numerical results, we quantify the parameters of the electron strahl and core part of the electron velocity distributions. The strahl fit parameters show that the density of the electron strahl is around 7% of the total electron density at a distance of 20 solar radii, the strahl bulk velocity and strahl temperature parallel to the background magnetic field stay approximately constant beyond a distance of 15 solar radii, and $β_{\parallel s}$ (i.e., the ratio between strahl parallel thermal pressure to the magnetic pressure) is approximately constant with heliocentric distance at a value of about 0.02. We compare our results with data measured by Parker Solar Probe. Furthermore, we provide theoretical evidence that the electron strahl is not scattered by the oblique fast-magnetosonic/whistler instability in the near-Sun environment.

preprint2022arXiv

The Stability of the Electron Strahl against the Oblique Fast-magnetosonic/Whistler Instability in the Inner Heliosphere

We analyze the micro-kinetic stability of the electron strahl in the solar wind depending on heliocentric distance. The oblique fast-magnetosonic/whistler (FM/W) instability has emerged in the literature as a key candidate mechanism for the effective scattering of the electron strahl into the electron halo population. Using data from Parker Solar Probe (PSP) and Helios, we compare the measured strahl properties with the analytical thresholds for the oblique FM/W instability in the low- and high-$β_{\parallel c}$ regimes, where $β_{\parallel c}$ is the ratio of the core parallel thermal pressure to the magnetic pressure. Our PSP and Helios data show that the electron strahl is on average stable against the oblique FM/W instability in the inner heliosphere. Our analysis suggests that the instability, if at all, can only be excited sporadically and on short timescales. We discuss the caveats of our analysis and potential alternative explanations for the observed scattering of the electron strahl in the solar wind. Furthermore, we recommend the numerical evaluation of the stability of individual distributions in the future to account for any uncertainties in the validity of the analytical expressions for the instability thresholds.

preprint2022arXiv

Whistler Waves As a Signature of Converging Magnetic Holes in Space Plasmas

Magnetic holes are plasma structures that trap a large number of particles in a magnetic field that is weaker than the field in its surroundings. The unprecedented high time-resolution observations by NASA's Magnetospheric Multi-Scale (MMS) mission enable us to study the particle dynamics in magnetic holes in the Earth's magnetosheath in great detail. We reveal the local generation mechanism of whistler waves by a combination of Landau-resonant and cyclotron-resonant wave-particle interactions of electrons in response to the large-scale evolution of a magnetic hole. As the magnetic hole converges, a pair of counter-streaming electron beams form near the hole's center as a consequence of the combined action of betatron and Fermi effects. The beams trigger the generation of slightly-oblique whistler waves. Our conceptual prediction is supported by a remarkable agreement between our observations and numerical predictions from the Arbitrary Linear Plasma Solver (ALPS). Our study shows that wave-particle interactions are fundamental to the evolution of magnetic holes in space and astrophysical plasmas.

preprint2021arXiv

Growth of Outward Propagating Fast-Magnetosonic/Whistler Waves in the Inner Heliosphere Observed by Parker Solar Probe

The solar wind in the inner heliosphere has been observed by Parker Solar Probe (PSP) to exhibit abundant wave activities. The cyclotron wave modes in the sense of ions or electrons are among the most crucial wave components. However, their origin and evolution in the inner heliosphere close to the Sun remain mysteries. Specifically, it remains unknown whether it is an emitted signal from the solar atmosphere or an eigenmode growing locally in the heliosphere due to plasma instability. To address and resolve this controversy, we must investigate the key quantity of the energy change rate of the wave mode. We develop a new technique to measure the energy change rate of plasma waves, and apply this technique to the wave electromagnetic fields measured by PSP. We provide the wave Poynting flux in the solar wind frame, identify the wave nature to be the outward propagating fast-magnetosonic/whistler wave mode instead of the sunward propagating waves. We provide the first evidence for growth of the fast-magnetosonic/whistler wave mode in the inner heliosphere based on the derived spectra of the real and imaginary parts of the wave frequencies. The energy change rate rises and stays at a positive level in the same wavenumber range as the bumps of the electromagnetic field power spectral densities, clearly manifesting that the observed fast-magnetosonic/whistler waves are locally growing to a large amplitude.

preprint2021arXiv

Stochastic Electron Acceleration by Temperature Anisotropy Instabilities Under Solar Flare Plasma Conditions

Using 2D particle-in-cell (PIC) plasma simulations we study electron acceleration by temperature anisotropy instabilities, assuming conditions typical of above-the-loop-top (ALT) sources in solar flares. We focus on the long-term effect of $T_{e,\perp} > T_{e,\parallel}$ instabilities by driving the anisotropy growth during the entire simulation time, through imposing a shearing or a compressing plasma velocity ($T_{e,\perp}$ and $T_{e,\parallel}$ are the temperatures perpendicular and parallel to the magnetic field). This magnetic growth makes $T_{e,\perp}/T_{e,\parallel}$ grow due to electron magnetic moment conservation, and amplifies the ratio $ω_{ce}/ω_{pe}$ from $\sim 0.53$ to $\sim 2$ ($ω_{ce}$ and $ω_{pe}$ are the electron cyclotron and plasma frequencies, respectively). In the regime $ω_{ce}/ω_{pe}\lesssim 1.2-1.7$ the instability is dominated by oblique, quasi-electrostatic (OQES) modes, and the acceleration is inefficient. When $ω_{ce}/ω_{pe}$ has grown to $ω_{ce}/ω_{pe}\gtrsim 1.2-1.7$, electrons are efficiently accelerated by the inelastic scattering provided by unstable parallel, electromagnetic z (PEMZ) modes. After $ω_{ce}/ω_{pe}$ reaches $\sim 2$, the electron energy spectra show nonthermal tails that differ between the shearing and compressing cases. In the shearing case, the tail resembles a power-law of index $α_s \sim$ 2.9 plus a high-energy bump reaching $\sim 300$ keV. In the compressing runs, $α_s \sim$ 3.7 with a spectral break above $\sim 500$ keV. This difference can be explained by the different temperature evolutions in these two types of simulations, suggesting a critical role played by the type of anisotropy driving, $ω_{ce}/ω_{pe}$ and the electron temperature in the efficiency of the acceleration.

preprint2020arXiv

Dependence of kinetic plasma waves on ion-to-electron mass ratio and light-to-Alfvén speed ratio

The magnetization $|Ω_{\mathrm e}|/ω_{\mathrm{e}}$ is an important parameter in plasma astrophysics, where $Ω_{\mathrm e}$ and $ω_{\mathrm{e}}$ are the electron gyro-frequency and electron plasma frequency, respectively. It only depends on the mass ratio $m_{\mathrm i}/m_{\mathrm e}$ and the light-to-Alfvén speed ratio $c/v_{\mathrm{Ai}}$, where $m_{\mathrm i}$ ($m_{\mathrm e}$) is the ion (electron) mass, $c$ is the speed of light, and $v_{\mathrm{Ai}}$ is the ion Alfvén speed. Nonlinear numerical plasma models such as particle-in-cell simulations must often assume unrealistic values for $m_{\mathrm i}/m_{\mathrm e}$ and for $c/v_{\mathrm{Ai}}$. Because linear theory yields exact results for parametric scalings of wave properties at small amplitudes, we use linear theory to investigate the dispersion relations of Alfvén/ion-cyclotron and fast-magnetosonic/whistler waves as prime examples for collective plasma behaviour depending on $m_{\mathrm i}/m_{\mathrm e}$ and $c/v_{\mathrm{Ai}}$. We analyse their dependence on $m_{\mathrm i}/m_{\mathrm e}$ and $c/v_{\mathrm{Ai}}$ in quasi-parallel and quasi-perpendicular directions of propagation with respect to the background magnetic field for a plasma with $β_j\sim1$, where $β_j$ is the ratio of the thermal to magnetic pressure for species $j$. Although their dispersion relations are largely independent of $c/v_{\mathrm{Ai}}$ for $c/v_{\mathrm{Ai}}\gtrsim 10$, the mass ratio $m_{\mathrm i}/m_{\mathrm e}$ has a strong effect at scales smaller than the ion inertial length. Moreover, we study the impact of relativistic electron effects on the dispersion relations. Based on our results, we recommend aiming for a more realistic value of $m_{\mathrm i}/m_{\mathrm e}$ than for a more realistic value of $c/v_{\mathrm{Ai}}$ in non-relativistic plasma simulations if such a choice is necessary, although $\dots$

preprint2020arXiv

Evolution of the Earth's Magnetosheath Turbulence: A statistical study based on MMS observations

Composed of shocked solar wind, the Earth's magnetosheath serves as a natural laboratory to study the transition of turbulence from low Alfv{é}n Mach number, $M_\mathrm{A}$, to high $M_\mathrm{A}$. The simultaneous observations of magnetic field and plasma moments with unprecedented high temporal resolution provided by NASA's \textit{Magnetospheric Multiscale} Mission enable us to study the magnetosheath turbulence at both magnetohydrodynamics (MHD) and sub-ion scales. Based on 1841 burst-mode segments of MMS-1 from 2015/09 to 2019/06, comprehensive patterns of the spatial evolution of magnetosheath turbulences are obtained: (1) from the sub-solar region to the flanks, $M_\mathrm{A}$ increases from $<$ 1 to $>$ 5. At MHD scales, the spectral indices of the magnetic-field and velocity spectra present a positive and negative correlation with $M_\mathrm{A}$. However, no obvious correlations between the spectral indices and $M_\mathrm{A}$ are found at sub-ion scales. (2) from the bow shock to the magnetopause, the turbulent sonic Mach number, $M_{\mathrm{turb}}$, generally decreases from $>$ 0.4 to $<$ 0.1. All spectra steepen at MHD scales and flatten at sub-ion scales, representing a positive/negative correlations with $M_\mathrm{turb}$. The break frequency increases by 0.1 Hz when approaching the magnetopause for the magnetic-field and velocity spectra, while it remains at 0.3 Hz for the density spectra. (3) In spite of some differences, similar results are found for the quasi-parallel and quasi-perpendicular magnetosheath. In addition, the spatial evolution of magnetosheath turbulence is found to be independent of the upstream solar wind conditions, e.g., the Z-component of the interplanetary magnetic field and the solar wind speed.

preprint2020arXiv

Polytropic Behavior of Solar Wind Protons Observed by Parker Solar Probe

A polytropic process describes the transition of a fluid from one state to another through a specific relationship between the fluid density and temperature. The value of the polytropic index that governs this relationship determines the heat transfer and the effective degrees of freedom during the process. In this study, we analyze solar wind proton plasma measurements, obtained by the Faraday cup instrument on-board Parker Solar Probe. We examine the large-scale variations of the proton plasma density and temperature within the inner heliosphere explored by the spacecraft. We also address a polytropic behavior in the density and temperature fluctuations in short-time intervals, which we analyze in order to derive the effective polytropic index of small time-scale processes. The large-scale variations of the solar wind proton density and temperature which are associated with the plasma expansion through the heliosphere, follow a polytropic model with a polytropic index ~5/3. On the other hand, the short time-scale fluctuations which may be associated with turbulence, follow a model with a larger polytropic index. We investigate possible correlations between the polytropic index of short time-scale fluctuations and the plasma speed, plasma beta, and the magnetic field direction. We discuss the scenario of mechanisms including energy transfer or mechanisms that restrict the particle effective degrees of freedom.

preprint2020arXiv

Statistics of Solar Wind Electron Breakpoint Energies Using Machine Learning Techniques

Solar wind electron velocity distributions at 1 au consist of a thermal "core" population and two suprathermal populations: "halo" and "strahl". The core and halo are quasi-isotropic, whereas the strahl typically travels radially outwards along the parallel and/or anti-parallel direction with respect to the interplanetary magnetic field. With Cluster-PEACE data, we analyse energy and pitch angle distributions and use machine learning techniques to provide robust classifications of these solar wind populations. Initially, we use unsupervised algorithms to classify halo and strahl differential energy flux distributions to allow us to calculate relative number densities, which are of the same order as previous results. Subsequently, we apply unsupervised algorithms to phase space density distributions over ten years to study the variation of halo and strahl breakpoint energies with solar wind parameters. In our statistical study, we find both halo and strahl suprathermal breakpoint energies display a significant increase with core temperature, with the halo exhibiting a more positive correlation than the strahl. We conclude low energy strahl electrons are scattering into the core at perpendicular pitch angles. This increases the number of Coulomb collisions and extends the perpendicular core population to higher energies, resulting in a larger difference between halo and strahl breakpoint energies at higher core temperatures. Statistically, the locations of both suprathermal breakpoint energies decrease with increasing solar wind speed. In the case of halo breakpoint energy, we observe two distinct profiles above and below 500 km/s. We relate this to the difference in origin of fast and slow solar wind.

preprint2020arXiv

Wave Composition, Propagation, and Polarization of MHD Turbulence within 0.3AU as Observed by PSP

Turbulence, a ubiquitous phenomenon in interplanetary space, is crucial for the energy conversion of space plasma at multiple scales. This work focuses on the propagation, polarization and wave composition properties of the solar wind turbulence within 0.3AU, and its variation with heliocentric distances at MHD scales (from 10s to 1000s in the spacecraft frame). We present the probability density function of propagation wavevectors (${\rm{PDF}}(k_\parallel,k_\perp)$) for solar wind turbulence winthin 0.3 AU for the first time: (1) wavevectors cluster quasi-(anti-)parallel to the local background magnetic field for $kd_{\rm i}<0.02$, where $d_{\rm i}$ is the ion inertial length; (2) wavevectors shift to quasi-perpendicular directions for $kd_{\rm i}>0.02$. Based on our wave composition diagnosis, we find that: the outward/anti-sunward Alfvén mode dominates over the whole range of scales and distances, the spectral energy density fraction of the inward/sunward fast mode decreases with distance, and the fractional energy densities of the inward and outward slow mode increase with distance. The outward fast mode and inward Alfvén mode represent minority populations throughout the explored range of distances and scales. On average, the degree of anisotropy of the magnetic fluctuations defined with respect to the minimum variation direction decreases with increasing scale, with no trend in distance at all scales. Our results provide comprehensive insight into the scenario of transport and transfer of the solar wind fluctuations/turbulence in the inner heliosphere.

preprint2019arXiv

The Impact of Turbulent Solar Wind Fluctuations on Solar Orbiter Plasma Proton Measurements

Solar Orbiter will observe the Sun and the inner heliosphere to study the connections between solar activity, coronal structure, and the origin of the solar wind. The plasma instruments on board Solar Orbiter will determine the three-dimensional velocity distribution functions of the plasma ions and electrons with high time resolution. The analysis of these distributions will determine the plasma bulk parameters, such as density, velocity, and temperature. This paper examines the effects of short-time-scale plasma variations on particle measurements and the estimated bulk parameters of the plasma. For the purpose of this study, we simulate the expected observations of solar wind protons, taking into account the performance of the Proton-Alpha Sensor (PAS) on board Solar Orbiter. We particularly examine the effects of Alfvénic and slow-mode-like fluctuations, commonly observed in the solar wind on timescales of milliseconds to hours, on the observations. We do this by constructing distribution functions from modeled observations and calculate their statistical moments in order to derive plasma bulk parameters. The comparison between the derived parameters with the known input allows us to estimate the expected accuracy of Solar Orbiter proton measurements in the solar wind under typical conditions. We find that the plasma fluctuations due to these turbulence effects have only minor effects on future SWA-PAS observations

preprint2019arXiv

The multi-scale nature of the solar wind

The solar wind is a magnetized plasma and as such exhibits collective plasma behavior associated with its characteristic spatial and temporal scales. The characteristic length scales include the size of the heliosphere, the collisional mean free paths of all species, their inertial lengths, their gyration radii, and their Debye lengths. The characteristic timescales include the expansion time, the collision times, and the periods associated with gyration, waves, and oscillations. We review the past and present research into the multi-scale nature of the solar wind based on in-situ spacecraft measurements and plasma theory. We emphasize that couplings of processes across scales are important for the global dynamics and thermodynamics of the solar wind. We describe methods to measure in-situ properties of particles and fields. We then discuss the role of expansion effects, non-equilibrium distribution functions, collisions, waves, turbulence, and kinetic microinstabilities for the multi-scale plasma evolution.

preprint2016arXiv

Collisionless Isotropization of the Solar-Wind Protons by Compressive Fluctuations and Plasma Instabilities

Compressive fluctuations are a minor yet significant component of astrophysical plasma turbulence. In the solar wind, long-wavelength compressive slow-mode fluctuations lead to changes in $β_{\parallel \mathrm p}\equiv 8πn_{\mathrm p}k_{\mathrm B}T_{\parallel \mathrm p}/B^2$ and in $R_{\mathrm p}\equiv T_{\perp \mathrm p}/T_{\parallel \mathrm p}$, where $T_{\perp \mathrm p}$ and $T_{\parallel \mathrm p}$ are the perpendicular and parallel temperatures of the protons, $B$ is the magnetic field strength, and $n_{\mathrm p}$ is the proton density. If the amplitude of the compressive fluctuations is large enough, $R_{\mathrm p}$ crosses one or more instability thresholds for anisotropy-driven microinstabilities. The enhanced field fluctuations from these microinstabilities scatter the protons so as to reduce the anisotropy of the pressure tensor. We propose that this scattering drives the average value of $R_{\mathrm p}$ away from the marginal stability boundary until the fluctuating value of $R_{\mathrm p}$ stops crossing the boundary. We model this "fluctuating-anisotropy effect" using linear Vlasov--Maxwell theory to describe the large-scale compressive fluctuations. We argue that this effect can explain why, in the nearly collisionless solar wind, the average value of $R_{\mathrm p}$ is close to unity.

preprint2016arXiv

Electrons under the dominant action of shock-electric fields

We consider a fast magnetosonic multifluid shock as a representation of the solar-wind termination shock. We assume the action of the transition happens in a three-step process: In the first step, the upstream supersonic solar-wind plasma is subject to a strong electric field that flashes up on a small distance scale $Δz\simeq U_1/ Ω_{\mathrm e}$ (first part of the transition layer), where $Ω_{\mathrm e}$ is the electron gyro-frequency and $U_1$ is the upstream speed. This electric field both decelerates the supersonic ion flow and accelerates the electrons up to high velocities. In this part of the transition region, the electric forces connected with the deceleration of the ion flow strongly dominate over the Lorentz forces. We, therefore, call this part the demagnetization region. In the second phase, Lorentz forces due to convected magnetic fields compete with the electric field, and the highly anisotropic and energetic electron distribution function is converted into a shell distribution with energetic shell electrons storing about 3/4 of the upstream ion kinetic energy. In the third phase, the plasma particles thermalize due to the relaxation of free energy by plasma instabilities. The first part of the transition region opens up a new thermodynamic degree of freedom never before taken into account for the electrons, since the electrons are usually considered to be enslaved to follow the behavior of the protons in all velocity moments like density, bulk velocity, and temperature. We show that electrons may be the downstream plasma fluid that dominates the downstream plasma pressure.

preprint2016arXiv

PIC Simulations of the Effect of Velocity Space Instabilities on Electron Viscosity and Thermal Conduction

In low-collisionality plasmas, velocity-space instabilities are a key mechanism providing an effective collisionality for the plasma. We use particle-in-cell (PIC) simulations to study the interplay between electron and ion-scale velocity-space instabilities and their effect on electron pressure anisotropy, viscous heating, and thermal conduction. The adiabatic invariance of the magnetic moment in low-collisionality plasmas leads to pressure anisotropy, $p_{\perp,j} > p_{||,j}$, if the magnetic field $\vec{B}$ is amplified ($p_{\perp,j}$ and $p_{||,j}$ denote the pressure of species $j$ [electron, ion] perpendicular and parallel to $\vec{B}$). If the resulting anisotropy is large enough, it can in turn trigger small-scale plasma instabilities. Our PIC simulations explore the nonlinear regime of the mirror, ion-cyclotron, and electron whistler instabilities, through continuous amplification of the magnetic field $|\vec{B}|$ by an imposed shear in the plasma. In the regime $1 \lesssim β_j \lesssim 20$ ($β_j \equiv 8πp_j/|\vec{B}|^2$), the saturated electron pressure anisotropy, $Δp_e/p_{||,e}$, is determined mainly by the (electron-lengthscale) whistler marginal stability condition, with a modest factor of $\sim 1.5-2$ decrease due to the trapping of electrons by the mirrors. We explicitly calculate the mean free path of the electrons and ions along the mean magnetic field and provide a simple physical prescription for the mean free path and thermal conductivity in low-collisionality $β_j \gtrsim 1$ plasmas. Our results imply that velocity-space instabilities likely decrease the thermal conductivity of plasma in the outer parts of massive, hot, galaxy clusters. We also discuss the implications of our results for electron heating and thermal conduction in low-collisionality accretion flows onto black holes, including Sgr A* in the Galactic Center.

preprint2015arXiv

A Modified Version of Taylor's Hypothesis for Solar Probe Plus Observations

The Solar Probe Plus (SPP) spacecraft will explore the near-Sun environment, reaching heliocentric distances less than $10 R_{\odot}$. Near Earth, spacecraft measurements of fluctuating velocities and magnetic fields taken in the time domain are translated into information about the spatial structure of the solar wind via Taylor's "frozen turbulence" hypothesis. Near the perihelion of SPP, however, the solar-wind speed is comparable to the Alfvén speed, and Taylor's hypothesis in its usual form does not apply. In this paper, we show that, under certain assumptions, a modified version of Taylor's hypothesis can be recovered in the near-Sun region. We consider only the transverse, non-compressive component of the fluctuations at length scales exceeding the proton gyroradius, and we describe these fluctuations using an approximate theoretical framework developed by Heinemann and Olbert. We show that fluctuations propagating away from the Sun in the plasma frame obey a relation analogous to Taylor's hypothesis when $V_{\rm sc,\perp} \gg z^-$ and $z^+ \gg z^-$, where $V_{\rm sc,\perp}$ is the component of the spacecraft velocity perpendicular to the mean magnetic field and $\bm{z}^+$ ($\bm{z}^-$) is the Elsasser variable corresponding to transverse, non-compressive fluctuations propagating away from (towards) the Sun in the plasma frame. Observations and simulations suggest that, in the near-Sun solar wind, the above inequalities are satisfied and $\bm{z}^+$ fluctuations account for most of the fluctuation energy. The modified form of Taylor's hypothesis that we derive may thus make it possible to characterize the spatial structure of the energetically dominant component of the turbulence encountered by SPP.

preprint2015arXiv

Deceleration of Alpha Particles in the Solar Wind by Instabilities and the Rotational Force: Implications for Heating, Azimuthal Flow, and the Parker Spiral Magnetic Field

Protons and alpha particles in the fast solar wind are only weakly collisional and exhibit a number of non-equilibrium features, including relative drifts between particle species. Two non-collisional mechanisms have been proposed for limiting differential flow between alpha particles and protons: plasma instabilities and the rotational force. Both mechanisms decelerate the alpha particles. In this paper, we derive an analytic expression for the rate $Q_{\mathrm{flow}}$ at which energy is released by alpha-particle deceleration, accounting for azimuthal flow and conservation of total momentum. We show that instabilities control the deceleration of alpha particles at $r< r_{\mathrm{crit}}$, and the rotational force controls the deceleration of alpha particles at $r> r_{\mathrm{crit}}$, where $r_{\mathrm{crit}} \simeq 2.5 \,\mathrm{AU}$ in the fast solar wind in the ecliptic plane. We find that $Q_{\mathrm{flow}}$ is positive at $r<r_{\mathrm{crit}}$ and $Q_{\mathrm{flow}} = 0$ at $r\geq r_{\mathrm{crit}}$, consistent with the previous finding that the rotational force does not lead to a release of energy. We compare the value of~$Q_{\mathrm{flow}}$ at $r< r_{\mathrm{crit}}$ with empirical heating rates for protons and alpha particles, denoted $Q_{\mathrm{p}}$ and $Q_α$, deduced from in-situ measurements of fast-wind streams from the \emph{Helios} and \emph{Ulysses} spacecraft. We find that $Q_{\mathrm{flow}}$ exceeds $Q_α$ at $r < 1\,\mathrm{AU}$, and that $Q_{\mathrm{flow}}/Q_{\rm p}$ decreases with increasing distance from the Sun from a value of about one at $r=0.29 - 0.42\,\mathrm{AU}$ to about 1/4 at 1 AU. We conclude that the continuous energy input from alpha-particle deceleration at $r< r_{\mathrm{crit}}$ makes an important contribution to the heating of the fast solar wind.

preprint2015arXiv

Effects of electron drift on the collisionless damping of kinetic Alfvén waves in the solar wind

The collisionless dissipation of anisotropic Alfvénic turbulence is a promising candidate to solve the solar wind heating problem. Extensive studies examined the kinetic properties of Alfvén waves in simple Maxwellian or bi-Maxwellian plasmas. However, the observed electron velocity distribution functions in the solar wind are more complex. In this study, we analyze the properties of kinetic Alfvén waves in a plasma with two drifting electron populations. We numerically solve the linearized Maxwell-Vlasov equations and find that the damping rate and the proton-electron energy partition for kinetic Alfvén waves are significantly modified in such plasmas, compared to plasmas without electron drifts. We suggest that electron drift is an important factor to take into account when considering the dissipation of Alfvénic turbulence in the solar wind or other $β\sim 1$ astrophysical plasmas.

preprint2015arXiv

On the Conservation of Cross Helicity and Wave Action in Solar-Wind Models with Non-WKB Alfven Wave Reflection

The interaction between Alfven-wave turbulence and the background solar wind affects the cross helicity in two ways. Non-WKB reflection converts outward-propagating Alfven waves into inward-propagating Alfven waves and vice versa, and the turbulence transfers momentum to the background flow. When both effects are accounted for, the total cross helicity is conserved. In the special case that the background density and flow speed are independent of time, the equations of cross-helicity conservation and total-energy conservation can be combined to recover a well-known equation derived by Heinemann and Olbert that has been interpreted as a non-WKB generalization of wave-action conservation. This latter equation (in contrast to cross-helicity and energy conservation) does not hold when the background varies in time.

preprint2015arXiv

PIC Simulations of Continuously Driven Mirror and Ion Cyclotron Instabilities in High Beta Astrophysical and Heliospheric Plasmas

We use particle-in-cell (PIC) simulations to study the nonlinear evolution of ion velocity space instabilities in an idealized problem in which a background velocity shear continuously amplifies the magnetic field. We simulate the astrophysically relevant regime where the shear timescale is long compared to the ion cyclotron period, and the plasma beta is ~ 1-100. The background field amplification in our calculation is meant to mimic processes such as turbulent fluctuations or MHD-scale instabilities. The field amplification continuously drives a pressure anisotropy with the perpendicular pressure larger than the parallel pressure, and the plasma becomes unstable to the mirror and ion cyclotron instabilities. In all cases, the nonlinear state is dominated by the mirror instability, not the ion cyclotron instability, and the plasma pressure anisotropy saturates near the threshold for the linear mirror instability. The magnetic field fluctuations initially undergo exponential growth but saturate in a secular phase in which the fluctuations grow on the same timescale as the background magnetic field (with δB ~ 0.3 <B> in the secular phase). At early times, the ion magnetic moment is well-conserved but once the fluctuation amplitudes exceed δB ~ 0.1 <B>, the magnetic moment is no longer conserved but instead changes on a timescale comparable to that of the mean magnetic field. We discuss the implications of our results for low-collisionality astrophysical plasmas, including the near-Earth solar wind and low-luminosity accretion disks around black holes.

preprint2015arXiv

The electron distribution function downstream of the solar-wind termination shock: Where are the hot electrons?

In the majority of the literature on plasma shock waves, electrons play the role of "ghost particles," since their contribution to mass and momentum flows is negligible, and they have been treated as only taking care of the electric plasma neutrality. In some more recent papers, however, electrons play a new important role in the shock dynamics and thermodynamics, especially at the solar-wind termination shock. They react on the shock electric field in a very specific way, leading to suprathermal nonequilibrium distributions of the downstream electrons, which can be represented by a kappa distribution function. In this paper, we discuss why this anticipated hot electron population has not been seen by the plasma detectors of the Voyager spacecraft downstream of the solar-wind termination shock. We show that hot nonequilibrium electrons induce a strong negative electric charge-up of any spacecraft cruising through this downstream plasma environment. This charge reduces electron fluxes at the spacecraft detectors to nondetectable intensities. Furthermore, we show that the Debye length $λ_{\mathrm D}^κ$ grows to values of about $λ_{\mathrm D}^κ/λ_{\mathrm D}\simeq 10^{6}$ compared to the classical value $λ_{\mathrm D}$ in this hot-electron environment. This unusual condition allows for the propagation of a certain type of electrostatic plasma waves that, at very large wavelengths, allow us to determine the effective temperature of the suprathermal electrons directly by means of the phase velocity of these waves. At moderate wavelengths, the electron-acoustic dispersion relation leads to nonpropagating oscillations with the ion-plasma frequency $ω_{\mathrm p}$ , instead of the traditional electron plasma frequency.

preprint2014arXiv

Magnetohydrodynamic Slow Mode with Drifting He$^{++}$: Implications for Coronal Seismology and the Solar Wind

The MHD slow mode wave has application to coronal seismology, MHD turbulence, and the solar wind where it can be produced by parametric instabilities. We consider analytically how a drifting ion species (e.g. He$^{++}$) affects the linear slow mode wave in a mainly electron-proton plasma, with potential consequences for the aforementioned applications. Our main conclusions are: 1. For wavevectors highly oblique to the magnetic field, we find solutions that are characterized by very small perturbations of total pressure. Thus, our results may help to distinguish the MHD slow mode from kinetic Alfvén waves and non-propagating pressure-balanced structures, which can also have very small total pressure perturbations. 2. For small ion concentrations, there are solutions that are similar to the usual slow mode in an electron-proton plasma, and solutions that are dominated by the drifting ions, but for small drifts the wave modes cannot be simply characterized. 3. Even with zero ion drift, the standard dispersion relation for the highly oblique slow mode cannot be used with the Alfvén speed computed using the summed proton and ion densities, and with the sound speed computed from the summed pressures and densities of all species. 4. The ions can drive a non-resonant instability under certain circumstances. For low plasma beta, the threshold drift can be less than that required to destabilize electromagnetic modes, but damping from the Landau resonance can eliminate this instability altogether, unless $T_{\mathrm e}/T_{\mathrm p}\gg1$.

preprint2014arXiv

Traveling solar-wind bulk-velocity fluctuations and their effects on electron heating in the inner heliosphere

Ambient plasma electrons undergo strong heating in regions associated with compressive traveling interplanetary solar-wind bulk-velocity jumps due to their specific interactions with the jump-inherent electric fields. After thermalization of this energy gain per shock passage through the operation of the Buneman instability, strong electron heating occurs that substantially influences the radial electron temperature profile. We describe the reduction of the jump amplitude due to energy expended by the traveling jump structure. We consider three effects; namely energy loss due to heating of electrons, energy loss due to work done against the pick-up-ion pressure gradient, and an energy gain due to nonlinear jump steepening. Taking these effects into account, we show that the decrease in jump amplitude with solar distance is more pronounced when the initial jump amplitude is higher in the inner solar system. Independent of the initial jump amplitude, it eventually decreases with increasing distance to a value of the order of $ΔU/U\simeq 0.1$ at the position of the heliospheric termination shock, where $ΔU$ is the jump amplitude, and $U$ is the average solar-wind bulk velocity.The electron temperature, on the other hand, is strongly correlated with the initial jump amplitude, leading to electron temperatures between 6000 K and 20 000 K at distances beyond 50 AU. We compare our results with in-situ measurements of the electron-core temperature from the Ulysses spacecraft in the plane of the ecliptic for $1.5\, \mathrm{AU}\leq r \leq 5\,\mathrm{AU}$, where $r$ is the distance from the Sun. We find a very good agreement between our results and these observations, which corroborates our extrapolated predictions beyond $r=5\,\mathrm{AU}$.

preprint2013arXiv

A Parallel-Propagating Alfvénic Ion-Beam Instability in the High-Beta Solar Wind

We investigate the conditions under which parallel-propagating Alfvén/ion-cyclotron waves are driven unstable by an isotropic ($T_{\perp α} = T_{\parallelα}$) population of alpha particles drifting parallel to the magnetic field at an average speed $U_α$ with respect to the protons. We derive an approximate analytic condition for the minimum value of $U_α$ needed to excite this instability and refine this result using numerical solutions to the hot-plasma dispersion relation. When the alpha-particle number density is $\simeq 5%$ of the proton number density and the two species have similar thermal speeds, the instability requires that $β_{\rm p} \gtrsim 1$, where $β_{\rm p}$ is the ratio of the proton pressure to the magnetic pressure. For $1\lesssim β_{\mathrm p}\lesssim 12$, the minimum $U_α$ needed to excite this instability ranges from $0.7v_{\mathrm A}$ to $0.9v_{\mathrm A}$, where $v_{\mathrm A}$ is the Alfvén speed. This threshold is smaller than the threshold of $\simeq 1.2v_{\mathrm A}$ for the parallel magnetosonic instability, which was previously thought to have the lowest threshold of the alpha-particle beam instabilities at $β_{\mathrm p}\gtrsim 0.5$. We discuss the role of the parallel Alfvénic drift instability for the evolution of the alpha-particle drift speed in the solar wind. We also analyze measurements from the \emph{Wind} spacecraft's Faraday cups and show that the $U_α$ values measured in solar-wind streams with $T_{\perp α}\approx T_{\parallelα}$ are approximately bounded from above by the threshold of the parallel Alfvénic instability.

preprint2013arXiv

Instabilities Driven by the Drift and Temperature Anisotropy of Alpha Particles in the Solar Wind

We investigate the conditions under which parallel-propagating Alfvén/ion-cyclotron (A/IC) waves and fast-magnetosonic/whistler (FM/W) waves are driven unstable by the differential flow and temperature anisotropy of alpha particles in the solar wind. We focus on the limit in which $w_{\parallel α} \gtrsim 0.25 v_{\mathrm A}$, where $w_{\parallel α} $ is the parallel alpha-particle thermal speed and $v_{\mathrm A}$ is the Alfvén speed. We derive analytic expressions for the instability thresholds of these waves, which show, e.g., how the minimum unstable alpha-particle beam speed depends upon $w_{\parallel α}/v_{\mathrm A}$, the degree of alpha-particle temperature anisotropy, and the alpha-to-proton temperature ratio. We validate our analytical results using numerical solutions to the full hot-plasma dispersion relation. Consistent with previous work, we find that temperature anisotropy allows A/IC waves and FM/W waves to become unstable at significantly lower values of the alpha-particle beam speed $U_α$ than in the isotropic-temperature case. Likewise, differential flow lowers the minimum temperature anisotropy needed to excite A/IC or FM/W waves relative to the case in which $U_α=0$. We discuss the relevance of our results to alpha particles in the solar wind near 1 AU.

preprint2013arXiv

Limits on Alpha Particle Temperature Anisotropy and Differential Flow from Kinetic Instabilities: Solar Wind Observations

Previous studies have shown that the observed temperature anisotropies of protons and alpha particles in the solar wind are constrained by theoretical thresholds for pressure-anisotropy-driven instabilities such as the Alfvén/ion-cyclotron (A/IC) and fast-magnetosonic/whistler (FM/W) instabilities. In this letter, we use a long period of in-situ measurements provided by the {\em Wind} spacecraft's Faraday cups to investigate the combined constraint on the alpha-proton differential flow velocity and the alpha-particle temperature anisotropy due to A/IC and FM/W instabilities. We show that the majority of the data are constrained to lie within the region of parameter space in which A/IC and FM/W waves are either stable or have extremely low growth rates. In the minority of observed cases in which the growth rate of the A/IC (FM/W) instability is comparatively large, we find relatively higher values of $T_{\perpα}/T_{\perp p}$ ($T_{\parallelα}/T_{\parallel p}$) when alpha-proton differential flow velocity is small, where $T_{\perpα}$ and $T_{\perp p}$ ($T_{\parallelα}$ and $T_{\parallel p}$) are the perpendicular (parallel) temperatures of alpha particles and protons. We conjecture that this observed feature might arise from preferential alpha-particle heating which can drive the alpha particles beyond the instability thresholds.

preprint2013arXiv

The Dispersion Relations and Instability Thresholds of Oblique Plasma Modes in the Presence of an Ion Beam

An ion beam can destabilize Alfvén/ion-cyclotron waves and magnetosonic/whistler waves if the beam speed is sufficiently large. Numerical solutions of the hot-plasma dispersion relation have previously shown that the minimum beam speed required to excite such instabilities is significantly smaller for oblique modes with $\vec k \times \vec B_0\neq 0$ than for parallel-propagating modes with $\vec k \times \vec B_0 = 0$, where $\vec k$ is the wavevector and $\vec B_0$ is the background magnetic field. In this paper, we explain this difference within the framework of quasilinear theory, focusing on low-$β$ plasmas. We begin by deriving, in the cold-plasma approximation, the dispersion relation and polarization properties of both oblique and parallel-propagating waves in the presence of an ion beam. We then show how the instability thresholds of the different wave branches can be deduced from the wave--particle resonance condition, the conservation of particle energy in the wave frame, the sign (positive or negative) of the wave energy, and the wave polarization. We also provide a graphical description of the different conditions under which Landau resonance and cyclotron resonance destabilize Alfvén/ion-cyclotron waves in the presence of an ion beam. We draw upon our results to discuss the types of instabilities that may limit the differential flow of alpha particles in the solar wind.

preprint2012arXiv

Compressive high-frequency waves riding on an Alfvén/ion-cyclotron wave in a multi-fluid plasma

In this paper, we study the weakly-compressive high-frequency plasma waves which are superposed on a large-amplitude Alfvén wave in a multi-fluid plasma consisting of protons, electrons, and alpha particles. For these waves, the plasma environment is inhomogenous due to the presence of the low-frequency Alfvén wave with a large amplitude, a situation that may apply to space plasmas such as the solar corona and solar wind. The dispersion relation of the plasma waves is determined from a linear stability analysis using a new eigenvalue method that is employed to solve the set of differential wave equations which describe the propagation of plasma waves along the direction of the constant component of the Alfvén wave magnetic field. This approach also allows one to consider weak compressive effects. In the presence of the background Alfvén wave, the dispersion branches obtained differ significantly from the situation of a uniform plasma. Due to compressibility, acoustic waves are excited and couplings between various modes occur, and even an instability of the compressive mode. In a kinetic treatment, these plasma waves would be natural candidates for Landau-resonant wave-particle interactions, and may thus via their damping lead to particle heating.

preprint2012arXiv

Kinetic cascade beyond magnetohydrodynamics of solar wind turbulence in two-dimensional hybrid simulations

The nature of solar wind turbulence in the dissipation range at scales much smaller than the large MHD scales remains under debate. Here a two-dimensional model based on the hybrid code abbreviated as A.I.K.E.F. is presented, which treats massive ions as particles obeying the kinetic Vlasov equation and massless electrons as a neutralizing fluid. Up to a certain wavenumber in the MHD regime, the numerical system is initialized by assuming a superposition of isotropic Alfvén waves with amplitudes that follow the empirically confirmed spectral law of Kolmogorov. Then turbulence develops and energy cascades into the dispersive spectral range, where also dissipative effects occur. Under typical solar wind conditions, weak turbulence develops as a superposition of normal modes in the kinetic regime. Spectral analysis in the direction parallel to the background magnetic field reveals a cascade of left-handed Alfvén/ion-cyclotron waves up to wave vectors where their resonant absorption sets in, as well as a continuing cascade of right-handed fast-mode and whistler waves. Perpendicular to the background field, a broad turbulent spectrum is found to be built up of fluctuations having a strong compressive component. Ion-Bernstein waves seem to be possible normal modes in this propagation direction for lower driving amplitudes. Also signatures of short-scale pressure-balanced structures (very oblique slow-mode waves) are found.

preprint2012arXiv

On nonlinear Alfvén-cyclotron waves in multi-species plasma

Large-amplitude Alfvén waves are ubiquitous in space plasmas and a main component of magnetohydrodynamic (MHD) turbulence in the heliosphere. As pump waves they are prone to parametric instability by which they can generate cyclotron and acoustic daughter waves. Here we revisit a related process within the framework of the multi-fluid equations for a plasma consisting of many species. The nonlinear coupling of the Alfvén wave to acoustic waves is studied, and a set of compressive and coupled wave equations for the transverse magnetic field and longitudinal electric field is derived for waves propagating along the mean-field direction. It turns out that slightly compressive Alfvén waves exert, through induced gyro-radius and kinetic-energy modulations, an electromotive force on the particles in association with a longitudinal electric field, which has a potential that is given by the gradient of the transverse kinetic energy of the particles gyrating about the mean field. This in turn drives electric fluctuations (sound and ion-acoustic waves) along the mean magnetic field, which can nonlinearly react back on the transverse magnetic field. Mutually coupled Alfvén-cyclotron-acoustic waves are thus excited, a nonlinear process that can drive a cascade of wave energy in the plasma and may generate compressive microturbulence. These driven electric fluctuations might have consequences for the dissipation of MHD turbulence and, thus, for the heating and acceleration of particles in the solar wind.

preprint2012arXiv

Parametric decay of oblique Alfvén waves in two-dimensional hybrid simulations

Certain types of plasma waves are known to become parametrically unstable under specific plasma conditions, in which the pump wave will decay into several daughter waves with different wavenumbers and frequencies. In the past, the related plasma instabilities have been treated analytically for various parameter regimes and by use of various numerical methods, yet the oblique propagation with respect to the background magnetic field has rarely been dealt with in two dimensions, mainly because of the high computational demand. Here we present a hybrid-simulation study of the parametric decay of a moderately oblique Alfvén wave having elliptical polarization. It is found that such a compressive wave can decay into waves with higher and lower wavenumbers than the pump.

preprint2011arXiv

Apparent temperature anisotropies due to wave activity in the solar wind

The fast solar wind is a collisionless plasma permeated by plasma waves on many different scales. A plasma wave represents the natural interplay between the periodic changes of the electromagnetic field and the associated coherent motions of the plasma particles. In this paper, a model velocity distribution function is derived for a plasma in a single, coherent, large-amplitude wave. This model allows one to study the kinetic effects of wave motions on particle distributions. They are by in-situ spacecraft measured by counting, over a certain sampling time, the particles coming from various directions and having different energies. We compare our results with the measurements by the Helios spacecraft, and thus find that by assuming high wave activity we are able to explain key observed features of the measured distributions within the framework of our model. We also address the recent discussions on nonresonant wave--particle interactions and apparent heating. The applied time-averaging procedure leads to an apparent ion temperature anisotropy which is connected but not identical to the intrinsic temperature of the underlying distribution function.