BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

G. Lapenta

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

ResearcherUnclaimed source record
physics.plasm-phphysics.space-phastro-ph.SRastro-ph.EPphysics.comp-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

12works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

12 published item(s)

preprint2022arXiv

Mixing the solar wind proton and electron scales. Theory and 2D-PIC simulations of firehose instability

Firehose-like instabilities (FIs) are cited in multiple astrophysical applications. Of particular interest are the kinetic manifestations in weakly-collisional or even collisionless plasmas, where these instabilities are expected to contribute to the evolution of macroscopic parameters. Relatively recent studies have initiated a realistic description of FIs, as induced by the interplay of both species, electrons and protons, dominant in the solar wind plasma. This work complements the current knowledge with new insights from linear theory and the first disclosures from 2D PIC simulations, identifying the fastest growing modes near the instability thresholds and their long-run consequences on the anisotropic distributions. Thus, unlike previous setups, these conditions are favorable to those aperiodic branches that propagate obliquely to the uniform magnetic field, with (maximum) growth rates higher than periodic, quasi-parallel modes. Theoretical predictions are, in general, confirmed by the simulations. The aperiodic electron FI (a-EFI) remains unaffected by the proton anisotropy, and saturates rapidly at low-level fluctuations. Regarding the firehose instability at proton scales, we see a stronger competition between the periodic and aperiodic branches. For the parameters chosen in our analysis, the a-PFI is excited before than the p-PFI, with the latter reaching a significantly higher fluctuation power. However, both branches are significantly enhanced by the presence of anisotropic electrons. The interplay between EFIs and PFIs also produces a more pronounced proton isotropization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-beam Energy Moments of Multibeam Particle Velocity Distributions

High resolution electron and ion velocity distributions, f(v), which consist of N effectively disjoint beams, have been measured by NASA's Magnetospheric Multi-Scale Mission (MMS) observatories and in reconnection simulations. Commonly used standard velocity moments generally assume a single mean-flow-velocity for the entire distribution, which can lead to counterintuitive results for a multibeam f(v). An example is the (false) standard thermal energy moment of a pair of equal and opposite cold particle beams, which is nonzero even though each beam has zero thermal energy. By contrast, a multibeam moment of two or more beams has no false thermal energy. A multibeam moment is obtained by taking a standard moment of each beam and then summing over beams. In this paper we will generalize these notions, explore their consequences and apply them to an f(v) which is sum of tri-Maxwellians. Both standard and multibeam energy moments have coherent and incoherent forms. Examples of incoherent moments are the thermal energy density, the pressure and the thermal energy flux (enthalpy flux plus heat flux). Corresponding coherent moments are the bulk kinetic energy density, the RAM pressure and the bulk kinetic energy flux. The false part of an incoherent moment is defined as the difference between the standard incoherent moment and the corresponding multibeam moment. The sum of a pair of corresponding coherent and incoherent moments will be called the undecomposed moment. Undecomposed moments are independent of whether the sum is standard or multibeam and therefore have advantages when studying moments of measured f(v).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Particle-In-Cell simulations of the parallel proton firehose instability influenced by the electron temperature anisotropy in solar wind conditions

In situ observations of the solar wind show a limited level of particle temperature anisotropy with respect to the interplanetary magnetic field direction. Kinetic electromagnetic instabilities are efficient to prevent the excessive growth of the anisotropy of particle velocity distribution functions. Among them, the firehose instabilities are often considered to prevent the increase of the parallel temperature and hence to shape the velocity distribution functions of electrons and protons in the solar wind. We present a non-linear modeling of the parallel firehose instability, retaining a kinetic description for both the electrons and protons. One-dimensional (1D) fully kinetic Particle-In-Cell simulations using the Energy Conserving semi-implicit method (ECsim) are performed to clarify the role of the electron temperature anisotropy in the development of the parallel proton firehose instability. We found that in the presence of an electron temperature anisotropy, such that the temperature parallel to the background magnetic field is higher than the temperature in the perpendicular direction, the onset of the parallel proton firehose instability occurs earlier and its growth rate is faster. The enhanced wave fluctuations contribute to the particle scattering reducing the temperature anisotropy to a stable, nearly isotropic state. The simulation results compare well with linear theory. A test case of 1D simulations at oblique angles with respect to the magnetic field is also considered, as a first step to study the cumulative effect of protons and electrons on the full spectrum of instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Local regimes of turbulence in 3D magnetic reconnection

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 much 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 form 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 pressure). The first regime develops in the reconnection outflows where beta is high and 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, remain laminar even though the electromagnetic fields are turbulent.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Particle Control in Phase Space by Global K-Means Clustering

