Direct-collapse black holes (DCBHs) are currently one of the leading contenders for the origins of the first quasars in the universe, over 300 of which have now been found at $z >$ 6. But the birth of a DCBH in an atomically-cooling halo does not by itself guarantee it will become a quasar by $z \sim$ 7, the halo must also be located in cold accretion flows or later merge with a series of other gas-rich halos capable of fueling the BH's rapid growth. Here, we present near infrared luminosities for DCBHs born in cold accretion flows in which they are destined to grow to 10$^9$ M$_{\odot}$ by $z \sim$ 7. Our observables, which are derived from cosmological simulations with radiation hydrodynamics with Enzo, reveal that DCBHs could be found by the James Webb Space Telescope at $z \lesssim$ 20 and strongly-lensed DCBHs might be found in future wide-field surveys by Euclid and the Wide-Field Infrared Space Telescope at $z \lesssim$ 15.

Very massive 140-260 Msun stars can die as highly-energetic pair-instability supernovae (PI SNe) with energies of up to 100 times those of core-collapse SNe that can completely destroy the star, leaving no compact remnant behind. These explosions can synthesize $0.1-30$ Msun of radioactive Ni56, which can cause them to rebrighten at later times when photons due to Ni56 decay diffuse out of the ejecta. However, heat from the decay of such large masses of Ni56 could also drive important dynamical effects deep in the ejecta that are capable of mixing elements and affecting the observational signatures of these events. We have now investigated the dynamical effect of Ni56 heating on PI SN ejecta with high-resolution two-dimensional hydrodynamic simulations performed with the CASTRO code. We find that expansion of the hot Ni56 bubble forms a shell at the base of the silicon layer of the ejecta about 200 days after the explosion but that no hydrodynamical instabilities develop that would mix Ni56 with the Si/O-rich ejecta. However, while the dynamical effects of Ni56 heating may be weak they could affect the observational signatures of some PI SNe by diverting decay energy into internal expansion of the ejecta at the expense of rebrightening at later times.

The leading contenders for the seeds of the first quasars are direct collapse black holes (DCBHs) formed during catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 20. The discovery of the Ly$\alpha$ emitter CR7 at $z =$ 6.6 was initially held to be the first detection of a DCBH, although this interpretation has since been challenged on the grounds of Spitzer IRAC and Very Large Telescope X-Shooter data. Here we determine if radio flux from a DCBH in CR7 could be detected and discriminated from competing sources of radio emission in the halo such as young supernovae and H II regions. We find that a DCBH would emit a flux of 10 - 200 nJy at 1.0 GHz, far greater than the sub-nJy signal expected for young supernovae but on par with continuum emission from star-forming regions. However, radio emission from a DCBH in CR7 could be distinguished from free-free emission from H II regions by its spectral evolution with frequency and could be detected by the Square Kilometer Array in the coming decade.

A rapidly spinning magnetar in a young supernova (SN) can produce a superluminous transient by converting a fraction of its rotational energy into radiation. Here, we present the first three-dimensional hydrodynamical simulations ever performed of a magnetar-powered SN in the circumstellar medium formed by the ejection of the outer layers of the star prior to the blast. We find that hydrodynamical instabilities form on two scales in the ejecta, not just one as in ordinary core-collapse SNe: in the hot bubble energized by the magnetar and in the forward shock of the SN as it plows up ambient gas. Pressure from the bubble also makes the instabilities behind the forward shock more violent and causes more mixing in the explosion than in normal SNe, with important consequences for the light curves and spectra of the event that cannot be captured by one-dimensional models. We also find that the magnetar can accelerate Ca and Si to velocities of $\sim $ 12000 km/s and account for their broadened emission lines in observations. Our simulations also reveal that energy from even weak magnetars can accelerate iron-group elements deep in the ejecta to $5000-7000$ km/s and explain the high-velocity Fe observed at early times in some core-collapse SNe such as SN 1987A.

Supermassive primordial stars forming during catastrophic baryon collapse in atomically-cooling halos at $z \sim$ 15 - 20 may be the origin of the first quasars in the universe. However, no simulation to date has followed the evolution of these halos at resolutions that are high enough or for times that are long enough to determine if collapse actually produces SMSs. Here we report new cosmological simulations of baryon collapse in atomically-cooled halos for times that are long enough for SMSs to form and die as direct-collapse black holes (DCBHs). We find that the high infall rates required to build up such stars do persist until the end of their lives and could fuel the rapid growth of their BHs thereafter. Our simulations also demonstrate that binary and even small multiples of SMSs can form in low-spin and high-spin halos, respectively. This discovery raises the exciting prospect of detecting gravitational waves from DCBH mergers with LISA and tidal disruption events in the near infrared with the {\em James Webb Space Telescope} and ground-based telescopes in the coming decade.

Supermassive primordial stars in hot, atomically-cooling haloes at $z \sim$ 15 - 20 may have given birth to the first quasars in the universe. Most simulations of these rapidly accreting stars suggest that they are red, cool hypergiants, but more recent models indicate that some may have been bluer and hotter, with surface temperatures of 20,000 - 40,000 K. These stars have spectral features that are quite distinct from those of cooler stars and may have different detection limits in the near infrared (NIR) today. Here, we present spectra and AB magnitudes for hot, blue supermassive primordial stars calculated with the TLUSTY and CLOUDY codes. We find that photometric detections of these stars by the James Webb Space Telescope (JWST) will be limited to $z \lesssim$ 10 - 12, lower redshifts than those at which red stars can be found, because of quenching by their accretion envelopes. With moderate gravitational lensing, Euclid and the Wide-Field Infrared Space Telescope (WFIRST) could detect blue supermassive stars out to similar redshifts in wide-field surveys.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

A major scientific goal of JWST is to probe the epoch of re-ionization of the Universe at z above 6, and up to 20 and beyond. At these redshifts, galaxies are just beginning to form and the observable objects are early black holes, supernovae, and cosmic infrared background. The JWST has the necessary sensitivity to observe these targets individually, but a public deep and wide science enabling survey in the wavelength range from 2-5 $\mu$m is needed to discover these black holes and supernovae and to cover the area large enough for cosmic infrared background to be reliably studied. This enabling survey will also discover a large number of other transients and enable sciences such as supernova cosmology up to z $\sim$ 5, star formation history at high redshift through supernova explosions, faint stellar objects in the Milky Way, and galaxy evolution up to z approaching 10. The results of this survey will also serve as an invaluable target feeder for the upcoming era of ELT and SKA.

The collapse of supermassive primordial stars in hot, atomically-cooled halos may have given birth to the first quasars at $z \sim$ 15 - 20. Recent numerical simulations of these rapidly accreting stars reveal that they are cool, red hypergiants shrouded by dense envelopes of pristine atomically-cooled gas at 6,000 - 8,000 K, with luminosities $L$ $\gtrsim$ 10$^{10}$ L$_{\odot}$. Could such luminous but cool objects be detected as the first stage of quasar formation in future near infrared (NIR) surveys? We have now calculated the spectra of supermassive primordial stars in their birth envelopes with the Cloudy code. We find that some of these stars will be visible to JWST at $z \lesssim$ 20 and that with modest gravitational lensing Euclid and WFIRST could detect them out to $z \sim$ 10 - 12. Rather than obscuring the star, its accretion envelope enhances its visibility in the NIR today by reprocessing its short-wavelength flux into photons that are just redward of the Lyman limit in the rest frame of the star.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

Supernovae (SNe) could be powerful probes of the properties of stars and galaxies at high redshifts in future surveys. Wide fields and longer exposure times are required to offset diminishing star formation rates and lower fluxes to detect useful numbers of events at high redshift. In principle, the Large Synoptic Survey Telescope (LSST) could discover large numbers of early SNe because of its wide fields but only at lower redshifts because of its AB mag limit of ~ 24. But gravitational lensing by galaxy clusters and massive galaxies could boost flux from ancient SNe and allow LSST to detect them at earlier times. Here, we calculate detection rates for lensed SNe at z ~ 5 - 7 for LSST. We find that the LSST Wide Deep Fast survey could detect up to 120 lensed Population (Pop) I and II SNe but no lensed Pop III SNe. Deep-drilling programs in a single 10 square degree FoV could detect Pop I and II core-collapse SNe and Pop III pair-instability SNe at AB magnitudes of 27 - 28 and 26, respectively. An alternative deep survey over 80 nights with a one-year cadence could find ~ 8 Pop III SNe.

