Recently, magnetic reconnection, where a relativistically hot plasma is confined by a strong magnetic field, has received great attention in relation to astrophysical objects (e.g., pulsar magnetosphere and magnetar). However, reconnection with relativistic high-speed drift current in the plasma sheet has not been investigated yet. Thus, from both the theoretical and computational points of view, we studied the growth rate of relativistic reconnection for the high-speed drift current in a relativistically hot plasma sheet. Consequently, we argue that, contrary to the conventional understanding of the fast energy dissipation by the relativistic reconnection, the growth of magnetic reconnection is suppressed.

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

Fully kinetic two-dimensional particle-in-cell simulations are used to study electron acceleration at high-Mach-number nonrelativistic perpendicular shocks. SNR shocks are mediated by the Weibel instability which is excited because of an interaction between shock-reflected and upstream ions. Nonlinear evolution of the Weibel instability leads to the formation of current sheets. At the turbulent shock ramp the current sheets decay through magnetic reconnection. The number of reconnection sites strongly depends on the ion-to-electron mass ratio and the Alfv\'enic Mach number of the simulated shock. Electron acceleration is observed at locations where magnetic reconnection operates. For the highest mass ratios almost all electrons are involved in magnetic reconnection, which makes the magnetic reconnection the dominant acceleration process for electrons at these shocks. We discuss the relevance of our results for 3D systems with realistic ion-to-electron mass ratio.

The first-order Fermi acceleration of electrons requires an injection of electrons into a mildly relativistic energy range. However, the mechanism of injection has remained a puzzle both in theory and observation. We present direct evidence for a novel stochastic shock drift acceleration theory for the injection obtained with Magnetospheric Multiscale (MMS) observations at Earth's bow shock. The theoretical model can explain electron acceleration to mildly relativistic energies at high-speed astrophysical shocks, which may provide a solution to the long-standing issue of electron injection.

Particle acceleration in collisionless plasma systems is a central question in astroplasma and astroparticle physics. The structure of the acceleration regions, electron-ion energy equilibration, preacceleration of particles at shocks to permit further energization by diffusive shock acceleration, require knowledge of the distribution function of particles besides the structure and dynamic of electromagnetic fields, and hence a kinetic description is desirable. Particle-in-cell simulations offer an appropriate, if computationally expensive method of essentially conducting numerical experiments that explore kinetic phenomena in collisionless plasma. We review recent results of PIC simulations of astrophysical plasma systems, particle acceleration, and the instabilities that shape them.

We explore electron pre-acceleration at high Mach-number nonrelativistic perpendicular shocks at, e.g., young supernova remnants, which are a prerequisite of further acceleration to very high energies via diffusive shock acceleration. Using fully kinetic particle-in-cell simulations of shocks and electron dynamics in them, we investigate the influence of shock-surfing acceleration at the shock foot on the nonthermal population of electrons downstream of the shock. The shock-surfing acceleration is followed by further energization at the shock ramp where the Weibel instability spawns a type of second-order Fermi acceleration. The combination of these two processes leads to the formation of a nonthermal electron population, but the importance of shock-surfing acceleration becomes smaller for larger ion-to-electron mass ratio in the simulation. We discuss the resulting electron spectra and the relevance of our results to the physics of systems with real ion-to-electron mass ratio and fully three-dimensional behavior.

We investigated electromagnetic precursor wave emission in relativistic shocks by using two-dimensional particle-in-cell simulations. We found that the wave amplitude is significantly enhanced by a positive feedback process associated with ion-electron coupling through the wakefields for high magnetization. The wakefields collapse during the nonlinear process of the parametric decay instability in the near-upstream region, where nonthermal electrons and ions are generated. The intense coherent emission and the particle acceleration may opperate in high-energy astrophysical objects.

We developed an experimental platform for studying magnetic reconnection in an external magnetic field with simultaneous measurements of plasma imaging, flow velocity, and magnetic-field variation. Here, we investigate the stagnation and acceleration in counter-streaming plasmas generated by high-power laser beams. A plasma flow perpendicular to the initial flow directions is measured with laser Thomson scattering. The flow is, interestingly, accelerated toward the high-density region, which is opposite to the direction of the acceleration by pressure gradients. This acceleration is possibly interpreted by the interaction of two magnetic field loops initially generated by Biermann battery effect, resulting in a magnetic reconnection forming a single field loop and additional acceleration by a magnetic tension force.

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.