We devise and explore an iterative optimization procedure for controlling particle populations in particle-in-cell (PIC) codes via merging and splitting of computational macro-particles. Our approach, is to compute an optimal representation of the global particle phase space structure while decreasing or increasing the entire particle population, based on k-means clustering of the data. In essence the procedure amounts to merging or splitting particles by statistical means, throughout the entire simulation volume in question, while minimizing a 6-dimensional total distance measure to preserve the physics. Particle merging is by far the most demanding procedure when considering conservation laws of physics; it amounts to lossy compression of particle phase space data. We demonstrate that our k-means approach conserves energy and momentum to high accuracy, even for high compression ratios, $\mathcal{R} \approx 3$ --- \emph{i.e.}, $N_{f} \lesssim 0.33N_{i}$. Interestingly, we find that an accurate particle splitting step can be performed using k-means as well; this from an argument of symmetry. The split solution, using k-means, places splitted particles optimally, to obtain maximal spanning on the phase space manifold. Implementation and testing is done using an electromagnetic PIC code, the \ppcode. Nonetheless, the k-means framework is general; it is not limited to Vlasov-Maxwell type PIC codes. We discuss advantages and drawbacks of this optimal phase space reconstruction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Where should MMS look for electron diffusion regions?

A great possible achievement for the MMS mission would be crossing electron diffusion regions (EDR). EDR are regions in proximity of reconnection sites where electrons decouple from field lines, breaking the frozen in condition. Decades of research on reconnection have produced a widely shared map of where EDRs are. We expect reconnection to take place around a so called x-point formed by the intersection of the separatrices dividing inflowing from outflowing plasma. The EDR forms around this x-point as a small electron scale box nested inside a larger ion diffusion region. But this point of view is based on a 2D mentality. We have recently proposed that once the problem is considered in full 3D, secondary reconnection events can form [Lapenta et al., Nature Physics, 11, 690, 2015] in the outflow regions even far downstream from the primary reconnection site. We revisit here this new idea confirming that even using additional indicators of reconnection and even considering longer periods and wider distances the conclusion remains true: secondary reconnection sites form downstream of a reconnection outflow causing a sort of chain reaction of cascading reconnection sites. If we are right, MMS will have an interesting journey even when not crossing necessarily the primary site. The chances are greatly increased that even if missing a primary site during an orbit, MMS could stumble instead on one of these secondary sites.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Plasma physical parameters along Coronal Mass Ejection-driven shocks: I observations

In this work UV and white light (WL) coronagraphic data are combined to derive the full set of plasma physical parameters along the front of a shock driven by a Coronal Mass Ejection. Pre-shock plasma density, shock compression ratio, speed and inclination angle are estimated from WL data, while pre-shock plasma temperature and outflow velocity are derived from UV data. The Rankine-Hugoniot (RH) equations for the general case of an oblique shock are then applied at three points along the front located between $2.2-2.6$ R$_\odot$ at the shock nose and at the two flanks. Stronger field deflection (by $\sim 46^\circ$), plasma compression (factor $\sim 2.7$) and heating (factor $\sim 12$) occur at the nose, while heating at the flanks is more moderate (factor $1.5-3.0$). Starting from a pre-shock corona where protons and electrons have about the same temperature ($T_p \sim T_e \sim 1.5 \cdot 10^6$ K), temperature increases derived with RH equations could better represent the protons heating (by dissipation across the shock), while the temperature increase implied by adiabatic compression (factor $\sim 2$ at the nose, $\sim 1.2-1.5$ at the flanks) could be more representative of electrons heating: the transit of the shock causes a decoupling between electron and proton temperatures. Derived magnetic field vector rotations imply a draping of field lines around the expanding flux rope. The shock turns out to be super-critical (sub-critical) at the nose (at the flanks), where derived post-shock plasma parameters can be very well approximated with those derived by assuming a parallel (perpendicular) shock.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Signatures of Secondary Collisionless Magnetic Reconnection Driven by Kink Instability of a Flux Rope