Supermassive black holes (SMBHs) are the central engines of luminous quasars and are found in most massive galaxies today. But the recent discoveries of ULAS J1120+0641, a $2 \times 10^9$ M$_{\odot}$ BH at $z =$ 7.1, and ULAS J1342+0928, a $8.0 \times 10^{8}$ M$_{\odot}$ BH at $z =$ 7.5, now push the era of quasar formation up to just 690 Myr after the Big Bang. Here we report new cosmological simulations of SMBHs with X-rays fully coupled to primordial chemistry and hydrodynamics that show that J1120 and J1342 can form from direct collapse black holes (DCBHs) if their growth is fed by cold, dense accretion streams, like those thought to fuel rapid star formation in some galaxies at later epochs. Our models reproduce all of the observed properties of J1120: its mass, luminosity, and H II region as well as star formation rates and metallicities in its host galaxy. They also reproduce the dynamical mass of the innermost 1.5 kpc of its emission region recently measured by ALMA and J-band magnitudes that are in good agreement with those found by the VISTA Hemisphere Survey.

Supermassive stars born from pristine gas in atomically-cooled haloes are thought to be the progenitors of supermassive black holes at high redshifts. However, the way they accrete their mass is still an unsolved problem. In particular, for accretion to proceed, a large amount of angular momentum has to be extracted from the collapsing gas. Here, we investigate the constraints stellar evolution imposes on this angular momentum problem. We present an evolution model of a supermassive Population III star including simultaneously accretion and rotation. We find that, for supermassive stars to form by accretion, the accreted angular momentum has to be about 1% of the Keplerian angular momentum. This tight constraint comes from the $\Omega\Gamma$-limit, at which the combination of radiation pressure and centrifugal force cancels gravity. It implies that supermassive stars are slow rotators, with a surface velocity less than 10-20% of their first critical velocity, at which the centrifugal force alone cancels gravity. At such low velocities, the deformation of the star due to rotation is negligible.

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

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.

Some ancient, dim, metal-poor stars may have formed in the ashes of the first supernovae (SNe). If their chemical abundances can be reconciled with the elemental yields of specific Population III (Pop III) explosions, they could reveal the properties of primordial stars. But multidimensional simulations of such explosions are required to predict their yields because dynamical instabilities can dredge material up from deep in the ejecta that would otherwise be predicted to fall back on to the central remnant and be lost in one-dimensional (1D) models. We have performed two-dimensional (2D) numerical simulations of two low-energy Pop III SNe, a 12.4 Msun explosion and a 60 Msun explosion, and find that they produce elemental yields that are a good fit to those measured in the most iron-poor star discovered to date, SMSS J031300.36-670839.3 (J031300). Fallback on to the compact remnant in these weak explosions accounts for the lack of measurable iron in J031300 and its low iron-group abundances in general. Our 2D explosions produce higher abundances of heavy elements (atomic number Z > 20) than their 1D counterparts due to dredge-up by fluid instabilities. Since almost no Ni is ejected by these weak SNe, their low luminosities will prevent their detection in the near-infrared with the James Webb Space Telescope and future 30-m telescopes on the ground. The only evidence that they ever occurred will be in the fossil abundance record.

Nascent neutron stars with millisecond periods and magnetic fields in excess of $10^{16}$ Gauss can drive highly energetic and asymmetric explosions known as magnetar-powered supernovae. These exotic explosions are one theoretical interpretation for supernovae Ic-BL which are sometimes associated with long gamma-ray bursts. Twisted magnetic field lines extract the rotational energy of the neutron star and release it as a disk wind or a jet with energies greater than 10$^{52}$ erg over $\sim 20$ sec. What fractions of the energy of the central engine go into the wind and the jet remain unclear. We have performed two-dimensional hydrodynamical simulations of magnetar-powered supernovae (SNe) driven by disk winds and jets with the CASTRO code to investigate the effect of the central engine on nucleosynthetic yields, mixing, and light curves. We find that these explosions synthesize less than 0.05 Msun of Ni and that this mass is not very sensitive to central engine type. The morphology of the explosion can provide a powerful diagnostic of the properties of the central engine. In the absence of a circumstellar medium these events are not very luminous, with peak bolometric magnitudes $M_b \sim -16.5 $ due to low Ni production.

Supermassive primordial stars are suspected to be the progenitors of the most massive quasars at z~6. Previous studies of such stars were either unable to resolve hydrodynamical timescales or considered stars in isolation, not in the extreme accretion flows in which they actually form. Therefore, they could not self-consistently predict their final masses at collapse, or those of the resulting supermassive black hole seeds, but rather invoked comparison to simple polytropic models. Here, we systematically examine the birth, evolution and collapse of accreting non-rotating supermassive stars under accretion rates of 0.01-10 solar masses per year, using the stellar evolution code KEPLER. Our approach includes post-Newtonian corrections to the stellar structure and an adaptive nuclear network, and can transition to following the hydrodynamic evolution of supermassive stars after they encounter the general relativistic instability. We find that this instability triggers the collapse of the star at masses of 150,000-330,000 solar masses for accretion rates of 0.1-10 solar masses per year, and that the final mass of the star scales roughly logarithmically with the rate. The structure of the star, and thus its stability against collapse, is sensitive to the treatment of convection, and the heat content of the outer accreted envelope. Comparison with other codes suggests differences here may lead to small deviations in the evolutionary state of the star as a function of time, that worsen with accretion rate. Since the general relativistic instability leads to the immediate death of these stars, our models place an upper limit on the masses of the first quasars at birth.

Supermassive primordial stars forming in atomically-cooled halos at $z \sim15-20$ are currently thought to be the progenitors of the earliest quasars in the Universe. In this picture, the star evolves under accretion rates of $0.1 - 1$ $M_\odot$ yr$^{-1}$ until the general relativistic instability triggers its collapse to a black hole at masses of $\sim10^5$ $M_\odot$. However, the ability of the accretion flow to sustain such high rates depends crucially on the photospheric properties of the accreting star, because its ionising radiation could reduce or even halt accretion. Here we present new models of supermassive Population III protostars accreting at rates $0.001 - 10$ $M_\odot$ yr$^{-1}$, computed with the GENEVA stellar evolution code including general relativistic corrections to the internal structure. We use the polytropic stability criterion to estimate the mass at which the collapse occurs, which has been shown to give a lower limit of the actual mass at collapse in recent hydrodynamic simulations. We find that at accretion rates higher than $0.001$ $M_\odot$ yr$^{-1}$ the stars evolve as red, cool supergiants with surface temperatures below $10^4$ K towards masses $>10^5$ $M_\odot$, and become blue and hot, with surface temperatures above $10^5$ K, only for rates $\lesssim0.001$ $M_\odot$ yr$^{-1}$. Compared to previous studies, our results extend the range of masses and accretion rates at which the ionising feedback remains weak, reinforcing the case for direct collapse as the origin of the first quasars.

Direct collapse black holes forming in pristine, atomically-cooling haloes at $z \approx 10-20$ may act as the seeds of supermassive black holes (BH) at high redshifts. In order to create a massive BH seed, the host halo needs to be prevented from forming stars. H$_2$ therefore needs to be irradiated by a large flux of Lyman-Werner (LW) UV photons in order to suppress H$_2$ cooling. A key uncertainty in this scenario is the escape fraction of LW radiation from first galaxies, the dominant source of UV photons at this epoch. To better constrain this escape fraction, we have performed radiation-hydrodynamical simulations of the growth of HII regions and their associated photodissociation regions in the first galaxies using the ZEUS-MP code. We find that the LW escape fraction crucially depends on the propagation of the ionisation front (I-front). For an R-type I-front overrunning the halo, the LW escape fraction is always larger than 95%. If the halo recombines later from the outside--in, due to a softened and weakened spectrum, the LW escape fraction in the rest-frame of the halo (the near-field) drops to zero. A detailed and careful analysis is required to analyse slowly moving, D-type I-fronts, where the escape fraction depends on the microphysics and can be as small as 3% in the near-field and 61% in the far-field or as large as 100% in both the near-field and the far-field.

Supernovae are important probes of the properties of stars at high redshifts because they can be detected at early epochs and their masses can be inferred from their light curves. Finding the first cosmic explosions in the universe will only be possible with the James Webb Space Telescope, the Wide-Field Infrared Survey Telescope and the next generation of extremely large telescopes. But strong gravitational lensing by massive clusters, like those in the Cluster Lensing and Supernova Survey with Hubble (CLASH), could reveal such events now by magnifying their flux by factors of 10 or more. We find that CLASH will likely discover at least 2 - 3 core-collapse supernovae at 5 < z < 12 and perhaps as many as ten. Future surveys of cluster lenses similar in scope to CLASH by the James Webb Space Telescope might find hundreds of these events out to z ~ 15 - 17. Besides revealing the masses of early stars, these ancient supernovae will also constrain cosmic star formation rates in the era of first galaxy formation.