Electron injection at high Mach-number nonrelativistic perpendicular shocks is studied here for parameters that are applicable to young SNR shocks. Using high-resolution large-scale two-dimensional fully kinetic particle-in-cell (PIC) simulations and tracing individual particles we in detail analyze the shock surfing acceleration (SSA) of electrons at the leading edge of the shock foot. The central question is to what degree the process can be captured in 2D3V simulations. We find that the energy gain in SSA always arises from the electrostatic field of a Buneman wave. Electron energization is more efficient in the out-of-plane orientation of the large-scale magnetic field because both the phase speed and the amplitude of the waves are higher than for the in-plane scenario. Also, a larger number of electrons is trapped by the waves compared to the in-plane configuration. We conclude that significant modifications of the simulation parameters are needed to reach the same level of SSA efficiency as in simulations with out-of-plane magnetic field or 3D simulations.

We discuss a mechanism for intense electromagnetic wave emission at an astrophysical relativistic shock in a magnetized collisionless plasma. At the magnetized shock, the particle reflection by a compressed magnetic field of the shock produces a ring-like distribution in momentum, which gives rise to plasma instabilities. Intense and coherent high-frequency electromagnetic waves will be emitted if the synchrotron maser instability (SMI) is excited, whereas non-propagating magnetic fluctuations will be generated when the Weibel instability (WI) is the dominant mode. The problem is of great astrophysical interest because if intense radiation is emitted, the interaction with the upstream medium induces a large-amplitude electrostatic field (or Wakefield), which may play a role for the acceleration of ultra-high-energy cosmic rays. We review our recent effort to measure the efficiency of the electromagnetic wave emission using fully self-consistent, two-dimensional (2D) particle-in-cell (PIC) simulations for pair plasmas. We found that the emission efficiency in 2D was systematically lower than one dimensional (1D) PIC simulation results. However, the power remains finite even when the WI is active to generate large-amplitude magnetic fluctuations. Astrophysical implications of the present results are briefly discussed.

The plasma heating during collisionless magnetic reconnection is investigated using particle-in-cell simulations. We analyze the time evolution of the plasma temperature associated with the motion of the reconnecting flux tube, where the plasma temperature is defined as second-order moment of velocity distribution function in the simulation frame/in the center of the flux tube frame, and we show that the plasma heating during magnetic reconnection can be separated into two distinct stages: the nonadiabatic heating stage, in which the magnetic field lines are just reconnecting in the X-type diffusion region, and the adiabatic heating stage, in which the flux tube is shrinking after two flux tubes merge. During the adiabatic heating stage, the plasma temperature $T$ can be approximated by $TV^{\gamma-1}=const.$, where $\gamma=5/3$ is the specific heat, and $V$ is the volume of the flux tube. In the nonadiabatic heating stage, we found numerically that the ratio of the increment of the ion temperature to that of the electron temperature can be approximated by $\Delta T_i/\Delta T_e \approx (m_i/m_e )^{1/4}$, where $m_i$ and $m_e$ are the ion and electron mass, respectively. We also present a theoretical model based on a magnetic-diffusion-dominated reconnection to explain the simulation result.

We investigated the precursor wave emission efficiency in magnetized purely perpendicular relativistic shocks in pair plasmas. We extended our previous study to include the dependence of upstream magnetic field orientations. We performed two-dimensional particle-in-cell simulations and focused on two magnetic field orientations: the magnetic field to be in the simulation plane (i.e., in-plane configuration) and perpendicular to the simulation plane (i.e., out-of-plane configuration). Our simulations in the in-plane configuration demonstrated that not only extraordinary but also ordinary mode waves are excited. We quantified the emission efficiency as a function of the magnetization parameter $\sigma_e$ and found that the large-amplitude precursor waves are emitted for a wide range of $\sigma_e$. We found that especially at low $\sigma_e$, the magnetic field generated by Weibel instability amplifies the ordinary mode wave power. The amplitude is large enough to perturb the upstream plasma, and transverse density filaments are generated as in the case of the out-of-plane configuration investigated in the previous study. We confirmed that our previous conclusion holds regardless of upstream magnetic field orientations with respect to the two-dimensional simulation plane. We discuss the precursor wave emission in three dimensions and the feasibility of wakefield acceleration in relativistic shocks based on our results.

How electrons get accelerated to relativistic energies in a high-Mach-number quasi-perpendicular shock is presented by means of ab initio particle-in-cell simulations in three dimensions. We found that coherent electrostatic Buneman waves and ion-Weibel magnetic turbulence coexist in a strong-shock structure whereby particles gain energy during shock-surfing and subsequent stochastic drift accelerations. Energetic electrons that initially experienced the surfing acceleration undergo pitch-angle diffusion by interacting with magnetic turbulence and continuous acceleration during confinement in the shock transition region. The ion-Weibel turbulence is the key to the efficient nonthermal electron acceleration.

