Collisionless shocks can be produced as a result of strong magnetic fields in a plasma flow, and therefore are common in many astrophysical systems. The Weibel instability is one candidate mechanism for the generation of sufficiently strong fields to create a collisionless shock. Despite their crucial role in astrophysical systems, observation of the magnetic fields produced by Weibel instabilities in experiments has been challenging. Using a proton probe to directly image electromagnetic fields, we present evidence of Weibel-generated magnetic fields that grow in opposing, initially unmagnetized plasma flows from laser-driven laboratory experiments. Three-dimensional particle-in-cell simulations reveal that the instability efficiently extracts energy from the plasma flows, and that the self-generated magnetic energy reaches a few percent of the total energy in the system. This result demonstrates an experimental platform suitable for the investigation of a wide range of astrophysical phenomena, including collisionless shock formation in supernova remnants, large-scale magnetic field amplification, and the radiation signature from gamma-ray bursts.

A two-dimensional particle-in-cell simulation is performed to investigate weakly magnetized perpendicular shocks with a magnetization parameter of 6 x 10^-5, which is equivalent to a high Alfv\'en Mach number M_A of ~130. It is shown that current filaments form in the foot region of the shock due to the ion-beam--Weibel instability (or the ion filamentation instability) and that they generate a strong magnetic field there. In the downstream region, these current filaments also generate a tangled magnetic field that is typically 15 times stronger than the upstream magnetic field. The thermal energies of electrons and ions in the downstream region are not in equipartition and their temperature ratio is T_e / T_i ~ 0.3 - 0.4. Efficient electron acceleration was not observed in our simulation, although a fraction of the ions are accelerated slightly on reflection at the shock. The simulation results agree very well with the Rankine-Hugoniot relations. It is also shown that electrons and ions are heated in the foot region by the Buneman instability (for electrons) and the ion-acoustic instability (for both electrons and ions). However, the growth rate of the Buneman instability is significantly reduced due to the relatively high temperature of the reflected ions. For the same reason, ion-ion streaming instability does not grow in the foot region.

A two-dimensional electromagnetic particle-in-cell simulation with the realistic ion-to-electron mass ratio of 1836 is carried out to investigate the electrostatic collisionless shocks in relatively high-speed (~3000 km s^-1) plasma flows and also the influence of both electrostatic and electromagnetic instabilities, which can develop around the shocks, on the shock dynamics. It is shown that the electrostatic ion-ion instability can develop in front of the shocks, where the plasma is under counter-streaming condition, with highly oblique wave vectors as was shown previously. The electrostatic potential generated by the electrostatic ion-ion instability propagating obliquely to the shock surface becomes comparable with the shock potential and finally the shock structure is destroyed. It is also shown that in front of the shock the beam-Weibel instability gradually grows as well, consequently suggesting that the magnetic field generated by the beam-Weibel instability becomes important in long-term evolution of the shock and the Weibel-mediated shock forms long after the electrostatic shock vanished. It is also observed that the secondary electrostatic shock forms in the reflected ions in front of the primary electrostatic shock.

Laboratory spectroscopy of non-thermal equilibrium plasmas photoionized by intense radiation is a key to understanding compact objects, such as black holes, based on astronomical observations. This paper describes an experiment to study photoionizing plasmas in laboratory under well-defined and genuine conditions. Photoionized plasma is here generated using a 0.5-keV Planckian x-ray source created by means of a laser-driven implosion. The measured x-ray spectrum from the photoionized silicon plasma resembles those observed from the binary stars Cygnus X-3 and Vela X-1 with the Chandra x-ray satellite. This demonstrates that an extreme radiation field was produced in the laboratory, however, the theoretical interpretation of the laboratory spectrum significantly contradicts the generally accepted explanations in x-ray astronomy. This model experiment offers a novel test bed for validation and verification of computational codes used in x-ray astronomy.