While the Lyman-$\alpha$ ($\mathrm{Ly}\alpha$) emission line serves as an important tool in the study of galaxies at $z\lesssim 6$, finding Ly$\alpha$ emitters (LAE) at significantly higher redshifts has been more challenging, probably because of the increasing neutrality of the intergalactic medium above $z\sim6$. Galaxies with extremely high rest-frame Ly$\alpha$ equivalent widths, EW(Ly$\alpha$) $\gtrsim 150$ \AA{}, at $z>6$ are good candidates for Ly$\alpha$ follow-up observations, and can stand out in multiband imaging surveys because of their unusual colors. We have conducted a photometric search for such objects in the Cluster Lensing And Supernova survey with Hubble (CLASH), and report here the identification of three likely gravitationally-lensed images of a single LAE candidate at $z\sim6.3$, behind the galaxy cluster Abell 2261($z = 0.225$). In the process, we also measured with Keck/MOSFIRE the first spectroscopic redshift of a multiply-imaged galaxy behind Abell 2261, at $z = 3.337$. This allows us to calibrate the lensing model, which in turn is used to study the properties of the candidate LAE. Population III galaxy spectral energy distribution (SED) model fits to the CLASH broadband photometry of the possible LAE provide a slightly better fit than Population I/II models. The best fitted model suggests intrinsic EW(Ly$\alpha$) $\approx 160$ \AA{} after absorption in the interstellar and intergalactic medium. Future spectroscopic observations will examine this prediction as well as shed more light on the morphology of this object, which indicates it may be a merger of two smaller galaxies.

With new observational facilities becoming available soon, discovering and characterising supernovae from the first stars will open up alternative observational windows to the end of the cosmic dark ages. Based on a semi-analytical merger tree model of early star formation we constrain Population III supernova rates. We find that our method reproduces the Population III supernova rates of large-scale cosmological simulations very well. Our computationally efficient model allows us to survey a large parameter space and to explore a wide range of different scenarios for Population III star formation. Our calculations show that observations of the first supernovae can be used to differentiate between cold and warm dark matter models and to constrain the corresponding particle mass of the latter. Our predictions can also be used to optimize survey strategies with the goal to maximize supernova detection rates.

CR7 is the brightest Lyman-$\alpha$ emitter observed at $z>6$, which shows very strong Lyman-$\alpha$ and HeII 1640\AA\ line luminosities, but no metal line emission. Previous studies suggest that CR7 hosts either young primordial stars with a total stellar mass of $\sim 10^7\,\mathrm{M}_\odot$ or a black hole of $\gtrsim 10^6\,\mathrm{M}_\odot$. Here, we explore different formation scenarios for CR7 with a semianalytical model, based on the random sampling of dark matter merger trees. We are unable to reproduce the observational constraints with a primordial stellar source, given our model assumptions, due to the short stellar lifetimes and the early metal enrichment. Black holes that are the remnants of the first stars are either not massive enough, or reside in metal-polluted haloes, ruling out this possible explanation of CR7. Our models instead suggest that direct collapse black holes, which form in metal-free haloes exposed to large Lyman-Werner fluxes, are more likely the origin of CR7. However, this result is derived under optimistic assumptions and future observations are necessary to further constrain the nature of CR7.

Supernovae are important probes of the properties of stars at high redshifts because they can be detected at early epochs and their masses can be inferred from their light curves. Direct detection of the first cosmic explosions in the universe will only be possible with JWST, WFIRST and the next generation of extremely large telescopes. But strong gravitational lensing by massive clusters, like those in the Frontier Fields, could reveal supernovae at slightly lower redshifts now by magnifying their flux by factors of 10 or more. We find that Frontier Fields will likely discover dozens of core-collapse supernovae at 5 $ < z <$ 12. Future surveys of cluster lenses similar in scope to Frontier Fields by JWST might find hundreds of these events out to $z \sim$ 15 - 17. Besides revealing the masses of early stars, these ancient supernovae could also constrain cosmic star formation rates in the era of first galaxy formation.

The first supernovae will soon be visible at the edge of the observable universe, revealing the birthplaces of Population III stars. With upcoming near-infrared missions, a broad analysis of the detectability of high-$z$ supernovae is paramount. We combine cosmological and radiation transport simulations, instrument specifications, and survey strategies to create synthetic observations of primeval core-collapse, Type IIn and pair-instability supernovae with the James Webb Space Telescope ($JWST$). We show that a dedicated observational campaign with the $JWST$ can detect up to $\sim 15$ pair-instability explosions, $\sim 300$ core-collapse supernovae, but less than one Type IIn explosion per year, depending on the Population III star formation history. Our synthetic survey also shows that $\approx 1-2 \times10^2$ supernovae detections, depending on the accuracy of the classification, are sufficient to discriminate between a Salpeter and flat mass distribution for high redshift stars with a confidence level greater than 99.5 per cent. We discuss how the purity of the sample affects our results and how supervised learning methods may help to discriminate between CC and PI SNe.

Direct collapse black holes (DCBHs) formed from the collapse of atomically-cooled primordial gas in the early Universe are strong candidates for the seeds of supermassive BHs. DCBHs are thought to form in atomic cooling haloes in the presence of a strong molecule-dissociating, Lyman-Werner (LW) radiation field. Given that star forming galaxies are likely to be the source of the LW radiation in this scenario, ionizing radiation from these galaxies may accompany the LW radiation. We present cosmological simulations resolving the collapse of primordial gas into an atomic cooling halo, including the effects of both LW and ionizing radiation. We find that in cases where the gas is not self-shielded from the ionizing radiation, the collapse can be delayed by ~ 25 Myr. When the ionized gas does collapse, the free electrons that are present catalyze H2 formation. In turn, H2 cooling becomes efficient in the center of the halo, and DCBH formation is prevented. We emphasize, however, that in many cases the gas collapsing into atomic cooling haloes at high redshift is self-shielding to ionizing radiation. Therefore, it is only in a fraction of such haloes in which DCBH formation is prevented due to reionization.

Pair Instability Supernovae have been suggested as candidates for some Super Luminous Supernovae, such as SN 2007bi, and as one of the dominant types of explosion occurring in the early Universe from massive, zero-metallicity Population III stars. The progenitors of such events can be rapidly rotating, therefore exhibiting different evolutionary properties due to the effects of rotationally-induced mixing and mass-loss. Proper identification of such events requires rigorous radiation hydrodynamics and radiative transfer calculations that capture not only the behavior of the light curve but also the spectral evolution of these events. We present radiation hydrodynamics and radiation transport calculations for 90-300 Msun rotating pair-instability supernovae covering both the shock break-out and late light curve phases. We also investigate cases of different initial metallicity and rotation rate to determine the impact of these parameters on the detailed spectral characteristics of these events. In agreement with recent results on non-rotating pair instability supernovae, we find that for a range of progenitor masses and rotation rates these events have intrinsically red colors in contradiction with observations of super-luminous supernovae. The spectroscopic properties of rotating pair instability supernovae are similar to those of non-rotating events with stripped hydrogen and helium envelopes. We find that the progenitor metallicity and rotation rate properties are erased after the explosion and cannot be identified in the resulting model spectra. It is the combined effects of pre-supernova mass-loss and the basic properties of the supernova ejecta such as mass, temperature and velocity that have the most direct impact in the model spectra of pair instability supernovae.

Population III supernovae have been of growing interest of late for their potential to directly probe the properties of the first stars, particularly the most energetic events that are visible near the edge of the observable universe. But until now, hypernovae, the unusually energetic Type Ib/c supernovae that are sometimes associated with gamma-ray bursts, have been overlooked as cosmic beacons at the highest redshifts. In this, the latest of a series of studies on Population III supernovae, we present numerical simulations of 25 - 50 M$_{\odot}$ hypernovae and their light curves done with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible at z = 10 - 15 to the James Webb Space Telescope (JWST) and z = 4 - 5 to the Wide-Field Infrared Survey Telescope (WFIRST), tracing star formation rates in the first galaxies and at the end of cosmological reionization. If, however, the hypernova crashes into a dense shell ejected by its progenitor, it is expected that a superluminous event will occur that may be seen at z ~ 20, in the first generation of stars.

