We calculate a new equation of state for baryons at sub-nuclear densities for the use in core-collapse simulations of massive stars. The formulation is the nuclear statistical equilibrium description and the liquid drop approximation of nuclei. The model free energy to minimize is calculated by relativistic mean field theory for nucleons and the mass formula for nuclei with atomic number up to ~ 1000. We have also taken into account the pasta phase. We find that the free energy and other thermodynamical quantities are not very different from those given in the standard EOSs that adopt the single nucleus approximation. On the other hand, the average mass is systematically different, which may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores.

Using three-dimensional (special relativistic) magnetohydrodynamics simulations, the amplification of magnetic field behind strong shock wave is studied. In supernova remnants and gamma-ray bursts, strong shock waves propagate through an inhomogeneous density field. When the shock wave hit a density bump or density dent, the Richtmyer-Meshkov instability is induced that cause a deformation of the shock front. The deformed shock leaves vorticity behind the shock wave that amplifies the magnetic field due to the stretching of field lines.

We investigate a proto-neutron star kick velocity estimated from kinetic momentum of a flow around the proto-neutron star after the standing accretion shock instability grows. In this study, ten different types of random perturbations are imposed on the initial flow for each neutrino luminosity. We found that the kick velocities of proto-neutron star are widely distributed from 40 km s−1 to 180 km s−1 when the shock wave reaches 2000 km away from the center of the star. The average value of kick velocity is 115 km s−1, whose value is smaller than the observational ones. The kick velocities do not depend on the neutrino luminosity.

We investigate the effects of neutrino-nucleus interactions on the production of Fluorine during normal supernovae and hypernovae, and discuss stellar mass, metallicity and explosion energy dependence of [F/Fe,Ne,O]. We find the clear trend of [F/Fe,O,Ne] with stellar mass and explosion energy, while no clear trend with metallicity. This trend of [F/O] can be used to constrain the contributed stellar mass by comparing with the observational abundance.

We investigate the lowest mass stars that produce Type-II supernovae, motivated by recent results showing that a large fraction of type-II supernova progenitors for which there are direct detections display unexpectedly low luminosity (for a review see e.g. Smartt 2009). There are three potential evolutionary channels leading to this fate. Alongside the standard ‘massive star’ Fe-core collapse scenario we investigate the likelihood of electron capture supernovae (EC-SNe) from super-AGB (S-AGB) stars in their thermal pulse phase, from failed massive stars for which neon burning and other advanced burning stages fail to prevent the star from contracting to the critical densities required to initiate rapid electron-capture reactions and thus the star's collapse. We find it indeed possible that both of these relatively exotic evolutionary channels may be realised but it is currently unclear for what proportion of stars. Ultimately, the supernova light curves, explosion energies, remnant properties (see e.g. Knigge et al. 2011) and ejecta composition are the quantities desired to establish the role that these stars at the lower edge of the massive star mass range play.

We have obtained optical integral field spectroscopy of the explosion sites of more than 25 nearby type-IIP/IIL/Ib/Ic supernovae using UH88/SNIFS, and additionally Gemini/GMOS IFU. This technique enables us to obtain both spatial and spectral information of the immediate environment of the supernovae. Using strong line method we measured the metallicity of the star cluster present at the explosion site, presumably the coeval parent stellar population of the supernova progenitor, and comparison with simple stellar population models gives age estimate of the cluster. With this method we were able to put constraints on the metallicity and age of the progenitor star. The age, i.e. lifetime, of the progenitor corresponds to the initial mass of the star. By far this is the most direct measurement of supernova progenitor metallicity and, if the cluster-progenitor association is confirmed, provides reliable determination of the initial mass of supernova progenitor stars.

The Gemini Observatories primarily operate a multi-instrument queue, with observers selecting observations that are best suited to weather and seeing conditions. The Target of Opportunity (ToO) observing mode is intended to allow observation of targets that cannot be specified in advance but which have a well defined external trigger such as distant supernovae or Gamma Ray bursts. In addition, the instrument and configuration best suited to observe the ToO may depend on properties of the event, such as brightness and redshift which again are impossible to know in advance. Queue observing naturally lends itself to Target of Opportunity (ToO) support since the time required to switch between programs and instruments is very short, and the staff observer is trained to operate all the available instruments and modes. Gemini Observatory has supported pre-approved ToO programs since beginning queue operations, and has implemented a rapid (less than 15 minutes response time) ToO mode since 2005. ToOs comprise a significant fraction of the queue (20–25% of the highest ranking band) nowadays. We discuss the ToO procedures, the statistics of rapid ToOs observing at Gemini North Observatory, the science related to GRBs and supernovae that this important mode has enabled.