Supercritical accretion flows inevitably produce radiation-pressure driven outflows, which will Compton up-scatter soft photons from the underlying accretion flow, thereby making hard emission. We perform two dimensional radiation hydrodynamic simulations of supercritical accretion flows and outflows, incorporating such Compton scattering effects, and demonstrate that there appears a new hard spectral state at higher photon luminosities than that of the slim-disk state. In this state, as the photon luminosity increases, the photon index decreases and the fraction of the hard emission increases. The Compton $y$-parameter is of the order of unity (and thus the photon index will be $\sim 2$) when the apparent photon luminosity is ${\sim}30L_{\rm E}$ (with $L_{\rm E}$ being the Eddington luminosity) for nearly face-on sources. This explains the observed spectral hardening of the ULX NGC1313 X-2 in its brightening phase and thus supports the model of supercritical accretion onto stellar mass black holes in this ULX.

We show that the Weibel-mediated collisionless shocks are driven at non-relativistic propagation speed (0.1c < V < 0.45c) in unmagnetized electron-ion plasmas by performing two-dimensional particle-in-cell simulations. It is shown that the profiles of the number density and the mean velocity in the vicinity of the shock transition region, which are normalized by the respective upstream values, are almost independent of the upstream bulk velocity, i.e., the shock velocity. In particular, the width of the shock transition region is ~100 ion inertial length independent of the shock velocity. For these shocks the energy density of the magnetic field generated by the Weibel-type instability within the shock transition region reaches typically 1-2% of the upstream bulk kinetic energy density. This mechanism probably explains the robust formation of collisionless shocks, for example, driven by young supernova remnants, with no assumption of external magnetic field in the universe.

We investigate numerically the hydrodynamic instability of an ionization front (IF) accelerating into a molecular cloud, with imposed initial perturbations of different amplitudes. When the initial amplitude is small, the imposed perturbation is completely stabilized and does not grow. When the initial perturbation amplitude is large enough, roughly the ratio of the initial amplitude to wavelength is greater than 0.02, portions of the IF temporarily separate from the molecular cloud surface, locally decreasing the ablation pressure. This causes the appearance of a large, warm HI region and triggers nonlinear dynamics of the IF. The local difference of the ablation pressure and acceleration enhances the appearance and growth of a multimode perturbation. The stabilization usually seen at the IF in the linear regimes does not work due to the mismatch of the modes of the perturbations at the cloud surface and in density in HII region above the cloud surface. Molecular pillars are observed in the late stages of the large amplitude perturbation case. The velocity gradient in the pillars is in reasonably good agreement with that observed in the Eagle Nebula. The initial perturbation is imposed in three different ways: in density, in incident photon number flux, and in the surface shape. All cases show both stabilization for a small initial perturbation and large growth of the second harmonic by increasing amplitude of the initial perturbation above a critical value.

We report the Compton scattering emission from the Poynting flux acceleration of electron- positron plasma simulated by the 2-1/2 dimensional particle-in-cell(PIC) code. We show these and other remarkable properties of Poynting flux acceleration and Compton spectral output, and discuss the agreement with the observed spectra of GRBs and XRFs.

We investigate the stability of differentially rotating proto-neutron stars (PNSs) with a toroidal magnetic field. Stability criteria for nonaxisymmetric MHD instabilities are derived using a local linear analysis. PNSs are expected to have much stronger radial shear in the rotation velocity compared to normal stars. We find that nonaxisymmetric magnetorotational instability (NMRI) with a large azimuthal wavenumber $m$ is dominant over the kink mode ($m=1$) in differentially rotating PNSs. The growth rate of the NMRI is of the order of the angular velocity $\Omega$ which is faster than that of the kink-type instability by several orders of magnitude. The stability criteria are analogous to those of the axisymmetric magnetorotational instability with a poloidal field, although the effects of leptonic gradients are considered in our analysis. The NMRI can grow even in convectively stable layers if the wavevectors of unstable modes are parallel to the restoring force by the Brunt-V\"ais\"al\"a oscillation. The nonlinear evolution of NMRI could amplify the magnetic fields and drive MHD turbulence in PNSs, which may lead to enhancement of the neutrino luminosity.