Gamma-ray bursts (GRBs) are the ultimate cosmic lighthouses, capable of illuminating the universe at its earliest epochs. Could such events probe the properties of the first stars at z $\sim$ 20, the end of the cosmic Dark Ages? Previous studies of Population III GRBs only considered explosions in the diffuse relic H II regions of their progenitors, or bursts that are far more more energetic than those observed to date. But the processes that produce GRBs at the highest redshifts likely reset their local environments, creating much more complicated structures than those in which relativistic jets have been modeled so far. These structures can greatly affect the luminosity of the afterglow, and hence the redshift at which it can be detected. We have now simulated Population III GRB afterglows in H II regions, winds, and dense shells ejected by the star during the processes that produce the burst. Our model, which has been used in previous work, has been extended to include contributions from reverse shocks, inverse Compton cooling and the effects of sphericity and beaming in the blast wave, and is valid in a variety of circumjet density profiles. We find that GRBs with E$_{\mathrm{iso},\gamma} =$ 10$^{ 51}$ - 10$^{53}$ erg will be visible at z $\gtrsim$ 20 to the next generation of near infrared and radio observatories. In many cases, the environment of the burst, and hence progenitor type, can be inferred from the afterglow light curve. Although some Population III GRBs are visible to Swift and the Very Large Array now, the optimal strategy for their detection will be future missions like EXIST and JANUS, which have large survey areas and onboard X-ray and infrared telescopes that can track their near infrared flux from the moment of the burst, thereby identifying its redshift.

Population III stars can regulate star formation in the primordial Universe in several ways. They can ionize nearby halos, and even if their ionizing photons are trapped by their own halos, their Lyman-Werner (LW) photons can still escape and destroy H$_2$ in other halos, preventing them from cooling and forming stars. LW escape fractions are thus a key parameter in cosmological simulations of early reionization and star formation but have not yet been parametrized for realistic halos by halo or stellar mass. To do so, we perform radiation hydrodynamical simulations of LW UV escape from 9--120 M$_{\odot}$ Pop III stars in $10^5$ to $10^7$ M$_{\odot}$ halos with ZEUS-MP. We find that photons in the LW lines (i.e. those responsible for destroying H$_{2}$ in nearby systems) have escape fractions ranging from 0% to 85%. No LW photons escape the most massive halo in our sample, even from the most massive star. Escape fractions for photons elsewhere in the 11.18--13.6~eV energy range, which can be redshifted into the LW lines at cosmological distances, are generally much higher, being above 60% for all but the least massive stars in the most massive halos. We find that shielding of H$_2$ by neutral hydrogen, which has been neglected in most studies to date, produces escape fractions that are up to a factor of three smaller than those predicted by H$_2$ self-shielding alone.

The era of the universe's first (Population III) stars is essentially unconstrained by observation. Ultra-luminous and massive stars from this time altered the chemistry of the cosmos, provided the radiative scaffolding to support the formation of the first protogalaxies, and facilitated the creation and growth of now-supermassive black holes. Unfortunately, because these stars lie literally at the edge of the observable universe, they will remain beyond the reach of even the next generation of telescopes such as the James Webb Space Telescope and the Thirty-Meter Telescope. In this paper, we provide a primer to supernovae modeling and the first stars to make our discussion accessible to those new to or outside our field. We review recent work of the Los Alamos Supernova Light Curve Project and Brigham Young University to explore the possibility of probing this era through observations of the spectacular deaths of the first stars. We find that many such brilliant supernova explosions will be observable as far back as $\sim 99$% of the universe's current age, tracing primordial star formation rates and the locations of their protogalaxies on the sky. The observation of Population III supernovae will be among the most spectacular discoveries in observational astronomy in the coming decade.

Observational evidence suggests that some very massive stars in the local Universe may die as pair-instability supernovae. We present 2D simulations of the pair-instability supernova of a non-zero metallicity star. We find that very little mixing occurs in this explosion because metals in the stellar envelope drive strong winds that strip the hydrogen envelope from the star prior to death. Consequently, a reverse shock cannot form and trigger fluid instabilities during the supernova. Only weak mixing driven by nuclear burning occurs in the earliest stages of the supernova, and it is too weak to affect the observational signatures of the explosion.

The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500 M$_{\odot}$ Z $\sim$ 0.1 - 0.4 Z$_{\odot}$ stars. These stars lose up to 80% of their mass to strong line-driven winds and explode as bare He cores. We find that their light curves and spectra are quite different from those of Population III pair-instability explosions, which therefore cannot be used as templates for low-redshift events. Although non-zero metallicity pair-instability supernovae are generally dimmer than their Population III counterparts, in some cases they will be bright enough to be detected at the earliest epochs at which they can occur, the formation of the first galaxies at $z \sim$ 10 - 15. Others can masquerade as dim, short duration supernovae that are only visible in the local universe and that under the right conditions could be hidden in a wide variety of supernova classes. We also report for the first time that some pair-instability explosions can create black holes with masses of $\sim$ 100 M$_{\odot}$.

Population III stars that die as pair-instability supernovae are usually thought to fall in the mass range of 140 - 260 M$_{\odot}$. But several lines of work have now shown that rotation can build up the He cores needed to encounter the pair instability at stellar masses as low as 90 $_{\odot}$. Depending on the slope of the initial mass function of Population III stars, there could be 4 - 5 times as many stars from 90 - 140 $_{\odot}$ in the primordial universe than in the usually accepted range. We present numerical simulations of the pair-instability explosions of such stars performed with the MESA, FLASH and RAGE codes. We find that they will be visible to supernova factories such as Pan-STARRS and LSST in the optical out to z $\sim$ 1 - 2 and to JWST and the 30 m-class telescopes in the NIR out to $z \sim$ 7 - 10. Such explosions will thus probe the stellar populations of the first galaxies and cosmic star formation rates in the era of cosmological reionization. These supernovae are also easily distinguished from more massive pair-instability explosions, underscoring the fact that there is far greater variety to the light curves of these events than previously understood.

Numerical studies of primordial star formation suggest that the first stars in the universe may have been very massive. Stellar models indicate that non-rotating Population III stars with initial masses of 140-260 Msun die as highly energetic pair-instability supernovae. We present new two-dimensional simulations of primordial pair-instability supernovae done with the CASTRO code. Our simulations begin at earlier times than previous multidimensional models, at the onset of core collapse, to capture any dynamical instabilities that may be seeded by collapse and explosive burning. Such instabilities could enhance explosive yields by mixing hot ash with fuel, thereby accelerating nuclear burning, and affect the spectra of the supernova by dredging up heavy elements from greater depths in the star at early times. Our grid of models includes both blue supergiants and red supergiants over the range in progenitor mass expected for these events. We find that fluid instabilities driven by oxygen and helium burning arise at the upper and lower boundaries of the oxygen shell $\sim$ 20 - 100 seconds after core bounce. Instabilities driven by burning freeze out after the SN shock exits the helium core. As the shock later propagates through the hydrogen envelope, a strong reverse shock forms that drives the growth of Rayleigh--Taylor instabilities. In red supergiant progenitors, the amplitudes of these instabilities are sufficient to mix the supernova ejecta.

Massive stars that end their lives with helium cores in the range of 35 to 65 Msun are known to produce repeated thermonuclear outbursts due to a recurring pair-instability. In some of these events, solar masses of material are ejected in repeated outbursts of several times 10$^{50}$ erg each. Collisions between these shells can sometimes produce very luminous transients that are visible from the edge of the observable universe. Previous 1D studies of these events produce thin, high-density shells as one ejection plows into another. Here, in the first multidimensional simulations of these collisions, we show that the development of a Rayleigh-Taylor instability truncates the growth of the high density spike and drives mixing between the shells. The progenitor is a 110 Msun solar-metallicity star that was shown in earlier work to produce a superluminous supernova. The light curve of this more realistic model has a peak luminosity and duration that are similar to those of 1D models but a structure that is smoother.

