Engineering topics

Whalen, Daniel J.

Publications and source records attributed to Whalen, Daniel J..

At least 19 records

DOE OSTI2023

Multidimensional Radiation Hydrodynamics Simulations of Pulsational Pair-instability Supernovae

Stars with masses of 80–130 M ⊙ can encounter pulsational pair-instability at the end of their lives, which triggers consecutive episodes of explosive burning that eject multiple massive shells. Collisions between these shells produce bright transients known as pulsational pair-instability supernovae (PPI SNe) that may explain some extreme supernovae. In this paper, we present the first 2D and 3D radiation hydrodynamics simulations of PPI SNe with the CASTRO code. Radiative cooling causes the collided shells to evolve into thin, dense structures with hot spots that can enhance the peak luminosity of the SN by factors of 2–3. The light curve peaks at 1.9–2.1 × 10 43 erg s -1 for 50 days and then plateaus at 2–3 × 10 42 erg s -1 for 200 days, depending on the viewing angle. The presence of 12 C and 16 O and the absence of 28 Si and 56 Fe in its spectra can uniquely identify this transient as a PPI SN in follow-up observations. Our models suggest that multidimensional radiation hydrodynamics is required to model the evolution and light curves of all shell-collision SNe, such as Type IIne, not just PPI SNe.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Radiation Hydrodynamical Simulations of the First Quasars

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 × 10 9 ${M}_{\odot }$ black hole (BH) at z = 7.1, and ULAS J1342+0928, a 8.0 × 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 which show that J1120 and J1342 can form from direct collapse black holes 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

On the Maximum Mass of Accreting Primordial Supermassive Stars

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 M {sub ⊙} yr{sup −1} 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 M {sub ⊙} for accretion rates of 0.1–10 M {sub ⊙} yr{sup −1}, 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Magnetar-powered Supernovae in Two Dimensions. II. Broad-line Supernovae Ic

