No abstract

A study of collisionless external shocks in gamma-ray bursts is presented. The shock structure, electromagnetic fields, and process of electron acceleration are assessed by performing a self-consistent 3D particle-in-cell (PIC) simulation. In accordance with hydrodynamic shock systems, the shock consists of a reverse shock (RS) and forward shock (FS) separated by a contact discontinuity (CD). The development and structure are controlled by the ion Weibel instability. The ion filaments are sources of strong transverse electromagnetic fields at both sides of the double shock structure over a length of 30 -100 ion skin depths. Electrons are heated up to a maximum energy ǫ ele ≈ √ ǫ b , where ǫ is the energy normalized to the total incoming energy. Jet electrons are trapped in the RS transition region due to the presence of an ambipolar electric field and reflection by the strong transverse magnetic fields in the shocked region. In a process similar to shock surfing acceleration (SSA) for ions, electrons experience drift motion and acceleration by ion filament transverse electric fields in the plane perpendicular to the shock propagation direction. Ultimately accelerated jet electrons are convected back into the upstream.

The analysis of plasma beam and shear problems in magnetic fields is usually based on a hydromagnetic fluid model. In a low-density collisionless plasma, however, the kinetic effects of the plasma, such as finite Larmor radius effects, are not yet clearly understood. In this paper, the kinetic equilibria of plasma shears in a uniform and fixed magnetic field, with full ion motion, are discussed by solving the Vlasov equation with a given electric field and drift velocity. In this model, the ion density profile through the plasma shear layer is quite different from the one predicted by a hydromagnetic model. As a result of a complicated ion gyromotion through the shear layer, single- and double-humped ion density profiles are obtained. The dependence on the temperature and the strength of the shear will be discussed. The results show a significant difference between positive and negative shears.

We use the nonstandard-finite-difference time-domain (NS-FDTD) method to investigate the interaction of light with the complicated microstructures in the Morpho butterfly scales, which produce the well-known brilliant blue coloring. The NS-FDTD algorithm is particularly suitable to analyze such complex structures because the calculation can be performed in a short time with high accuracy on a relatively coarse numerical grid. We analyze (1) the microstructure obtained directly by binarizing an electron microgram of the cross section of a scale, (2) the reflection and diffraction properties of three model structures--flat, alternating, and tree-shaped alternating multilayers, and (3) an array of alternating multilayers with random noise superposed on the height of the structures. We found that the actual microstructure well reproduced the reflection spectrum in a blue region by integrating the reflection intensities over all the reflection angles. Under normal incidence, the flat multilayer mainly stresses on multilayer interference except for shorter wavelengths, while alternating multilayer rather enhances the effect of diffraction grating due to longitudinally repeating structure by strongly suppressing the reflection toward the normal direction. In the array of alternating multilayers, the reflection into larger angles is considerably suppressed and the spectral shape becomes different from that expected for a single alternating multilayer. This suppression mainly comes from the scattering of reflected light by adjacent structures, which is particularly prominent for the TM mode. Thus a clear difference between the TE and TM modes is observed with respect to the origin of spectral shape, though the obtained spectra are similar to each other. Finally, the polarization dependence of the reflection and the importance of the alternating multilayer are discussed.

The course of non-thermal electron ejection in relativistic unmagnetized electron-ion shocks is investigated by performing self-consistent particle-in-cell simulations. The shocks are excited through the injection of relativistic jet into ambient plasma, leading to two distinct shocks (named as the trailing shock and leading shock) and a contact discontinuity. The Weibel-like instabilities heat the electrons up to approximately half of ion kinetic energy. The double layers formed in the trailing and leading edges then accelerated the electrons by the ion kinetic energy. The electron distribution function in the leading edge shows a clear non-thermal power-law tail which contains ∼ 1% of electrons and ∼ 8% of electron energy. Its power-law index is -2.6. The acceleration efficiency is ∼ 23% by number and ∼ 50% by energy and the power-law index is -1.8 for electron distribution function in the trailing edge. The effect of the dimensionality is examined by comparing results of 3D simulation with 2D ones. It exhibits that the electron acceleration is more efficient in 2D.