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

Giovanni Lapenta

Giovanni 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-phMachine Learningastro-phastro-ph.GAastro-ph.HEBiological PhysicsComputational Engineering, Finance, and Scienceeess.SPphysics.data-anphysics.ed-ph

Catalog footprint

What is connected

40works
14topics
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

40 published item(s)

preprint2022arXiv

Anisotropic electron heating in turbulence-driven magnetic reconnection in the near-Sun solar wind

We perform a high-resolution two-dimensional fully-kinetic numerical simulation of a turbulent plasma system with observation-driven conditions, in order to investigate the interplay between turbulence, magnetic reconnection, and particle heating from ion to sub-electron scales in the near-Sun solar wind. We find that the power spectra of the turbulent plasma and electromagnetic fluctuations show multiple power-law intervals down to scales smaller than the electron gyroradius. Magnetic reconnection is observed to occur in correspondence of current sheets with a thickness of the order of the electron inertial length, which form and shrink due to interacting ion-scale vortexes. In some cases, both ion and electron outflows are observed (the classic reconnection scenario), while in others -- typically for the shortest current sheets -- only electron jets are presents ("electron-only reconnection"). At the onset of reconnection, the electron temperature starts to increase and a strong parallel temperature anisotropy develops. This suggests that in strong turbulence electron-scale coherent structures may play a significant role for electron heating, as impulsive and localized phenomena such as magnetic reconnection may transfer energy from the electromagnetic fields to particles more efficiently than damping mechanisms related to interactions with wave-like fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Do we need to consider electron kinetic effects to properly model a planetary magnetosphere: the case of Mercury

The magnetosphere of Mercury is studied using an implicit full particle in cell simulation (PIC). We use a hybrid simulation where ions are full particles and electrons are considered as a fluid to start a full PIC simulation where electrons are also particles and follow their distribution function. This approach allows us to estimate the changes introduced by the electron kinetic physics. We find that the overall macroscopic state of the magnetosphere of Mercury is little affected but several physical processes are significantly modified in the full PIC simulation: the foreshock region is more active with more intense shock reformation, the Kelvin-Helmholtz rippling effects on the nightside magnetopause are sharper, and the magnetotail current sheet becomes thinner than those predicted by the hybrid simulation. The greatest effect of the electron physics, comes from the processes of particle energization. Both species, not just the electrons, are found to gain more energy when kinetic electron processes are taken into account. The region with the most energetic plasma is found on the dusk side of the tail where magnetic flux ropes are formed due to reconnection. We find that the ion and electron energization is associated with the regions of reconnection and the development of kinetic instabilities caused by counter-streaming electron populations. The resulting electron distributions are highly non Maxwellian, a process that neither MHD nor hybrid models can describe.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Identification of high order closure terms from fully kinetic simulations using machine learning

Simulations of large-scale plasma systems are typically based on a fluid approximation approach. These models construct a moment-based system of equations that approximate the particle-based physics as a fluid, but as a result lack the small-scale physical processes available to fully kinetic models. Traditionally, empirical closure relations are used to close the moment-based system of equations, which typically approximate the pressure tensor or heat flux. The more accurate the closure relation, the stronger the simulation approaches kinetic-based results. In this paper, new closure terms are constructed using machine learning techniques. Two different machine learning models, a multi-layer perceptron and a gradient boosting regressor, synthesize a local closure relation for the pressure tensor and heat flux vector from fully kinetic simulations of a 2D magnetic reconnection problem. The models are compared to an existing closure relation for the pressure tensor, and the applicability of the models is discussed. The initial results show that the models can capture the diagonal components of the pressure tensor accurately, and show promising results for the heat flux, opening the way for new experiments in multi-scale modeling. We find that the sampling of the points used to train both models play a capital role in their accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Kinetic heating by Alfvén waves in magnetic shears

With first-principles kinetic simulations, we show that a large-scale Alfvén wave (AW) propagating in an inhomogeneous background decays into kinetic Alfvén waves (KAWs), triggering ion and electron energization. We demonstrate that the two species can access unequal amounts of the initial AW energy, experiencing differential heating. During the decay process, the electric field carried by KAWs produces non-Maxwellian features in the particle VDFs, in accordance with space observations. The process we present solely requires the interaction of a large-scale AW with a magnetic shear and may be relevant for several astrophysical and laboratory plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Automatic classification of plasma regions in near-Earth space with supervised machine learning: application to Magnetospheric Multi Scale 2016-2019 observations