The kinetic features of secondary magnetic reconnection in a single flux rope undergoing internal kink instability are studied by means of three-dimensional Particle-in-Cell simulations. Several signatures of secondary magnetic reconnection are identified in the plane perpendicular to the flux rope: a quadrupolar electron and ion density structure and a bipolar Hall magnetic field develop in proximity of the reconnection region. The most intense electric fields form perpendicularly to the local magnetic field, and a reconnection electric field is identified in the plane perpendicular to the flux rope. An electron current develops along the reconnection line in the opposite direction of the electron current supporting the flux rope magnetic field structure. Along the reconnection line, several bipolar structures of the electric field parallel to the magnetic field occur making the magnetic reconnection region turbulent. The reported signatures of secondary magnetic reconnection can help to localize magnetic reconnection events in space, astrophysical and fusion plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear evolution of the magnetized Kelvin-Helmholtz instability: from fluid to kinetic modeling

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spectroscopic indication of suprathermal ions in the solar corona

Using spectroscopic data, we support the possibility of suprathermal distributions of coronal ions by fitting the equivalent kappa functions to their emission line profiles. We fit different kappa and Gaussian model functions to line profiles of the strong Fe xv line at 284.16 Angstroms, across two large-field spectroscopic rasters taken in a solar active region. Both single- and double-component Gaussian models are applied, as well as two kappa models, one with a free width parameter allowing for and the other with a constrained width that precludes "microturbulence". We then compare the goodness of fit of the computed best fits for each model. The kappa distribution is a generalization, or superset, of the Maxwellian, so they are able to fit line profiles more precisely than a Gaussian. In most of the data, the best-fit kappa model produces much lower residuals across the profile than any single Gaussian and sometimes double Gaussian. Most importantly, the distribution of estimated kappa values is found to lie mostly in the low-kappa range, implying ion populations far from Maxwellian. Even when the width is removed as a free parameter of fit, the kappa model is still able to fit the data credibly, again with low best-fit values of kappa. We find the shape of the Fe xv line, in the vast majority of the data analyzed, to be indicative of a highly suprathermal ion population.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Multi Level Multi Domain Method for Particle In Cell Plasma Simulations

A novel adaptive technique for electromagnetic Particle In Cell (PIC) plasma simulations is presented here. Two main issues are identified in designing adaptive techniques for PIC simulation: first, the choice of the size of the particle shape function in progressively refined grids, with the need to avoid the exertion of self-forces on particles, and, second, the necessity to comply with the strict stability constraints of the explicit PIC algorithm. The adaptive implementation presented responds to these demands with the introduction of a Multi Level Multi Domain (MLMD) system (where a cloud of self-similar domains is fully simulated with both fields and particles) and the use of an Implicit Moment PIC method as baseline algorithm for the adaptive evolution. Information is exchanged between the levels with the projection of the field information from the refined to the coarser levels and the interpolation of the boundary conditions for the refined levels from the coarser level fields. Particles are bound to their level of origin and are prevented from transitioning to coarser levels, but are repopulated at the refined grid boundaries with a splitting technique. The presented algorithm is tested against a series of simulation challenges.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bipolar Electric Field Signatures of Reconnection Separatrices for a Hydrogen Plasma at Realistic Guide Fields

In preparation for the MMS mission we ask the question: how common are bipolar signatures linked to the presence of electron holes along separatrices emanating from reconnection regions? To answer this question, we conduct massively parallel simulations for realistic conditions and for the hydrogen mass ratio in boxes larger than considered in similar previous studies. The magnetic field configuration includes both a field reversal and a out of plane guide field, as typical of many space situations. The guide field is varied in strength from low values (typical of the Earth magnetotail) to high values comparable to the in plane reconnecting field (as in the magnetopause). In all cases, along the separatrices a strong electron flow is observed, sufficient to lead to the onset of streaming instabilities and to form bipolar parallel electric field signatures. The presence of bipolar structures at all guide fields allows the control of the MMS mission to consider the presence of bipolar signatures as a general flag of the presence of a nearby reconnection site both in the nightside and in the dayside of the magnetosphere.

0 citations0 reviewsNo code linkedOpen access