The formation of supermassive Population III stars with masses $\gtrsim$ 10,000 Msun in primeval galaxies in strong UV backgrounds at $z \sim$ 15 may be the most viable pathway to the formation of supermassive black holes by $z \sim$ 7. Most of these stars are expected to live for short times and then directly collapse to black holes, with little or no mass loss over their lives. But we have now discovered that non-rotating primordial stars with masses close to 55,000 Msun can instead die as highly energetic thermonuclear supernovae powered by explosive helium burning, releasing up to 10$ ^{55}$ erg, or about 10,000 times the energy of a Type Ia supernova. The explosion is triggered by the general relativistic contribution of thermal photons to gravity in the core of the star, which causes the core to contract and explosively burn. The energy release completely unbinds the star, leaving no compact remnant, and about half of the mass of the star is ejected into the early cosmos in the form of heavy elements. The explosion would be visible in the near infrared at $z \lesssim$ 20 to {\it Euclid} and the Wide-Field Infrared Survey Telescope (WFIRST), perhaps signaling the birth of supermassive black hole seeds and the first quasars.

Population III supernovae have been the focus of growing attention because of their potential to directly probe the properties of the first stars, particularly the most energetic events that can be seen at the edge of the observable universe. But until now pair-pulsation supernovae, in which explosive thermonuclear burning in massive stars fails to unbind them but can eject their outer layers into space, have been overlooked as cosmic beacons at the earliest redshifts. These shells can later collide and, like Type IIn supernovae, produce superluminous events in the UV at high redshifts that could be detected in the near infrared today. We present numerical simulations of a 110 M$_{\odot}$ pair-pulsation explosion done with the Los Alamos radiation hydrodynamics code RAGE. We find that collisions between consecutive pair pulsations are visible in the near infrared out to z $\sim$ 15 - 20 and can probe the earliest stellar populations at cosmic dawn.

We present a fully self-consistent simulation of a synthetic survey of the furthermost cosmic explosions. The appearance of the first generation of stars (Population III) in the Universe represents a critical point during cosmic evolution, signaling the end of the dark ages, a period of absence of light sources. Despite their importance, there is no confirmed detection of Population III stars so far. A fraction of these primordial stars are expected to die as pair-instability supernovae (PISNe), and should be bright enough to be observed up to a few hundred million years after the big bang. While the quest for Population III stars continues, detailed theoretical models and computer simulations serve as a testbed for their observability. With the upcoming near-infrared missions, estimates of the feasibility of detecting PISNe are not only timely but imperative. To address this problem, we combine state-of-the-art cosmological and radiative simulations into a complete and self-consistent framework, which includes detailed features of the observational process. We show that a dedicated observational strategy using $\lesssim 8$ per cent of total allocation time of the James Webb Space Telescope mission can provide us up to $\sim 9-15$ detectable PISNe per year.

The existence of supermassive black holes as early as z ~ 7 is one of the great unsolved problems in cosmological structure formation. One leading theory argues that they are born during catastrophic baryon collapse in z ~ 15 protogalaxies in strong Lyman-Werner UV backgrounds. Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to directly collapse to 10^4 - 10^5 solar-mass black holes. We have now discovered that some of these fragments can instead become supermassive stars that eventually explode as pair-instability supernovae with energies of ~ 10^55 erg, the most energetic explosions in the universe. We have calculated light curves and spectra for supermassive Pop III PI SNe with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible in NIR all-sky surveys by Euclid out to z ~ 10 - 15 and by WFIRST and WISH out to z ~ 15 - 20, perhaps revealing the birthplaces of the first quasars.

One of the leading contenders for the origin of supermassive black holes at $z \gtrsim$ 7 is catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 15. In this scenario, a few protogalaxies form in the presence of strong Lyman-Werner UV backgrounds that quench H$_2$ formation in their constituent halos, preventing them from forming stars or blowing heavy elements into the intergalactic medium prior to formation. At masses of 10$^ 8$ \Ms\ and virial temperatures of 10$^4$ K, gas in these halos rapidly cools by H lines, in some cases forming 10$^4$ - 10$^6$ \Ms\ Pop III stars and, a short time later, the seeds of supermassive black holes. Instead of collapsing directly to black holes some of these stars died in the most energetic thermonuclear explosions in the universe. We have modeled the explosions of such stars in the dense cores of line-cooled protogalaxies in the presence of cosmological flows. In stark contrast to the explosions in diffuse regions in previous simulations, these SNe briefly engulf the protogalaxy but then collapse back into its dark matter potential. Fallback drives turbulence that efficiently distributes metals throughout the interior of the halo and fuels the rapid growth of nascent black holes at its center. The accompanying starburst and x-ray emission from these line-cooled galaxies easily distinguish them from more slowly evolving neighbors and might reveal the birthplaces of supermassive black holes on the sky.

The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova explosions, which may soon be detected by a new generation of NIR observatories such as JWST and WFIRST. We present light curves and spectra for Pop III pair-instability supernovae calculated with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan-STARRS will find them at z ~ 7 - 8. The discovery of these ancient explosions will probe the first stellar populations and reveal the existence of primitive galaxies that might not otherwise have been detected.

Supermassive primordial stars are expected to form in a small fraction of massive protogalaxies in the early universe, and are generally conceived of as the progenitors of the seeds of supermassive black holes (BHs). Supermassive stars with masses of ~55,000 M_Sun, however, have been found to explode and completely disrupt in a supernova (SN) with an energy of up to ~10^55 erg instead of collapsing to a BH. Such events, ~10,000 times more energetic than typical SNe today, would be among the biggest explosions in the history of the universe. Here we present a simulation of such a SN in two stages. Using the RAGE radiation hydrodynamics code we first evolve the explosion from an early stage through the breakout of the shock from the surface of the star until the blast wave has propagated out to several parsecs from the explosion site, which lies deep within an atomic cooling dark matter (DM) halo at z ~ 15. Then, using the GADGET cosmological hydrodynamics code we evolve the explosion out to several kiloparsecs from the explosion site, far into the low-density intergalactic medium. The host DM halo, with a total mass of 4 x 10^7 M_Sun, much more massive than typical primordial star-forming halos, is completely evacuated of high density gas after < 10 Myr, although dense metal-enriched gas recollapses into the halo, where it will likely form second-generation stars with metallicities of ~ 0.05 Z_Sun after > 70 Myr. The chemical signature of supermassive star explosions may be found in such long-lived second-generation stars today.

The first primitive galaxies formed from accretion and mergers by z ~ 15, and were primarily responsible for cosmological reionization and the chemical enrichment of the early cosmos. But a few of these galaxies may have formed in the presence of strong Lyman-Werner UV fluxes that sterilized them of H_2, preventing them from forming stars or expelling heavy elements into the IGM prior to assembly. At masses of 10^8 Ms and virial temperatures of 10^4 K, these halos began to rapidly cool by atomic lines, perhaps forming 10^4 - 10^6 Ms Pop III stars and, later, the seeds of supermassive black holes. We have modeled the explosion of a supermassive Pop III star in the dense core of a line-cooled protogalaxy with the ZEUS-MP code. We find that the supernova (SN) expands to a radius of ~ 1 kpc, briefly engulfing the entire galaxy, but then collapses back into the potential well of the dark matter. Fallback fully mixes the interior of the protogalaxy with metals, igniting a violent starburst and fueling the rapid growth of a massive black hole at its center. The starburst would populate the protogalaxy with stars in greater numbers and at higher metallicities than in more slowly-evolving, nearby halos. The SN remnant becomes a strong synchrotron source that can be observed with eVLA and eMERLIN and has a unique signature that easily distinguishes it from less energetic SN remnants. Such explosions, and their attendant starbursts, may well have marked the birthplaces of supermassive black holes on the sky.

Recent observations of quasars powered by supermassive black holes (SMBHs) out to z > 7 constrain both the initial seed masses and the growth of the most massive black holes (BHs) in the early universe. Here we elucidate the implications of the radiative feedback from early generations of stars and from BH accretion for popular models for the formation and growth of seed BHs. We show that by properly accounting for (1) the limited role of mergers in growing seed BHs as inferred from cosmological simulations of early star formation and radiative feedback, (2) the sub-Eddington accretion rates of BHs expected at the earliest times, and (3) the large radiative efficiencies (e_rad) of the most massive BHs inferred from observations of active galactic nuclei at high redshift (e_rad > 0.1), we are led to the conclusion that the initial BH seeds may have been as massive as > 10^5 M_Sun. This presents a strong challenge to the Population III seed model, which calls for seed masses of ~ 100 M_Sun and, even with constant Eddington-limited accretion, requires e_rad < 0.09 to explain the highest-z SMBHs in today's standard LambdaCDM cosmological model. It is, however, consistent with the prediction of the direct collapse scenario of SMBH seed formation, in which a supermassive primordial star forms in a region of the universe with a high molecule-dissociating background radiation field, and collapses directly into a 10^4--10^6 M_Sun seed BH. These results corroborate recent cosmological simulations and observational campaigns which suggest that these massive BHs were the seeds of a large fraction of the SMBHs residing in the centers of galaxies today.

