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.

A preliminary numerical experiment is conducted for laboratory experiments on the generation of magnetized collisionless shocks with high-power lasers by using one-dimensional particle-in-cell simulation. The present study deals with the interaction between a moving Aluminum plasma and a Nitrogen plasma at rest. In the numerical experiment, the Nitrogen plasma is unmagnetized or magnetized by a weak external magnetic field. Since the previous study suggested the generation of spontaneous magnetic field in the piston (Aluminum) plasma due to the Biermann battery, the effect of the magnetic field is of interest. Sharp jumps of electron density and magnetic field are observed around the interface between the two plasmas as long as one of the two plasmas is magnetized, which indicates the formation of tangential electron-magneto-hydro-dynamic discontinuity. When the Aluminum plasma is magnetized, strong compression of both density and magnetic field takes place in the pure Aluminum plasma during the gyration of Nitrogen ions in the Aluminum plasma region. The formation of a shock downstream is indicated from the shock jump condition. The result suggests that the spontaneous magnetic field in the piston (Aluminum) plasma plays an essential role in the formation of a perpendicular collisionless shock.

A simple form of the Boris solver in particle-in-cell (PIC) simulation is proposed. It employs an exact solution of the Lorentz-force part, and it is equivalent to the Boris solver with a gyrophase correction. As a favorable property for stable schemes, this form preserves a volume in the phase space. Numerical tests of the Boris solvers are conducted by test-particle simulations and by PIC simulations. The proposed form provides better accuracy than the popular form, while it only requires few additional computation time.

The structure of the diffusion regions in antiparallel magnetic reconnection is investigated by means of a theory and a Vlasov simulation. The magnetic diffusion is considered as relaxation to the frozen-in state, which depends on a reference velocity field. A field-aligned component of the frozen-in condition is proposed to evaluate a diffusion-like process. Diffusion signatures with respect to ion and electron bulk flows indicate the ion and electron diffusion regions near the reconnection site. The electron diffusion region resembles the energy dissipation region. These results are favorable to a previous expectation that an electron-scale dissipation region is surrounded by an ion-scale Hall-physics region.

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.

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.

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.

A full particle simulation study is carried out on a perpendicular collisionless shock with a relatively low Alfven Mach number (M_A=5). In the present study, we have performed a two-dimensional (2D) electromagnetic full particle simulation with a "shock-rest-frame model". The simulation domain is taken to be larger than the ion inertial length in order to include full kinetics of both electrons and ions. The present simulation result has confirmed the transition of shock structures from the cyclic self-reformation to the quasi-stationary shock front. During the transition, electrons and ions are thermalized in the direction parallel to the shock magnetic field. Ions are thermalized by low-frequency electromagnetic waves (or rippled structures) excited by strong ion temperature anisotropy at the shock foot, while electrons are thermalized by high-frequency electromagnetic waves (or whistler mode waves) excited by electron temperature anisotropy at the shock overshoot. Ion acoustic waves are also excited at the shock overshoot where the electron parallel temperature becomes higher than the ion parallel temperature. We expect that ion acoustic waves are responsible for parallel diffusion of both electrons and ions, and that a cross-scale coupling between an ion-scale mesoscopic instability and an electron-scale microscopic instability is important for structures and dynamics of a collisionless perpendicular shock.

Nonlinear evolution of circularly polarized Alfv\'en waves are discussed by using the recently developed Vlasov-MHD code, which is a generalized Landau-fluid model. The numerical results indicate that as far as the nonlinearity in the system is not so large, the Vlasov-MHD model can validly solve time evolution of the Alfv\'enic turbulence both in the linear and nonlinear stages. The present Vlasov-MHD model is proper to discuss the solar coronal heating and solar wind acceleration by Alfve\'n waves propagating from the photosphere.

A full particle simulation study is carried out on the electron acceleration at a collisionless, relatively low Alfven Mach number (M_A=5), perpendicular shock. Recent self-consistent hybrid shock simulations have demonstrated that the shock front of perpendicular shocks has a dynamic rippled character along the shock surface of low-Mach-number perpendicular shocks. In this paper, the effect of the rippling of perpendicular shocks on the electron acceleration is examined by means of large-scale (ion-scale) two-dimensional full particle simulations. It has been shown that a large-amplitude electric field is excited at the shock front in association with the ion-scale rippling, and that reflected ions are accelerated upstream at a localized region where the shock-normal electric field of the rippled structure is polarized upstream. The current-driven instability caused by the highly-accelerated reflected ions has a high growth rate to large-amplitude electrostatic waves. Energetic electrons are then generated by the large-amplitude electrostatic waves via electron surfing acceleration at the leading edge of the shock transition region. The present result suggests that the electron surfing acceleration is also a common feature at low-Mach-number perpendicular collisionless shocks.

A two-dimensional (2D) shock-rest-frame model for particle simulations is developed. Then full kinetic dynamics of a perpendicular collisionless shock is examined by means of a 2D full particle simulation. We found that in the 2D simulation there are fewer nonthermal electrons due to surfing acceleration which was seen in the previous 1D simulations of a high Mach number perpendicular shock in a low-beta and weakly magnetized plasma. This is because the particle motion along the ambient magnetic field disturbs the formation of coherent electrostatic solitary structures which is necessary for electron surfing acceleration.

The full kinetic dynamics of a perpendicular collisionless shock is studied by means of a one-dimensional electromagnetic full particle simulation. The present simulation domain is taken in the shock rest frame in contrast to the previous full particle simulations of shocks. Preliminary results show that the downstream state falls into a unique cyclic reformation state for a given set of upstream parameters through the self-consistent kinetic processes.