We investigate hydrodynamic instability of accelerating ionization fronts using two dimensional hydrodynamic simulations that include detailed energy deposition and release due to the absorption of UV radiation, recombination of hydrogen, radiative molecular cooling, and magnetic pressure. We consider linear perturbation growth and find that the stabilization mechanism associated with non-accelerated fronts remains a significant factor even when acceleration is present. Conversely, if recombination in the ionized region is turned off,Rayleigh-Taylor (RT) instability becomes effective, and the classical RT growth rate recovered.

We investigate the dynamics and morphology of jets propagating into the interstellar medium using 2D relativistic hydrodynamic simulations. The calculations are performed assuming axisymmetric geometry and trace a long propagation of jets. The jets are assumed to be 'light' with the density ratio between the beam to the ambient gas much less than unity. We examine the mechanism for the appearance of vortices at the head of jets in the hot spot. Such vortices are known as a trigger of a deceleration phase which appears after a short phase in which one dimensional analysis is dominant. We find that an oblique shock at the boundary rim near the end of the beam strongly affects the flow structure in and around the hot spot. Weakly shocked gas passed through this oblique shock and becomes a trigger for the generation of vortices. We also find the parameter dependence of these effects for the propagation and dynamics of the jets. The jet with slower propagation velocity is weakly pinched, has large vortices and shows very complex structure at the head of the jets and extended synchrotron emissivity.

We perform 2-dimensional relativistic hydrodynamical simulations in the context of collapsar model. Calculations of explosive nucleosynthesis are also accomplished. We investigate the influence of the structure of the progenitor and energy deposition rate on the resulting explosive nucleosynthesis, assuming that 56Ni is mainly synthesized in the jet launched by the neutrino heating. We show the amount of 56Ni is very sensitive to the energy deposition rate. Thus we conclude that it is quite natural not to detect an underlying supernova in some X-ray afterglows as in GRB 010921. We also point out the possibility that the relative abundance of the elements with intermediate mass number such as Si and S in the X-ray afterglow of GRB 011211 may be naturally explained if the energy deposition rate at the central engine is relatively long because little amount of 56Ni should be synthesized under such an environment. If this discussion is true, there should be correlation between the line features in the X-ray afterglow and duration of the GRB. It should be noted that the duration of GRB 011211 is 270 seconds, making it the longest burst ever observed by Beppo-SAX although it suffers from the effect of red-shift (z_host=2.14), which supports our conclusion. Our results also suggest that the type I collapsar model in which the energy deposition rate is relatively low (\dot{E} \sim 10^{51} erg/s) might have difficulty in reproducing the observed amount of 56Ni in a hypernova such as SN 1998bw. This means that the mechanism of the central engine of a hypernova accompanying GRB may be constrained by the discussion of explosive nucleosynthesis.

In this paper we present a numerical study of plasma jets produced by intense laser matter interactions. Through this study we hope to better understand astrophysical jets and their recent experimental simulations in the laboratory. We paid special attention to radiation cooling and the interaction of the jet with ambient gas. Four cases are presented in this paper; two of them deal with the propagation of jets in vacuum, while in the other two the propagation takes place in the ambient gas. Available experimental results are reproduced to good accuracy in the vacuum case. For jets in ambient gas, we find that the existence of the surrounding gas confines the jet into a narrow cylindrical shape so that both the density and temperature of the jet remain high enough for effective radiation cooling. As a result, a collimated plasma jet is formed in these cases. The dimensionless parameters characterizing the laboratory jets and protostellar jets have overlapping domains. We also discuss the cooling lengths for our model and compare them with the corresponding values in the astrophysical jets. A plasma jet in the ambient gas experiment is proposed which is within the reach of present day technology, and can be relevant to astrophysical phenomena.