The proper classification of plasma regions in near-Earth space is crucial to perform unambiguous statistical studies of fundamental plasma processes such as shocks, magnetic reconnection, waves and turbulence, jets and their combinations. The majority of available studies have been performed by using human-driven methods, such as visual data selection or the application of predefined thresholds to different observable plasma quantities. While human-driven methods have allowed performing many statistical studies, these methods are often time-consuming and can introduce important biases. On the other hand, the recent availability of large, high-quality spacecraft databases, together with major advances in machine-learning algorithms, can now allow meaningful applications of machine learning to in-situ plasma data. In this study, we apply the fully convolutional neural network (FCN) deep machine-leaning algorithm to the recent Magnetospheric Multi Scale (MMS) mission data in order to classify ten key plasma regions in near-Earth space for the period 2016-2019. For this purpose, we use available intervals of time series for each such plasma region, which were labeled by using human-driven selective downlink applied to MMS burst data. We discuss several quantitative parameters to assess the accuracy of both methods. Our results indicate that the FCN method is reliable to accurately classify labeled time series data since it takes into account the dynamical features of the plasma data in each region. We also present good accuracy of the FCN method when applied to unlabeled MMS data. Finally, we show how this method used on MMS data can be extended to data from the Cluster mission, indicating that such method can be successfully applied to any in situ spacecraft plasma database.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Domain of Influence analysis: implications for Data Assimilation in space weather forecasting

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamic Time Warping as a New Evaluation for Dst Forecast with Machine Learning

Models based on neural networks and machine learning are seeing a rise in popularity in space physics. In particular, the forecasting of geomagnetic indices with neural network models is becoming a popular field of study. These models are evaluated with metrics such as the root-mean-square error (RMSE) and Pearson correlation coefficient. However, these classical metrics sometimes fail to capture crucial behavior. To show where the classical metrics are lacking, we trained a neural network, using a long short-term memory network, to make a forecast of the disturbance storm time index at origin time $t$ with a forecasting horizon of 1 up to 6 hours, trained on OMNIWeb data. Inspection of the model's results with the correlation coefficient and RMSE indicated a performance comparable to the latest publications. However, visual inspection showed that the predictions made by the neural network were behaving similarly to the persistence model. In this work, a new method is proposed to measure whether two time series are shifted in time with respect to each other, such as the persistence model output versus the observation. The new measure, based on Dynamical Time Warping, is capable of identifying results made by the persistence model and shows promising results in confirming the visual observations of the neural network's output. Finally, different methodologies for training the neural network are explored in order to remove the persistence behavior from the results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tokamak disruption prediction using different machine learning techniques

Disruption prediction and mitigation is of key importance in the development of sustainable tokamakreactors. Machine learning has become a key tool in this endeavour. In this paper multiple machinelearning models will be tested and compared. A particular focus has been placed on their portability.This describes how easily the models can be used with data from new devices. The methods used inthis paper are support vector machine, 2-tiered support vector machine, random forest, gradient boostedtrees and long-short term memory. The results show that the support vector machine performanceis marginally better among the standard models, while the gradient boosted trees performed the worst.The portable variant of each model had lower performance. Random forest obtained the highest portableperformance. Results also suggest that disruptions can be detected as early as 600ms before the event.An analysis of the computational cost showed all models run in less than 1ms, allowing sufficient timefor disruption mitigation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Visualizing and Interpreting Unsupervised Solar Wind Classifications

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Characterizing magnetic reconnection regions using Gaussian mixture models on particle velocity distributions

We present a method based on unsupervised machine learning to identify regions of interest using particle velocity distributions as a signature pattern. An automatic density estimation technique is applied to particle distributions provided by PIC simulations to study magnetic reconnection. The key components of the method involve: i) a Gaussian mixture model determining the presence of a given number of subpopulations within an overall population, and ii) a model selection technique with Bayesian Information Criterion to estimate the appropriate number of subpopulations. Thus, this method identifies automatically the presence of complex distributions, such as beams or other non-Maxwellian features, and can be used as a detection algorithm able to identify reconnection regions. The approach is demonstrated for specific double Harris sheet simulations but it can in principle be applied to any other type of simulation and observational data on the particle distribution function.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Energy exchanges in reconnection outflows