The detection of Pop III supernovae could directly probe the primordial IMF for the first time, unveiling the properties of the first galaxies, early chemical enrichment and reionization, and the seeds of supermassive black holes. Growing evidence that some Pop III stars were less massive than 100 solar masses may complicate prospects for their detection, because even though they would have been more plentiful they would have died as core-collapse supernovae, with far less luminosity than pair-instability explosions. This picture greatly improves if the SN shock collides with a dense circumstellar shell ejected during a prior violent LBV type eruption. Such collisions can turn even dim SNe into extremely bright ones whose luminosities can rival those of pair-instability SNe. We present simulations of Pop III Type IIn SN light curves and spectra performed with the Los Alamos RAGE and SPECTRUM codes. Taking into account Lyman-alpha absorption in the early universe and cosmological redshifting, we find that 40 solar mass Pop III Type IIn SNe will be visible out to z ~ 20 with JWST and out to z ~ 7 with WFIRST. Thus, even low mass Pop III SNe can be used to probe the primeval universe.

Understanding the properties of Pop III stars is prerequisite to elucidating the nature of primeval galaxies, the chemical enrichment and reionization of the early IGM, and the origin of supermassive black holes. While the primordial IMF remains unknown, recent evidence from numerical simulations and stellar archaeology suggests that some Pop III stars may have had lower masses than previously thought, 15 - 50 \Ms in addition to 50 - 500 \Ms. The detection of Pop III supernovae by JWST, WFIRST or the TMT could directly probe the primordial IMF for the first time. We present numerical simulations of 15 - 40 \Ms Pop III core-collapse SNe done with the Los Alamos radiation hydrodynamics code RAGE. We find that they will be visible in the earliest galaxies out to z ~ 10 - 15, tracing their star formation rates and in some cases revealing their positions on the sky. Since the central engines of Pop III and solar-metallicity core-collapse SNe are quite similar, future detection of any Type II supernovae by next-generation NIR instruments will in general be limited to this epoch.

We have entered the era of explosive transient astronomy, in which upcoming real-time surveys like the Large Synoptic Survey Telescope (LSST), the Palomar Transient Factory (PTF) and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) will detect supernovae in unprecedented numbers. Future telescopes such as the James Webb Space Telescope may discover supernovae from the earliest stars in the universe and reveal their masses. The observational signatures of these astrophysical transients are the key to unveiling their central engines, the environments in which they occur, and to what precision they will pinpoint cosmic acceleration and the nature of dark energy. We present a new method for modeling supernova light curves and spectra with the radiation hydrodynamics code RAGE coupled with detailed monochromatic opacities in the SPECTRUM code. We include a suite of tests that demonstrate how the improved physics is indispensable to modeling shock breakout and light curves.

Primordial stars are key to primeval structure formation as the first stellar components of primeval galaxies, the sources of cosmic chemical enrichment and likely cosmic reionization, and they possibly gave rise to the supermassive black holes residing at the centres of galaxies today. While the direct detection of individual Pop III stars will likely remain beyond reach for decades to come, we show their supernova remnants may soon be detectable in the radio. We calculate radio synchrotron signatures between 0.5 - 35 GHz from hydrodynamical computations of the supernova remnants of Pop III stars in minihaloes. We find that hypernovae yield the brightest systems, with observed radio fluxes as high as 1 - 10 muJy. Less energetic Type II supernovae yield remnants about a factor of 30 dimmer and pair-instability supernova remnants are dimmer by a factor of more than 10,000. Because of the high gas densities of the progenitor environments, synchrotron losses severely limit the maximum emission frequencies, producing a distinctive peaked radio spectrum distinguishable from normal galactic supernova remnant spectra. Hypernovae radio remnants should be detectable by existing radio facilities like eVLA and eMERLIN while Type II supernova remnants will require the Square Kilometre Array. The number counts of hypernova remnants at z > 20 with fluxes above 1 muJy are expected to be one per hundred square degree field, increasing to a few per square degree if they form down to z = 10. The detection of a z > 20 Type II supernova remnant brighter than 1 nJy would require a 100 - 200 square degree field, although only a 1 - 2 square degree field for those forming down to z = 10. Hypernova and Type II supernova remnants are easily separated from one another by their light curves, which will enable future surveys to use them to constrain the initial mass function of Pop III stars.

The first stars ended the cosmic Dark Ages and created the first heavy elements necessary for the formation of planets and life. The properties of these stars remain uncertain, and it may be decades before individual Pop III stars are directly observed. Their masses, however, can be inferred from their supernova explosions, which may soon be found in both deep-field surveys by JWST and in all-sky surveys by WFIRST. We have performed radiation hydrodynamical simulations of the near infrared signals of Pop III pair-instability supernovae in realistic circumstellar environments with Lyman absorption by the neutral intergalactic medium. We find that JWST and WFIRST will detect these explosions out to z ~ 30 and 20, respectively, unveiling the first generation of stars in the universe.

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.

Pop III stars are the key to the character of primeval galaxies, the first heavy elements, the onset of cosmological reionization, and the seeds of supermassive black holes. Unfortunately, in spite of their increasing sophistication, numerical models of Pop III star formation cannot yet predict the masses of the first stars. Because they lie at the edge of the observable universe, individual Pop III stars will also remain beyond the reach of telescopes for the foreseeable future, and so their properties remain unknown. However, it will soon be possible to constrain their masses by the direct detection of their supernovae and by reconciling their nucleosynthetic yields to the chemical abundances measured in ancient metal-poor stars in the Galactic halo, some of which may be bear the ashes of the first stars. Here, I review current problems on the simulation frontier in Pop III star formation and discuss the best prospects for constraining their properties observationally in the near term.

The existence of 10$^9$ M$_{\odot}$ black holes (BH) in massive galaxies by $z \ sim$ 7 is one of the great unsolved mysteries in cosmological structure formation. One theory argues that they originate from the black holes of Pop III stars at $z \sim$ 20 and then accrete at the Eddington limit down to the epoch of reionization, which requires that they have constant access to rich supplies of fuel. Because early numerical simulations suggested that Pop III stars were $\gtrsim$ 100 M$_{\odot}$, the supermassive black hole seeds considered up to now were 100 - 300 M$_{\odot}$. However, there is a growing numerical and observational consensus that some Pop III stars were tens of solar masses, not hundreds, and that 20 - 40 M$_{\odot}$ black holes may have been much more plentiful at high redshift. However, we find that natal kicks imparted to 20 - 40 M$_{\odot}$ Pop III BHs during formation eject them from their halos and hence their fuel supply, precluding them from Eddington-limit growth. Consequently, supermassive black holes are far less likely to form from low-mass Pop III stars than from very massive ones.

One of the most promising explanations for the origin of the billion solar mass black holes (BHs) inferred to power quasars at redshifts z > 6 is that supermassive stars (SMSs) with masses > 10,000 solar masses collapse to form the seed BHs from which they grow. Here we review recent theoretical advances which provide support for this scenario. Firstly, given sufficiently high accretion rates of gas into the cores of primordial protogalaxies, it appears that neither the high energy radiation emitted from the stellar surface nor the limited lifetime of SMSs can prevent their growth to masses of up to > 100,000 solar masses. Secondly, recent cosmological simulations suggest that the high fluxes of molecule-dissociating radiation which may be required in order to achieve such high accretion rates may be more common in the early universe than previously thought. We conclude that the majority of supermassive BHs may originate from SMSs at high redshifts.

The collapse of baryons into extremely massive stars with masses exceeding 10^4 M_Sun in a small fraction of protogalaxies at z > 10 is a promising candidate for the origin of supermassive black holes, some of which grow to a billion solar masses by z ~ 7. We determine the maximum masses such stars can attain by accreting primordial gas. We find that at relatively low accretion rates the strong ionizing radiation of these stars limits their masses to M_* ~ 10^3 M_Sun (dM_acc/dt / 10^-3 M_Sun yr^-1)^8/7, where dM_acc/dt is the rate at which the star gains mass. However, at the higher central infall rates usually found in numerical simulations of protogalactic collapse (>~ 0.1 M_Sun yr^-1), the lifetime of the star instead limits its final mass to >~ 10^6 M_Sun. Furthermore, for the spherical accretion rates at which the star can grow, its ionizing radiation is confined deep within the protogalaxy, so the evolution of the star is decoupled from that of its host galaxy. Lyman alpha emission from the surrounding H II region is trapped in these heavy accretion flows and likely reprocessed into strong Balmer series emission, which may be observable by the James Webb Space Telescope. This, along with strong He II 1640 Angstrom and continuum emission, are likely to be the key observational signatures of the progenitors of supermassive black holes at high redshift.

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.

