A shock tube problem is solved numerically by using one-dimensional full particle-in-cell simulations under the condition that a relatively tenuous and weakly magnetized plasma is continuously pushed by a relatively dense and strongly magnetized plasma having supersonic relative velocity. A forward and a reverse shock and a contact discontinuity are self-consistently reproduced. The spatial width of the contact discontinuity increases as the angle between the discontinuity normal and ambient magnetic field decreases. The inner structure of the discontinuity shows different profiles between magnetic field and plasma density, or pressure, which is caused by a non-MHD effect of the local plasma. The region between the two shocks is turbulent. The fluctuations in the relatively dense plasma are compressible and propagating away from the contact discontinuity, although the fluctuations in the relatively tenuous plasma contain both compressible and incompressible components. The source of the compressible fluctuations in the relatively dense plasma is in the relatively tenuous plasma. Only compressible fast mode fluctuations generated in the relatively tenuous plasma are transmitted through the contact discontinuity and propagate in the relatively dense plasma. These fast mode fluctuations are steepened when passing the contact discontinuity. This wave steepening and probably other effects may cause the broadening of the wave spectrum in the very local interstellar medium plasma. The results are discussed in the context of the heliospheric boundary region or heliopause.

A generation mechanism of super Alfv\'{e}nic (SPA) waves in multi-ion species plasma is proposed and the associated heavy ion acceleration process is discussed. The SPA waves are thought to play important roles in particle acceleration since they have large wave electric field because of their high phase velocity. It is demonstrated by using full particle-in-cell simulation that large amplitude proton cyclotron waves, excited due to proton temperature anisotropy, nonlinearly destabilize SPA waves through parametric decay instability in a three component plasma composed of electrons, protons, and $\alpha$ particles. At the same time, $\alpha$ cyclotron waves get excited via another decay instability. A pre-accelerated $\alpha$ particle resonates simultaneously with the two daughter waves, the SPA waves and the $\alpha$ cyclotron waves, and it is further accelerated perpendicular to the ambient magnetic field. The process may work in astrophysical environments where sufficiently large temperature anisotropy of lower mass ions occurs.

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.

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.

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.

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 review is given for cyclotron resonant interactions in space plasmas. After giving a simple formulation for the test particle approach, illustrative examples for resonant interactions are given. It is shown that for obliquely propagating whistler waves, not only fundamental cyclotron resonance, but also other resonances, such as transit-time resonance, anomalous cyclotron resonance, higher-harmonic cyclotron resonance, and even subharmonic resonance can come into play. A few recent topics of cyclotron resonant interactions, such as electron injection in shocks, cyclotron resonant heating of solar wind heavy ions, and relativistic modifications, are also reviewed.

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.

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.

A new particle acceleration process in a developing Alfv\'{e}n turbulence in the course of successive parametric instabilities of a relativistic pair plasma is investigated by utilyzing one-dimensional electromagnetic full particle code. Coherent wave-particle interactions result in efficient particle acceleration leading to a power-law like energy distribution function. In the simulation high energy particles having large relativistic masses are preferentially accelerated as the turbulence spectrum evolves in time. Main acceleration mechanism is simultaneous relativistic resonance between a particle and two different waves. An analytical expression of maximum attainable energy in such wave-particle interactions is derived.

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.