Search

The interaction between a thin foil target and a circularly polarized laser light injected along an external magnetic field is investigated numerically by particle-in-cell simulations. A standing wave appears at the front surface of the target, overlapping the injected and partially reflected waves. Hot electrons are efficiently generated at the standing wave due to the relativistic two-wave resonant acceleration if the magnetic field amplitude of the standing wave is larger than the ambient field. A bifurcation occurs in the gyration motion of electrons, allowing all electrons with non-relativistic velocities to acquire relativistic energy through the cyclotron resonance. The optimal conditions for the highest energy and the most significant fraction of hot electrons are derived precisely through a simple analysis of test-particle trajectories in the standing wave. Since the number of hot electrons increases drastically by many orders of magnitude compared to the conventional unmagnetized cases, this acceleration could be a great advantage in laser-driven ion acceleration and its applications., Comment: 22 pages, 15 figures, 4 tables, accepted for publication in PRE

Fast radio bursts (FRBs) are millisecond transient astrophysical phenomena and bright at radio frequencies. The emission mechanism, however, remains unsolved yet. One scenario is a coherent emission associated with the magnetar flares and resulting relativistic shock waves. Here, we report unprecedentedly large-scale simulations of relativistic magnetized ion-electron shocks, showing that strongly linear-polarized electromagnetic waves are excited. The kinetic energy conversion to the emission is so efficient that the wave amplitude is responsible for the brightness. We also find a polarization angle swing reflecting shock front modulation, implicating the polarization property of some repeating FRBs. The results support the shock scenario as an origin of the FRBs., Comment: accepted to PRL

Magnetic reconnection is an important source of energetic particles in systems ranging from astrophysics to the laboratory. The large separation of spatiotemporal scales involved makes it critical to determine the minimum physical model containing the necessary physics for modeling particle acceleration. By resolving the energy gain from ideal and nonideal magnetohydrodynamic electric fields self-consistently in kinetic particle-in-cell simulations of reconnection, we conclusively show the dominant role of the nonideal field for the early stage of energization known as injection. The importance of the nonideal field increases with magnetization, guide field, and in three-dimensions, indicating its general importance for reconnection in natural astrophysical systems. We obtain the statistical properties of the injection process from the simulations, paving the way for the development of extended MHD models capable of accurately modeling particle acceleration in large-scale systems. The novel analysis method developed in this study can be applied broadly to give new insight into a wide range of processes in plasma physics.

Particle acceleration at magnetized purely perpendicular relativistic shocks in electron-ion plasmas are studied by means of two-dimensional particle-in-cell simulations. Magnetized shocks with the upstream bulk Lorentz factor $\gamma_1 \gg 1$ are known to emit intense electromagnetic waves from the shock front, which induce electrostatic plasma waves (wakefield) and transverse filamentary structures in the upstream region via the stimulated/induced Raman scattering and the filamentation instability, respectively. The wakefield and filaments inject a fraction of incoming particles into a particle acceleration process, in which particles are once decoupled from the upstream bulk flow by the wakefield, and are piked up again by the flow. The picked-up particles are accelerated by the motional electric field. The maximum attainable Lorentz factor is estimated as $\gamma_{max,e} \sim \alpha\gamma_1^3$ for electrons and $\gamma_{max,i} \sim (1+m_e\gamma_1/m_i)\gamma_1^2$ for ions, where $\alpha \sim 10$ is determined from our simulation results. $\alpha$ can increase up to $\gamma_1$ for weakly magnetized shock if $\gamma_1$ is sufficiently large. This result indicates that highly relativistic astrophysical shocks such as external shocks of gamma-ray bursts can be an efficient particle accelerator., Comment: 14 pages, 13 figures, accepted to ApJ

Using large-scale fully-kinetic two-dimensional particle-in-cell simulations, we investigate the effects of shock rippling on electron acceleration at low-Mach-number shocks propagating in high-$\beta$ plasmas, in application to merger shocks in galaxy clusters. We find that the electron acceleration rate increases considerably when the rippling modes appear. The main acceleration mechanism is stochastic shock-drift acceleration, in which electrons are confined at the shock by pitch-angle scattering off turbulence and gain energy from the motional electric field. The presence of multi-scale magnetic turbulence at the shock transition and the region immediately behind the main shock overshoot is essential for electron energization. Wide-energy non-thermal electron distributions are formed both upstream and downstream of the shock. The maximum energy of the electrons is sufficient for their injection into diffusive shock acceleration. We show for the first time that the downstream electron spectrum has a~power-law form with index $p\approx 2.5$, in agreement with observations., Comment: 15 pages, 14 figures, to be published in ApJ