This paper presents a series of stratified shearing-box simulations of collisionless accretion disks in the recently developed framework of kinetic magnetohydrodynamics (MHD), which can handle finite non-gyrotropy of a pressure tensor. Although a fully kinetic simulation predicted a more efficient angular-momentum transport in collisionless disks than in the standard MHD regime, the enhanced transport has not been observed in past kinetic MHD approaches to gyrotropic pressure anisotropy. For the purpose of investigating this missing link between the fully kinetic and MHD treatments, this paper pays attention to the role of non-gyrotropic pressure, and makes a first attempt to incorporate certain collisionless effects into disk-scale, stratified disk simulations. When the timescale of gyrotropization was longer than, or comparable to, the disk rotation frequency of the orbit, we found that the finite non-gyrotropy selectively remaining in the vicinity of current sheets contributes to suppressing magnetic reconnection in the shearing-box system. This leads to increases both in the saturated amplitude of the MHD turbulence driven by magnetorotational instabilities and in the resultant efficiency of angular-momentum transport. Our results seem favorable for fast advection of magnetic fields toward the rotation axis of a central object, which is required to launch an ultra-relativistic jet from a black-hole accretion system in, for example, a magnetically arrested disk state.

Pulsar Wind Nebulae, Blazars, Gamma Ray Bursts and Magnetars all contain regions where the electromagnetic energy density greatly exceeds the plasma energy density. These sources exhibit dramatic flaring activity where the electromagnetic energy distributed over large volumes, appears to be converted efficiently into high energy particles and gamma-rays. We call this general process magnetoluminescence. Global requirements on the underlying, extreme particle acceleration processes are described and the likely importance of relativistic beaming in enhancing the observed radiation from a flare is emphasized. Recent research on fluid descriptions of unstable electromagnetic configurations are summarized and progress on the associated kinetic simulations that are needed to account for the acceleration and radiation is discussed. Future observational, simulation and experimental opportunities are briefly summarized.

We investigated the efficiency of coherent upstream large-amplitude electromagnetic wave emission via synchrotron maser instability at relativistic magnetized shocks by using two-dimensional particle-in-cell simulations. We considered the purely perpendicular shock in an electron-positron plasma. The coherent wave emission efficiency was measured as a function of the magnetization parameter {\sigma}, which is defined by the ratio of the Poynting flux to the kinetic energy flux. The wave amplitude was systematically smaller than that observed in one-dimensional simulations. However, it continued to persist, even at a considerably low magnetization rate, where the Weibel instability dominated the shock transition. The emitted electromagnetic waves were sufficiently strong to disturb the upstream medium, and transverse filamentary density structures of substantial amplitude were produced. Based on this result, we discuss the possibility of the wakefield acceleration model for the production of non-thermal electrons in a relativistic magnetized ion-electron shock.

A new type of instability that is expected to drive magnetohydrodynamic (MHD) turbulence from a purely toroidal magnetic field in an accretion disk is presented. It is already known that in a differentially rotating system, the uniform toroidal magnetic field is unstable due to a magnetorotational instability (MRI) under a non-axisymmetric and vertical perturbation, while it is stable under a purely vertical perturbation. Contrary to the previous study, this paper proposes an unstable mode completely confined to the equatorial plane, driven by the expansive nature of the magnetic pressure gradient force under a non-uniform toroidal field. The basic nature of this growing eigenmode, to which we give a name "magneto-gradient driven instability", is studied using linear analysis, and the corresponding nonlinear evolution is then investigated using two-dimensional ideal MHD simulations. Although a single localized magnetic field channel alone cannot provide sufficient Maxwell stress to contribute significantly to the angular momentum transport, we find that the mode coupling between neighboring toroidal fields under multiple localized magnetic field channels drastically generates a highly turbulent state and leads to the enhanced transport of angular momentum, comparable to the efficiency seen in previous studies on MRIs. This horizontally confined mode may play an important role in the saturation of an MRI through complementray growth with the toroidal MHDs and coupling with magnetic reconnection.

