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.

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.

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.

This White Paper outlines the importance of addressing the fundamental science theme <<How are charged particles energized in space plasmas>> through a future ESA mission. The White Paper presents five compelling science questions related to particle energization by shocks, reconnection,waves and turbulence, jets and their combinations. Answering these questions requires resolving scale coupling, nonlinearity and nonstationarity, which cannot be done with existing multi-point observations. In situ measurements from a multi-point, multi-scale L-class plasma observatory consisting of at least 7 spacecraft covering fluid, ion and electron scales are needed. The plasma observatory will enable a paradigm shift in our comprehension of particle energization and space plasma physics in general, with very important impact on solar and astrophysical plasmas. It will be the next logical step following Cluster, THEMIS and MMS for the very large and active European space plasmas community. Being one of the cornerstone missions of the future ESA Voyage 2035-2050 science program, it would further strengthen the European scientific and technical leadership in this important field.

In collision-poor plasmas from space, e.g., solar wind or stellar outflows, the heat-flux carried by the strahl or beaming electrons is expected to be regulated by the self-generated instabilities. Recently, simultaneous field and particle observations have indeed revealed enhanced whistler-like fluctuations in the presence of counter-beaming populations of electrons, connecting these fluctuations to the whistler heat-flux instability (WHFI). This instability is predicted only for limited conditions of electron beam-plasmas, and was not captured in numerical simulations yet. In this letter we report the first simulations of WHFI in particle-in-cell (PIC) setups, realistic for the solar wind conditions, and without temperature gradients or anisotropies to trigger the instability in the initiation phase. The velocity distributions have a complex reaction to the enhanced whistler fluctuations conditioning the instability saturation by a decrease of the relative drifts combined with induced (effective) temperature anisotropies (heating the core electrons and pitch-angle and energy scattering the strahl). These results are in good agreement with a recent quasilinear approach, and support therefore a largely accepted belief that WHFI saturates at moderate amplitudes. In anti-sunward direction the strahl becomes skewed with a pitch-angle distribution decreasing in width as electron energy increases, that seems to be characteristic to self-generated whistlers and not to small-scale turbulence.

The statistical properties of ions in two-dimensional fully developed turbulence have been compared between two different numerical algorithms. In particular, we compare Hybrid Particle In Cell (hybrid PIC with fluid electrons) and full PIC simulations, focusing on particle diffusion and acceleration phenomena. To investigate several heliospheric plasma conditions, a series of numerical simulations has been performed by varying the plasma $\beta$ - the ratio between kinetic and magnetic pressure. These numerical studies allow the exploration of different scenarios, going from the solar corona (low $\beta$) to the solar wind ($\beta \sim 1$), as well as the Earth's magnetosheath (high $\beta$). It has been found that the two approaches compare pretty well, especially for the spectral properties of the magnetic field and the ion diffusion statistics. Small differences among the models have been found regarding the electric field behaviour at sub-ion scales and the acceleration statistics, due evidently to the more consistent treatment of the plasma in the full PIC approach.

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 new diagnostics, the topographical fluctuations analysis (TFA) to study the spectrum of fluctuations in the various regions around a reconnection site. We find that fluctuations belong to two very different regimes. The first regime is better known, it develops in the reconnection outflows 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, new, regime: it develops in the inflow and in the region around the separatrix surfaces, including the reconnection diffusion region itself. In this new regime the plasma remains laminar but the electromagnetic fields fluctuates strongly. We present an analogy with the smooth continuous motion of the bow of a violin producing the vibrations of the strings to emit music.