Reconnection outflows are highly energetic directed flows that interact with the ambient plasma or with flows from other reconnection regions. Under these conditions the flow becomes highly unstable and chaotic, as any flow jets interacting with a medium. We report here massively parallel simulations of the two cases of interaction between outflow jets and between a single outflow with an ambient plasma. We find in both case the development of a chaotic magnetic field, subject to secondary reconnection events that further complicate the topology of the field lines. The focus of the present analysis is on the energy balance. We compute each energy channel (electromagnetic, bulk, thermal, for each species) and find where the most energy is exchanged and in what form. The main finding is that the largest energy exchange is not at the reconnection site proper but in the regions where the outflowing jets are destabilized.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the ions acceleration via collisionless magnetic reconnection in laboratory plasmas

This work presents an analysis of the ion outflow from magnetic reconnection throughout fully kinetic simulations with typical laboratory plasmas values. A symmetric initial configuration for the density and magnetic field is considered across the current sheet. After analyzing the behavior of a set of nine simulations with a reduced mass ratio and with a permuted value of three initial electron temperature and magnetic field intensity, the best ion acceleration scenario is further studied with a realistic mass ratio in terms of the ion dynamics and energy budget. Interestingly, a series of shock waves structures are observed in the outflow, resembling the shock discontinuities found in recent magnetohydrodynamic (MHD) simulations. An analysis of the ion outflow at several distances from the reconnection point is presented, in light of possible laboratory applications. The analysis suggests that magnetic reconnection could be used as a tool for plasma acceleration, with applications ranging from electric propulsion to production of ion thermal beams.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Energetics of kinetic reconnection in a three-dimensional null points cluster

We performed three-dimensional Particle-in-Cell simulations of magnetic reconnection with multiple magnetic null points. Magnetic field energy conversion into kinetic energy was about five times higher than in traditional Harris sheet configuration. More than 85% of initial magnetic field energy was transferred to particle energy during 25 reversed ion cyclofrequencies. Magnetic reconnection in the cluster of null points evolved in three phases. During the first phase, ion beams were excited, that then gave part of their energy back to magnetic field in the second phase. In the third phase, magnetic reconnection occurs in many small patches around the current channels formed along the stripes of low magnetic field. Magnetic reconnection in null points presents essentially three-dimensional features, with no two dimensional symmetries or current sheets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Energy dissipation in magnetic null points at kinetic scales