The magnetohydrodynamic linear stability with the localized bulk flow oriented parallel to the neutral sheet is investigated, by including the Hall effect and the guide magnetic field. We observe three different unstable modes: a "streaming tearing" mode at a slow flow speed, a "streaming sausage" mode at a medium flow speed, and a "streaming kink" mode at a fast flow speed. The streaming tearing and sausage modes have a standard tearing mode-like structure with symmetric density fluctuations in the neutral sheet, while the kink mode has an asymmetric fluctuation. The growth rate of the streaming tearing mode decreases with increasing magnetic Reynolds number, while the growth rates of the sausage and kink modes do not depend strongly on the Reynolds number. The sausage and kink modes can be unstable for not only super-Alfv\'enic flow but also sub-Alfv\'enic flow when the lobe density is low. The wavelengths of these unstable modes are of the same order of magnitude as the thickness of the plasma sheet. Their maximum growth rates are higher than that of a standard tearing mode, and under a strong guide magnetic field, the growth rates of the sausage and kink modes are enhanced, while under a weak guide magnetic field, they are suppressed. For a thin plasma sheet with the Hall effect, the fluctuations of the streaming modes can exist over the plasma sheet. These unstable modes may be regarded as being one of the processes generating Alfv\'enic turbulence in the plasma sheet during magnetic reconnection.

Angular momentum transport and particle acceleration during the magnetorotational instability (MRI) in a collisionless accretion disk are investigated using three-dimensional particle-in-cell (PIC) simulation. We show that the kinetic MRI can provide not only high energy particle acceleration but also enhancement of angular momentum transport. We find that the plasma pressure anisotropy inside the channel flow with $p_{\|} > p_{\perp}$ induced by active magnetic reconnection suppresses the onset of subsequent reconnection, which in turn leads to high magnetic field saturation and enhancement of Maxwell stress tensor of angular momentum transport. Meanwhile, during the quiescent stage of reconnection the plasma isotropization progresses in the channel flow, and the anisotropic plasma with $p_{\perp} > p_{\|}$ due to the dynamo action of MRI outside the channel flow contributes to rapid reconnection and strong particle acceleration. This efficient particle acceleration and enhanced angular momentum transport in a collisionless accretion disk may explain the origin of high energy particles observed around massive black holes.

An evolution of a magnetic reconnection in a collisionless accretion disk is investigated using a 2.5 dimensional hybrid code simulation. In astrophysical disks, magnetorotational instability (MRI) is considered to play an important role by generating turbulence in the disk and contributes to an effective angular momentum transport through a turbulent viscosity. Magnetic reconnection, on the other hand, also plays an important role on the evolution of the disk through a dissipation of a magnetic field enhanced by a dynamo effect of MRI. In this study, we developed a hybrid code to calculate an evolution of a differentially rotating system. With this code, we first confirmed a linear growth of MRI. We also investigated a behavior of a particular structure of a current sheet, which would exist in the turbulence in the disk. From the calculation of the magnetic reconnection, we found an asymmetric structure in the out-of-plane magnetic field during the evolution of reconnection, which can be understood by a coupling of the Hall effect and the differential rotation. We also found a migration of X-point whose direction is determined only by an initial sign of J_0 \times \Omega_0, where J_0 is the initial current density in the neutral sheet and \Omega_0 is the rotational vector of the background Keplerian rotation. Associated with the migration of X-point, we also found a significant enhancement of the perpendicular magnetic field compared to an ordinary MRI. MRI-Magnetic reconnection coupling and the resulting magnetic field enhancement can be an effective process to sustain a strong turbulence in the accretion disk and to a transport of angular momentum.

Particle acceleration in the magnetic reconnection of electron-positron plasmas is studied by using a particle-in-cell simulation. It is found that a significantly large number of nonthermal particles are generated by the inductive electric fields around an X-type neutral line when the reconnection outflow velocity, which is known to be an Alfv\'{e}n velocity, is on the order of the speed of light. In such a relativistic reconnection regime, we also find that electrons and positrons form a power-law-like energy distribution through their drift along the reconnection electric field under the relativistic Speiser motion. A brief discussion of the relevance of these results to the current sheet structure, which has an antiparallel magnetic field in astrophysical sources of synchrotron radiation, is presented.

Through the enhancement of transport, turbulence is expected to contribute to the fast reconnection. However the effects of turbulence are not so straightforward. In addition to the enhancement of transport, turbulence under some environment shows effects that suppress the transport. In the presence of turbulent cross helicity, such a dynamic balance between the transport enhancement and suppression occurs. As this result of dynamic balance, the region of effective enhanced magnetic diffusivity is confined to a narrow region, leading to the fast reconnection. In order to confirm this idea, a self-consistent turbulence model for the magnetic reconnection is proposed. With the aid of numerical simulations where turbulence effects are incorporated in a consistent manner through the turbulence model, the dynamic balance in the turbulence magnetic reconnection is confirmed.