Nascent neutron stars (NSs) with millisecond periods and magnetic fields in excess of 10{sup 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 NS and release it as a disk wind or a jet with energies greater than 10{sup 52} erg over ∼20 s. What fraction 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 M{sub ⊙} of {sup 56}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 of M{sub b}∼−16.5 due to low {sup 56}Ni production.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Low-energy Population III supernovae and the origin of extremely metal-poor stars

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 needed 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 M ⊙ explosion and a 60 M ⊙ 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 56 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 remain in the fossil abundance record

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Finding the first cosmic explosions. IV. 90–140 M{sub ⊙} pair-instability supernovae

Population III stars that die as pair-instability supernovae are usually thought to fall in the mass range of 140–260 M{sub ⊙}. 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 M{sub ⊙}. 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 M{sub ⊙} in the primordial universe as 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 ∼ 1–2 and to James Webb Space Telescope and the 30 m class telescopes in the near-infrared (NIR) out to z ∼ 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

EMISSION FROM PAIR-INSTABILITY SUPERNOVAE WITH ROTATION

Pair-instability supernovae (PISNe) have been suggested as candidates for some superluminous 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 M {sub ☉} rotating PISNe covering both the shock breakout 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 PISNe, we find that for a range of progenitor masses and rotation rates these events have intrinsically red colors in contradiction with observations of superluminous supernovae. The spectroscopic properties of rotating PISNe 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 PISNe.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

POPULATION III HYPERNOVAE

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. 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 {sub ☉} 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 and z = 4-5 to the Wide-Field Infrared Survey Telescope, 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Pair-instability supernovae in the local universe

The discovery of 150-300 M {sub ☉} 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 {sub ☉} Z ∼ 0.1-0.4 Z {sub ☉} 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 ∼ 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 ∼100 M {sub ☉}.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The general relativistic instability supernova of a supermassive population III star

The formation of supermassive Population III stars with masses ≳10,000 M{sub ☉} in primeval galaxies in strong ultraviolet backgrounds at z ∼ 15 may be the most viable pathway to the formation of supermassive black holes by z ∼ 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. However, we have now discovered that non-rotating primordial stars with masses close to 55,000 M{sub ☉} can instead die as highly energetic thermonuclear supernovae powered by explosive helium burning, releasing up to 10{sup 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 ≲ 20 to Euclid and the Wide-Field Infrared Survey Telescope, perhaps signaling the birth of supermassive black hole seeds and the first quasars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The first gamma-ray bursts in the universe

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 ∼ 20, the end of the cosmic Dark Ages? Previous studies of Population III (Pop III) GRBs only considered explosions in the diffuse relic H II regions of their progenitors or bursts that are far more energetic than those observed to date. However, 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 Pop III GRB afterglows in H II regions, winds, and dense shells ejected by the star during the processes that produce the burst. We find that GRBs with E {sub iso,γ} = 10{sup 51}-10{sup 53} erg will be visible at z ≳ 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 Pop III GRBs are visible to Swift and the Very Large Array now, the optimal strategy for their detection will be future missions like the proposed EXIST and JANUS missions with 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 their redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Finding the first cosmic explosions. III. Pulsational pair-instability supernovae

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 pulsational pair-instability 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 {sub ☉} pulsational pair-instability explosion done with the Los Alamos radiation hydrodynamics code Radiation Adaptive Grid Eulerian. We find that collisions between consecutive pulsations are visible in the near infrared out to z ∼ 15-20 and can probe the earliest stellar populations at cosmic dawn.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Supermassive population III supernovae and the birth of the first quasars

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 that form 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{sup 4}-10{sup 5} M {sub ☉} black holes. We have now discovered that some of these fragments can instead become supermassive stars that eventually explode as thermonuclear supernovae (SNe) with energies of ∼10{sup 55} erg, the most energetic explosions in the universe. We have calculated light curves and spectra for supermassive Pop III SNe with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible in near-infrared 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

FINDING THE FIRST COSMIC EXPLOSIONS. I. PAIR-INSTABILITY SUPERNOVAE

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 (Pop 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 near-IR (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 intergalactic medium 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

THE BIGGEST EXPLOSIONS IN THE UNIVERSE. II

One of the leading contenders for the origin of supermassive black holes (SMBHs) at z ∼> 7 is catastrophic baryon collapse in atomically cooled halos at z ∼ 15. In this scenario, a few protogalaxies form in the presence of strong Lyman-Werner UV backgrounds that quench H{sub 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{sup 8} M{sub ☉} and virial temperatures of 10{sup 4} K, gas in these halos rapidly cools by H lines, in some cases forming 10{sup 4}-10{sup 6} M{sub ☉} Population III stars and, a short time later, the seeds of SMBHs. Instead of collapsing directly to black holes (BHs), 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 supernovae 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 BHs 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 SMBHs on the sky.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

THE BIGGEST EXPLOSIONS IN THE UNIVERSE

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{sub ☉}, however, have been found to explode and completely disrupt in a supernova (SN) with an energy of up to ∼10{sup 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 × 10{sup 7} M{sub ☉}, 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{sub ☉} after ∼> 70 Myr. The chemical signature of supermassive star explosions may be found in such long-lived second-generation stars today.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

THE SUPERNOVA THAT DESTROYED A PROTOGALAXY: PROMPT CHEMICAL ENRICHMENT AND SUPERMASSIVE BLACK HOLE GROWTH

The first primitive galaxies formed from accretion and mergers by z {approx} 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{sub 2}, preventing them from forming stars or expelling heavy elements into the intergalactic medium prior to assembly. At masses of 10{sup 8} M{sub Sun} and virial temperatures of 10{sup 4} K, these halos began to rapidly cool by atomic lines, perhaps forming 10{sup 4}-10{sup 6} M{sub Sun} 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 {approx}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.

79 ASTRONOMY AND ASTROPHYSICS↗