Progenitors Of Long-duration Gamma-ray Bursts Research Articles

Nuclear burning in collapsar accretion discs

ABSTRACT The core collapse of massive, rapidly-rotating stars are thought to be the progenitors of long-duration gamma-ray bursts (GRB) and their associated hyperenergetic supernovae (SNe). At early times after the collapse, relatively low angular momentum material from the infalling stellar envelope will circularize into an accretion disc located just outside the black hole horizon, resulting in high accretion rates necessary to power a GRB jet. Temperatures in the disc mid-plane at these small radii are sufficiently high to dissociate nuclei, while outflows from the disc can be neutron-rich and may synthesize r-process nuclei. However, at later times, and for high progenitor angular momentum, the outer layers of the stellar envelope can circularize at larger radii ≳ 107 cm, where nuclear reactions can take place in the disc mid-plane (e.g. 4He + 16O → 20Ne + γ). Here we explore the effects of nuclear burning on collapsar accretion discs and their outflows by means of hydrodynamical α-viscosity torus simulations coupled to a 19-isotope nuclear reaction network, which are designed to mimic the late infall epochs in collapsar evolution when the viscous time of the torus has become comparable to the envelope fall-back time. Our results address several key questions, such as the conditions for quiescent burning and accretion versus detonation and the generation of 56Ni in disc outflows, which we show could contribute significantly to powering GRB SNe. Being located in the slowest, innermost layers of the ejecta, the latter could provide the radioactive heating source necessary to make the spectral signatures of r-process elements visible in late-time GRB-SNe spectra.

THE ARDUOUS JOURNEY TO BLACK HOLE FORMATION IN POTENTIAL GAMMA-RAY BURST PROGENITORS

We present a quantitative study on the properties at death of fast-rotating massive stars evolved at low-metallicity, objects that are proposed as likely progenitors of long-duration gamma-ray bursts (LGRBs). We perform 1D+rotation stellar-collapse simulations on the progenitor models of Woosley & Heger (2006) and critically assess their potential for the formation of a black hole and a Keplerian disk (namely a collapsar) or a proto-magnetar. We note that theoretical uncertainties in the treatment of magnetic fields and the approximate handling of rotation compromises the accuracy of stellar-evolution models. We find that only the fastest rotating progenitors achieve sufficient compactness for black-hole formation while the bulk of models possess a core density structure typical of garden-variety core-collapse supernova (SN) progenitors evolved without rotation and at solar metallicity. Of the models that do have sufficient compactness for black-hole formation, most of them also retain a large amount of angular momentum in the core, making them prone to a magneto-rotational explosion, therefore preferentially leaving behind a proto-magnetar. A large progenitor angular-momentum budget is often the sole criterion invoked in the community today to assess the suitability for producing a collapsar. This simplification ignores equally important considerations such as the core compactness, which conditions black-hole formation, the core angular momentum, which may foster a magneto-rotational explosion preventing black-hole formation, or the metallicity and the residual envelope mass which must be compatible with inferences from observed LGRB/SNe. Our study suggests that black-hole formation is non trivial, that there is room for accommodating both collapsars and proto-magnetars as LGRB progenitors, although proto-magnetars seem much more easily produced by current stellar-evolutionary models.

SIMULATIONS OF ACCRETION POWERED SUPERNOVAE IN THE PROGENITORS OF GAMMA-RAY BURSTS

Observational evidence suggests a link between long duration gamma ray bursts (LGRBs) and Type Ic supernovae. Here, we propose a potential mechanism for Type Ic supernovae in LGRB progenitors powered solely by accretion energy. We present spherically-symmetric hydrodynamic simulations of the long-term accretion of a rotating gamma-ray burst progenitor star, a "collapsar," onto the central compact object, which we take to be a black hole. The simulations were carried out with the adaptive mesh refinement code FLASH in one spatial dimension and with rotation, an explicit shear viscosity, and convection in the mixing length theory approximation. Once the accretion flow becomes rotationally supported outside of the black hole, an accretion shock forms and traverses the stellar envelope. Energy is carried from the central geometrically thick accretion disk to the stellar envelope by convection. Energy losses through neutrino emission and nuclear photodisintegration are calculated but do not seem important following the rapid early drop of the accretion rate following circularization. We find that the shock velocity, energy, and unbound mass are sensitive to convective efficiency, effective viscosity, and initial stellar angular momentum. Our simulations show that given the appropriate combinations of stellar and physical parameters, explosions with energies ~5x10^50 ergs, velocities 3000 km/s, and unbound material masses >6 solar masses are possible in a rapidly rotating 16 solar mass main sequence progenitor star. Further work is needed to constrain the values of these parameters, to identify the likely outcomes in more plausible and massive LRGB progenitors, and to explore nucleosynthetic implications.