Electron acceleration associated with various plasma kinetic instabilities in a nonrelativistic, very-high-Alfv\'en Mach-number ($M_A \sim 45$) shock is revealed by means of a two-dimensional fully kinetic PIC simulation. Electromagnetic (ion Weibel) and electrostatic (ion-acoustic and Buneman) instabilities are strongly activated at the same time in different regions of the two-dimensional shock structure. Relativistic electrons are quickly produced predominantly by the shock surfing mechanism with the Buneman instability at the leading edge of the foot. The energy spectrum has a high-energy tail exceeding the upstream ion kinetic energy accompanying the main thermal population. This gives a favorable condition for the ion acoustic instability at the shock front, which in turn results in additional energization. The large-amplitude ion Weibel instability generates current sheets in the foot, implying another dissipation mechanism via magnetic reconnection in a three-dimensional shock structure in the very-high-$M_A$ regime.

We study the formation of slow-mode shocks in collisionless magnetic reconnection by using one- and two-dimensional collisionless MHD codes based on the double adiabatic approximation and the Landau closure model. We bridge the gap between the Petschek-type MHD reconnection model accompanied by a pair of slow shocks and the observational evidence of the rare occasion of in-situ slow shock observation. Our results showed that once magnetic reconnection takes place, a firehose-sense pressure anisotropy arises in the downstream region, and the generated slow shocks are quite weak comparing with those in an isotropic MHD. In spite of the weakness of the shocks, however, the resultant reconnection rate is 10-30% higher than that in an isotropic case. This result implies that the slow shock does not necessarily play an important role in the energy conversion in the reconnection system, and is consistent with the satellite observation in the Earth's magnetosphere.

The stability of the cosmic ray modified shock (CRMS) is studied by means of numerical simulations. Owing to the nonlinear feedback of cosmic-ray (CR) acceleration, a downstream state of the modified shock can no longer be uniquely determined for given upstream parameters. It is known that up to three distinct solutions exist, which are characterized by the CR production efficiency as "efficient", "intermediate" and "inefficient" branches. The stability of these solutions is investigated by performing direct time-dependent simulations of a two-fluid model. It is found that both the efficient and inefficient branches are stable even against a large-amplitude perturbation, while the intermediate one is always unstable and evolves into the inefficient state as a result of nonlinear time development. This bistable feature is robust in a wide range of parameters and does not depend on the injection model. Fully nonlinear time evolution of a hydrodynamic shock with injection results in the least efficient state in terms of the CR production, consistent with the bistable feature. This suggests that the CR production efficiency at supernova remnant shocks may be lower than previously discussed in the framework of the nonlinear shock acceleration theory considering the efficient solution of the CRMS.

Particle acceleration during the magnetorotational instability (MRI) in a collisionless accretion disk was investigated by using a particle-in-cell (PIC) simulation. We discuss the important role that magnetic reconnection plays not only on the saturation of MRI but also on the relativistic particle generation. The plasma pressure anisotropy of $p_{\perp} > p_{\para}$ induced by the action of MRI dynamo leads to rapid growth in magnetic reconnection, resulting in the fast generation of nonthermal particles with a hard power-law spectrum. This efficient particle acceleration mechanism involved in a collisionless accretion disk may be a possible model to explain the origin of high energy particles observed around massive black holes.

We report simulation results for turbulent magnetic reconnection obtained using a newly developed Reynolds-averaged magnetohydrodynamics model. We find that the initial Harris current sheet develops in three ways, depending on the strength of turbulence: laminar reconnection, turbulent reconnection, and turbulent diffusion. The turbulent reconnection explosively converts the magnetic field energy into both kinetic and thermal energy of plasmas, and generates open fast reconnection jets. This fast turbulent reconnection is achieved by the localization of turbulent diffusion. Additionally, localized structure forms through the interaction of the mean field and turbulence.

Electron accelerations at high Mach number collision-less shocks are investigated by means of two-dimensional electromagnetic Particle-in-Cell simulations with various Alfven Mach numbers, ion-to-electron mass ratios, and the upstream electron beta_e (the ratio of the thermal pressure to the magnetic pressure). We found electrons are effectively accelerated at a super-high Mach number shock (MA~30) with a mass ratio of M/m=100 and beta_e=0.5. The electron shock surfing acceleration is an effective mechanism for accelerating the particles toward the relativistic regime even in two dimensions with the large mass ratio. Buneman instability excited at the leading edge of the foot in the super-high Mach number shock results in a coherent electrostatic potential structure. While multi-dimensionality allows the electrons to escape from the trapping region, they can interact with the strong electrostatic field several times. Simulation runs in various parameter regimes indicate that the electron shock surfing acceleration is an effective mechanism for producing relativistic particles in extremely-high Mach number shocks in supernova remnants, provided that the upstream electron temperature is reasonably low.