The observational signatures of the first cosmic explosions and their chemical imprint on second-generation stars both crucially depend on how heavy elements mix within the star at the earliest stages of the blast. We present numerical simulations of the early evolution of Population III pair-instability supernovae with the new adaptive mesh refinement code CASTRO. In stark contrast to 15 - 40 Msun core-collapse primordial supernovae, we find no mixing in most 150 - 250 Msun pair-instability supernovae out to times well after breakout from the surface of the star. This may be the key to determining the mass of the progenitor of a primeval supernova, because vigorous mixing will cause emission lines from heavy metals such as Fe and Ni to appear much sooner in the light curves of core-collapse supernovae than in those of pair-instability explosions. Our results also imply that unlike low-mass Pop III supernovae, whose collective metal yields can be directly compared to the chemical abundances of extremely metal-poor stars, further detailed numerical simulations will be required to determine the nucleosynthetic imprint of very massive Pop III stars on their direct descendants.

Massive Population III stars from 140 - 260 solar masses ended their lives as pair-instability supernovae (PISNe), the most energetic thermonuclear explosions in the universe. Detection of these explosions could directly constrain the primordial IMF for the first time, which is key to the formation of the first galaxies, early cosmological reionization, and the chemical enrichment of the primeval IGM. We present radiation hydrodynamical calculations of Pop III PISN light curves and spectra performed with the RAGE code. We find that the initial radiation pulse due to shock breakout from the surface of the star, although attenuated by the Lyman-alpha forest, will still be visible by JWST at z ~ 10 - 15, and possibly out to z ~ 20 with strong gravitational lensing. We have also studied metal mixing at early stages of the explosion prior to breakout from the surface of the star with the CASTRO AMR code and find vigorous mixing in primordial core-collapse explosions but very little in PISNe. This implies that the key to determining progenitor masses of the first cosmic explosions is early spectroscopy just after shock breakout, and that multidimensional mixing is crucial to accurate low-mass Pop III SNe light curves and spectra.

Observations of the explosions of Population III (Pop III) stars have the potential to teach us much about the formation and evolution of these zero-metallicity objects. To realize this potential, we must tie observed emission to an explosion model, which requires accurate light curve and spectra calculations. Here, we discuss many of the pitfalls and problems involved in such models, presenting some preliminary results from radiation-hydrodynamics simulations.

We present an overview of recent work that focuses on understanding the radiative feedback processes that are potentially important during Population III star formation. Specifically, we examine the effect of the Lyman-Werner (photodissociating) background on the early stages of primordial star formation, which serves to delay the onset of star formation in a given halo but never suppresses it entirely. We also examine the effect that both photodissociating and ionizing radiation in I-fronts from nearby stellar systems have on the formation of primordial protostellar clouds. Depending on the strength of the incoming radiation field and the central density of the halos, Pop III star formation can be suppressed, unaffected, or even enhanced. Understanding these and other effects is crucial to modeling Population III star formation and to building the earliest generations of galaxies in the Universe.

The inclusion of rotationally-induced mixing in stellar evolution can alter the structure and composition of presupernova stars. We survey the effects of progenitor rotation on nucleosynthetic yields in Population III and II supernovae using the new adaptive mesh refinement (AMR) code CASTRO. We examine spherical explosions in 15, 25 and 40 solar mass stars at Z = 0 and 10^-4 solar metallicity with three explosion energies and two rotation rates. Rotation in the Z = 0 models resulted in primary nitrogen production and a stronger hydrogen burning shell which led all models to die as red supergiants. On the other hand, the Z=10^-4 solar metallicity models that included rotation ended their lives as compact blue stars. Because of their extended structure, the hydrodynamics favors more mixing and less fallback in the metal free stars than the Z = 10^-4 models. As expected, higher energy explosions produce more enrichment and less fallback than do lower energy explosions, and less massive stars produce more enrichment and leave behind smaller remnants than do more massive stars. We compare our nucleosynthetic yields to the chemical abundances in the three most iron-poor stars yet found and reproduce the abundance pattern of one, HE 0557-4840, with a zero metallicity 15 solar mass, 2.4 x 10^51 erg supernova. A Salpeter IMF averaged integration of our yields for Z=0 models with explosion energies of 2.4x10^51 ergs or less is in good agreement with the abundances observed in larger samples of extremely metal-poor stars, provided 15 solar mass stars are included. Since the abundance patterns of extremely metal-poor stars likely arise from a representative sample of progenitors, our yields suggest that low-mass supernovae contributed the bulk of the metals to the early universe.

Ionizing UV radiation and supernova flows amidst clustered minihalos at high redshift regulated the rise of the first stellar populations in the universe. Previous studies have addressed the effects of very massive primordial stars on the collapse of nearby halos into new stars, but the absence of the odd-even nucleosynthetic signature of pair-instability supernovae in ancient metal-poor stars suggests that Population III stars may have been less than 100 M$_{\odot}$. We extend our earlier survey of local UV feedback on star formation to 25 - 80 M$_{\odot}$ stars and include kinetic feedback by supernovae for 25 - 40 M$_{\odot}$ stars. We find radiative feedback to be relatively uniform over this mass range, primarily because the larger fluxes of more massive stars are offset by their shorter lifetimes. Our models demonstrate that prior to the rise of global UV backgrounds, Lyman-Werner photons from nearby stars cannot prevent halos from forming new stars. These calculations also reveal that violent dynamical instabilities can erupt in the UV radiation front enveloping a primordial halo but that they ultimately have no effect on the formation of a star. Finally, our simulations suggest that relic H II regions surrounding partially evaporated halos may expel Lyman-Werner backgrounds at lower redshifts, allowing stars to form that were previously suppressed. We provide fits to radiative and kinetic feedback on star formation for use in both semianalytic models and numerical simulations.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

The cosmic dark ages are the mysterious epoch during which the pristine gas began to condense and ultimately form the first stars. Although these beginnings have long been a topic of theoretical interest, technology has only recently allowed the beginnings of observational insight into this epoch. Many questions surround the formation of stars in metal-free gas and the history of the build-up of metals in the intergalactic medium: (1) What were the properties of the first stellar and galactic sources to form in pristine (metal-free) gas? (2) When did the epoch of Population III (metal-free) star formation take place and how long did it last? (3) Was the stellar initial mass function dramatically different for the first stars and galaxies? These questions are all active areas of theoretical research. However, new observational constraints via the direct detection of Population III star formation are vital to making progress in answering the broader questions surrounding how galaxies formed and how the cosmological properties of the universe have affected the objects it contains.

Current numerical studies suggest that the first galaxies formed a few stars at a time and were enriched only gradually by the first heavy elements. However, the large box sizes in these models cannot resolve primordial supernova explosions or the mixing of their metals with ambient gas, which could result in intervening, prompt generations of low-mass stars. We present multiscale 1D models of Population III supernovae in cosmological halos that evolve the blast from its earliest stages as a free expansion. We find that if the star ionizes the halo, the ejecta strongly interacts with the dense shell swept up by the H II region, potentially cooling and fragmenting it into clumps that are unstable to gravitational collapse. If the star fails to ionize the halo, the explosion propagates metals out to 20 - 40 pc and then collapses, enriching tens of thousands of solar masses of primordial gas, in contrast to previous models that suggest such explosions 'fizzle'. Rapid formation of low-mass stars trapped in the gravitational potential well of the halo appears to be inevitable in such circumstances. Consequently, it is possible that far more stars were swept up into the first galaxies, at earlier times and with distinct chemical signatures, than in present models. Upcoming measurements by JWST and ALMA may discriminate between these two paradigms.