Metallicity properties of the simulated host galaxies of long gamma-ray bursts and the fundamental metallicity relation

By combining high-resolution N-body simulations with semi-analytical models of galaxy formation, we study the implications of the collapsar model for long-duration gamma-ray bursts (LGRBs) on the metallicity properties of the host galaxies. The cosmological model that we use reproduces the fundamental metallicity relation – the metallicity decreases with increasing star formation rate for galaxies of a given stellar mass. This was recently discovered for the Sloan Digital Sky Survey galaxies. We select host galaxies that house pockets of gas particles, that are young and that have different thresholds for their metallicities; these can be sites of LRGB events, according to the collapsar model. The simulated samples are compared with 18 observed LGRB hosts with the aim of discovering whether the metallicity is a primary parameter. We find that a threshold of metallicity for the LGRB progenitors, within the model galaxies, is not necessary to reproduce the observed distribution of host metallicities. The low metallicities of most LGRB hosts are consistent with the expectation that GRBs trace star formation. The star formation rate appears to be the primary parameter tracing the occurrence of a burst event. Finally, we show that only a few LGRBs are observed in massive, highly extinct galaxies, despite the fact that these galaxies are expected to produce many such events. We identify these missing events with the fraction of dark LGRBs.

On the presence and absence of disks around O-type stars

(abridged) As the favoured progenitors of long-duration gamma-ray bursts, massive stars may represent our best signposts of individual objects in the early Universe, but special conditions seem required to make these bursters, which might originate from the progenitor's rapid rotation and associated asymmetry. To obtain empirical constraints on the interplay between stellar rotation and wind asymmetry, we perform linear Halpha spectropolarimetry on a sample of 18 spectroscopically peculiar massive O stars, including OVz, Of?p, Oe, and Onfp stars, supplemented by an earlier sample of 20 O supergiants. Despite their rapid rotation (with vsin(i) up to 400 km/s) most O-type stars are found to be spherically symmetric, but with notable exceptions amongst specific object classes. We divide the peculiar O stars into four distinct categories: Groups III and IV include the Oe stars and Onfp stars, which are on the high-end tail of the O star rotation distribution and have in the past been claimed to be embedded in disks. Here we report the detection of a classical depolarization ``line effect'' in the Oe star HD 45314, but the overall incidence of line effects amongst Oe stars is significantly lower (1 out of 6) than amongst Be stars. The chance that the Oe and Be datasets are drawn from the same parent population is negligible (with 95% confidence). This implies there is as yet no evidence for a disk hypothesis in Oe stars, providing relevant constraints on the physical mechanism that is responsible for the Be phenomenon. Finally, we find that 3 out of 4 of the group IV Onfp stars show evidence for complex polarization effects, which are likely related to rapid rotation, and we speculate on the evolutionary links to B[e] stars.

Using Spatial Distributions to Constrain Progenitors of Supernovae and Gamma‐Ray Bursts

We carry out a comprehensive theoretical examination of the relationship between the spatial distribution of optical transients and the properties of their progenitor stars. By constructing analytic models of star-forming galaxies and the evolution of stellar populations within them, we are able to place constraints on candidate progenitors for core-collapse supernovae (SNe), long-duration gamma ray bursts, and supernovae Ia. In particular we first construct models of spiral galaxies that reproduce observations of core-collapse SNe, and we use these models to constrain the minimum mass for SNe Ic progenitors to approximately 25 solar masses. Secondly, we lay out the parameters of a dwarf irregular galaxy model, which we use to show that the progenitors of long-duration gamma-ray bursts are likely to have masses above approximately 43 solar masses. Finally, we introduce a new method for constraining the time scale associated with SNe Ia and apply it to our spiral galaxy models to show how observations can better be analyzed to discriminate between the leading progenitor models for these objects.

