Acta Polytechnica doi:10.14311/AP.2013.53.0573 Acta Polytechnica 53(Supplement):573–578, 2013 © Czech Technical University in Prague, 2013 available online at http://ojs.cvut.cz/ojs/index.php/ap THE FIRST STARS Daniel J. Whalen∗ McWilliams Fellow, Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213 ∗ corresponding author: dwhalen@lanl.gov Abstract. 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 also lie at the edge of the observable universe, individual Pop III stars will remain beyond the reach of observatories for decades to come, and so their properties are unknown. However, it will soon be possible to constrain their masses by 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 bear the ashes of the first stars. Here, I review the state of the art in numerical simulations of primordial stars and attempts to directly and indirectly constrain their properties. Keywords: early universe, galaxies: high redshift, stars: early type, supernovae: general, radiative transfer, hydrodynamics, shocks. 1. The Simulation Frontier Unlike with star formation in the Galaxy today, there is little disagreement over the initial conditions of star formation in the primordial universe. The original numerical simulations suggested that Pop III stars formed in small pregalactic structures known as cos- mological minihalos at z ∼ 20 ÷ 30, or ∼ 200 Myr after the Big Bang [1–4, 24]. These models predicted that the stars formed in isolation, one per halo, and that they were 100 ÷ 500 M�. Pop III stars profoundly transformed the halos that gave birth to them, ex- pelling their baryons in supersonic ionized flows and later exploding as supernovae (e.g. [21, 38, 39, 44]). Radiation fronts from these stars also engulfed nearby halos, either promoting or suppressing star formation in them, thereby regulating the rise of the first stellar populations [11, 15, 30, 33, 34, 40, 41]. The original estimates of Pop III stellar masses were not obtained by modeling the actual formation and evolution of the stars. They were derived by compar- ing infall rates at the center of the halo at early stages of collapse to Kelvin–Helmholz contraction times to place upper limits on the final mass of the star. Later simulation campaigns in the same vein revealed a much broader range of final masses for Pop III stars, 30 ÷ 300 M� [27], and that they could form as bina- ries in a fraction of the halos [35]. Heroic numerical efforts have only recently achieved the formation of a hydrostatic protostar at the center of the halo [51] and the collapse of the central flow into an accretion disk [6, 9, 10, 31, 32]. In particular, the disk calculations indicate that they are unstable to fragmentation, raising the possi- bility that Pop III stars may have only been tens of solar masses, not hundreds, and that they may have Figure 1. The formation and fragmentation of a Pop III protostellar disk in the Arepo code [10]. formed in small swarms of up to a dozen at the centers of primeval halos. Computer models of ionizing UV breakout in the final stages of Pop III protostellar disks have also found that the I-front of the nascent star exits the disk in bipolar outflows that terminate accretion onto the star and mostly evaporate the disk by the time the star reaches ∼ 40 M� [13]. This result reinforces the sentiments of some in the community that while the Pop III IMF was top-heavy, primordial stars may only have been 10 ÷ 40 �. 573 http://dx.doi.org/10.14311/AP.2013.53.0573 http://ojs.cvut.cz/ojs/index.php/ap Daniel J. Whalen Acta Polytechnica 1.1. High-mass or low-mass Pop III stars? In spite of their increasing sophistication, these sim- ulations should be taken to be very preliminary for several reasons. First, Pop III accretion disks form in smoothed-particle hydrodynamics (SPH) models but not in adaptive mesh refinement (AMR) simula- tions, although the AMR models have not evolved the collapse of the halo to the times achieved by SPH calculations. This raises the question of whether one technique better captures the transport of angular momentum out of the center of the cloud than the other, and whether accretion is ultimately spherical or through a disk. Second, the stability of the disk itself remains an open question because although the simu- lations can now fully resolve the disk they do not yet incorporate all of its relevant physics. In particular, they lack high-order radiation transport, which regu- lates the thermal state of the disk and its tendency to fragment. Furthermore, the role of primordial mag- netic fields in the formation and evolution of the disk is not well understood [36, 48]. The use of sink particles to represent disk fragments in the original SPH simulations of Pop III protostellar disk formation called into question the longevity of the fragments. Once they are created in the simulations they are never destroyed, unlike real fragments that could be torn apart by gravitational torques and vis- cous forces [25]. More recent moving mesh simulations performed with the Arepo code that do not rely on sink particles find that the fragments persist but only evolve the disk for 10÷20 yr [10]. Perhaps most impor- tantly, no simulation has followed the disks for more than a few centuries, far short of the time required to assemble a massive star. Thus, it remains unclear if the fragments in the disk remain distinct objects or merge with the largest one at the center, building it up into a very massive star over time through protracted, clumpy accretion. 1.2. Accretion cutoff and the final masses of Pop III stars This latter point directly impacts estimates of final masses for Pop III stars inferred from numerical sim- ulations that attempt to model ho ionizing UV from the star reverses infall and evaporates the accretion disk. At the heart of such models is a simple recipe for the evolution of the protostar that provides a pre- scription for its radius and luminosity as a function of time and acquired mass. The Hosokawa et al. [13] 2D calculations take the growth of the protostar to be relatively steady, in which case it contracts and settles onto the main sequence at ∼ 30 M�. At this point the star becomes extremely luminous in ionizing UV radiation that halts accretion onto the star in a few hundred kyr at a final mass of ∼ 40 M�. If accretion instead turns out to be clumpy, the protostar could remain puffy and cool and reach much larger masses before burning off the disk. A finer point is that all current accretion cutoff simulations evolve both radiation and hydrodynamics on the Courant time, a practice which is known to lead to serious inaccuracies in I-front propagation in density gradients [43]. Such coarse time steps may result in premature I-front breakout and accretion cut- off, and hence underestimates of the final mass of the star. Three dimensional simulations with more accu- rate radiation–matter coupling schemes, both steady and clumpy accretion scenarios, more realistic pre- scriptions for protostellar evolution based on nucle- osynthesis codes such as KEPLER [37, 49] and a variety of halo environments may better constrain the Pop III IMF. However, in judging the power of such simulations to model the masses of the first stars, it should be remembered that no simulations realisti- cally bridge the gap in time between the formation and fragmentation of a protostellar disk and its photo- evaporation up to a Myr later. We note in passing that fragments can also stop accreting if they are ejected from the disk by 3-body gravitational effects [9, 17]. These fragments could become very low-mass Pop III stars (∼ 1 M�); if so, some of them may live today. 2. Constraining the Pop III IMF with stellar archaeology Unfortunately, because they lie at the edge of the observable universe, individual Pop III stars will re- main beyond the reach of direct detection for decades to come, even with their enormous luminosities [29] and the advent of the next generation of near in- frared (NIR) observatories such as the James Webb Space Telescope (JWST) and the Thirty-Meter Tele- scope (TMT). However, there have been attempts to indirectly constrain the masses of Pop III stars by comparing the cumulative elemental yield of their su- pernovae to the fossil chemical abundances found in ancient metal-poor stars in the Galactic halo, some of which may be second-generation stars. Stellar evo- lution models indicate that 15 ÷ 40 M� primordial stars died in core collapse (CC) supernovae (SNe) and that 40 ÷ 140 M� stars collapsed to black holes, perhaps with violent pulsational mass loss prior to death [12]. Pop III stars between 140 and 260 M� can die in pair-instability (PI) SNe, with energies up to 100 times those of Type Ia SNe that completely unbind the star and leave no compact remnant (Chat- zopoulos & Wheeler 2012 have recently discovered that rotating Pop III stars down to 65 M� can also die as PI SNe). These explosions were the first great nucleosynthetic engines of the universe, expelling up to half the mass of the progenitor in heavy elements into the early IGM. Primordial stars above 260 M� collapsed directly to black holes, with no mass loss. Joggerst et al. [18] recently calculated the chemical imprint of low-mass Pop III SNe on later generations of stars by modeling mixing and fallback onto the central black hole in 15 ÷ 40 M� Pop III core collapse explosions with the CASTRO AMR code. As shown 574 vol. 53 supplement/2013 The First Stars Figure 2. Comparing Pop III SN yields to the chemical abundances of three of the most metal-poor stars (left panel) and the extremely metal-poor (EMP) stars in the Cayrel et al. [5] and Lai et al. [22] surveys (right panel). In the left panel, the abundances in HE0557-4840 agree well with the yields from SN model z15G in the Joggerst et al. [18] study. In the right panel we show that higher explosion energy rotating Z = 0 stars reproduce EMP abundances well. The existence of 15 M� Pop III stars is required to produce this good agreement with observations. in Fig. 2, a simple power-law IMF average of the elemental yields of these explosions is in good agree- ment with the fossil abundances in a sample of 130 extremely metal poor stars with Z < 10−4 Z� [5, 22]. Although these results suggest that low-mass Pop III stars shouldered the bulk of the chemical enrichment of the early IGM, 40 ÷ 60 M� hypernova explosions, whose energies are intermediate to those of CC and PI SNe, may also have contributed metals at high redshifts [16]. To date, the telltale odd-even nucleosynthetic signa- ture of PI SNe has not been found in the fossil abun- dance record, leading some to assert that Pop III stars could not have been very massive. However, the odd- even effect may have been masked by observational bias in previous surveys [19]. Reconciling Pop III SN yields to the elemental patterns in metal-poor stars is still in its infancy for several reasons. First, only small numbers of extremely metal-poor stars have been dis- covered to date, and larger sample sizes would better constrain early SN yields. Second, measurements of some elements in low-metallicity stars are challeng- ing and in the past have been subject to systematic error. Finally, there are many intervening hydrody- namical processes between the expulsion of the first metals and their uptake into new stars that are not yet understood. 3. Finding the first cosmic explosions Detection of Pop III SNe would unambiguously probe the masses of primordial stars for the first time. Since these explosions are 100,000 times brighter than either their progenitors or the primitive galaxies that host them, they could be found by JWST or the Wide- Field Infrared Survey Telescope (WFIRST). However, unlike the Type Ia SNe used to constrain cosmic accel- eration, light from primeval supernovae must traverse the vast cosmic web of neutral hydrogen that filled the universe prior to the epoch of reionization. Lyman absorption by this hydrogen removes or scatters most of the light from ancient supernovae out of our line of sight, obscuring them. Whalen et al. [45–47] have calculated JWST NIR light curves for Pop III pair-instability SNe with the Los Alamos RAGE and SPECTRUM codes [7, 8], which are shown for z = 10, 15, 20 and 30 in Fig. 3. These simulations include radiation hydrodynamical calculations of the SN light curve and spectra in the local frame, cosmological redshifting, and Lyman ab- sorption by intergalactic hydrogen. JWST detection limits at 2 ÷ 4 µm are AB magnitude 31 ÷ 32, so it is clear that JWST will be able to detect the first cosmic explosions in the universe if they are PI SNe (and even perform spectrometry on them). Even given JWST’s very narrow fields of view at high redshifts, recent calculations indicate that at least a few PI SNe should be present in any given JWST survey [14]. Also, be- cause WFIRST detection limits will be AB magnitude 26.5 at 2.2 µm, it is clear from Fig. 3 that Pop III PI SNe will be visible to WFIRST out to z ∼ 15 ÷ 20. Since it is an all-sky survey, and because this redshift range may favor the formation of very massive Pop III stars because of the rise of Lyman-Werner UV backgrounds [28], WFIRST will detect much larger numbers of Pop III SNe. 575 Daniel J. Whalen Acta Polytechnica 0 200 400 600 800 1000 time (days) 36 34 32 30 28 26 24 M a g ( F 2 0 0 ) 150 M O • 175 M O • 200 M O • 225 M O • 250 M O • z = 10 0 200 400 600 800 1000 time (days) 36 34 32 30 28 26 24 M a g ( F 2 0 0 ) 150 M O • 175 M O • 200 M O • 225 M O • 250 M O • z = 15 0 200 400 600 800 1000 time (days) 36 34 32 30 28 26 24 M a g ( F 2 7 7 ) 150 M O • 175 M O • 200 M O • 225 M O • 250 M O • z = 20 0 200 400 600 800 1000 time (days) 36 34 32 30 28 26 24 M a g ( F 3 5 6 ) 150 M O • 175 M O • 200 M O • 225 M O • 250 M O • z = 30 Figure 3. Pop III PI SN light curves for the JWST NIRCam. Clockwise from the upper left panel, the redshifts are z = 10, 15, 20 and 30. The optimum filter for each redshift is noted on the y-axis labels and the times on the x-axes are for the observer frame. F200, F277 and F356 are at 2.0, 2.77 and 3.56 µm, respectively. Could Pop III SNe be detected at later stages by other means? Whalen et al. [40], found that most of the energy of Pop III CC SNe is eventually radiated away as H and He lines as the remnant sweeps up and shocks surrounding gas. At later epochs this en- ergy would instead be lost to fine structure cooling by metals. In both cases the emission is too dim, redshifted and drawn out over time to be detected by any upcoming instruments. However, PI SNe deposit up to half of their energy into the cosmic microwave background (CMB) by inverse Compton scattering at z ∼ 20 [21, 40] and could impose excess power on the CMB at small scales [26]. The resolution of current ground-based CMB telescopes such as the At- acama Cosmology Telescope and South Pole Telescope approaches that required to directly image Sunyaev- Zeldovich (SZ) fluctuations from individual Pop III PI SN remnants, so future observatories may detect them. Unlike PI SNe, CC SNe deposit little of their energy into the CMB at z ∼ 20, and even less at lower redshifts because the density of CMB photons falls with cosmological expansion. The extreme NIR luminosities of primordial PI SNe could contribute to a NIR background excess, as has been suggested for Pop III stars themselves, i.e. [20]. New calculations reveal that enough synchrotron emis- sion from CC SN remnants is redshifted into the 21 cm band above z ∼ 10 to be directly detected by the Square Kilometer Array (SKA) [23]. Somewhat more energetic hypernovae could be detected by existing fa- cilities such as the Extended Very-Large Array (eVLA) and eMERLIN. PI SN remnants generally expand into ambient media that are too diffuse to generate a detectable synchrotron signal. Pop III SN event rates make it unlikely that they will be found in absorption at 21 cm at z > 10. The detection of the first cosmic explosions will be one of the most spectacular discoveries in extragalac- tic astronomy in the coming decade, opening our first observational window on the era of first light and the 576 vol. 53 supplement/2013 The First Stars end of the cosmic Dark Ages at z ∼ 30. They will unveil the nature of primordial stars and constrain sce- narios for early cosmological reionization, the process whereby the universe was gradually transformed from a cold, dark, featureless void into the vast, hot, ionized expanse of galaxies we observe today. At somewhat lower redshifts (z ∼ 10 ÷ 15), detections of Pop III supernovae will probe the era of primitive galaxy for- mation, marking the positions of nascent galaxies on the sky that might otherwise have eluded detection by JWST. Finally, finding the first supernovae could also reveal the masses of the seeds of the supermassive black holes lurking at the centers of massive galaxies today. Acknowledgements DJW was supported by the Bruce and Astrid McWilliams Center for Cosmology at Carnegie Mellon University. Work at LANL was done under the auspices of the National Nu- clear Security Administration of the U.S. Dept of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396. All numerical simulations were performed on Institutional Computing (IC) and Yellow network platforms at LANL (Conejo, Lobo and Yellowrail). References [1] Abel, T., Bryan, G. & Norman, M. L.: 2000, 540, 39 [2] Abel, T., Bryan, G. & Norman, M. L.: 2002, Science, 295, 93 [3] Bromm, V., Coppi, P. S. and Larson, R. B.: 1999, ApJL, 527, L5 [4] Bromm, V., Coppi, P. S. and Larson, R. B.: 2002, ApJ, 564, 23 [5] Cayrel, R. et al.: 2004, A&A, 416, 1117 [6] Clark, P. C. et al.: 2011, Science, 331, 1040 [7] Frey, L. H. et al.: 2012, ApJS, submitted, arXiv:1203.5832 [8] Gittings, M. et al.: 2008, Comp. Sci. & Disc., 1, 015005 [9] Greif, T. H. et al.: 2011, ApJ, 737, 75 [10] Greif, T. H. et al.: 2012, MNRAS, 424, 399 [11] Hasegawa, K. et al.: 2009, MNRAS, 445 [12] Heger, A. & Woosley, S. E.: 2002, ApJ, 567, 532 [13] Hosokawa, T. et al.: 2011, Science, 334, 1250 [14] Hummel, J. et al.: 2012, ApJ, submitted, arXiv:1112.5207 [15] Iliev, I. et al.: 2005, MNRAS, 361, 405 [16] Iwamoto, N. et al.: 2005, Science, 309, 451 [17] Johnson, J. L. & Khochfar, S.: 2011, MNRAS, 413, 1184 [18] Joggerst, C. C. et al.: 2010, ApJ, 709, 11 [19] Karlsson, T., Johnson, J. L. & Bromm, V.: 2008, ApJ, 679, 6 [20] Kashlinsky, A. et al.: 2005, Nature, 438, 45 [21] Kitayama K. & Yoshida, N.: 2005, ApJ, 630, 675 [22] Lai, D. K. et al.: 2008, ApJ, 681, 1524 [23] Meiksen, A. & Whalen, D. J.: 2012, MNRAS, submitted, arXiv:1209:1915 [24] Nakamura, F. & Umemura, M.: 2001, ApJ, 548, 19 [25] Norman, M. L.: 2010, in First Stars and Galaxies: Challenges in the Coming Decade, AIP Conf. Ser. 1294, 17 [26] Oh, S. P., Cooray, A. & Kamionkowski, M.: 2003, MNRAS, 342, L20 [27] O’Shea, B. L. & Norman, M. L.: 2007, ApJ, 654, 66 [28] O’Shea, B. L. & Norman, M. L.: 2008, ApJ, 673, 14 [29] Schaerer, D.: 2002, A&A, 382, 28 [30] Shapiro, P. R. et al.: 2004, MNRAS, 348, 753 [31] Smith, R. J. et al.: 2011, MNRAS, 414, 3633 [32] Stacy, A. et al.: 2010, MNRAS, 403, 45 [33] Susa, H. & Umemura, M.: 2006, ApJL, 645, L93 [34] Susa, H. et al.: 2009, ApJ, 702, 480 [35] Turk, M. J. et al.: 2009, Science, 325, 601 [36] Turk, M. J. et al.: 2012, ApJ, 745, 154 [37] Weaver, T. A., Zimmerman, G. B. & Woosley, S. E.: 1978, ApJ, 225, 1021 [38] Whalen, D. J., Abel, T. & Norman, M. L.: 2004, ApJ, 610, 14 [39] Whalen, D. J. et al.: 2008a, ApJ, 682, 49 [40] Whalen, D. J. et al.: 2008b, ApJ, 679, 925 [41] Whalen, D. J. et al.: 2010, ApJ, 712, 101 [42] Whalen, D. J. & Fryer, C. L.: 2012, ApJL, ApJL, 756, L19 [43] Whalen, D. J. & Norman, M. L.: 2006, ApJ, 162, 281 [44] Whalen, D. J. & Norman, M. L.: 2008, ApJ, 673, 664 [45] Whalen, D. J. et al.: 2012a, ApJ, submitted [46] Whalen, D. J. et al.: 2012b, ApJ, submitted [47] Whalen, D. J. et al.: 2012c, ApJ, submitted [48] Widrow, L. M. et al.: 2012, SSR, 166, 37 [49] Woosley, S. E., Heger, A. & Weaver, T. A.: 2002, Rev. Mod. Phys., 74, 1015 [50] Yorke, H. & Welz, A.: 1996, A&A, 315, 555 [51] Yoshida, N., Omukai, K & Hernquist, L.: 2008, Science, 321 577 Daniel J. Whalen Acta Polytechnica Discussion Wolfgang Kundt — You spoke of a mass range (> 40 M�) for primordial stars for which they would col- lapse to black holes. Where are the nearest of them today? Daniel Whalen — This is an excellent question. We recently published a letter (Whalen & Fryer 2012) in which we found that most 20 ÷ 40 M� Pop III black holes would be ejected from the cosmological halos that gave birth to them at velocities of 500 ÷ 1000 km s−1 by natal kicks due to asymmetries in the core-collapse engine. Such velocities are far above the escape velocity of any halo they would encounter for over a Hubble time, so there is a good chance that these black holes would be exiled to the voids between galaxies today. Black holes above 40 M� are unlikely to be born with kicks and remain in the halo, accreting and growing over cosmic time. These black holes are much more likely to reside in the galaxies into which their host halos were taken, a few of which could become the supermassive black holes found in the SDSS quasars today. Maurice Van Putten — What fraction of the gas in a cosmological halo ends up in primordial stars? Daniel Whalen — It is currently thought that the minimum halo mass for forming a Pop III star is ∼ 105 M� and that 1–10 stars are formed with masses of 30÷300 M�. Thus, a conservative estimate is that 0.1 ÷ 1 % of the baryons in the halo are converted into stars, and that the rest are evicted from the halo by strong ionized flows over the life of the stars. 578 Acta Polytechnica 53(Supplement):573–578, 2013 1 The Simulation Frontier 1.1 High-mass or low-mass Pop III stars? 1.2 Accretion cutoff and the final masses of Pop III stars 2 Constraining the Pop III IMF with stellar archaeology 3 Finding the first cosmic explosions Acknowledgements References