-
Anisotropic Heating and Parallel Heat Flux in Electron-only Magnetic Reconnection with Intense Guide Fields
Authors:
Jincai Ren,
Giuseppe Arrò,
Maria Elena Innocenti,
Giovanni Lapenta
Abstract:
Electron-only reconnection is a process recently observed in the Earth's magnetosheath, where magnetic reconnection occurs at electron kinetic scales with only electrons involved. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales. This paper investigates electron-only reconnection under different inte…
▽ More
Electron-only reconnection is a process recently observed in the Earth's magnetosheath, where magnetic reconnection occurs at electron kinetic scales with only electrons involved. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales. This paper investigates electron-only reconnection under different intensities of strong guide fields via two-dimensional fully kinetic Particle-in-Cell (PIC) simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. The quadrupolar structures of the out-of-plane magnetic field and asymmetric structures of density are found to become narrower in space as the guide field increases. Electron velocities are considerably larger than ion velocities, and the motion of both species is notably affected by the strength of guide fields. Similarly to previous experiments studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a six-polar to an eight-polar structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and strong parallel electron heat flux is observed.
△ Less
Submitted 27 June, 2024; v1 submitted 15 June, 2024;
originally announced June 2024.
-
Applications of Fast Magnetic Reconnection Models to the Atmospheres of the Sun and Protoplanetary Disks
Authors:
Fulvia Pucci,
Alkendra Singh,
Uma Gorti,
Marco Velli,
Neal Turner,
Disha Varshney,
Maria Elena Innocenti
Abstract:
Partially-ionized plasmas consist of charged and neutral particles whose mutual collisions modify magnetic reconnection compared with the fully-ionized case. The collisions alter the rate and locations of the magnetic dissipation heating and the distribution of energies among the particles accelerated into the non-thermal tail. We examine the collisional regimes for the onset of fast reconnection…
▽ More
Partially-ionized plasmas consist of charged and neutral particles whose mutual collisions modify magnetic reconnection compared with the fully-ionized case. The collisions alter the rate and locations of the magnetic dissipation heating and the distribution of energies among the particles accelerated into the non-thermal tail. We examine the collisional regimes for the onset of fast reconnection in two environments: the partially-ionized layers of the solar atmosphere and the protoplanetary disks that are the birthplaces for planets around young stars. In both these environments, magnetic nulls readily develop into resistive current sheets in the regime where the charged and neutral particles are fully coupled by collisions, but the current sheets quickly break down under the ideal tearing instability. The current sheets collapse repeatedly, forming magnetic islands at successively smaller scales, till they enter a collisionally-decoupled regime where the magnetic energy is rapidly turned into heat and charged-particle kinetic energy. Small-scale, decoupled fast reconnection in the solar atmosphere may lead to preferential heating and energization of ions and electrons that escape into the corona. In protoplanetary disks such reconnection causes localized heating in the atmospheric layers that produce much of the infrared atomic and molecular line emission observed with the Spitzer and James Webb Space Telescopes.
△ Less
Submitted 14 May, 2024;
originally announced May 2024.
-
Fully Kinetic Simulations of Proton-Beam-Driven Instabilities from Parker Solar Probe Observations
Authors:
Luca Pezzini,
Andrei N. Zhukov,
Fabio Bacchini,
Giuseppe Arrò,
Rodrigo A. López,
Alfredo Micera,
Maria Elena Innocenti,
Giovanni Lapenta
Abstract:
The expanding solar wind plasma ubiquitously exhibits anisotropic non-thermal particle velocity distributions. Typically, proton Velocity Distribution Functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes ex…
▽ More
The expanding solar wind plasma ubiquitously exhibits anisotropic non-thermal particle velocity distributions. Typically, proton Velocity Distribution Functions (VDFs) show the presence of a core and a field-aligned beam. Novel observations made by Parker Solar Probe (PSP) in the innermost heliosphere have revealed new complex features in the proton VDFs, namely anisotropic beams that sometimes experience perpendicular diffusion. This phenomenon gives rise to VDFs that resemble a "hammerhead". In this study, we use a 2.5D fully kinetic simulation to investigate the stability of proton VDFs with anisotropic beams observed by PSP. Our setup consists of a core and an anisotropic beam populations that drift with respect to each other. This configuration triggers a proton-beam instability from which nearly parallel fast magnetosonic modes develop. Our results demonstrate that before this instability reaches saturation, the waves resonantly interact with the beam protons, causing significant perpendicular heating at the expense of the parallel temperature. Furthermore, the proton perpendicular heating induces a hammerhead-like shape in the resulting VDF. Our results suggest that this mechanism probably contribute to producing the observed hammerhead distributions.
△ Less
Submitted 13 May, 2024;
originally announced May 2024.
-
Electron-driven instabilities in the solar wind
Authors:
Daniel Verscharen,
Benjamin D. G. Chandran,
Elisabetta Boella,
Jasper Halekas,
Maria Elena Innocenti,
Vamsee K. Jagarlamudi,
Alfredo Micera,
Viviane Pierrard,
Stepan Stverak,
Ivan Y. Vasko,
Marco Velli,
Phyllis L. Whittlesey
Abstract:
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-insta…
▽ More
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.
△ Less
Submitted 21 June, 2022;
originally announced June 2022.
-
Unsupervised classification of simulated magnetospheric regions
Authors:
Maria Elena Innocenti,
Jorge Amaya,
Joachim Raeder,
Romain Dupuis,
Banafsheh Ferdousi,
Giovanni Lapenta
Abstract:
In magnetospheric missions, burst mode data sampling should be triggered in the presence of processes of scientific or operational interest. We present an unsupervised classification method for magnetospheric regions, that could constitute the first-step of a multi-step method for the automatic identification of magnetospheric processes of interest. Our method is based on Self Organizing Maps (SOM…
▽ More
In magnetospheric missions, burst mode data sampling should be triggered in the presence of processes of scientific or operational interest. We present an unsupervised classification method for magnetospheric regions, that could constitute the first-step of a multi-step method for the automatic identification of magnetospheric processes of interest. Our method is based on Self Organizing Maps (SOMs), and we test it preliminarily on data points from global magnetospheric simulations obtained with the OpenGGCM-CTIM-RCM code. The dimensionality of the data is reduced with Principal Component Analysis before classification. The classification relies exclusively on local plasma properties at the selected data points, without information on their neighborhood or on their temporal evolution. We classify the SOM nodes into an automatically selected number of classes, and we obtain clusters that map to well defined magnetospheric regions. We validate our classification results by plotting the classified data in the simulated space and by comparing with K-means classification. For the sake of result interpretability, we examine the SOM feature maps (magnetospheric variables are called features in the context of classification), and we use them to unlock information on the clusters. We repeat the classification experiments using different sets of features, we quantitatively compare different classification results, and we obtain insights on which magnetospheric variables make more effective features for unsupervised classification.
△ Less
Submitted 10 September, 2021;
originally announced September 2021.
-
On the role of solar wind expansion as a source of whistler waves: scattering of suprathermal electrons and heat flux regulation in the inner heliosphere
Authors:
A. Micera,
A. N. Zhukov,
R. A. López,
E. Boella,
A. Tenerani,
M. Velli,
G. Lapenta,
M. E. Innocenti
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 Parker Solar Probe observations during its first perihelion at 0.166 au, consisting of a dense core and a…
▽ More
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 Parker Solar Probe observations during its first perihelion at 0.166 au, consisting of a dense core and an anti-sunward 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. 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.
△ Less
Submitted 1 July, 2021; v1 submitted 30 June, 2021;
originally announced June 2021.
-
Particle-In-Cell simulation of whistler heat flux instabilities in the solar wind: heat flux regulation and electron halo formation
Authors:
A. Micera,
A. N. Zhukov,
R. A. López,
M. E. Innocenti,
M. Lazar,
E. Boella,
G. Lapenta
Abstract:
We present results of 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 Parke…
▽ More
We present results of 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, 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, towards 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 velocity distribution function 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.
△ Less
Submitted 22 October, 2020; v1 submitted 21 October, 2020;
originally announced October 2020.
-
Domain of Influence analysis: implications for Data Assimilation in space weather forecasting
Authors:
Dimitrios Millas,
Maria Elena Innocenti,
Brecht Laperre,
Joachim Raeder,
Stefaan Poedts,
Giovanni Lapenta
Abstract:
Solar activity, ranging from the background solar wind to energetic coronal mass ejections (CMEs), is the main driver of the conditions in the interplanetary space and in the terrestrial space environment, known as space weather. A better understanding of the Sun-Earth connection carries enormous potential to mitigate negative space weather effects with economic and social benefits. Effective spac…
▽ More
Solar activity, ranging from the background solar wind to energetic coronal mass ejections (CMEs), is the main driver of the conditions in the interplanetary space and in the terrestrial space environment, known as space weather. A better understanding of the Sun-Earth connection carries enormous potential to mitigate negative space weather effects with economic and social benefits. Effective space weather forecasting relies on data and models. In this paper, we discuss some of the most used space weather models, and propose suitable locations for data gathering with space weather purposes. We report on the application of \textit{Representer analysis (RA)} and \textit{Domain of Influence (DOI) analysis} to three models simulating different stages of the Sun-Earth connection: the OpenGGCM and Tsyganenko models, focusing on solar wind - magnetosphere interaction, and the PLUTO model, used to simulate CME propagation in interplanetary space. Our analysis is promising for space weather purposes for several reasons. First, we obtain quantitative information about the most useful locations of observation points, such as solar wind monitors. For example, we find that the absolute values of the DOI are extremely low in the magnetospheric plasma sheet. Since knowledge of that particular sub-system is crucial for space weather, enhanced monitoring of the region would be most beneficial. Second, we are able to better characterize the models. Although the current analysis focuses on spatial rather than temporal correlations, we find that time-independent models are less useful for Data Assimilation activities than time-dependent models. Third, we take the first steps towards the ambitious goal of identifying the most relevant heliospheric parameters for modelling CME propagation in the heliosphere, their arrival time, and their geoeffectiveness at Earth.
△ Less
Submitted 9 September, 2020;
originally announced September 2020.
-
Collisionless heat flux regulation via electron firehose instability in presence of a core and suprathermal population in the expanding solar wind
Authors:
Maria Elena Innocenti,
Elisabetta Boella,
Anna Tenerani,
Marco Velli
Abstract:
The evolution of the electron heat flux in the solar wind is regulated by the interplay between several effects: solar wind expansion, that can potentially drive velocity-space instabilties, turbulence and wave-particle interactions, and, possibly, collisions. Here we address the respective role played by the solar wind expansion and the electron firehose instability, developing in the presence of…
▽ More
The evolution of the electron heat flux in the solar wind is regulated by the interplay between several effects: solar wind expansion, that can potentially drive velocity-space instabilties, turbulence and wave-particle interactions, and, possibly, collisions. Here we address the respective role played by the solar wind expansion and the electron firehose instability, developing in the presence of multiple electron populations, in regulating the heat flux. We carry out fully kinetic, Expanding Box Model simulations and separately analyze the enthalpy, bulk and velocity distribution function skewness contributions for each of the electron species. We observe that the key factor determining electron energy flux evolution is the reduction of the drift velocity of the electron populations in the rest frame of the solar wind. In our simulations, redistribution of the electron thermal energy from the parallel to the perpendicular direction after the onset of the electron firehose instability is observed. However, this process seems to impact energy flux evolution only minimally. Hence, reduction of the electron species drift velocity in the solar wind frame appears to directly correlate with efficiency for heat flux instabilities
△ Less
Submitted 14 July, 2020;
originally announced July 2020.
-
Oblique tearing mode instability: guide field and Hall effect
Authors:
Chen Shi,
Marco Velli,
Fulvia Pucci,
Anna Tenerani,
Maria Elena Innocenti
Abstract:
The tearing mode instability is one important mechanism that may explain the triggering of fast magnetic reconnection in astrophysical plasmas such as the solar corona and the Earth's magnetosphere. In this paper, the linear stability analysis of the tearing mode is carried out for a current sheet in the presence of a guide field, including the Hall effect. We show that the presence of a strong gu…
▽ More
The tearing mode instability is one important mechanism that may explain the triggering of fast magnetic reconnection in astrophysical plasmas such as the solar corona and the Earth's magnetosphere. In this paper, the linear stability analysis of the tearing mode is carried out for a current sheet in the presence of a guide field, including the Hall effect. We show that the presence of a strong guide field does not modify the most unstable mode in the two-dimensional wave vector space orthogonal to the current gradient direction, which remains the fastest growing parallel mode. With the Hall effect, the inclusion of a guide field turns the non-dispersive propagation along the guide field direction to a dispersive one. The oblique modes have a wave-like structure along the normal direction of the current sheet and a strong guide field suppresses this structure while making the eigen-functions asymmetric.
△ Less
Submitted 4 September, 2020; v1 submitted 1 July, 2020;
originally announced July 2020.
-
Visualizing and Interpreting Unsupervised Solar Wind Classifications
Authors:
Jorge Amaya,
Romain Dupuis,
Maria Elena Innocenti,
Giovanni Lapenta
Abstract:
One of the goals of machine learning is to eliminate tedious and arduous repetitive work. The manual and semi-automatic classification of millions of hours of solar wind data from multiple missions can be replaced by automatic algorithms that can discover, in mountains of multi-dimensional data, the real differences in the solar wind properties. In this paper we present how unsupervised clustering…
▽ More
One of the goals of machine learning is to eliminate tedious and arduous repetitive work. The manual and semi-automatic classification of millions of hours of solar wind data from multiple missions can be replaced by automatic algorithms that can discover, in mountains of multi-dimensional data, the real differences in the solar wind properties. In this paper we present how unsupervised clustering techniques can be used to segregate different types of solar wind. We propose the use of advanced data reduction methods to pre-process the data, and we introduce the use of Self-Organizing Maps to visualize and interpret 14 years of ACE data. Finally, we show how these techniques can potentially be used to uncover hidden information, and how they compare with previous manual and automatic categorizations.
△ Less
Submitted 18 August, 2020; v1 submitted 28 April, 2020;
originally announced April 2020.
-
On the Electron Agyrotropy during Rapid Asymmetric Magnetic Island Coalescence in Presence of a Guide Field
Authors:
Emanuele Cazzola,
Maria Elena Innocenti,
Martin V. Goldman,
David L. Newman,
Stefano Markidis,
Giovanni Lapenta
Abstract:
We present an analysis of the properties of the electron velocity distribution during island coalescence in asymmetric reconnection with and without guide field. In a previous study, three main domains were identified, in the case without guide field, as X-, D- and M-regions featuring different reconnection evolutions {Cazzola et al. 2015). These regions are also identified here in the case with g…
▽ More
We present an analysis of the properties of the electron velocity distribution during island coalescence in asymmetric reconnection with and without guide field. In a previous study, three main domains were identified, in the case without guide field, as X-, D- and M-regions featuring different reconnection evolutions {Cazzola et al. 2015). These regions are also identified here in the case with guide field. We study the departure from isotropic and gyrotropic behavior by means of different robust detection algorithms proposed in the literature. While in the case without guide field these metrics show an overall agreement, when the guide field is present a discrepancy in the agyrotropy within some relevant regions is observed, such as at the separatrices and inside magnetic islands. Moreover, in light of the new observations from the Multiscale MagnetoSpheric mission, an analysis of the electron velocity phase-space in these domains is presented.
△ Less
Submitted 2 August, 2016;
originally announced August 2016.
-
Magnetic null points in kinetic simulations of space plasmas
Authors:
Vyacheslav Olshevsky,
Jan Deca,
Andrey Divin,
Ivy Bo Peng,
Stefano Markidis,
Maria Elena Innocenti,
Emanuele Cazzola,
Giovanni Lapenta
Abstract:
We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic Particle-in-Cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interactin…
▽ More
We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic Particle-in-Cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interacting with the solar wind; and a relaxing turbulent configuration with multiple null points. Spiral nulls are more likely created in space plasmas: in all our simulations except lunar magnetic anomaly and quadrupolar mini-magnetosphere the number of spiral nulls prevails over the number of radial nulls by a factor of 3-9. We show that often magnetic nulls do not indicate the regions of intensive energy dissipation. Energy dissipation events caused by topological bifurcations at radial nulls are rather rare and short-lived. The so-called X-lines formed by the radial nulls in the Harris current sheet and lunar magnetic anomaly simulations are rather stable and do not exhibit any energy dissipation. Energy dissipation is more powerful in the vicinity of spiral nulls enclosed by magnetic flux ropes with strong currents at their axes (their cross-sections resemble 2D magnetic islands). These null lines reminiscent of Z-pinches efficiently dissipate magnetic energy due to secondary instabilities such as the two-stream or kinking instability, accompanied by changes in magnetic topology. Current enhancements accompanied by spiral nulls may signal magnetic energy conversion sites in the observational data.
△ Less
Submitted 7 December, 2015;
originally announced December 2015.
-
On the electron dynamics during island coalescence in asymmetric magnetic reconnection
Authors:
Emanuele Cazzola,
Maria Elena Innocenti,
Stefano Markidis,
Martin V. Goldman,
David L. Newman,
Giovanni Lapenta
Abstract:
We present an analysis of the electron dynamics during rapid island merging in asymmetric magnetic reconnection. We consider a doubly periodic system with two asymmetric transitions. The upper layer is an asymmetric Harris sheet initially perturbed to promote a single reconnection site. The lower layer is a tangential discontinuity that promotes the formation of many X-points, separated by rapidly…
▽ More
We present an analysis of the electron dynamics during rapid island merging in asymmetric magnetic reconnection. We consider a doubly periodic system with two asymmetric transitions. The upper layer is an asymmetric Harris sheet initially perturbed to promote a single reconnection site. The lower layer is a tangential discontinuity that promotes the formation of many X-points, separated by rapidly merging islands. Across both layers the magnetic field and the density have a strong jump, but the pressure is held constant. Our analysis focuses on the consequences of electron energization during island coalescence. We focus first on the parallel and perpendicular components of the electron temperature to establish the presence of possible anisotropies and non-gyrotropies. Thanks to the direct comparison between the two different layers simulated, we can distinguish three main types of behavior characteristic of three different regions of interest. The first type represents the regions where traditional asymmetric reconnections take place without involving island merging. The second type of regions instead show reconnection events between two merging islands. Finally, the third regions identifies the regions between two diverging island and where typical signature of reconnection is not observed. Electrons in these latter regions additionally show a flat-top distribution resulting from the saturation of a two-stream instability generated by the two interacting electron beams from the two nearest reconnection points. Finally, the analysis of agyrotropy shows the presence of a distinct double structure laying all over the lower side facing the higher magnetic field region. The distinguishing features found for the three types of regions investigated provide clear indicators to the recently launched MMS NASA mission for investigating magnetopause reconnection involving multiple islands.
△ Less
Submitted 18 November, 2015;
originally announced November 2015.