We use kinetic particle-in-cell and magnetohydrodynamic simulations supported by an observational dataset to investigate magnetic reconnection in clusters of null points in space plasma. The magnetic configuration under investigation is driven by fast adiabatic flux rope compression that dissipates almost half of the initial magnetic field energy. In this phase powerful currents are excited producing secondary instabilities, and the system is brought into a state of `intermittent turbulence' within a few ion gyro-periods. Reconnection events are distributed all over the simulation domain and energy dissipation is rather volume-filling. Numerous spiral null points interconnected via their spines form null lines embedded into magnetic flux ropes; null point pairs demonstrate the signatures of torsional spine reconnection. However, energy dissipation mainly happens in the shear layers formed by adjacent flux ropes with oppositely directed currents. In these regions radial null pairs are spontaneously emerging and vanishing, associated with electron streams and small-scale current sheets. The number of spiral nulls in the simulation outweighs the number of radial nulls by a factor of 5\---10, in accordance with Cluster observations in the Earth's magnetosheath. Twisted magnetic fields with embedded spiral null points might indicate the regions of major energy dissipation for future space missions such as Magnetospheric Multiscale Mission (MMS).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetic null points in kinetic simulations of space plasmas

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the electron dynamics during island coalescence in asymmetric magnetic reconnection

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Role of Z-pinches in magnetic reconnection in space plasmas

A widely accepted scenario of magnetic reconnection in collisionless space plasmas is the breakage of magnetic field lines in X-points. In laboratory, reconnection is commonly studied in pinches, current channels embedded into twisted magnetic fields. No model of magnetic reconnection in space plasmas considers both null-points and pinches as peers. We have performed a particle-in-cell simulation of magnetic reconnection in a three-dimensional configuration where null-points are present initially, and Z-pinches are formed during the simulation along the lines of spiral null-points. The non-spiral null-points are more stable than spiral ones, and no substantial energy dissipation is associated with them. On the contrary, turbulent magnetic reconnection in the pinches causes the magnetic energy to decay at a rate of ~1.5% per ion gyro period. Dissipation in similar structures is a likely scenario in space plasmas with large fraction of spiral null-points.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Understanding space weather to shield society: A global road map for 2015-2025 commissioned by COSPAR and ILWS

There is a growing appreciation that the environmental conditions that we call space weather impact the technological infrastructure that powers the coupled economies around the world. With that comes the need to better shield society against space weather by improving forecasts, environmental specifications, and infrastructure design. [...] advanced understanding of space weather requires a coordinated international approach to effectively provide awareness of the processes within the Sun-Earth system through observation-driven models. This roadmap prioritizes the scientific focus areas and research infrastructure that are needed to significantly advance our understanding of space weather of all intensities and of its implications for society. Advancement of the existing system observatory through the addition of small to moderate state-of-the-art capabilities designed to fill observational gaps will enable significant advances. Such a strategy requires urgent action: key instrumentation needs to be sustained, and action needs to be taken before core capabilities are lost in the aging ensemble. We recommend advances through priority focus (1) on observation-based modeling throughout the Sun-Earth system, (2) on forecasts more than 12 hrs ahead of the magnetic structure of incoming coronal mass ejections, (3) on understanding the geospace response to variable solar-wind stresses that lead to intense geomagnetically-induced currents and ionospheric and radiation storms, and (4) on developing a comprehensive specification of space climate, including the characterization of extreme space storms to guide resilient and robust engineering of technological infrastructures. The roadmap clusters its implementation recommendations by formulating three action pathways, and outlines needed instrumentation and research programs and infrastructure for each of these. [...]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Web-based description of the space radiation environment using the Bethe-Bloch model

Space weather is a rapidly growing area not only in scientific and engineering applications but also in physics education and in the interest of the public. We focus especially on space radiation and its impact on space exploration. The topic is highly interdisciplinary bringing together fundamental concepts of nuclear physics with aspects of radiation protection and space science. We present a new approach to presenting the topic by developing a web-based tool that combines some of the fundamental concepts from these two fields in a single tool that can be developed in the context of advanced secondary or undergraduate university education. We present DREADCode, an outreach or teaching tool to asses rapidly the current conditions of the radiation field in space. DREADCode uses the available data feeds from a number of ongoing space missions to produce a first order approximation of the dose an astronaut would receive during a mission of exploration in deep space. DREADcode is based on a intuitive GUI interface available online from the European Space Weather Portal. The core of the radiation transport computation to produce the radiation dose from the observed fluence of radiation observed by the spacecraft fleet considered is based on a relatively simple approximation: the Bethe-Block equation. DREADCode assumes also a simplified geometry and material configuration for the shields used to compute the dose. The approach is approximate and it sacrifices some important physics on the altar of a rapid execution time allowing a real time operation scenario. There is no intention here to produce an operational tool for use in the space science and engineering. Rather we present an educational tool at undergraduate level that uses modern web-based and programming methods to learn some of the most important concepts in the application of radiation protection to space weather problems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electromagnetic Energy Conversion in Downstream Fronts from 3D Kinetic Reconnection

The electromagnetic energy equation is analyzed term by term in a 3D simulation of kinetic reconnection previously reported by \citet{vapirev2013formation}. The evolution presents the usual 2D-like topological structures caused by an initial perturbation independent of the third dimension. However, downstream of the reconnection site, where the jetting plasma encounters the yet unperturbed pre-existing plasma, a downstream front (DF) is formed and made unstable by the strong density gradient and the unfavorable local acceleration field. The energy exchange between plasma and fields is most intense at the instability, reaching several $pW/m^3$, alternating between load (energy going from fields to particles) and generator (energy going from particles to fields) regions. Energy exchange is instead purely that of a load at the reconnection site itself in a region focused around the x-line and elongated along the separatrix surfaces. Poynting fluxes are generated at all energy exchange regions and travel away from the reconnection site transporting an energy signal of the order of about $\mathbf S \approx 10^{-3} W/m^2$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High Energy Electron Confinement in a Magnetic Cusp Configuration

We report experimental results validating the concept that plasma confinement is enhanced in a magnetic cusp configuration when beta (plasma pressure/magnetic field pressure) is order of unity. This enhancement is required for a fusion power reactor based on cusp confinement to be feasible. The magnetic cusp configuration possesses a critical advantage: the plasma is stable to large scale perturbations. However, early work indicated that plasma loss rates in a reactor based on a cusp configuration were too large for net power production. Grad and others theorized that at high beta a sharp boundary would form between the plasma and the magnetic field, leading to substantially smaller loss rates. The current experiment validates this theoretical conjecture for the first time and represents critical progress toward the Polywell fusion concept which combines a high beta cusp configuration with an electrostatic fusion for a compact, economical, power-producing nuclear fusion reactor.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kinetic plasma turbulence during the nonlinear stage of the Kelvin-Helmholtz instability

Using a full kinetic, implicit particle-in-cell code, iPiC3D, we studied the properties of plasma kinetic turbulence, such as would be found at the interface between the solar wind and the Earth magnetosphere at low latitude during northwards periods. In this case, in the presence of a magnetic field B oriented mostly perpendicular to the velocity shear, turbulence is fed by the disruption of a Kelvin-Helmholtz vortex chain via secondary instabilities, vortex pairing and non-linear interactions. We found that the magnetic energy spectral cascade between ion and electron inertial scales, $d_i$ and $d_e$, is in agreement with satellite observations and other previous numerical simulations; however, in our case the spectrum ends with a peak beyond $d_e$ due to the occurrence of the lower hybrid drift instability. The electric energy spectrum is influenced by effects of secondary instabilities: anomalous resistivity, fed by the development of the lower hybrid drift instability, steepens the spectral decay and, depending on the alignment or anti-alignment of B and the shear vorticity, peaks due to ion-Bernstein waves may dominate the spectrum around $d_i$. These waves are generated by counter-streaming flow structures, through flux freezing also responsible for reconnection of the in-plane component of the magnetic field, which then generates electron pressure anisotropy and flattening of the field-aligned component of the electron distribution function.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kinetic Simulations of Plasmoid Chain Dynamics

The dynamics of a plasmoid chain is studied with three dimensional Particle-in-Cell simulations. The evolution of the system with and without a uniform guide field, whose strength is 1/3 the asymptotic magnetic field, is investigated. The plasmoid chain forms by spontaneous magnetic reconnection: the tearing instability rapidly disrupts the initial current sheet generating several small-scale plasmoids, that rapidly grow in size coalescing and kinking. The plasmoid kink is mainly driven by the coalescence process. It is found that the presence of guide field strongly influences the evolution of the plasmoid chain. Without a guide field, a main reconnection site dominates and smaller reconnection regions are included in larger ones, leading to an hierarchical structure of the plasmoid-dominated current sheet. On the contrary in presence of a guide field, plasmoids have approximately the same size and the hierarchical structure does not emerge, a strong core magnetic field develops in the center of the plasmoid in the direction of the existing guide field, and bump-on-tail instability, leading to the formation of electron holes, is detected in proximity of the plasmoids.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collisionless magnetic reconnection in a plasmoid chain

The kinetic features of plasmoid chain formation and evolution are investigated by two dimensional Particle-in-Cell simulations. Magnetic reconnection is initiated in multiple X points by the tearing instability. Plasmoids form and grow in size by continuously coalescing. Each chain plasmoid exhibits a strong out-of plane core magnetic field and an out-of-plane electron current that drives the coalescing process. The disappearance of the X points in the coalescence process are due to anti-reconnection, a magnetic reconnection where the plasma inflow and outflow are reversed with respect to the original reconnection flow pattern. Anti-reconnection is characterized by the Hall magnetic field quadrupole signature. Two new kinetic features, not reported by previous studies of plasmoid chain evolution, are here revealed. First, intense electric fields develop in-plane normally to the separatrices and drive the ion dynamics in the plasmoids. Second, several bipolar electric field structures are localized in proximity of the plasmoid chain. The analysis of the electron distribution function and phase space reveals the presence of counter-streaming electron beams, unstable to the two stream instability, and phase space electron holes along the reconnection separatrices.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rethinking Electrostatic Solvers in Particle Simulations for the Exascale Era

In preparation to the exascale era, an alternative approach to calculate the electrostatic forces in Particle Mesh (PM) methods is proposed. While the traditional techniques are based on the calculation of the electrostatic potential by solving the Poisson equation, in the new approach the electric field is calculated by solving the Ampere's law. When the Ampere's law is discretized explicitly in time, the electric field values on the mesh are simply updated from the previous values. In this way, the electrostatic solver becomes an embarrassingly parallel problem, making the algorithm extremely scalable and suitable for exascale computing platforms. An implementation of a one dimensional PM code is presented to show that the proposed method produces correct results, and it is a very promising algorithm for exascale PM simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Three dimensional density cavities in guide field collisionless magnetic reconnection

Particle-in-Cell simulations of collisionless magnetic reconnection with a guide field reveal for the first time the three dimensional features of the low density regions along the magnetic reconnection separatrices, the so-called "cavities". It is found that structures with further lower density develop within the cavities. Because their appearance is similar to the rib shape, these formations are here called "low density ribs". Their location remains approximately fixed in time and their density progressively decreases, as electron currents along the cavities evacuate them. They develop along the magnetic field lines and are supported by a strong perpendicular electric field that oscillates in space. In addition, bipolar parallel electric field structures form as isolated spheres between the cavities and the outflow plasma, along the direction of the low density ribs and of magnetic field lines.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Achieving Fast Reconnection in Resistive MHD Models via Turbulent Means

Astrophysical fluids are generally turbulent and this preexisting turbulence must be taken into account for the models of magnetic reconnection which are attepmted to be applied to astrophysical, solar or heliospheric environments. In addition, reconnection itself induces turbulence which provides an important feedback on the reconnection process. In this paper we discuss both theoretical model and numerical evidence that magnetic reconnection gets fast in the approximation of resistive MHD. We consider the relation between the Lazarian & Vishniac turbulent reconnection theory and Lapenta's numerical experiments testifying of the spontaneous onset of turbulent reconnection in systems which are initially laminar.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Development and performance analysis of a UPC Particle-in-Cell code

The development and the implementation of a Particle-in-Cell code written in the Unified Parallel C (UPC) language for plasma simulations with application to astrophysics and fusion nuclear energy machines are presented. A simple one dimensional electrostatic Particle-in-Cell code has been developed first to investigate the implementation details in the UPC language, and second to study the UPC performance on parallel computers. The initial simulations of plasmas with the UPC Particle-in-Cell code and a study of parallel speed-up of the UPC code up to 128 cores are shown.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic simulations of magnetic reconnection in presence of a background O+ population

Particle-in-Cell simulations of magnetic reconnection with an H+ current sheet and a mixed background plasma of H+ and O+ ions are completed using physical mass ratios. Four main results are shown. First, the O+ presence slightly decreases the reconnection rate and the magnetic reconnection evolution depends mainly on the lighter H+ ion species in the presented simulations. Second, the Hall magnetic field is characterized by a two-scale structure in presence of O+ ions: it reaches sharp peak values in a small area in proximity of the neutral line, and then decreases slowly over a large region. Third, the two background species initially separate in the outflow region because H+ and O+ ions are accelerated by different mechanisms occurring on different time scales and with different strengths. Fourth, the effect of a guide field on the O+ dynamics is studied: the O+ presence does not change the reconnected flux and all the characteristic features of guide field magnetic reconnection are still present. Moreover, the guide field introduces an O+ circulation pattern between separatrices that enhances high O+ density areas and depletes low O+ density regions in proximity of the reconnection fronts. The importance and the validity of these results are finally discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measurement of the Electric Current in a Kpc-Scale Jet

We present radio emission, polarization, and Faraday rotation maps of the radio jet of the galaxy 3C303. From this data we derive the magnetoplasma and electrodynamic parameters of this 50 kpc long jet. For a {$\sim 2$} kpc segment of this jet we obtain for the first time a direct determination of a {\it galactic}-scale electric current ($\sim 3 \times 10^{18}$ A), and its direction $-$ {\it positive} away from the AGN. Our analysis strongly supports a model where the jet energy flow is mainly electromagnetic.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self-consistent Seeding of the Interchange Instability in Dipolarization Fronts

We report a 3D magnetohydrodynamics simulation that studies the formation of dipolarization fronts during magnetotail reconnection. The crucial new feature uncovered in the present 3D simulation is that the process of reconnection produces flux ropes developing within the reconnection region. These flux ropes are unstable to the kink mode and introduce a spontaneous structure in the dawn-dusk direction. The dipolarization fronts forming downstream of reconnection are strongly affected by the kinking ropes. At the fronts, a density gradient is present with opposite direction to that of the acceleration field and leads to an interchange instability. We present evidence for a causal link where the perturbations of the kinking flux ropes with their natural and well defined scales drive and select the scales for the interchange mode in the dipolarization fronts. The results of the simulation are validated against measured structures observed by the Themis mission.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Space Weather Prediction with Exascale Computing

Space weather refers to conditions on the Sun, in the interplanetary space and in the Earth space environment that can influence the performance and reliability of space-borne and ground-based technological systems and can endanger human life or health. Adverse conditions in the space environment can cause disruption of satellite operations, communications, navigation, and electric power distribution grids, leading to a variety of socioeconomic losses. The conditions in space are also linked to the Earth climate. The activity of the Sun affects the total amount of heat and light reaching the Earth and the amount of cosmic rays arriving in the atmosphere, a phenomenon linked with the amount of cloud cover and precipitation. Given these great impacts on society, space weather is attracting a growing attention and is the subject of international efforts worldwide. We focus here on the steps necessary for achieving a true physics-based ability to predict the arrival and consequences of major space weather storms. Great disturbances in the space environment are common but their precise arrival and impact on human activities varies greatly. Simulating such a system is a grand- challenge, requiring computing resources at the limit of what is possible not only with current technology but also with the foreseeable future generations of super computers

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spontaneous transition to a fast 3D turbulent reconnection regime

We show how the conversion of magnetic field energy via magnetic reconnection can progress in a fully three-dimensional, fast, volume-filling regime. An initial configuration representative of many laboratory, space and astrophysical plasmas spontaneously evolves from the well-known regime of slow, resistive reconnection to a new regime that allows to explain the rates of energy transfer observed in jets emitted from accretion disks, in stellar/solar flare processes as well as in laboratory plasmas. This process does not require any pre-existing turbulence seed which often is not observed in the host systems prior to the onset of the energy conversion. The dynamics critically depends on the interplay of perturbations developing along the magnetic field lines and across them, a process possible only in three-dimensions. The simulations presented here are the first able to show this transition in a fully three-dimensional configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Energy Conserving Particle-in-Cell Method

A new Particle-in-Cell (PIC) method, that conserves energy exactly, is presented. The particle equations of motion and the Maxwell's equations are differenced implicitly in time by the midpoint rule and solved concurrently by a Jacobian-free Newton Krylov (JFNK) solver. Several tests show that the finite grid instability is eliminated in energy conserving PIC simulations, and the method correctly describes the two-stream and Weibel instabilities, conserving exactly the total energy. The computational time of the energy conserving PIC method increases linearly with the number of particles, and it is rather insensitive to the number of grid points and time step. The kinetic enslavement technique can be effectively used to reduce the problem matrix size and the number of JFNK solver iterations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Investigation of Dynamics of Self-Similarly Evolving Magnetic Clouds

Magnetic clouds (MCs) are "magnetized plasma clouds" moving in the solar wind. MCs transport magnetic flux and helicity away from the Sun. These structures are not stationary but feature temporal evolution. Commonly, simplified MC models are considered. The goal of the present study is to investigate the dynamics of more general, radially expanding MCs. They are considered as cylindrically symmetric magnetic structures with low plasma β. In order to study MC`evolution the self-similar approach method and a numerical approach are used. It is shown that the forces are balanced in the considered self-similarly evolving, cylindrically symmetric magnetic structures. Explicit analytical expressions for magnetic field, plasma velocity, density and pressure within MCs are derived. These solutions are characterized by conserved values of magnetic flux and helicity. We also investigate the dynamics of self-similarly evolving MCs by means of the numerical code "Graale". In addition, their expansion in a medium with higher density and higher plasma β is studied. It is shown that the physical parameters of the MCs maintain their self-similar character throughout their evolution. Conclusions. A comparison of the different self-similar and numerical solutions allows us to conclude that the evolving MCs are quite adequately described by our self-similar solutions - they retain their self-similar, coherent nature for quite a long time and over large distances from the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Jet deflection by very weak guide fields during magnetic reconnection

