A new look at the structure of the electron diffusion region in collisionless magnetic reconnection is presented. The research is based on a particle-in-cell simulation of asymmetric magnetic reconnection, which include a temperature gradient across the current layer in addition to density and magnetic field gradient. We find that none of X-point, flow stagnation point, and local current density peak coincide. Current and energy balance analyses around the flow stagnation point and current density peak show consistently that current dissipation is associated with the divergence of nongyrotropic electron pressure. Furthermore, the same pressure terms, when combined with shear-type gradients of the electron flow velocity, also serve to maintain local thermal energy against convective losses. These effects are similar to those found also in symmetric magnetic reconnection. In addition, we find here significant effects related to the convection of current, which we can relate to a generalized diamagnetic drift by the nongyrotropic pressure divergence. Therefore, only part of the pressure force serves to dissipate the current density. However, the prior conclusion that the role of the reconnection electric field is to maintain the current density, which was obtained for a symmetric system, applies here as well. Finally, we discuss related features of electron distribution function in the EDR.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

During magnetically dominated relativistic reconnection, inflowing plasma depletes the initial relativistic pressure at the x-line and collisionless plasma heating inside the diffusion region is insufficient to overcome this loss. The resulting pressure drop causes a collapse at the x-line, essentially a localization mechanism of the diffusion region necessary for fast reconnection. The extension of this low-pressure region further explains the bursty nature of anti-parallel reconnection because a once opened outflow exhaust can also collapse, which repeatedly triggers secondary tearing islands. However, a stable single x-line reconnection can be achieved when an external guide field exists, since the reconnecting magnetic field component rotates out of the reconnection plane at outflows, providing additional magnetic pressure to sustain the opened exhausts.

The spreading of the X-line out of the reconnection plane under a strong guide field is investigated using large-scale three-dimensional (3D) particle-in-cell (PIC) simulations in asymmetric magnetic reconnection. A simulation with a thick, ion-scale equilibrium current sheet (CS) reveals that the X-line spreads at the ambient ion/electron drift speeds, significantly slower than the Alfv\'en speed based on the guide field $V_{Ag}$. Additional simulations with a thinner, sub-ion-scale CS show that the X-line spreads at $V_{Ag}$ (Alfv\'enic spreading), much higher than the ambient species drifts. An Alfv\'enic signal consistent with kinetic Alfv\'en waves develops and propagates, leading to CS thinning and extending, which then results in reconnection onset. The continuous onset of reconnection in the signal propagation direction manifests as Alfv\'enic X-line spreading. The strong dependence on the CS thickness of the spreading speeds, and the X-line orientation are consistent with the collisionless tearing instability. Our simulations indicate that when the collisionless tearing growth is sufficiently strong in a thinner CS such that $\gamma/\Omega_{ci}\gtrsim\mathcal{O}(1)$, Alfv\'enic X-line spreading can take place. Our results compare favorably with a number of numerical simulations and recent magnetopause observations. A key implications is that the magnetopause CS is typically too thick for Alfv\'enic X-line spreading to effectively take place.

Using 3D particle-in-cell (PIC) simulations, we study magnetic reconnection with the x-line being spatially confined in the current direction. We include thick current layers to prevent reconnection at two ends of a thin current sheet that has a thickness on an ion inertial (di) scale. The reconnection rate and outflow speed drop significantly when the extent of the thin current sheet in the current direction is < O(10 di). When the thin current sheet extent is long enough, we find it consists of two distinct regions; an inactive region (on the ion-drifting side) exists adjacent to the active region where reconnection proceeds normally as in a 2D case. The extent of this inactive region is ~ O(10 di), and it suppresses reconnection when the thin current sheet extent is comparable or shorter. The time-scale of current sheet thinning toward fast reconnection can be translated into the spatial-scale of this inactive region; because electron drifts inside the ion diffusion region transport the reconnected magnetic flux, that drives outflows and furthers the current sheet thinning, away from this region. This is a consequence of the Hall effect in 3D. While this inactive region may explain the shortest possible azimuthal extent of dipolarizing flux bundles at Earth, it may also explain the dawn-dusk asymmetry observed at the magnetotail of Mercury, that has a global dawn-dusk extent much shorter than that of Earth.