We compare relativistic particle integrators commonly used in plasma physics showing several test cases relevant for astrophysics. Three explicit particle pushers are considered, namely the Boris, Vay, and Higuera-Cary schemes. We also present a new relativistic fully implicit particle integrator that is energy conserving. Furthermore, a method based on the relativistic guiding center approximation is included. The algorithms are described such that they can be readily implemented in magnetohydrodynamics codes or Particle-in-Cell codes. Our comparison focuses on the strengths and key features of the particle integrators. We test the conservation of invariants of motion, and the accuracy of particle drift dynamics in highly relativistic, mildly relativistic, and non-relativistic settings. The methods are compared in idealized test cases, i.e., without considering feedback on the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields, $\mathbf{E}\times\mathbf{B}$-fields, force-free fields, and setups relevant for high-energy astrophysics, e.g., a magnetic mirror, a magnetic dipole, and a magnetic null. These tests have direct relevance for particle acceleration in shocks and in magnetic reconnection.

In the Earth magnetosphere, sunspots and magnetic cusp fusion devices, the boundary between the plasma and the magnetic field is marked by a diamagnetic current layer with a rapid change in plasma pressure and magnetic field strength. First principles numerical simulations were conducted to investigate this boundary layer with a spatial resolution beyond electron gyroradius while incorporating a global equilibrium structure. The boundary layer thickness is discovered to be on the order of electron gyroradius scale due to a self-consistent electric field suppressing ion gyromotion at the boundary. Formed at the scale of the electron gyroradius, the electric field plays a critical role in determining equilibrium structure and plasma transport. The discovery highlights the necessity to incorporate electron gyroradius scale physics in studies aimed at advancing our understanding of fusion devices, the magnetosphere and sunspots.

The properties of the turbulence which develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the the dissipation mechanisms in the turbulent outflows similar to the one observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.

The collision of magnetic reconnection jets is studied by means of a three dimensional numerical simulation at kinetic scale, in the presence of a strong guide field. We show that turbulence develops due to the jets collision producing several current sheets in reconnection outflows, aligned with the guide field direction. The turbulence is mainly two-dimensional, with stronger gradients in the plane perpendicular to the guide field and a low wave-like activity in the parallel direction. First, we provide a numerical method to isolate the central turbulent region. Second, we analyze spatial second-order structure function and prove that turbulence is confined in this region. Finally, we compute local magnetic and electric frequency spectra, finding a trend in the sub-ion range that differs from typical cases for which the Taylor hypothesis is valid, as well as wave activity in the range between ion and electron cyclotron frequencies. Our results are relevant to understand observations of reconnection jets collisions in space plasmas.

We present the first detection of solar flare emission at middle-ultraviolet wavelengths around 2000 A by the channel 2 of the Large-Yield RAdiometer (LYRA) onboard the PROBA2 mission. The flare (SOL20170906) was also observed in the channel 1 of LYRA centered at the H I Lyman-{\alpha} line at 1216 A, showing a clear non-thermal profile in both channels. The flare radiation in channel 2 is consistent with the hydrogen Balmer continuum emission produced by an optically thin chromospheric slab heated up to 10000 K. Simultaneous observations in channels 1 and 2 allow the separation of the line emission (primarily from the Lyman-{\alpha} line) from the Balmer continuum emission. Together with the recent detection of the Balmer continuum emission in the near-ultraviolet by IRIS, the LYRA observations strengthen the interpretation of broadband flare emission as the hydrogen recombination continua originating in the chromosphere.

Reconnection outflows are regions of intense recent scrutiny, from in situ observations and from simulations. These regions are host to a variety of instabilities and intense energy exchanges, often even superior to the main reconnection site. We report here a number of results drawn from investigation of simulations. First, the outflows are observed to become unstable to drift instabilities. Second, these instabilities lead to the formation of secondary reconnection sites. Third, the secondary processes are responsible for large energy exchanges and particle energization. Finally, the particle distribution function are modified to become non-Maxwellian and include multiple interpenetrating populations.