We perform a two-dimensional simulation by using an electromagnetic hybrid code to study the formation of slow-mode shocks in collisionless magnetic reconnection in low beta plasmas, and we focus on the relation between the formation of slow shocks and the ion temperature anisotropy enhanced at the shock downstream region. It is known that as magnetic reconnection develops, the parallel temperature along the magnetic field becomes large in association with the anisotropic PSBL (plasma sheet boundary layer) ion beams, and this temperature anisotropy has a tendency to suppress the formation of slow shocks. Based on our simulation result, we found that the slow shock formation is suppressed due to the large temperature anisotropy near the X-type region, but the ion temperature anisotropy relaxes with increasing the distance from the magnetic neutral point. As a result, two pairs of current structures, which are the strong evidence of dissipation of magnetic field in slow shocks, are formed at the distance x > 115 ion inertial lengths from the neutral point.

A nonthermal particle acceleration mechanism involving the interaction of a charged particle with multiple magnetic islands is proposed. The original Fermi acceleration model, which assumes randomly distributed magnetic clouds moving at random velocity $V_c$ in the interstellar medium, is known to be of second-order acceleration of $O(V_c/c)^2$ owing to the combination of head-on and head-tail collisions. In this letter, we reconsider the original Fermi model by introducing multiple magnetic islands during reconnect ion instead of magnetic clouds. We discuss that the energetic particles have a tendency to be distributed outside the magnetic islands, and they mainly interact with reconnection outflow jets. As a result, the acceleration efficiency becomes first-order of $O(V_A/c)$, where $V_A$ and $c$ are the Alfv\'en velocity and the speed of light, respectively.

Roles of turbulence in the context of magnetic reconnection are investigated with special emphasis on the mutual interaction between flow (large-scale inhomogeneous structure) and turbulence. In order to evaluate the effective transport due to turbulence, in addition to the {\it intensity} information of turbulence represented by the turbulent energy, the {\it structure} information represented by pseudoscalar statistical quantities (helicities) is important. On the basis of the evolution equation, mechanisms that provide turbulence with cross helicity are presented. Magnetic-flux freezing in highly turbulent media is considered with special emphasis on the spatial distribution of the turbulent cross helicity. The cross-helicity effects in the context of magnetic reconnection are also investigated. It is shown that the large-scale flow and magnetic-field configurations favorable for the cross-helicity generation is compatible with the fast reconnection. In this sense, turbulence and large-scale structures promote magnetic reconnection mediated by the turbulent cross helicity.

Electron acceleration in collisionless shocks with arbitrary magnetic field orientations is discussed. It is shown that the injection of thermal electrons into diffusive shock acceleration process is achieved by an electron beam with a loss-cone in velocity space that is reflected back upstream from the shock through shock drift acceleration mechanism. The electron beam is able to excite whistler waves which can scatter the energetic electrons themselves when the Alfven Mach number of the shock is sufficiently high. A critical Mach number for the electron injection is obtained as a function of upstream parameters. The application to supernova remnant shocks is discussed.

A new rapid energization process within a supernova shock transition region (STR) is reported by utilizing numerical simulation. Although the scale of a STR as a main dissipation region is only several hundreds of thousands km, several interesting structures are found relating to generation of a root of the energetic particles. The nonlinear evolution of plasma instabilities lead to a dynamical change in the ion phase space distribution which associates with change of the field properties. As a result, different types of large-amplitude field structures appear. One is the leading wave packet and another is a series of magnetic solitary humps. Each field structure has a microscopic scale (~ the ion inertia length). Through the multiple nonlinear scattering between these large-amplitude field structures, electrons are accelerated directly. Within a STR, quick thermalization realizes energy equipartition between the ion and electron, hot electrons play an important role in keeping these large-amplitude field structures on the ion-acoustic mode. The hot electron shows non-Maxwellian distribution and could be the seed of further non-thermal population. The "shock system", where fresh incoming and reflected ions are supplied constantly, play an essential role in our result. With a perpendicular shock geometry, the maximum energy of the electron is estimated by equating a width of the STR to a length of the Larmor radius of the energetic electron. Under some realistic condition of M_A = 170 and \omega_{pe}/\Omega_{ce} = 120, maximum energy is estimated to ~ 10 MeV at an instant only within the STR.