Simulations of antiparallel reconnection have shown collimated electron jets outflowing from the x-point, and associated highly elongated "outer electron diffusion regions." New PIC simulations with ion/electron mass ratios as large as 1836 show that jets are deflected towards the magnetic separatrix by out-of-plane guide fields, Bg, as small as 0.05 times the asymptotic reconnecting field, B0. The outer electron diffusion region is distorted and broken up, but the diffusion rate is unchanged. These results are interpreted in terms of electron dynamics and are compared to recent measurements of reconnection jets in the magnetosheath.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Analysis of the energy release for different magnetic reconnection regimes within the solar environment

A 2.5-dimensional magnetohydrodynamics simulation analysis of the energy release for three different reconnection regimes is presented. The system under investigation consists in a current-sheet located in a medium with a strong density variation along the current layer: such system is modeled as it were located in the high chromosphere/low solar corona as in the case of pre- flare and coronal mass ejection (CME) configurations or in the aftermath of such explosive phenomena. By triggering different magnetic-reconnection dynamics, that is from a laminar slow evolution to a spontaneous non-steady turbulent reconnection [1,2,3], we observe a rather different efficiency and temporal behavior with regard to the energy fluxes associated with each of these reconnection-driven evolutions. These discrepancies are fundamental key-properties to create realistic models of the triggering mechanisms and initial evolution of all those phenomena requiring fast (and high power) magnetic reconnection events within the solar environment. 1. G. Lapenta, Phys. Rev. Lett. 100, 235001 (2008). 2. L. Bettarini, and G. Lapenta, ApJ Submitted (2009). 3. M. Skender, and G. Lapenta, Phys. Plasmas submitted (2009).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Spontaneous non-steady magnetic reconnection within the solar environment