MMS observations recently confirmed that crescent-shaped electron velocity distributions in the plane perpendicular to the magnetic field occur in the electron diffusion region near reconnection sites at Earth's magnetopause. In this paper, we re-examine the origin of the crescent-shaped distributions in the light of our new finding that ions and electrons are drifting in opposite directions when displayed in magnetopause boundary-normal coordinates. Therefore, ExB drifts cannot cause the crescent shapes. We performed a high-resolution multi-scale simulation capturing sub-electron skin depth scales. The results suggest that the crescent-shaped distributions are caused by meandering orbits without necessarily requiring any additional processes found at the magnetopause such as the highly asymmetric magnetopause ambipolar electric field. We use an adiabatic Hamiltonian model of particle motion to confirm that conservation of canonical momentum in the presence of magnetic field gradients causes the formation of crescent shapes without invoking asymmetries or the presence of an ExB drift. An important consequence of this finding is that we expect crescent-shaped distributions also to be observed in the magnetotail, a prediction that MMS will soon be able to test.

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.

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.

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.

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.

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.

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.

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.

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.

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).

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. [...]

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.

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.

Reconnection is one of the key processes in astrophysical and laboratory plasmas: it is the opposite of a dynamo. Looking at energy, a dynamo transforms kinetic energy in magnetic energy while reconnection takes magnetic energy and returns is to its kinetic form. Most plasma processes at their core involve first storing magnetic energy accumulated over time and then releasing it suddenly. We focus here on this release. A key concept in analysing reconnection is that of the separatrix, a surface (line in 2D) that separates the fresh unperturbed plasma embedded in magnetic field lines not yet reconnected with the hotter exhaust embedded in reconnected field lines. In kinetic physics, the separatrices become a layer where many key processes develop. We present here new results relative to the processes at the separatrices that regulate the plasma flow, the energisation of the species, the electromagnetic fields and the instabilities developing at the separatrices.

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.

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$.

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.

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.

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.

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.

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.

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

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.

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.

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.

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 {\beta}. 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 {\beta} 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.

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).

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).

Based on a recent theory (Coppi, Nuclear Fusion, 42, 1, 2002) of spontaneous toroidal rotation in tokamaks (Lee et al, Phys Rev Lett, 91, 205003, 2003) and in astrophysical accretion disks, we propose that an analogous process could be at play also in the Earth space environment. We use fully kinetic PIC simulations to study the evolution of drift instabilities and we show that indeed a macroscopic velocity shear is generated spontaneously in the plasma. As in tokamaks, the microscopic fluctuations remain limited to the edge of the plasma channel but the momentum spreads over the whole macroscopic system.

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.

The results of kinetic simulations of magnetic reconnection in Harris current sheets are analyzed. A range of guide fields is considered to study reconnection in plasmas characterized by different $\beta$ values, $\beta > m_e/m_i$. Both an implicit Particle-in-Cell (PIC) simulation method and a parallel explicit PIC code are used. Simulations with mass ratios up to the physical value are performed. The simulations show that the reconnection rate decreases with the guide field and depends weakly on the mass ratio. The off-diagonal components of the electron pressure tensor break the frozen-in condition, even in low $\beta$ plasmas. In high $\beta$ plasmas, evidence is presented that whistler waves play a key role in the enhanced reconnection, while in low $\beta$ plasmas the kinetic Alfv\'en waves are important. The in-plane and the out-of-plane ion and electron motion are also considered, showing that they are influenced by the mass ratio and the plasma $\beta$.

We discuss electron acceleration and heating during collisionless magnetic reconnection by using the results of implicit kinetic simulations of Harris current sheets. We consider and compare electron dynamics in plasmas with different \beta values and perform simulations up to the physical mass ratio. We analyze the typical trajectory of electrons passing through the reconnection region, we study the electron velocity, focusing on the out-of-plane velocity, and we discuss the electron heating along the in-plane and out-of-plane directions.

A new class of soliton-like solutions is derived for the Grad-Shafranov (GS) equations. A mathematical analogy between the GS equation for MHD equilibria and the cubic Schr\"odinger (CS) equation for non-linear wave propagation forms the basis to derive the new class of solutions. The soliton-like solutions are considered for their possible relevance to astrophysics and solar physics problems. We discuss how a soliton-like solution can be generated by a repetitive process of magnetic arcade stretching and plasmoid formation induced by the differential rotation of the solar photosphere or of an accretion disk.