The linear and non-linear evolution of a relativistic current sheet of pair ($e^{\pm}$) plasmas is investigated by three-dimensional particle-in-cell simulations. In a Harris configuration, it is obtained that the magnetic energy is fast dissipated by the relativistic drift kink instability (RDKI). However, when a current-aligned magnetic field (the so-called "guide field") is introduced, the RDKI is stabilized by the magnetic tension force and it separates into two obliquely-propagating modes, which we call the relativistic drift-kink-tearing instability (RDKTI). These two waves deform the current sheet so that they trigger relativistic magnetic reconnection at a crossover thinning point. Since relativistic reconnection produces a lot of non-thermal particles, the guide field is of critical importance to study the energetics of a relativistic current sheet.

Electron acceleration mechanism at high Mach number collisionless shocks propagating in a weakly magnetized medium is investigated by a self-consistent two-dimensional particle-in-cell simulation. Simulation results show that strong electrostatic waves are excited via the electron-ion electrostatic two-stream instability at the leading edge of the shock transition region as in the case of earlier one-dimensional simulations. We observe strong electron acceleration that is associated with the turbulent electrostatic waves in the shock transition region. The electron energy spectrum in the shock transition region exhibits a clear power-law distribution with spectral index of $2.0 {\rm -} 2.5$. By analyzing the trajectories of accelerated electrons, we find that the acceleration mechanism is very similar to shock surfing acceleration of ions. In contrast to the ion shock surfing, however, the energetic electrons are reflected by electron-scale electrostatic fluctuations in the shock transition region, but not by the ion-scale cross-shock electrostatic potential. The reflected electrons are then accelerated by the convective electric field in front of the shock. We conclude that the multidimensional effects as well as the self-consistent shock structure are essential for the strong electron acceleration at high Mach number shocks.

We study the role of the guide field in relativistic magnetic reconnection in a Harris current sheet of pair ($e^{\pm}$) plasmas, using linear theories and particle-in-cell (PIC) simulations. Two-dimensional PIC simulations exhibit the guide field dependence to the linear instabilities; the tearing or reconnection modes are relatively insensitive, while the relativistic drift-kink instability (RDKI), the fastest mode in a relativistic current sheet, is stabilized by the guide field. Particle acceleration in the nonlinear stage is also investigated. A three-dimensional PIC simulation demonstrates that the current sheet is unstable to the RDKI, although small reconnection occurs in the deformed current sheet. Another three-dimensional PIC simulation with a guide field demonstrates a completely different scenario. Secondary magnetic reconnection is triggered by nonlinear coupling of oblique instabilities, which we call the relativistic drift-sausage tearing instability. Therefore, particle acceleration by relativistic guide field reconnection occurs in three-dimensional configuration. Based on the plasma theories, we discuss an important role of the guide field: to enable non-thermal particle acceleration by magnetic reconnection.

We study linear and nonlinear development of relativistic and ultrarelativistic current sheets of pair plasmas with antiparallel magnetic fields. Two types of two-dimensional problems are investigated by particle-in-cell simulations. First, we present the development of relativistic magnetic reconnection, whose outflow speed is an order of the light speed c. It is demonstrated that particles are strongly accelerated in and around the reconnection region, and that most of magnetic energy is converted into "nonthermal" part of plasma kinetic energy. Second, we present another two-dimensional problem of a current sheet in a cross-field plane. In this case, the relativistic drift kink instability (RDKI) occurs. Particle acceleration also takes place, but the RDKI fast dissipates the magnetic energy into plasma heat. We discuss the mechanism of particle acceleration and the theory of the RDKI in detail. It is important that properties of these two processes are similar in the relativistic regime of T > mc^2, as long as we consider the kinetics. Comparison of the two processes indicates that magnetic dissipation by the RDKI is more favorable process in the relativistic current sheet. Therefore the striped pulsar wind scenario should be reconsidered by the RDKI.