Particle acceleration and heating at mildly relativistic magnetized shocks in electron-ion plasma are investigated with unprecedentedly high-resolution two-dimensional particle-in-cell simulations that include ion-scale shock rippling. Electrons are super-adiabatically heated at the shock, and most of the energy transfer from protons to electrons takes place at or downstream of the shock. We are the first to demonstrate that shock rippling is crucial for the energization of electrons at the shock. They remain well below equipartition with the protons. The downstream electron spectra are approximately thermal with a limited supra-thermal power-law component. Our results are discussed in the context of wakefield acceleration and the modelling of electromagnetic radiation from blazar cores., Comment: 11 pages, 11 figures; This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Mildly relativistic shocks in magnetized electron-ion plasmas are investigated with 2D kinetic particle-in-cell simulations of unprecedentedly high resolution and large scale for conditions that may be found at internal shocks in blazar cores. Ion-scale effects cause corrugations along the shock surface whose properties somewhat depend on the configuration of the mean perpendicular magnetic field, that is either in or out of the simulation plane. We show that the synchrotron maser instability persists to operate in mildly relativistic shocks in agreement with theoretical predictions and produces coherent emission of upstream-propagating electromagnetic waves. Shock front ripples are excited in both mean-field configurations and they engender effective wave amplification. The interaction of these waves with upstream plasma generates electrostatic wakefields., Comment: 14 pages, 19 figures; This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Microinstabilities and waves excited at moderate-Mach-number perpendicular shocks in the near-Sun solar wind are investigated by full particle-in-cell (PIC) simulations. By analyzing the dispersion relation of fluctuating field components directly issued from the shock simulation, we obtain key findings concerning wave excitations at the shock front: (1) at the leading edge of the foot, two types of electrostatic (ES) waves are observed. The relative drift of the reflected ions versus the electrons triggers an electron cyclotron drift instability (ECDI) which excites the first ES wave. Because the bulk velocity of gyro-reflected ions shifts to the direction of the shock front, the resulting ES wave propagates oblique to the shock normal. Immediately, a fraction of incident electrons are accelerated by this ES wave and a ring-like velocity distribution is generated. They can couple with the hot Maxwellian core and excite the second ES wave around the upper hybrid frequency. (2) from the middle of the foot all the way to the ramp, electrons can couple with both incident and reflected ions. ES waves excited by ECDI in different directions propagate across each other. Electromagnetic (EM) waves (X mode) emitted toward upstream are observed in both regions. They are probably induced by a small fraction of relativistic electrons. Results shed new insight on the mechanism for the occurrence of ES wave excitations and possible EM wave emissions at young CME-driven shocks in the near-Sun solar wind.

Ion temperature anisotropy is a common feature for (quasi-)perpendicular collisionless shocks. By using two-dimensional full particle simulations, it is shown, that the ion temperature component perpendicular to the shock magnetic field at the shock foot region is proportional to the square of the Alfven Mach number divided by the plasma beta. This result is also explained by a simple analytical argument, in which the reflected ions get energy from upstream plasma flow. By comparing our analytic and numerical results, it is also confirmed that the fraction of the reflected ions hardly depends on the plasma beta and the Alfven Mach number when the square of the Alfven Mach number divided by the plasma beta is larger than about 20., Comment: 6 pages, 3 figures, 1 table. AIP Advances, in press

Large-scale two-dimensional (2D) full particle-in-cell simulations are carried out for studying the relationship between the dynamics of a perpendicular shock and microinstabilities generated at the shock foot. The structure and dynamics of collisionless shocks are generally determined by Alfven Mach number and plasma beta, while microinstabilities at the shock foot are controlled by the ratio of the upstream bulk velocity to the electron thermal velocity and the ratio of the plasma-to-cyclotron frequency. With a fixed Alfven Mach number and plasma beta, the ratio of the upstream bulk velocity to the electron thermal velocity is given as a function of the ion-to-electron mass ratio. The present 2D full PIC simulations with a relatively low Alfven Mach number (M_A ~ 6) show that the modified two-stream instability is dominant with higher ion-to-electron mass ratios. It is also confirmed that waves propagating downstream are more enhanced at the shock foot near the shock ramp as the mass ratio becomes higher. The result suggests that these waves play a role in the modification of the dynamics of collisionless shocks through the interaction with shock front ripples., Comment: 12 pages,7 figures, 1 table, Physics of Plasmas, in press; the paper with full resolution images is http://www.phys.aoyama.ac.jp/~ryo/papers/microinsta_PoP2.pdf arXiv admin note: text overlap with arXiv:1204.2539