We demonstrate the dragging of the magnetic field by the super-Alfvenic shear flows out of the reconnection plane can strongly localize the reconnection x-line in collisionless plasmas, reversing the current direction at the x-line. Reconnection with this new morphology, which is impossible in resistive-magnetohydrodynamic (MHD), is enabled by electron inertia. Surprisingly, the quasi-steady reconnection rate remains of order 0.1 even though the aspect ratio of the local x-line geometry is larger than unity. We explain this by examining the transport of the reconnected magnetic flux and the opening angle made by the upstream magnetic field, concluding that the reconnection rate is still limited by the constraint imposed at the inflow region. This study further suggests the nearly universal fast rate value of order 0.1 cannot be explained by the physics of tearing modes, nor can it be explained by a universal localization mechanism.

The orientation and stability of the reconnection x-line in asymmetric geometry is studied using three-dimensional (3D) particle-in-cell simulations. We initiate reconnection at the center of a large simulation domain to minimize the boundary effect. The resulting x-line has sufficient freedom to develop along an optimal orientation, and it remains laminar. Companion 2D simulations indicate that this x-line orientation maximizes the reconnection rate. The divergence of the non-gyrotropic pressure tensor breaks the frozen-in condition, consistent with its 2D counterpart. We then design 3D simulations with one dimension being short to fix the x-line orientation, but long enough to allow the growth of the fastest growing oblique tearing modes. This numerical experiment suggests that reconnection tends to radiate secondary oblique tearing modes if it is externally (globally) forced to proceed along an orientation not favored by the local physics. The development of oblique structure easily leads to turbulence inside small periodic systems.

Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earth's magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large electric-field amplitudes (up to 100 mV/m). These wave structures are associated with spatially oscillatory dissipation, which appears as alternatingly positive and negative values of J dot E (dissipation). For small guide magnetic fields the wave structures occur in the electron stagnation region at the magnetosphere edge of the EDR. For larger guide fields the structures also occur near the reconnection x-line. This difference is explained in terms of channels for the out-of-plane current (agyrotropic electrons at the stagnation point and guide-field-aligned electrons at the x-line).

A prediction of the steady-state reconnection electric field in asymmetric reconnection is obtained by maximizing the reconnection rate as a function of the opening angle made by the upstream magnetic field on the weak magnetic field (magnetosheath) side. The prediction is within a factor of two of the widely examined asymmetric reconnection model [Cassak and Shay, Phys. Plasmas 14, 102114, 2007] in the collisionless limit, and they scale the same over a wide parameter regime. The previous model had the effective aspect ratio of the diffusion region as a free parameter, which simulations and observations suggest is on the order of 0.1, but the present model has no free parameters. In conjunction with the symmetric case [Liu et al., Phys. Rev. Lett. 118, 085101, 2017], this work further suggests that this nearly universal number 0.1, essentially the normalized fast reconnection rate, is a geometrical factor arising from maximizing the reconnection rate within magnetohydrodynamic (MHD)-scale constraints.

Simulations suggest collisionless steady-state magnetic reconnection of Harris-type current sheets proceeds with a rate of order 0.1, independent of dissipation mechanism. We argue this long-standing puzzle is a result of constraints at the magnetohydrodynamic (MHD) scale. We perform a scaling analysis of the reconnection rate as a function of the opening angle made by the upstream magnetic fields, finding a maximum reconnection rate close to 0.2. The predictions compare favorably to particle-in-cell simulations of relativistic electron-positron and non-relativistic electron-proton reconnection. The fact that simulated reconnection rates are close to the predicted maximum suggests reconnection proceeds near the most efficient state allowed at the MHD-scale. The rate near the maximum is relatively insensitive to the opening angle, potentially explaining why reconnection has a similar fast rate in differing models.

There is growing evidence that when magnetic reconnection occurs in high Lundquist number plasmas such as in the Solar Corona or the Earth's Magnetosphere it does so within a fragmented, rather than a smooth current layer. Within the extent of these fragmented current regions the associated magnetic flux transfer and energy release occurs simultaneously in many different places. This investigation focusses on how best to quantify the rate at which reconnection occurs in such layers. An analytical theory is developed which describes the manner in which new connections form within fragmented current layers in the absence of magnetic nulls. It is shown that the collective rate at which new connections form can be characterized by two measures; a total rate which measures the true rate at which new connections are formed and a net rate which measures the net change of connection associated with the largest value of the integral of $E_{\|}$ through all of the non-ideal regions. Two simple analytical models are presented which demonstrate how each should be applied and what they quantify.