Understanding the intrinsic cosmic long gamma-ray burst (GRB) rate is essential in many aspects of astrophysics and cosmology, such as revealing the connection between GRBs, supernovae (SNe), and stellar evolution. Swift, a multi-wavelength space telescope, is quickly expanding the GRB category by observing hundreds of GRBs and their redshifts. However, it remains difficult to determine the intrinsic GRB rate due to the complex trigger algorithm adopted by Swift. Current studies of the GRB rate usually approximate the Swift trigger algorithm by a single detection threshold. Nevertheless, unlike the previously flown GRB instruments, Swift has over 500 trigger criteria based on count rates and additional thresholds for localization. To investigate possible systematic biases and further explore the intrinsic GRB rate as a function of redshift and the GRB luminosity function, we adopt a Monte Carlo approach by simulating all trigger criteria used by Swift. A precise estimation of the intrinsic GRB rate is important to reveal the GRB origins and their relation to the black-hole forming SNe. Additionally, the GRB rate at high redshifts provides a strong probe of the star formation history in the early universe, which is hard to measure directly through other methods.

The Ultra-Fast Flash Observatory (UFFO) is a space mission to detect the early moments of an explosion from Gamma-ray bursts (GRBs), thus enhancing our understanding of the GRB mechanism. It consists of the UFFO Burst & Trigger telescope (UBAT) for the recognition of GRB positions using hard X-ray from GRBs. It also contains the Slewing Mirror Telescope (SMT) for the fast detection of UV-optical photons from GRBs. It is designed to begin the UV-optical observations in less than a few seconds after the trigger. The UBAT is based on a coded-mask X-ray camera with a wide field of view (FOV) and is composed of the coded mask, a hopper and a detector module. The SMT has a fast rotatable mirror which allows a fast UV-optical detection after the trigger. The telescope is a modified Ritchey-Chrétien telescope with the aperture size of 10 cm diameter, and an image intensifier readout by CCD. The UFFO pathfinder is scheduled to launch into orbit on 2012 June by the Lomonosov spacecraft. It is a scaled-down version of UFFO in order to make the first systematic study of early UV/optical light curves, including the rise phase of GRBs. We expect UBAT to trigger ~44 GRBs/yr and expect SMT to detect ~10 GRBs/yr.

By monitoring WC7 and the O5.5 binary WR 140 with the Suzaku telescope, we demonstrate a new method to measure the mass loss rates of both stars. By using the absorption column density, we found a mass-loss rate for the WC7 component: ṀWC7 ≈ 1.2 × 10−5 M⊙ yr−1. We also measured the mass-loss rate of the companion O component using a luminosity variation in phases: ṀO5.5 ≈ 5 × 10−7 M⊙ yr−1.

Long gamma-ray bursts (GRBs) can be linked to the massive stars and their host galaxies are assumed to be the star-forming galaxies within small dark matter halos. We apply a galaxy evolution model, in which the star formation process inside the virialized dark matter halo at a given redshift is achieved. The star formation rates (SFRs) in the GRB host galaxies at different redshifts can be derived from our model. The related stellar masses, luminosities, and metalicities of these GRB host galaxies are estimated. We also calculate the X-ray and optical absorption of GRB afterglow emission. At higher redshift, the SFR of host galaxy is stronger, and the absorption in the X-ray and optical bands of GRB afterglow is stronger, when the dust and metal components are locally released, surrounding the GRB environment. These model predictions are compared with some observational data as well.

Photospheric thermal radiation components from gamma-ray burst (GRB) jets are estimated based on relativistic hydrodynamic simulations of jet propagation. The light curves and spectra are derived, considering viewing angle effects. The light curves exhibit several seconds time variability and the luminosity is as large as that of GRB prompt emission. For observers at a viewing angle of several degrees the spectra below the peak energy are much softer than that of Planck distribution and close to typical GRB spectrum. Whereas the spectra for observers at small viewing angle are hard and close to Planck distribution. Numerical Amati and Yonetoku relations are reproduced.

We present results of the simulation of a magneto-rotational supernova explosion. We show that, due to the differential rotation of the collapsing iron core, the magnetic field increases with time. The magnetic field transfers angular momentum and a MHD shock wave forms. This shock wave produces the supernova explosion. The explosion energy computed in our simulations is 0.5-2.5 ċ 1051erg. We used two different equations of state for the simulations. The results are rather similar.

