We present new two-dimensional (2D) axisymmetric neutrino radiation/hydrodynamic models of core-collapse supernova (CCSN) cores. We use the CASTRO code, which incorporates truly multi-dimensional, multi-group, flux-limited diffusion (MGFLD) neutrino
transport, including all relevant $mathcal{O}(v/c)$ terms. Our main motivation for carrying out this study is to compare with recent 2D models produced by other groups who have obtained explosions for some progenitor stars and with recent 2D VULCAN results that did not incorporate $mathcal{O}(v/c)$ terms. We follow the evolution of 12, 15, 20, and 25 solar-mass progenitors to approximately 600 milliseconds after bounce and do not obtain an explosion in any of these models. Though the reason for the qualitative disagreement among the groups engaged in CCSN modeling remains unclear, we speculate that the simplifying ``ray-by-ray approach employed by all other groups may be compromising their results. We show that ``ray-by-ray calculations greatly exaggerate the angular and temporal variations of the neutrino fluxes, which we argue are better captured by our multi-dimensional MGFLD approach. On the other hand, our 2D models also make approximations, making it difficult to draw definitive conclusions concerning the root of the differences between groups. We discuss some of the diagnostics often employed in the analyses of CCSN simulations and highlight the intimate relationship between the various explosion conditions that have been proposed. Finally, we explore the ingredients that may be missing in current calculations that may be important in reproducing the properties of the average CCSNe, should the delayed neutrino-heating mechanism be the correct mechanism of explosion.
Starting as highly relativistic collimated jets, gamma-ray burst outflows gradually decelerate and become non-relativistic spherical blast waves. Although detailed analytical solutions describing the afterglow emission received by an on-axis observer
during both the early and late phases of the outflow evolution exist, a calculation of the received flux during the intermediate phase and for an off-axis observer requires either a more simplified analytical model or direct numerical simulations of the outflow dynamics. In this paper we present light curves for off-axis observers covering the long-term evolution of the blast wave calculated from a high resolution two-dimensional relativistic hydrodynamics simulation using a synchrotron radiation model. We compare our results to earlier analytical work and calculate the consequence of the observer angle with respect to the jet axis both for the detection of orphan afterglows and for jet break fits to the observational data. We find that observable jet breaks can be delayed for up to several weeks for off-axis observers, potentially leading to overestimation of the beaming corrected total energy. When using our off-axis light curves to create synthetic Swift X-ray data, we find that jet breaks are likely to remain hidden in the data. We also confirm earlier results in the literature finding that only a very small number of local Type Ibc supernovae can harbor an orphan afterglow.
Direct multi-dimensional numerical simulation is the most reliable approach for calculating the fluid dynamics and observational signatures of relativistic jets in gamma-ray bursts (GRBs). We present a two-dimensional relativistic hydrodynamic simula
tion of a GRB outflow during the afterglow phase, which uses the fifth-order weighted essentially non-oscillatory scheme and adaptive mesh refinement. Initially, the jet has a Lorentz factor of 20. We have followed its evolution up to 150 years. Using the hydrodynamic data, we calculate synchrotron radiation based upon standard afterglow models and compare our results with previous analytic work. We find that the sideways expansion of a relativistic GRB jet is a very slow process and previous analytic works have overestimated its rate. In our computed lightcurves, a very sharp jet break is seen and the post-break lightcurves are steeper than analytic predictions. We find that the jet break in GRB afterglow lightcurves is mainly caused by the missing flux when the edge of the jet is observed. The outflow becomes nonrelativistic at the end of the Blandford-McKee phase. But it is still highly nonspherical, and it takes a rather long time for it to become a spherical Sedov-von Neumann-Taylor blast wave. We find that the late-time afterglows become increasingly flatter over time. But we disagree with the common notion that there is a sudden flattening in lightcurves due to the transition into the Sedov-von Neumann-Taylor solution. We have also found that there is a bump in lightcurves at very late times ($sim 1000$ days) due to radiation from the counter jet. We speculate that such a counter jet bump might have already been observed in GRB 980703.
Magnetic field strengths inferred for relativistic outflows including gamma-ray bursts (GRB) and active galactic nuclei (AGN) are larger than naively expected by orders of magnitude. We present three-dimensional relativistic magnetohydrodynamics (MHD
) simulations demonstrating amplification and saturation of magnetic field by a macroscopic turbulent dynamo triggered by the Kelvin-Helmholtz shear instability. We find rapid growth of electromagnetic energy due to the stretching and folding of field lines in the turbulent velocity field resulting from non-linear development of the instability. Using conditions relevant for GRB internal shocks and late phases of GRB afterglow, we obtain amplification of the electromagnetic energy fraction to $epsilon_B sim 5 times 10^{-3}$. This value decays slowly after the shear is dissipated and appears to be largely independent of the initial field strength. The conditions required for operation of the dynamo are the presence of velocity shear and some seed magnetization both of which are expected to be commonplace. We also find that the turbulent kinetic energy spectrum for the case studied obeys Kolmogorovs 5/3 law and that the electromagnetic energy spectrum is essentially flat with the bulk of the electromagnetic energy at small scales.