A model of global magnetic reconnection rate in relativistic collisionless plasmas is developed and validated by the fully kinetic simulation. Through considering the force balance at the upstream and downstream of the diffusion region, we show that the global rate is bounded by a value $\sim 0.3$ even when the local rate goes up to $\sim O(1)$ and the local inflow speed approaches the speed of light in strongly magnetized plasmas. The derived model is general and can be applied to magnetic reconnection under widely different circumstances.

Using fully kinetic simulations, we study the suppression of asymmetric reconnection in the limit where the diamagnetic drift speed >> Alfven speed and the magnetic shear angle is moderate. We demonstrate that the slippage between electrons and the magnetic flux facilitates reconnection, and can even result in fast reconnection that lacks one of the outflow jets. Through comparing a case where the diamagnetic drift is supported by the temperature gradient with a companion case that has a density gradient instead, we identify a robust suppression mechanism. The drift of the x-line is slowed down locally by the asymmetric nature of the current sheet and the resulting tearing modes, then the x-line is run over and swallowed by the faster-moving following flux.

Results are presented of a first study of collisionless magnetic reconnection starting from a recently found exact nonlinear force-free Vlasov-Maxwell equilibrium. The initial state has a Harris sheet magnetic field profile in one direction and a non-uniform guide field in a second direction, resulting in a spatially constant magnetic field strength as well as a constant initial plasma density and plasma pressure. It is found that the reconnection process initially resembles guide field reconnection, but that a gradual transition to anti-parallel reconnection happens as the system evolves. The time evolution of a number of plasma parameters is investigated, and the results are compared with simulations starting from a Harris sheet equilibrium and a Harris sheet plus constant guide field equilibrium.

Using fully kinetic simulations, we study the x-line orientation of magnetic reconnection in an asymmetric configuration. A spatially localized perturbation is employed to induce a single x-line, that has sufficient freedom to choose its orientation in three-dimensional systems. The effect of ion to electron mass ratio is investigated, and the x-line appears to bisect the magnetic shear angle across the current sheet in the large mass ratio limit. The orientation can generally be deduced by scanning through corresponding 2D simulations to find the reconnection plane that maximizes the peak reconnection electric field. The deviation from the bisection angle in the lower mass ratio limit can be explained by the physics of tearing instability.

It was recently proposed that the electron-frame dissipation measure, the energy transfer from the electromagnetic field to plasmas in the electron's rest frame, identifies the dissipation region of collisionless magnetic reconnection [Zenitani et al. Phys. Rev. Lett. 106, 195003 (2011)]. The measure is further applied to the electron-scale structures of antiparallel reconnection, by using two-dimensional particle-in-cell (PIC) simulations. The size of the central dissipation region is controlled by the electron-ion mass ratio, suggesting that electron physics is essential. A narrow electron jet extends along the outflow direction until it reaches an electron shock. The jet region appears to be anti-dissipative. At the shock, electron heating is relevant to a magnetic cavity signature. The results are summarized to a unified picture of the single dissipation region in a Hall magnetic geometry.

We investigate the one-dimensional interaction of a relativistic jet and an external medium. Relativistic magnetohydrodynamic simulations show an anomalous boost of the jet fluid in the boundary layer, as previously reported. We describe the boost mechanism using an ideal relativistic fluid and magnetohydrodynamic theory. The kinetic model is also examined for further understanding. Simple scaling laws for the maximum Lorentz factor are derived, and verified by the simulations.

Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter $\sigma$ and approaches the light speed when $\sigma > O(100)$, leading to an enhanced reconnection rate. In all regimes, the divergence of pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region remains $\sim 0.1$ in both the non-relativistic and relativistic limits.

The fluid-scale evolution of relativistic magnetic reconnection is investigated by using two-fluid and magnetofluid simulation models. Relativistic two-fluid simulations demonstrate the meso-scale evolution beyond the kinetic scales, and exhibit quasi-steady Petschek-type reconnection. Resistive relativistic MHD simulations further show new shock structures in and around the downstream magnetic island (plasmoid). General discussions on these models are presented.

A new measure to identify a small-scale dissipation region in collisionless magnetic reconnection is proposed. The energy transfer from the electromagnetic field to plasmas in the electron's rest frame is formulated as a Lorentz-invariant scalar quantity. The measure is tested by two-dimensional particle-in-cell simulations in typical configurations: symmetric and asymmetric reconnection, with and without the guide field. The innermost region surrounding the reconnection site is accurately located in all cases. We further discuss implications for nonideal MHD dissipation.

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