Constraining GRB progenitor models by probing Wolf-Rayet wind geometries in the Large Magellanic Cloud

Context. The favoured progenitors of long-duration gamma-ray bursts (GRBs) are rapidly rotating Wolf-Rayet (WR) stars. However, most Galactic WR stars are slow rotators, as stellar winds are thought to remove angular momentum. This poses a serious challenge to the collapsar model. Recent observations indicate that GRBs occur predominately in low metallicity (Z) environments, which may resolve the problem: lower Z leads to less mass loss, which may inhibit angular momentum removal, allowing WR stars to remain rotating rapidly until collapse. Aims. We wish to determine whether low Z WR stars rotate on average more rapidly than Galactic WR stars. Methods. We perform a Very Large Telescope (VLT) linear spectropolarimetry survey of WR stars in the low Z environment of the Large Magellanic Cloud (LMC) and compare our results with the Galactic sample of Harries et al. (1998, MNRAS, 296, 1072). Results. We find that only 2 out of 13 (i.e. 15%) of LMC WR stars show line polarization effects, compared to a similarly low fraction of ∼15-20% for Galactic WR stars. Conclusions. The low incidence of line polarization effects in LMC WR stars suggests that the threshold metallicity where significant differences in WR rotational properties occur is below that of the LMC (Z ∼ 0.5 Z ⊙ ), possibly constraining GRB progenitor channels to this upper metallicity.

The host galaxies of long-duration gamma-ray bursts in a cosmological hierarchical scenario

We developed a Monte Carlo code to generate long-duration gamma-ray burst (LGRB) events within cosmological hydrodynamical simulations consistent with the concordance Λ cold dark matter model. As structure is assembled, LGRBs are generated in the substructure that formed galaxies today. We adopted the collapsar model so that LGRBs are produced by single, massive stars at the final stage of their evolution. We found that the observed properties of the LGRB host galaxies (HGs) are reproduced if LGRBs are also required to be generated by low-metallicity stars. The low-metallicity condition imposed on the progenitor stars of LGRBs selects a sample of HGs with mean gas abundances of 12 + log O/H ≈ 8.6. For z < 1 the simulated HGs of low-metallicity LGRB progenitors tend to be faint, slow rotators with high star formation efficiency, compared with the general galaxy population, in agreement with observations. At higher redshift, our results suggest that larger systems with high star formation activity could also contribute to the generation of LGRBs from low-metallicity progenitors since the fraction of low-metallicity gas available for star formation increases for all systems with look-back time. Under the hypothesis of our LGRB model, our results support the claim that LGRBs could be unbiased tracers of star formation at high redshifts.

The Observed Offset Distribution of Gamma-Ray Bursts from Their Host Galaxies: A Robust Clue to the Nature of the Progenitors

We present a comprehensive study to measure the locations of gamma-ray bursts (GRBs) relative to their host galaxies. In total, we find the offsets of 20 long-duration GRBs from their apparent host galaxy centers by utilizing ground-based images from Palomar and Keck and space-based images from the Hubble Space Telescope (HST). We discuss in detail how a host galaxy is assigned to an individual GRB and the robustness of the assignment process. The median projected angular (physical) offset is 017 (1.3 kpc). The median offset normalized by the individual host half-light radii is 0.98, suggesting a strong connection of GRB locations with the UV light of their hosts. This provides strong observational evidence for the connection of GRBs to star formation. We further compare the observed offset distribution with the predicted burst locations of leading stellar-mass progenitor models. In particular, we compare the observed offset distribution with an exponential disk, a model for the location of collapsars and promptly bursting binaries (e.g., helium star–black hole binaries). The statistical comparison shows good agreement, given the simplicity of the model, with the Kolmogorov-Smirnov probability that the observed offsets derive from the model distribution of PKS = 0.45. We also compare the observed GRB offsets with the expected offset distribution of delayed merging remnant progenitors (black hole–neutron star and neutron star–neutron star binaries). We find that delayed merging remnant progenitors, insofar as the predicted offset distributions from population synthesis studies are representative, can be ruled out at the 2 × 10-3 level. This is arguably the strongest observational constraint yet against delayed merging remnants as the progenitors of long-duration GRBs. In the course of this study, we have also discovered the putative host galaxies of GRB 990510 and GRB 990308 in archival HST data.