Gamma-Ray Bursts in Circumstellar Shells: A Possible Explanation for Flares

It is now generally accepted that long-duration gamma ray bursts (GRBs) are due to the collapse of massive rotating stars. The precise collapse process itself, however, is not yet fully understood. Strong winds, outbursts, and intense ionizing UV radiation from single stars or strongly interacting binaries are expected to destroy the molecular cloud cores that give birth to them and create highly complex circumburst environments for the explosion. Such environments might imprint features on GRB light curves that uniquely identify the nature of the progenitor and its collapse. We have performed numerical simulations of realistic environments for a variety of long-duration GRB progenitors with ZEUS-MP, and have developed an analytical method for calculating GRB light curves in these profiles. Though a full, three-dimensional, relativistic magnetohydrodynamical computational model is required to precisely describe the light curve from a GRB in complex environments, our method can provide a qualitative understanding of these phenomena. We find that, in the context of the standard afterglow model, massive shells around GRBs produce strong signatures in their light curves, and that this can distinguish them from those occurring in uniform media or steady winds. These features can constrain the mass of the shell and the properties of the wind before and after the ejection. Moreover, the interaction of the GRB with the circumburst shell is seen to produce features that are consistent with observed X-ray flares that are often attributed to delayed energy injection by the central engine. Our algorithm for computing light curves is also applicable to GRBs in a variety of environments such as those in high-redshift cosmological halos or protogalaxies, both of which will soon be targets of future surveys such as JANUS or Lobster.Comment: 12 pages, 5 figures, Accepted by Ap

Radiation Hydrodynamical Instabilities in Cosmological and Galactic Ionization Fronts

Ionization fronts, the sharp radiation fronts behind which H/He ionizing photons from massive stars and galaxies propagate through space, were ubiquitous in the universe from its earliest times. The cosmic dark ages ended with the formation of the first primeval stars and galaxies a few hundred Myr after the Big Bang. Numerical simulations suggest that stars in this era were very massive, 25 - 500 solar masses, with H II regions of up to 30,000 light-years in diameter. We present three-dimensional radiation hydrodynamical calculations that reveal that the I-fronts of the first stars and galaxies were prone to violent instabilities, enhancing the escape of UV photons into the early intergalactic medium (IGM) and forming clumpy media in which supernovae later exploded. The enrichment of such clumps with metals by the first supernovae may have led to the prompt formation of a second generation of low-mass stars, profoundly transforming the nature of the first protogalaxies. Cosmological radiation hydrodynamics is unique because ionizing photons coupled strongly to both gas flows and primordial chemistry at early epochs, introducing a hierarchy of disparate characteristic timescales whose relative magnitudes can vary greatly throughout a given calculation. We describe the adaptive multistep integration scheme we have developed for the self-consistent transport of both cosmological and galactic ionization fronts.Comment: 6 pages, 4 figures, accepted for proceedings of HEDLA2010, Caltech, March 15 - 18, 201

How the First Stars Regulated Star Formation. II. Enrichment by Nearby Supernovae

Metals from Population III (Pop III) supernovae led to the formation of less massive Pop II stars in the early universe, altering the course of evolution of primeval galaxies and cosmological reionization. There are a variety of scenarios in which heavy elements from the first supernovae were taken up into second-generation stars, but cosmological simulations only model them on the largest scales. We present small-scale, high-resolution simulations of the chemical enrichment of a primordial halo by a nearby supernova after partial evaporation by the progenitor star. We find that ejecta from the explosion crash into and mix violently with ablative flows driven off the halo by the star, creating dense, enriched clumps capable of collapsing into Pop II stars. Metals may mix less efficiently with the partially exposed core of the halo, so it might form either Pop III or Pop II stars. Both Pop II and III stars may thus form after the collision if the ejecta do not strip all the gas from the halo. The partial evaporation of the halo prior to the explosion is crucial to its later enrichment by the supernova.Comment: Accepted to Ap

The Molecular Hydrogen Deficit in Gamma-Ray Burst Afterglows

Recent analysis of five gamma-ray burst (GRB) afterglow spectra reveal the absence of molecular hydrogen absorption lines, a surprising result in light of their large neutral hydrogen column densities and the detection of H$_2$ in similar, more local star-forming regions like 30 Doradus in the Large Magellanic Cloud (LMC). Observational evidence further indicates that the bulk of the neutral hydrogen column in these sight lines lies 100 pc beyond the progenitor and that H$_2$ was absent prior to the burst, suggesting that direct flux from the star, FUV background fields, or both suppressed its formation. We present one-dimensional radiation hydrodynamical models of GRB host galaxy environments, including self-consistent radiative transfer of both ionizing and Lyman-Werner photons, nine-species primordial chemistry with dust formation of H$_2$, and dust extinction of UV photons. We find that a single GRB progenitor is sufficient to ionize neutral hydrogen to distances of 50 - 100 pc but that a galactic Lyman-Werner background is required to dissociate the molecular hydrogen in the ambient ISM. Intensities of 0.1 - 100 times the Galactic mean are necessary to destroy H$_2$ in the cloud, depending on its density and metallicity. The minimum radii at which neutral hydrogen will be found in afterglow spectra is insensitive to the mass of the progenitor or the initial mass function (IMF) of its cluster, if present.Comment: 12 pages, 7 figures, accepted for Ap

A Multistep Algorithm for the Radiation Hydrodynamical Transport of Cosmological Ionization Fronts and Ionized Flows

Radiation hydrodynamical transport of ionization fronts in the next generation of cosmological reionization simulations holds the promise of predicting UV escape fractions from first principles as well as investigating the role of photoionization in feedback processes and structure formation. We present a multistep integration scheme for radiative transfer and hydrodynamics for accurate propagation of I-fronts and ionized flows from a point source in cosmological simulations. The algorithm is a photon-conserving method which correctly tracks the position of I-fronts at much lower resolutions than non-conservative techniques. The method applies direct hierarchical updates to the ionic species, bypassing the need for the costly matrix solutions required by implicit methods while retaining sufficient accuracy to capture the true evolution of the fronts. We review the physics of ionization fronts in power-law density gradients, whose analytical solutions provide excellent validation tests for radiation coupling schemes. The advantages and potential drawbacks of direct and implicit schemes are also considered, with particular focus on problem timestepping which if not properly implemented can lead to morphologically plausible I-front behavior that nonetheless departs from theory. We also examine the effect of radiation pressure from very luminous central sources on the evolution of I-fronts and flows.Comment: 25 pages, 16 figures, accepted to ApJ. Minor revisions included. Full resolution PDF available at http://cosmos.ucsd.edu/~dwhalen/downloads/dwhalen_zeusmp_method.pd