This work presents a 2.5-dimensional simulation study of the instability of current-sheets located in a medium with a strong density variation along the current layer. The initial force-free configuration is observed to undergo a two-stage evolution consisting of an abrupt regime transition from a slow to a fast reconnection process leading the system to a final chaotic configuration. Yet, the onset of the fast phase is not determined by the presence of any anomalous enhancement in plasma's local resistivity, but rather is the result of a new mechanism discovered in Lapenta (2008)* and captured only when sufficient resolution is present. Hence, the effects of the global resistivity, the global viscosity and the plasma beta on the overall dynamics are considered. This mechanism allowing the transition from slow to fast reconnection provides a simple but effective model of several processes taking place within the solar atmosphere from the high chromosphere up to the low corona. In fact, the understanding of a spontaneous transition to a self-feeding fast reconnection regime as well as its macroscopic evolution is the first and fundamental step to produce realistic models of all those phenomena requiring fast (and high power) triggering events (* Lapenta G. 2008, Phys. Rev. Lett., 100, 235001).

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Blob formation and acceleration in the solar wind: role of converging flows and viscosity

The effect of viscosity and of converging flows on the formation of blobs in the slow solar wind is analysed by means of resistive MHD simulations. The regions above coronal streamers where blobs are formed (Sheeley et al., 1997) are simulated using a model previously proposed by Einaudi et al. (1999). The result of our investigation is twofold. First, we demonstrate a new mechanism for enhanced momentum transfer between a forming blob and the fast solar wind surrounding it. The effect is caused by the longer range of the electric field caused by the tearing instability forming the blob. The electric field reaches into the fast solar wind and interacts with it, causing a viscous drag that is global in nature rather than local across fluid layers as it is the case in normal uncharged fluids (like water). Second, the presence of a magnetic cusp at the tip of a coronal helmet streamer causes a converging of the flows on the two sides of the streamer and a direct push of the forming island by the fast solar wind, resulting in a more efficient momentum exchange.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Plasma sheet structure in the magnetotail: kinetic simulation and comparison with satellite observations

We use the results of a three-dimensional kinetic simulation of an Harris current sheet to propose an explanation and to reproduce the ISEE-1/2, Geotail, and Cluster observations of the magnetotail current sheet structure. Current sheet flapping, current density bifurcation, and reconnection are explained as the results of the kink and tearing instabilities, which dominate the current sheet evolution.

0 citations0 reviewsNo code linkedOpen access