A particle acceleration mechanism by radiation pressure of precursor waves in a relativistic shock is studied. For a relativistic, perpendicular shock with the upstream bulk Lorentz factor of $\gamma_1 \gg 1$, large amplitude electromagnetic (light) waves are known to be excited in the shock front due to the synchrotron maser instability, and those waves can propagate towards upstream as precursor waves. We find that non-thermal, high energy electrons and ions can be quickly produced by an action of electrostatic wakefields generated by the ponderomotive force of the precursor waves. The particles can be quickly accelerated up to $\epsilon_{\rm max}/\gamma_1 m_e c^2 \sim \gamma_1$ in the upstream coherent wakefield region, and they can be further accelerated during the nonlinear stage of the wakefield evolution. The maximum attainable energy is estimated by $\epsilon_{\rm max}/\gamma_1 m_e c^2 \sim L_{\rm sys}/(c/\omega_{pe})$, where $L_{\rm sys}$ and $c/\omega_{pe}$ are the size of an astrophysical object and the electron inertial length, respectively.

Electron injection process at high Mach number collisionless quasi-perpendicular shock waves is investigated by means of one-dimensional electromagnetic particle-in-cell simulations. We find that energetic electrons are generated through the following two steps: (1) electrons are accelerated nearly perpendicular to the local magnetic field by shock surfing acceleration at the leading edge of the shock transition region. (2) the preaccelerated electrons are further accelerated by shock drift acceleration. As a result, energetic electrons are preferentially reflected back to the upstream. Shock surfing acceleration provides sufficient energy required for the reflection. Therefore, it is important not only for the energization process by itself, but also for triggering the secondary acceleration process. We also present a theoretical model of the two-step acceleration mechanism based on the simulation results, which can predict the injection efficiency for subsequent diffusive shock acceleration process. We show that the injection efficiency obtained by the present model agrees well with the value obtained by Chandra X-ray observations of SN 1006. At typical supernova remnant shocks, energetic electrons injected by the present mechanism can self-generate upstream Alfven waves, which scatter the energetic electrons themselves.

We discuss that energetic electrons are generated near the X-type magnetic reconnection region due to a surfing acceleration mechanism. In a thin plasma sheet, the polarization electric fields pointing towards the neutral sheet are induced around the boundary between the lobe and plasma sheet in association with the Hall electric current. By using a particle-in-cell simulation, we demonstrate that the polarization electric fields are strongly enhanced in an externally driven reconnection system, and some electrons can be trapped by the electrostatic potential well of the polarization field. During the trapping phase, the electrons can gain their energies from the convection/inductive reconnection electric fields. We discuss that relativistic electrons with MeV energies are quickly generated in and around the X-type neutral region by utilizing the surfing acceleration.

Two-dimensional particle simulations of a relativistic Harris current sheet of pair plasmashave demonstrated that the system is unstable to the relativistic drift kink instability (RDKI) and that a new kind of acceleration process takes place in the deformed current sheet. This process contributes to the generation of non-thermal particles and contributes to the fast magnetic dissipation in the current sheet structure. The acceleration mechanism and a brief comparison with relativistic magnetic reconnection are presented.

On December 27, 2004, plasma particle detectors on the GEOTAIL spacecraft detected an extremely strong signal of hard X-ray photons from the giant flare of SGR1806-20, a magnetar candidate. While practically all gamma-ray detectors on any satellites were saturated during the first ~500 ms interval after the onset, one of the particle detectors on GEOTAIL was not saturated and provided unique measurements of the hard X-ray intensity and the profile for the first 600 ms interval with 5.48 ms time resolution. After ~50 ms from the initial rapid onset, the peak photon flux (integrated above ~50 keV) reached the order of 10^7 photons sec^{-1} cm^{-2}. Assuming a blackbody spectrum with kT=175 keV, we estimate the peak energy flux to be 21 erg sec^{-1} cm^{-2} and the fluence (for 0-600 ms) to be 2.4 erg cm^{-2}. The implied energy release comparable to the magnetic energy stored in a magnetar (~10^{47} erg) suggests an extremely efficient energy release mechanism.

We study the suprathermal electron acceleration mechanism in a perpendicular magnetosonic shock wave in a high Mach number regime by using a particle-in-cell simulation. We find that shock surfing/surftron acceleration producing the suprathermal electrons occurs in the shock transition region where a series of large amplitude electrostatic solitary waves (ESWs) are excited by Buneman instability under the interaction between the reflected ions and the incoming electrons. It is shown that the electrons are likely to be trapped by ESWs, and during the trapping phase they can be effectively accelerated by the shock motional/convection electric field. We discuss that suprathermal electrons can be accelerated up to $m_i c^2 (v_0/c)$, where $m_i c^2$ is the ion rest mass energy and $v_0$ is the shock upstream flow velocity. Furthermore, some of these suprathermal electrons may be effectively trapped for infinitely long time when Alfv\'en Mach number $M_A$ exceeds several 10, and they are accelerated up to the shock potential energy determined by the global shock size.