We present numerical simulations of primordial supernovae in cosmological minihalos at $z \sim$ 20. We consider Type II supernovae, hypernovae, and pair instability supernovae (PISN) in halos from 6.9 $\times$ 10$^5$ - 1.2 $\times$ 10$^7$ $\Ms$, those in which Population III stars are expected to form via H$_2$ cooling. The supernovae evolve along two evolutionary paths according to whether they explode in \ion{H}{2} regions or neutral halos. Those in \ion{H}{2} regions first expand adiabatically and then radiate strongly upon collision with baryons ejected from the halo during its photoevaporation by the progenitor. Explosions in neutral halos promptly emit most of their kinetic energy as x-rays, but retain enough momentum to seriously disrupt the halo. We find that the least energetic of the supernovae are capable of destroying halos $\lesssim$ 10$^7$ $\Ms$, while a single PISN can destroy even more massive halos. Blasts in \ion{H}{2} regions disperse heavy elements into the IGM, but neutral halos confine the explosion and its metals. In \ion{H}{2} regions, a prompt second generation of stars may form in the remnant at radii of 100 - 200 pc in the halo. Explosions confined by large halos instead recollapse, with infall rates in excess of 10$^{-2}$ $\Ms$ yr$^{-1}$ that heavily contaminate their interior. This fallback may either fuel massive black hole growth at very high redshifts or create the first globular cluster with a radius of 10 - 20 pc at the center of the halo. Our findings allow the possibility that the first primitive galaxies formed sooner, with greater numbers of stars and distinct chemical abundance patterns, than in current models.

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.

We present numerical simulations of how a 120 M$_\odot$ primordial star regulates star formation in nearby cosmological halos at $z \sim$ 20 by photoevaporation. Our models include nine-species primordial chemistry and self-consistent multifrequency conservative transfer of UV photons with all relevant radiative processes. Whether or not new stars form in halos clustered around a Population III star ultimately depends on their core densities and proximity to the star. Diffuse halos with central densities below 2 - 3 cm$^{-3}$ are completely ionized and evaporated anywhere in the cluster. Evolved halos with core densities above 2000 cm$^{-3}$ are impervious to both ionizing and Lyman-Werner flux at most distances from the star and collapse as quickly as they would in its absence. Star formation in halos of intermediate density can be either promoted or suppressed depending on how the I-front remnant shock compresses, deforms and enriches the core with H$_2$. We find that the 120 M$_\odot$ star photodissociates H$_2$ in most halos in the cluster but that catalysis by H- restores it a few hundred kyr after the death of the star, with little effect on star formation. Our models exhibit significant departures from previous one-dimensional spherically-symmetric simulations, which are prone to serious errors due to unphysical geometric focusing effects.

Radiative cooling by metals in shocked gas mediates the formation of ionization front instabilities in the galaxy today that are responsible for a variety of phenomena in the interstellar medium, from the morphologies of nebulae to triggered star formation in molecular clouds. An important question in early reionization and chemical enrichment of the intergalactic medium is whether such instabilities arose in the H II regions of the first stars and primeval galaxies, which were devoid of metals. We present three-dimensional numerical simulations that reveal both shadow and thin-shell instabilities readily formed in primordial gas. We find that the hard UV spectra of Population III stars broadened primordial ionization fronts, causing H2 formation capable of inciting violent thin- shell instabilities in D-type fronts, even in the presence of intense Lyman-Werner flux. The high post- front gas temperatures associated with He ionization sustained and exacerbated shadow instabilities, unaided by molecular hydrogen cooling. Our models indicate that metals eclipsed H2 cooling in I-front instabilities at modest concentrations, from 0.001- 0.01 solar. We conclude that ionization front instabilities were prominent in the H II regions of the first stars and galaxies, influencing the escape of ionizing radiation and metals into the early universe.

We present numerical simulations of the photoevaporation of cosmological halos clustered around a 120 M$_\odot$ primordial star, confining our study to structures capable of hosting Population III star formation. The calculations include self-consistent multifrequency conservative transfer of UV photons together with nine-species primordial chemistry and all relevant radiative processes. The ultimate fates of these halos varies with central density and proximity to the central source but generally fall into one of four categories. Diffuse halos with central densities below 2 - 3 cm$^{-3}$ are completely ionized and evaporated by the central star anywhere in the cluster. More evolved halo cores at densities above 2000 cm$^{-3}$ are impervious to both ionizing and Lyman-Werner flux at most distances from the star and collapse of their cores proceeds without delay. Radiative feedback in halos of intermediate density can be either positive or negative, depending on how the I-front remnant shock both compresses and deforms the core and enriches it with H$_2$. We find that the 120 M$_\odot$ star photodissociates H$_2$ in most halos within the cluster but that catalysis by H- rapidly restores molecular hydrogen within a few hundred Kyr after the death of the star, with little delay in star formation. Our models exhibit significant departures from previous one-dimensional spherically-symmetric simulations, which are prone to serious errors due to unphysical geometric focusing effects.

Ionization front instabilities have long been of interest for their suspected role in a variety of phenomena in the galaxy, from the formation of bright rims and 'elephant trunks' in nebulae to triggered star formation in molecular clouds. Numerical treatments of these instabilities have historically been limited in both dimensionality and input physics, leaving important questions about their true evolution unanswered. We present the first three-dimensional radiation hydrodynamical calculations of both R-type and D-type ionization front instabilities in galactic environments (i.e., solar metallicity gas). Consistent with linear stability analyses of planar D-type fronts, our models exhibit many short-wavelength perturbations growing at early times that later evolve into fewer large-wavelength structures. The simulations demonstrate that both self-consistent radiative transfer and three-dimensional flow introduce significant morphological differences to unstable modes when compared to earlier two-dimensional approximate models. We find that the amplitude of the instabilities in the nonlinear regime is primarily determined by the efficiency of cooling within the shocked neutral shell. Strong radiative cooling leads to long, extended structures with pronounced clumping while weaker cooling leads to saturated modes that devolve into turbulent flows. These results suggest that expanding H II regions may either promote or provide turbulent support against the formation of later generations of stars, with potential consequences for star formation rates in the galaxy today.

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

We present results from the analysis of the optical spectra of 47 radio-selected narrow-line Seyfert 1 galaxies (NLS1s). These objects are a subset of the First Bright Quasar Survey (FBQS) and were initially detected at 20 cm (flux density limit ~1 mJy) in the VLA FIRST Survey. We run Spearman rank correlation tests on several sets of parameters and conclude that, except for their radio properties, radio-selected NLS1 galaxies do not exhibit significant differences from traditional NLS1 galaxies. Our results are also in agreement with previous studies suggesting that NLS1 galaxies have small black hole masses that are accreting very close to the Eddington rate. We have found 16 new radio-loud NLS1 galaxies, which increases the number of known radio-loud NLS1 galaxies by a factor of ~5.

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.

There has been considerable theoretical debate over whether photoionization and supernova feedback from the first Population III stars facilitate or suppress the formation of the next generation of stars. We present results from an Eulerian adaptive mesh refinement simulation demonstrating the formation of a primordial star within a region ionized by an earlier nearby star. Despite the higher temperatures of the ionized gas and its flow out of the dark matter potential wells, this second star formed within 23 million years of its neighbor's death. The enhanced electron fraction within the HII region catalyzes rapid molecular hydrogen formation that leads to faster cooling in the subsequent star forming halos than in the first halos. This "second generation" primordial protostar has a much lower accretion rate because, unlike the first protostar, it forms in a rotationally supported disk of approx. 10-100 solar masses. This is primarily due to the much higher angular momentum of the halo in which the second star forms. In contrast to previously published scenarios, such configurations may allow binaries or multiple systems of lower mass stars to form. These first high resolution calculations offer insight into the impact of feedback upon subsequent populations of stars and clearly demonstrate how primordial chemistry promotes the formation of subsequent generations of stars even in the presence of the entropy injected by the first stars into the IGM.

We simulate the ionization environment of z ~ 20 luminous objects formed within the framework of the current CDM cosmology and compute their UV escape fraction. These objects are likely single very massive stars that are copious UV emitters. We present analytical estimates as well as one--dimensional radiation hydrodynamical calculations of the evolution of these first HII regions in the universe. The initially D--type ionization front evolves to become R--type within $\lesssim 10^5$ yrs at a distance $\sim1$ pc. This ionization front then completely overruns the halo, accelerating an expanding shell of gas outward to velocities in excess of 30 km s$^{-1}$, about ten times the escape velocity of the confining dark matter halo. We find that the evolution of the HII region depends only weakly on the assumed stellar ionizing luminosities. Consequently, most of the gas surrounding the first stars will leave the dark halo whether or not the stars produce supernovae. If they form the first massive seed black holes these are unlikely to accrete within a Hubble time after they formed until they are incorporated into larger dark matter halos that contain more gas. Because these I--fronts exit the halo on timescales much shorter than the stars' main sequence lifetimes their host halos have UV escape fractions of $\gtrsim 0.95$, fixing an important parameter for theoretical studies of cosmological hydrogen reionization.