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.

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.

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 propose a novel electron acceleration mechanism, which we call stochastic shock drift acceleration (SSDA), that extends the standard shock drift acceleration (SDA) for low-energy electrons at a quasi-perpendicular shock to include the effect of stochastic pitch-angle scattering. We demonstrate that the steady-state energy spectrum of electrons accelerated within the shock transition region becomes a power-law in the limit of strong scattering. The spectral index is independent of the pitch-angle scattering coefficient. On the other hand, the maximum energy attainable through the mechanism scales linearly with the pitch-angle scattering coefficient. These results have been confirmed by Monte Carlo simulations that include finite pitch-angle anisotropy. We find that the theory can reasonably well explain in-situ observations of quasi-perpendicular Earth's bow shock. Theoretical scaling law suggests that the maximum energy increases in proportion to the square of the shock speed, indicating that the thermal electrons may be accelerated up to mildly relativistic energies by the SSDA at quasi-perpendicular supernova remnant shocks. Therefore, the mechanism provides a plausible solution to the long-standing electron injection problem.

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.

A generalized fluid-particle hybrid model for collisionless plasmas under the assumption of quasi-neutrality is presented. The system consists of fluid ions and electrons as well as arbitrary numbers of species whose dynamics is governed by the Vlasov equation. The proposed model is thus a generalized version of the well-known standard hybrid plasma simulation model, in which the ions are fully kinetic whereas the electrons are assumed to be a fluid. Since the proposed model employs the exact form of the generalized Ohm's law, the mass and energy densities, as well as the charge-to-mass ratio of the kinetic species, are taken to be arbitrary. In the absence of the kinetic species, it reduces to the quasi-neutral two-fluid model (Amano, 2015, JCP). In the opposite situation where the mass and energy densities of the kinetic species are much larger than the fluid ions, it is nothing more than the standard hybrid model with finite electron inertia effect. If the kinetic species is an energetic particle (EP) population having negligible mass density but substantial energy density and the scale size is much larger than the ion and electron inertial lengths, it describes the self-consistent coupling between the magnetohydrodynamics (MHD) and the EP dynamics. The energetic-particle-magnetohydrodynamics (EP-MHD) hybrid model is thus a special case of the more general model described in this paper. Advantages of this approach over the existing models are discussed. A three-dimensional simulation code solving the proposed set of equations is described. The code combines the Particle-in-Cell scheme for solving the kinetic species and a Riemann-solver-based code for the two-fluid equations. Several benchmark simulation results are shown to confirm that the code successfully captures the dynamics of the EP population interacting self-consistently with the MHD fluid.

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.

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 multidimensional simulation code for relativistic two-fluid electrodynamics (RTFED) is described. The basic equations consist of the full set of Maxwell's equations coupled with relativistic hydrodynamic equations for separate two charged fluids, representing the dynamics of either an electron-positron or an electron-proton plasma. It can be recognized as an extension of conventional relativistic magnetohydrodynamics (RMHD). Finite resistivity may be introduced as a friction between the two species, which reduces to resistive RMHD in the long wavelength limit without suffering from a singularity at infinite conductivity. A numerical scheme based on HLL (Harten-Lax-Van Leer) Riemann solver is proposed that exactly preserves the two divergence constraints for Maxwell's equations simultaneously. Several benchmark problems demonstrate that it is capable of describing RMHD shocks/discontinuities at long wavelength limit, as well as dispersive characteristics due to the two-fluid effect appearing at small scales. This shows that the RTFED model is a promising tool for high energy astrophysics application.

We report the first extragalactic detection of CF+, the fluoromethylidynium ion, in the z=0.89 absorber toward PKS1830-211. We estimate an abundance of ~3E-10 relative to H2 and that ~1% of fluorine is captured in CF+. The absorption line profile of CF+ is found to be markedly different from that of other species observed within the same tuning, and is notably anti-correlated with CH3OH. On the other hand, the CF+ profile resembles that of [C I]. Our results are consistent with expected fluorine chemistry and point to chemical differentiation in the column of absorbing gas.

