Relativistic unmagnetized cloud-plasma interaction is analyzed by performing linear analysis and particle-in-cell simulation. This course consists of an electron-ion cloud injected into a stationary ambient plasma and has long been a favorite topic in laboratory and space plasmas. An oblique electromagnetic instability dominates the unstable spectrum. In the interaction with the generated electromagnetic fields, the cloud electrons are entirely mixed with the ambient ones and form a hot electron population. The velocity of the cloud ions, however, has not changed significantly from the initial bulk velocity. As this ion cloud propagates into the plasma, it derives an electrostatic field which can accelerate the electrons up to energy equipartition between electrons and ions. The electrostatic field is amplified at the expense of the kinetic energy of ions, and its spatial scale is in order of the electron skin depth. The electron acceleration in such an electrostatic field is, therefore, a likely process for pre-acceleration of electrons in unmagnetized plasmas.

We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ~10 Mo and a highly variable luminosity of 10^{38}-10^{39} erg/s, comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ~10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10^{-3} pc, which is ~100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation.

Direct collapse within dark matter (DM) halos is a promising path to form supermassive black hole (SMBH) seeds at high redshifts. The outer part of this collapse remains optically thin, and has been studied intensively using numerical simulations. However, the innermost region of the collapse is expected to become optically thick and requires us to follow the radiation field in order to understand its subsequent evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation, for isolated halos. For direct collapse within isolated DM halos, we find that (1) the photosphere forms at ~10^{-6} pc and rapidly expands outward. (2) A central core forms, with a mass of ~1 Mo, supported by thermal gas pressure gradients and rotation. (3) Growing thermal gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow, and another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere encompassing these cores is ~5 x 10^{37}-5 x 10^{38} erg/s, of order the Eddington luminosity. (6) Two variability timescales are associated with this process: a long one, which is related to the accretion flow within the central ~10^{-4}-10^{-3} pc, and ~0.1 yr, which is related to radiation diffusion. (7) Adiabatic models have been run for comparison and their evolution differs profoundly from that of the FLD models, by forming a central geometrically-thick disk. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, mainly because the radiation in the FLD is capable of escaping due to anisotropy in the optical depth and associated gradients.

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 $\sim1\%$ of electrons and $\sim8\%$ of electron energy. Its power-law index is -2.6. The acceleration efficiency is $\sim23\%$ by number and $\sim50\%$ 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.

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 $\epsilon_{\rm ele}\approx \sqrt{\epsilon_{\rm b}}$, where $\epsilon$ 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.