Abstract In this paper, we introduce and discuss an exactly energy-conserving particle-in-cell method for arbitrary curvilinear coordinates. The flexibility provided by curvilinear coordinates ...enables the study of plasmas in complex-shaped domains by aligning the grid to the given geometry or by focusing grid resolution on regions of interest without overresolving the surrounding, potentially uninteresting domain. We have achieved this through the introduction of the metric tensor, the Jacobian matrix, and contravariant operators combined with an energy-conserving fully implicit solver. We demonstrate the method’s capabilities using a Python implementation to study several one- and two-dimensional test cases: the electrostatic two-stream instability, the electromagnetic Weibel instability, and the geomagnetic environment modeling reconnection challenge. The test results confirm the capability of our new method to reproduce theoretical expectations (e.g., instability growth rates) and the corresponding results obtained with a Cartesian uniform grid when using curvilinear grids. Simultaneously, we show that the method conserves energy to machine precision in all cases.
We compare relativistic particle integrators commonly used in plasma physics, showing several test cases relevant for astrophysics. Three explicit particle pushers are considered, namely, the Boris, ...Vay, and Higuera-Cary schemes. We also present a new relativistic fully implicit particle integrator that is energy conserving. Furthermore, a method based on the relativistic guiding center approximation is included. The algorithms are described such that they can be readily implemented in magnetohydrodynamics codes or Particle-in-Cell codes. Our comparison focuses on the strengths and key features of the particle integrators. We test the conservation of invariants of motion and the accuracy of particle drift dynamics in highly relativistic, mildly relativistic, and non-relativistic settings. The methods are compared in idealized test cases, i.e., without considering feedback onto the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields, fields, force-free fields, and setups relevant for high-energy astrophysics, e.g., a magnetic mirror, a magnetic dipole, and a magnetic null. These tests have direct relevance for particle acceleration in shocks and in magnetic reconnection.
The Magnetosphere Multiscale Mission (MMS) will provide the first opportunity to probe electron-scale physics during magnetic reconnection in Earth’s magnetopause and magnetotail. This article will ...address only tail reconnection—as a
non
-steady-state process in which the first reconnected field lines advance away from the
x
-point in
flux pile-up
fronts
directed Earthward and anti-Earthward. An up-to-date microscopic physical picture of electron and ion-scale collisionless tail reconnection processes is presented based on 2-D Particle-In-Cell (PIC) simulations initiated from a Harris current sheet and on Cluster and Themis measurements of tail reconnection. The successes and limitations of simulations when compared to measured reconnection are addressed in detail. The main focus is on particle and field diffusion region signatures in the tail reconnection geometry. The interpretation of these signatures is vital to enable spacecraft to identify physically significant reconnection events, to trigger meaningful data transfer from MMS to Earth and to construct a useful overall physical picture of tail reconnection. New simulation results and theoretical interpretations are presented for energy transport of particles and fields, for the size and shape of electron and ion diffusion regions, for processes occurring near the fronts and for the
j
×
B
(Hall) electric field.
Context.
Recent satellite measurements in the turbulent magnetosheath of Earth have given evidence of an unusual reconnection mechanism that is driven exclusively by electrons. This newly observed ...process was called electron-only reconnection, and its interplay with plasma turbulence is a matter of great debate.
Aims.
By using 2D-3V hybrid Vlasov–Maxwell simulations of freely decaying plasma turbulence, we study the role of electron-only reconnection in the development of plasma turbulence. In particular, we search for possible differences with respect to the turbulence associated with standard ion-coupled reconnection.
Methods.
We analyzed the structure functions of the turbulent magnetic field and ion fluid velocity fluctuations to characterize the structure and the intermittency properties of the turbulent energy cascade.
Results.
We find that the statistical properties of turbulent fluctuations associated with electron-only reconnection are consistent with those of turbulent fluctuations associated with standard ion-coupled reconnection, and no peculiar signature related to electron-only reconnection is found in the turbulence statistics. This result suggests that the turbulent energy cascade in a collisionless magnetized plasma does not depend on the specific mechanism associated with magnetic reconnection. The properties of the dissipation range are discussed as well, and we claim that only electrons contribute to the dissipation of magnetic field energy at sub-ion scales.
The process of magnetic reconnection when studied in nature or when modeled in 3D simulations differs in one key way from the standard 2D paradigmatic cartoon: it is accompanied by many fluctuations ...in the electromagnetic fields and plasma properties. We developed a diagnostics to study the spectrum of fluctuations in the various regions around a reconnection site. We define the regions in terms of the local value of the flux function that determines the distance from the reconnection site, with positive values in the outflow and negative values in the inflow. We find that fluctuations belong to two very different regimes depending on the local plasma beta (defined as the ratio of plasma and magnetic pressures). The first regime develops in the reconnection outflows where beta is high and it is characterized by a strong link between plasma and electromagnetic fluctuations, leading to momentum and energy exchanges via anomalous viscosity and resistivity. But there is a second, low-beta regime: it develops in the inflow and in the region around the separatrix surfaces, including the reconnection electron diffusion region itself. It is remarkable that this low-beta plasma, where the magnetic pressure dominates, remains laminar even though the electromagnetic fields are turbulent.
While it is certain that the fast solar wind originates from coronal holes, where and how the slow solar wind (SSW) is formed remains an outstanding question in solar physics even in the post-SOHO ...era. The quest for the SSW origin forms a major objective for the planned future missions such as the Solar Orbiter and Solar Probe Plus. Nonetheless, results from spacecraft data, combined with theoretical modeling, have helped to investigate many aspects of the SSW. Fundamental physical properties of the coronal plasma have been derived from spectroscopic and imaging remote-sensing data and in situ data, and these results have provided crucial insights for a deeper understanding of the origin and acceleration of the SSW. Advanced models of the SSW in coronal streamers and other structures have been developed using 3D MHD and multi-fluid equations.
Abstract
The role of solar wind expansion in generating whistler waves is investigated using the EB-iPic3D code, which models solar wind expansion self-consistently within a fully kinetic ...semi-implicit approach. The simulation is initialized with an electron velocity distribution function modeled after observations of the Parker Solar Probe during its first perihelion at 0.166 au, consisting of a dense core and an antisunward strahl. This distribution function is initially stable with respect to kinetic instabilities. Expansion drives the solar wind into successive regimes where whistler heat flux instabilities are triggered. These instabilities produce sunward whistler waves initially characterized by predominantly oblique propagation with respect to the interplanetary magnetic field. The excited waves interact with the electrons via resonant scattering processes. As a consequence, the strahl pitch angle distribution broadens and its drift velocity reduces. The strahl electrons are scattered in the direction perpendicular to the magnetic field, and an electron halo is formed. At a later stage, resonant electron firehose instability is triggered and further affects the electron temperature anisotropy as the solar wind expands. Wave–particle interaction processes are accompanied by a substantial reduction of the solar wind heat flux. The simulated whistler waves are in qualitative agreement with observations in terms of wave frequencies, amplitudes, and propagation angles. Our work proposes an explanation for the observations of oblique and parallel whistler waves in the solar wind. We conclude that solar wind expansion has to be factored in when trying to explain kinetic processes at different heliocentric distances.
We investigate current‐driven modes developing at jet fronts during collisionless reconnection. Initial evolution of the reconnection is simulated using conventional 2‐D setup starting from the ...Harris equilibrium. Three‐dimensional PIC calculations are implemented at later stages, when fronts are fully formed. Intense currents and enhanced wave activity are generated at the fronts because of the interaction of the fast flow plasma and denser ambient current sheet plasma. The study reveals that the lower hybrid drift instability develops quickly in the 3‐D simulation. The instability produces strong localized perpendicular electric fields, which are several times larger than the convective electric field at the front, in agreement with Time History of Events and Macroscale Interactions during Substorms observations. The instability generates waves, which escape the front edge and propagate into the undisturbed plasma ahead of the front. The parallel electron pressure is substantially larger in the 3‐D simulation compared to that of the 2‐D. In a time ∼Ωci−1, the instability forms a layer, which contains a mixture of the jet plasma and current sheet plasma. The results confirm that the lower hybrid drift instability is important for the front evolution and electron energization.
Key Points
Lower hybrid drift instability is excited at the reconnection jet front edge
LHDI heats electrons predominantly parallel to the magnetic field
The front acts as a source of waves traveling ahead of the front
Abstract
We present results of a two-dimensional fully kinetic particle-in-cell simulation in order to shed light on the role of whistler waves in the scattering of strahl electrons and in the ...heat-flux regulation in the solar wind. We model the electron velocity distribution function as initially composed of core and strahl populations as typically encountered in the near-Sun solar wind as observed by Parker Solar Probe. We demonstrate that, as a consequence of the evolution of the electron velocity distribution function (VDF), two branches of the whistler heat-flux instability can be excited, which can drive whistler waves propagating in the direction oblique or parallel to the background magnetic field. First, oblique whistler waves induce pitch-angle scattering of strahl electrons, toward higher perpendicular velocities. This leads to the broadening of the strahl pitch-angle distribution and hence to the formation of a halo-like population at the expense of the strahl. Later on, the electron VDF experiences the effect of parallel whistler waves, which contributes to the redistribution of the particles scattered in the perpendicular direction into a more symmetric halo, in agreement with observations. Simulation results show a remarkable agreement with the linear theory of the oblique whistler heat-flux instability. The process is accompanied by a significant decrease of the heat flux carried by the strahl population.