In various astrophysics settings it is common to have a two-fluid relativistic plasma that interacts with the electromagnetic field. While it is common to ignore the displacement current in the ideal, classical magnetohydrodynamic limit, when the flows become relativistic this approximation is less than absolutely well-justified. In such a situation, it is more natural to consider a positively charged fluid made up of positrons or protons interacting with a negatively charged fluid made up of electrons. The two fluids interact collectively with the full set of Maxwell's equations. As a result, a solution strategy for that coupled system of equations is sought and found here. Our strategy extends to higher orders, providing increasing accuracy. Three important innovations are reported here. In our first innovation, the magnetic field within each zone is reconstructed in a divergence-free fashion while the electric field within each zone is reconstructed in a form that is consistent with Gauss' law. In our second innovation, a multidimensionally upwinded strategy is presented which ensures that the magnetic field can be updated via a discrete interpretation of Faraday's law and the electric field can be updated via a discrete interpretation of the generalized Ampere's law. Our third innovation consists of an efficient design for several popular IMEX schemes so that they provide strong coupling between the finite-volume-based fluid solver and the electromagnetic fields at high order. Several accuracy analyses are presented showing that our method meets its design accuracy in the MHD limit as well as in the limit of electromagnetic wave propagation. Several stringent test problems are also presented. We also present a relativistic version of the GEM problem, which shows that our algorithm can successfully adapt to challenging problems in high energy astrophysics.

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.

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.

The dynamics of a standing shock front in a Poynting-flux dominated relativistic flow is investigated by using a one-dimensional, relativistic, two-fluid simulation. An upstream flow containing a circularly polarized, sinusoidal magnetic shear wave is considered, mimicking a wave driven by an obliquely rotating pulsar. It is demonstrated that this wave is converted into large amplitude electromagnetic waves with superluminal phase speeds by interacting with the shock when the shock-frame frequency of the wave exceeds the proper plasma frequency. The superluminal waves propagate in the upstream, modify the shock structure substantially, and form a well-developed precursor region ahead of a subshock. Dissipation of Poynting flux occurs in the precursor as well as in the downstream region through a parametric instability driven by the superluminal waves. The Poynting flux remaining in the downstream region is carried entirely by the superluminal waves. The downstream plasma is therefore an essentially unmagnetized, relativistically hot plasma with a non-relativistic flow speed, as suggested by observations of pulsar wind nebulae.

We present an extension of the multi-moment advection scheme (Minoshima et al., 2011, J. Comput. Phys.) to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

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 present a new numerical scheme for solving the advection equation and its application to the Vlasov simulation. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. We have developed one- and two-dimensional schemes and show that they provide quite accurate solutions compared to other existing schemes with the same memory usage. The two-dimensional scheme can solve the solid body rotation problem of a gaussian profile with little numerical diffusion. This is a very important property for Vlasov simulations of magnetized plasma. The application of the scheme to the electromagnetic Vlasov simulation of collisionless shock waves is presented as a benchmark test.

Particle acceleration in the inner ~ 200 pc of the Galaxy is discussed, where diffuse TeV gamma-rays have been detected by the High Energy Stereoscopic System (HESS) observation. The diffuse gamma-ray emission has a strong correlation with molecular clouds with large velocity dispersion, indicating the presence of turbulence. It is argued that the turbulence may contribute to the acceleration of cosmic rays via stochastic acceleration. The stochastic acceleration may energize cosmic-ray protons up to ~ 100 TeV and electrons to ~ 1 TeV in a relatively tenuous medium. The difference in the efficiency between protons and electrons supports the hadronic scenario of the diffuse TeV gamma-ray emission.

We have compared the TeV gamma-rays with the new 12CO J=2-1 data toward HESS J1745-303 in the Galactic center and confirmed that the molecular gas MG358.9-0.5 toward (l, b)=(358.9, -0.5 at VLSR=-100-0 km s-1 shows a reasonable positional agreement with the primary peak (northern part) of the gamma-ray source. For the southern part of HESS J1745-303, we see no CO counterpart, whereas the HI gas in the Parkes 21 cm HI dataset shows a possible counterpart to the gamma-ray source. This HI gas may be optically thick as supported by the HI line shape similar to the optically thick 12CO. We estimate the total mass of interstellar protons including both the molecular and atomic gas to be 2x10^6 Mo and the cosmic-ray proton energy to be 6x10^{48} ergs in the hadronic scenario. We discuss possible origins of the cosmic-ray protons including the nearby SNR G359.1-0.5. The SNR may be able to explain the northern gamma-ray source but the southern source seems to be too far to be energized by the SNR. As an alternative, we argue that the second-order Fermi acceleration in the inter-clump space surrounded by randomly moving high-velocity clumps may offer a possible mechanism to accelerate protons. The large turbulent motion with velocity dispersion of ~15 km s-1 has energy density two orders of magnitude higher than in the solar vicinity and is viable as the energy source.

We present a new numerical scheme for solving the advection equation and its application to Vlasov simulations. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, for better conservation of the information entropy. We have developed one- and two-dimensional schemes and show that they provide quite accurate solutions within reasonable usage of computational resources compared to other existing schemes. The two-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile for more than hundred rotation periods with little numerical diffusion. This is crucially important for Vlasov simulations of magnetized plasmas. Applications of the one- and two-dimensional schemes to electrostatic and electromagnetic Vlasov simulations are presented with some benchmark tests.

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.

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.

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.