A full particle simulation study is carried out for studying microinstabilities generated at the shock front of perpendicular collisionless shocks. The structure and dynamics of shock waves are determined by Alfven Mach number and plasma beta, while microinstabilities are controlled by the ratio of the upstream bulk velocity to the electron thermal velocity and the plasma-to-cyclotron frequency. Thus, growth rates of microinstabilities are changed by the ion-to-electron mass ratio, even with the same Mach number and plasma beta. The present two-dimensional simulations show that the electron cyclotron drift instability is dominant for a lower mass ratio, and electrostatic electron cyclotron harmonic waves are excited. For a higher mass ratio, the modified two-stream instability is dominant and oblique electromagnetic whistler waves are excited, which can affect the structure and dynamics of collisionless shocks by modifying shock magnetic fields., Comment: 13 pages, 7 figures, Physics of Plasmas, in press; the paper with full resolution images is http://www.phys.aoyama.ac.jp/~ryo/papers/microinsta_PoP.pdf

An extreme case of electron shock drift acceleration in low Mach number collisionless shocks is investigated as a plausible mechanism of initial acceleration of relativistic electrons in large-scale shocks in galaxy clusters where upstream plasma temperature is of the order of 10 keV and a degree of magnetization is not too small. One-dimensional electromagnetic full particle simulations reveal that, even though a shock is rather moderate, a part of thermal incoming electrons are accelerated and reflected through relativistic shock drift acceleration and form a local nonthermal population just upstream of the shock. The accelerated electrons can self-generate local coherent waves and further be back-scattered toward the shock by those waves. This may be a scenario for the first stage of the electron shock acceleration occurring at the large-scale shocks in galaxy clusters such as CIZA J2242.8+5301 which has well defined radio relics., Comment: 26 pages, 10 figures, accepted for publication in ApJ

Efficiency of electron heating through microinstabilities generated in the transition region of a quasi-perpendicular shock for wide ange of Mach numbers is investigated by utilizing PIC (Particle-In-Cell) simulation and model analyses. In the model analyses saturation levels of effective electron temperature as a result of microinstabilities are estimated from an extended quasilinear (trapping) analysis for relatively low (high) Mach number shocks. Here, MTSI (modified two-stream instability) is assumed to become dominant in low Mach number regime, while BI (Buneman instability) to become dominant in high Mach number regime, respectively. It is revealed that Mach number dependence of the effective electron temperature in the MTSI dominant case is essentially different from that in the BI dominant case. The effective electron temperature through the MTSI does not depend much on the Mach number, although that through the BI increases with the Mach number as in the past studies. The results are confirmed to be consistent with the PIC simulations both in qualitative and quantitative levels. The model analyses predict that a critical Mach number above which steep rise of electron heating rate occurs may arise at the Mach number of a few tens., Comment: 9 pages, 5 figures, Phys. Plasmas in press

We discuss the parametric instabilities of large-amplitude parallel propagating Alfven waves using the 2-D PIC simulation code. First, we confirmed the results in the past study [Sakai et al, 2005] that the electrons are heated due to the modified two stream instability and that the ions are heated by the parallel propagating ion acoustic waves. However, although the past study argued that such parallel propagating longitudinal waves are excited by transverse modulation of parent Alfven wave, we consider these waves are more likely to be generated by the usual, parallel decay instability. Further, we performed other simulation runs with different polarization of the parent Alfven waves or the different ion thermal velocity. Numerical results suggest that the electron heating by the modified two stream instability due to the large amplitude Alfven waves is unimportant with most parameter sets., Comment: submitted to New journal of Physics

The linear dispersion relation of a finite amplitude, parallel, circularly polarized Alfv\'en wave in a relativistic electron-positron plasma is derived. In the nonrelativistic regime, the dispersion relation has two branches, one electromagnetic wave, with a low frequency cutoff at $\sqrt{1+2\omega_p^2/\Omega_p^2}$ (where $\omega_p=(4\pi n e^2/m)^{1/2}$ is the electron/positron plasma frequency), and an Alfv\'en wave, with high frequency cutoff at the positron gyrofrequency $\Omega_p$. There is only one forward propagating mode for a given frequency. However, due to relativistic effects, there is no low frequency cutoff for the electromagnetic branch, and there appears a critical wave number above which the Alfv\'en wave ceases to exist. This critical wave number is given by $ck_c/\Omega_p=a/\eta$, where $a=\omega_p^2/\Omega_p^2$ and $\eta$ is the ratio between the Alfv\'en wave magnetic field amplitude and the background magnetic field. In this case, for each frequency in the Alfv\'en branch, two additional forward propagating modes exist with equal frequency. A simple numerical example is studied: by numerically solving the coupled system of fluid and Maxwell equations, normal incidence of a finite amplitude Alfv\'en wave on an interface between two electron-positron plasmas of different densities is considered., Comment: 12th International Congress on Plasma Physics, 25-29 October 2004, Nice (France)