We investigated influence of stellar oscillations on the electrodynamics of pulsars as well as magnetars magnetosphere. Besides finding noticeable modification of electromagnetic field and charge density in the polar cap vicinity of oscillating neutron stars we proposed qualitative hypotheses explaining phenomena of part time pulsars as well as sporadic radio emission from generally radio-quiet magnetars with the help of stellar oscillations.

The Carnegie Supernova Project (CSP) has performed, during the period 2004-2009, the optical and NIR follow up of 253 supernovae (SNe) of all types. Among those, 124 were core collapse events, comprising 93 SNe of type II and 31 of types Ib/Ic/IIb. Our follow up consisted of photometric observations suitable to build detailed light curves and a considerable amount of optical spectroscopy.

The bulk of our observations is carried out at Las Campanas Observatory, while access to other facilities is also provided thanks to our strong collaboration with the Millennium Center for Supernova Studies (MCSS).

Our spectroscopic observations were primarily aimed at typing possible new SNe, and follow-up the evolution of CSP targets. One of the goals of the follow-up of type II SNe is the application of independent distance indicators such as the Standard Candle (SCM) and the Expanding Photosphere (EPM) methods. Moreover, through the study of the spectroscopic evolution of these objects, from as early as possible after explosion to the nebular phases, we hope to contribute to their further understanding. Specific analysis of particular objects is underway by members of the CSP and an extended collaboration.

We investigate the propagation of accretion-powered jets in various types of progenitor candidates of GRBs. We perform two dimensional axisymmetric simulations of relativistic hydrodynamics taking into account both the envelope collapse and the jet propagation. In our simulations, the accretion rate is estimated by the mass flux going through the inner boundary, and the jet is injected with a constant accretion-to-jet conversion efficiency η. By varying the efficiency η and opening angle θop for more than 30 models, we find that the jet can make a relativistic breakout from all types of progenitors for GRBs if a simple condition η ≳ 10−3 (θop/20°)2 is satisfied, that is consistent with analytical estimates, otherwise no explosion or some failed spherical explosions occur.

We have investigated the revival of a shock wave by nuclear burning reactions at the central region of core-collapse supernovae. For this purpose, we performed hydrodynamic simulations of core collapse and bounce for 15 M⊙ progenitor model, using ZEUS-MP code in axi-symmetric coordinates. Our numerical code is equipped with a simple nuclear reaction network including 13 α nuclei form 4He to 56Ni, and accounting for energy feedback from nuclear reactions as well as neutrino heating and cooling. We found that the energy released by nuclear reactions is significantly helpful in accelerating shock waves and is able to produce energetic explosion even if the input neutrino luminosity is low.

Some supernovae and gamma-ray bursts are thought to accompany a black hole formation. In the process of a black hole formation, a central core becomes hot and dense enough for hyperons and quarks to appear. In this study, we perform neutrino-radiation hydrodynamical simulations of a stellar core collapse and black hole formation taking into account such exotic components. In our computation, general relativity is fully considered under spherical symmetry. As a result, we find that the additional degrees of freedom soften the equation of state of matter and promote the black hole formation. Furthermore, their effects are detectable as a neutrino signal. We believe that the properties of hot and dense matter at extreme conditions are essential for the studies on the astrophysical black hole formation. This study will be hopefully a first step toward a physics of the central engine of gamma-ray bursts.

Some theoretical studies on the origin of long gamma-ray bursts (GRBs) using stellar evolution models suggest that a low metallicity environment may be a necessary condition for a GRB to occur. However, recent discoveries of high-metallicity host galaxies of some GRBs cast doubt on the requirement of low-metallicity in GRB occurrence. In this study, we predict the metallicity distribution of GRB host galaxies, assuming empirical formulations of galaxy properties. We take internal dispersion of metallicity within each galaxy into account. Assuming GRBs trace low-metallicity star formation 12+log(O/H) < 8.2, we find that ≳ 10% of GRB host galaxies may have Z > Z⊙, depending on the internal dispersion of metallicity within galaxies.

Recent PIC simulations of relativistic electron-positron (electron-ion) jets injected into a stationary medium show that particle acceleration occurs in the shocked regions. Simulations show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields and for particle acceleration. These magnetic fields contribute to the electron's transverse deflection behind the shock. The “jitter” radiation from deflected electrons in turbulent magnetic fields has properties different from synchrotron radiation calculated in a uniform magnetic field. This jitter radiation may be important for understanding the complex time evolution and/or spectral structure of gamma-ray bursts, relativistic jets in general, and supernova remnants. In order to calculate radiation from first principles and go beyond the standard synchrotron model, we have used PIC simulations. We present synthetic spectra to compare with the spectra obtained from Fermi observations.