Resistive relativistic magnetohydrodynamic (RRMHD) simulations are applied to investigate the system evolution of relativistic magnetic reconnection. A time-split Harten--Lan--van Leer method is employed. Under a localized resistivity, the system exhibits a fast reconnection jet with an Alfv\'{e}nic Lorentz factor inside a narrow Petschek-type exhaust. Various shock structures are resolved in and around the plasmoid such as the post-plasmoid vertical shocks and the "diamond-chain" structure due to multiple shock reflections. Under a uniform resistivity, Sweet--Parker-type reconnection slowly evolves. Under a current-dependent resistivity, plasmoids are repeatedly formed in an elongated current sheet. It is concluded that the resistivity model is of critical importance for RRMHD modeling of relativistic magnetic reconnection.

The nonlinear evolution of relativistic magnetic reconnection in sheared magnetic configuration (with a guide field) is investigated by using two-dimensional relativistic two-fluid simulations. Relativistic guide field reconnection features the charge separation and the guide field compression in and around the outflow channel. As the guide field increases, the composition of the outgoing energy changes from enthalpy-dominated to Poynting-dominated. The inertial effects of the two-fluid model play an important role to sustain magnetic reconnection. Implications for the single fluid magnetohydrodynamic approach and the physics models of relativistic reconnection are briefly addressed.

We investigate the large scale evolution of a relativistic magnetic reconnection in an electron-positron pair plasma by a relativistic two-fluid magnetohydrodynamic (MHD) code. We introduce an inter-species friction force as an effective resistivity to dissipate magnetic fields. We demonstrate that magnetic reconnection successfully occurs in our two-fluid system, and that it involves Petschek-type bifurcated current layers in later stage. We further observe a quasi-steady evolution thanks to an open boundary condition, and find that the Petschek-type structure is stable over the long time period. Simulation results and theoretical analyses exhibit that the Petschek outflow channel becomes narrower when the reconnection inflow contains more magnetic energy, as previously claimed. Meanwhile, we find that the reconnection rate goes up to ~1 in extreme cases, which is faster than previously thought. The role of the resistivity, implications for reconnection models in the magnetically dominated limit, and relevance to kinetic reconnection works are discussed.

We investigate properties of the reconnecting current layer in relativistic pair plasma reconnection. We found that the current layer self-regulates its thickness when the current layer runs out current carriers, and so relativistic reconnection retains a fast reconnection rate. Constructing a steady state Sweet-Parker model, we discuss conditions for the current sheet expansion. Based on the energy argument, we conclude that the incompressible assumption is invalid in relativistic Sweet-Parker reconnection. The guide field cases are more incompressible than the anti-parallel cases, and we find a more significant current sheet expansion.

The role of the Weibel instability is investigated for the first time in the context of the large-scale magnetic reconnection problem. A late-time evolution of magnetic reconnection in relativistic pair plasmas is demonstrated by particle-in-cell (PIC) simulations. In the outflow regions, powerful reconnection jet piles up the magnetic fields and then a tangential discontinuity appears there. Further downstream, it is found that the two-dimensional extension of the relativistic Weibel instability generates electro-magnetic fields, which are comparable to the anti-parallel or piled-up fields. In a microscopic viewpoint, the instability allows plasma's multiple interactions with the discontinuity. In a macroscopic viewpoint, the instability leads to rapid expansion of the current sheet and then the reconnection jet front further propagates into the downstream. Possible application to the three-dimensional case is briefly discussed.

We present INTEGRAL and XMM-Newton observations of the prompt gamma-ray emission and the X-ray afterglow of GRB030227, the first GRB for which the quick localization obtained with the INTEGRAL Burst Alert System (IBAS) has led to the discovery of X-ray and optical afterglows. GRB030227 had a duration of about 20 s and a peak flux of 1.1 photons cm^-2 s^-1 in the 20-200 keV energy range. The time averaged spectrum can be fit by a single power law with photon index about 2 and we find some evidence for a hard to soft spectral evolution. The X-ray afterglow has been detected starting only 8 hours after the prompt emission, with a 0.2-10 keV flux decreasing as t^-1 from 1.3x10e-12 to 5x10e-13 erg cm^-2 s^-1. The afterglow spectrum is well described by a power law with photon index 1.94+/-0.05 modified by a redshifted neutral absorber with column density of several 10e22 cm^-2. A possible emission line at 1.67 keV could be due to Fe for a redshift z=3, consistent with the value inferred from the absorption.