Source author record

Alexander Heger

Alexander Heger appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.HE astro-ph.CO astro-ph.GA astro-ph nucl-ex nucl-th physics.comp-ph hep-ph

Catalog footprint

What is connected

67works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2024arXiv

Spectacular nucleosynthesis from early massive stars

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys.

preprint2022arXiv

Long-term evolution of post-explosion Helium-star Companions of Type Iax Supernovae

Supernovae of Type Iax (SNe Iax) are an accepted faint subclass of hydrogen-free supernovae. Their origin, the nature of the progenitor systems, however, is an open question. Recent studies suggest that the weak deflagration explosion of a near-Chandrasekhar-mass white dwarf in a binary system with a helium star donor could be the origin of SNe Iax. In this scenario, the helium star donor is expected to survive the explosion. We use the one-dimensional stellar evolution codes \textsc{MESA} and \textsc{Kepler} to follow the post-impact evolution of the surviving helium companion stars. The stellar models are based on our previous hydrodynamical simulations of ejecta-donor interaction, and we explore the observational characteristics of these surviving helium companions. We find that the luminosities of the surviving helium companions increase significantly after the impact: They could vary from $2\mathord,500\,\mathrm{L_{\odot}}$ to $16\mathord,000\,\mathrm{L_{\odot}}$ for a Kelvin-Helmholtz timescale of about $10^{4}\,\mathrm{yr}$. After the star reaches thermal equilibrium, it evolves as an O-type hot subdwarf (sdO) star and continues its evolution along the evolutionary track of a normal sdO star with the same mass. Our results will help to identify the surviving helium companions of SNe Iax in future observations and to place new constraints on their progenitor models.

preprint2022arXiv

The Impact of the New $^{65\!}$As(p,$γ$)$^{66\!}$Se Reaction Rate on the Two-Proton Sequential Capture of $^{64}\!$Ge, Weak GeAs Cycles, and Type-I X-Ray Bursts such as the Clocked Burster GS 1826$-$24

We re-assess $^{65}$As(p,$γ$)$^{66}$Se reaction rates based on a set of proton thresholds of $^{66}$Se, $S_\mathrm{p}$($^{66}$Se), estimated from the experimental mirror nuclear masses, theoretical mirror displacement energies, and full $pf$-model space shell-model calculation. The self-consistent relativistic Hartree-Bogoliubov theory is employed to obtain the mirror displacement energies with much reduced uncertainty, and thus reducing the proton-threshold uncertainty up to 161 keV compared to the AME2020 evaluation. Using the simulation instantiated by the one-dimensional multi-zone hydrodynamic code, KEPLER, that closely reproduces the observed GS 1826$-$24 clocked bursts, the present forward and reverse $^{65}$As(p,$γ$)$^{66}$Se reaction rates based on a selected $S_\mathrm{p}$($^{66}$Se) = 2.469$\pm$0.054 MeV, and the latest $^{22}$Mg($α$,p)$^{25}$Al, $^{56}$Ni(p,$γ$)$^{57}$Cu(p,$γ$)$^{58}$Zn, $^{55}$Ni(p,$γ$)$^{56}$Cu, and $^{64}$Ge(p,$γ$)$^{65}$As reaction rates, we find that though the GeAs cycles is weakly established in the rapid-proton capture process path, the $^{65}$As(p,$γ$)$^{66}$Se reaction still strongly characterizes the burst tail end due to the two-proton sequential capture on $^{64}$Ge, not found by Cyburt et al. (2016) sensitivity study. The $^{65}$As(p,$γ$)$^{66}$Se reaction influences the abundances of nuclei $A$ = 64, 68, 72, 76, and 80 up to a factor of 1.4. The new $S_\mathrm{p}$($^{66}$Se) and the inclusion of the updated $^{22}$Mg($α$,p)$^{25}$Al reaction rate increases the production of $^{12}$C up to a factor of $4.5$ that is not observable and could be the main fuel for superburst. The waiting point status of and two-proton sequential capture on $^{64}$Ge, weak-cycle feature of GeAs at region heavier than $^{64}$Ge, and impact of other possible $S_\mathrm{p}$($^{66}$Se) are also discussed.

preprint2022arXiv

The Regulated NiCu Cycles with the new $^{57}$Cu(p,$γ$)$^{58}$Zn reaction rate and the Influence on Type-I X-Ray Bursts: GS 1826$-$24 Clocked Burster

During the X-ray bursts of GS 1826$-$24, "clocked burster", the nuclear reaction flow that surges through the rapid-proton capture process path has to pass through the NiCu cycles before reaching the ZnGa cycles that moderate the further extent of hydrogen burning in the region above germanium and selenium isotopes. The $^{57}$Cu(p,$γ$)$^{58}$Zn reaction located in the NiCu cycles plays an important role in influencing the burst light curves as found by Cyburt et al. (2016). We deduce the $^{57}$Cu(p,$γ$)$^{58}$Zn reaction rate based on the experimentally determined important nuclear structure information, isobaric-multiplet-mass equation, and large-scale shell model calculations. Based on the isobaric-multiplet-mass equation, we propose a possible order of $1^+_1$ and $2^+_3$ dominant resonance states and constrain the resonance energy of the $1^+_2$ state. The latter reduces the contribution of the $1^+_2$ dominant resonance state. The new reaction rate is up to a factor of four lower than the Forstner et al. (2001) rate recommended by JINA REACLIB v2.2 at the temperature regime sensitive to clocked bursts of GS 1826$-$24. Using the simulation from the one-dimensional implicit hydrodynamic code, KEPLER, to model the thermonuclear X-ray bursts of GS 1826$-$24 clocked burster, we find that the new $^{57}$Cu(p,$γ$)$^{58}$Zn coupled with the latest $^{56}$Ni(p,$γ$)$^{57}$Cu and $^{55}$Ni(p,$γ$)$^{56}$Cu reaction rates redistributes the reaction flow in the NiCu cycles and strongly influences the burst ash composition, whereas the $^{59}$Cu(p,$α$)$^{56}$Ni and $^{59}$Cu(p,$γ$)$^{60}$Zn reactions suppress the influence of the $^{57}$Cu(p,$γ$)$^{58}$Zn reaction and diminish the impact of nuclear reaction flow that by-passes the important $^{56}$Ni waiting point induced by the $^{55}$Ni(p,$γ$)$^{56}$Cu reaction on burst light curve.

preprint2022arXiv

The Role of the Hadron-Quark Phase Transition in Core-Collapse Supernovae

The hadron-quark phase transition in quantum chromodyanmics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F\_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors in the range of $14\texttt{-}100\,\text{M}_\odot$. We find explosions only for two low-compactness models ($14 \text{M}_\odot$ and $16\,\text{M}_\odot$) with the DD2F\_SF EoS, both with low explosion energies of $\mathord{\sim}10^{50}\,\mathrm{erg}$. These weak explosions are characterised by a neutrino signal with several mini-bursts in the explosion phase due to complex reverse shock dynamics, in addition to the typical second neutrino burst for phase-transition driven explosions. The nucleosynthesis shows significant overproduction of nuclei such as $^{90}\mathrm{Zr}$ for the $14\,\text{M}_\odot$ zero-metallicity model and $^{94}\mathrm{Zr}$ for the $16\,\text{M}_\odot$ solar-metallicity model, but the overproduction factors are not large enough to place constraints on the occurrence of such explosions. Several other low-compactness models using the DD2F\_SF EoS and two high-compactness models using the STOS EoS end up as failed explosions and emit a second neutrino burst. For the CMF EoS, the phase transition never leads to a second bounce and explosion. For all three EoS, inverted convection occurs deep in the core of the proto-compact star due to anomalous behaviour of thermodynamic derivatives in the mixed phase, which heats the core to entropies up to $4k_\text{B}/\text{baryon}$ and may have a distinctive gravitational wave signature, also for a second-order phase transition.

preprint2021arXiv

On the origin of nitrogen at low metallicity

Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment. We show that several heretofore unexplained features of the N/O distribution at low O/H can be explained by the N seen in metal-poor galaxies being mostly primary nitrogen that is returned to the ISM via pre-supernova winds from rapidly rotating massive stars ($M \gtrsim 10$ M$_\odot$, $v/v_{\rm crit} \gtrsim 0.4$). This mechanism naturally produces the observed N/O plateau at low O/H. We show that the large scatter in N/O at low O/H also arises naturally from variations in star-formation efficiency. By contrast, models in which the N and O come primarily from supernovae provide a very poor fit to the observed abundance distribution. We propose that the peculiar abundance patterns we observe at low O/H are a signature that dwarf galaxies retain little of their SN ejecta, leaving them with abundance patterns typical of winds.

preprint2020arXiv

A Minimum Dilution Scenario for Supernovae and Consequences for Extremely Metal-Poor Stars

To date no metal-free stars have been identified by direct observations. The most common method of constraining their properties is searching the spectra of the most metal-poor stars for the chemical elements created in the first stars and their supernova. In this approach, modelled supernova yields are compared to the observed abundance patterns in extremely metal-poor stars. The method typically only uses the abundance ratios, i.e., the yields are diluted to the observed level. Following the usual assumption of spherical symmetry we compute a simple lower limit of the mass a supernova can mix with and find that it is consistent with all published simulations of early chemical enrichment in the interstellar medium. For three different cases, we demonstrate that this dilution limit can change the conclusions from the abundance fitting. There is a large discrepancy between the dilution found in simulations of SN explosions in minihaloes and the dilution assumed in many abundance fits. Limiting the dilution can significantly alter the likelihood of which supernovae are possible progenitors of observed CEMP-no stars. In particular, some of the faint, very low-yield SNe, which have been suggested as models for the abundance pattern of SMSS0313-6708, cannot explain the measured metal abundances, as their predicted metal yields are too small by two orders of magnitude. Altogether, the new dilution model presented here emphasizes the need to better understand the mixing and dilution behaviour of aspherical SNe.

preprint2020arXiv

Helium and Nitrogen Enrichment in Massive Main Sequence Stars: Mechanisms and Implications for the Origin of WNL Stars

The evolutionary paths taken by massive stars with $M \gtrsim 60 \, \mathrm{M}_\odot$ remain substantially uncertain. They begin their lives as main sequence (MS) O-stars. Depending on their masses, rotation rates, and metallicities, they can then encounter a wide range of evolutionary states with an equally broad set of possible surface compositions and spectral classifications. We present a new grid of calculations for the evolution of such stars that covers a broad range in mass, M/M$_\odot = 60$ to $150$, rotation rate, $v \, / \, v_{\rm crit} = 0$ to $0.6$, metallicity, $[\mathrm{Fe}/\mathrm{H}] = -4$ to $0$, and $α$-element enhancement, $[α/\mathrm{Fe}] = 0$ to $0.4$. We show that rotating stars undergo rotationally-induced dredge-up of nucleosynthetic products, mostly He and N, to their surfaces while still on the MS. Non-rotating metal-rich stars also reveal the products of nucleosynthesis on their surfaces because even modest amounts of mass loss expose their "fossil" convective cores: regions that are no longer convective, but which were part of the convective core at an early stage in the star's evolution. Thus surface enhancement of He and N is expected for rotating stars at all metallicities, and for non-rotating stars if they are relatively metal-rich. We calculate a stellar atmosphere for a representative model from our grid, properly accounting for He- and N-enhancement, and show that the resulting spectrum provides a good match to observed WNL stars, strongly suggesting that the physical mechanisms we have identified are the ultimate cause of the WNL phase.

preprint2020arXiv

Multi-epoch X-ray burst modelling: MCMC with large grids of 1D simulations

Type-I X-ray bursts are recurring thermonuclear explosions on the surface of accreting neutron stars. Matching observed bursts to computational models can help to constrain system properties, such as the neutron star mass and radius, crustal heating rates, and the accreted fuel composition, but systematic parameter studies to date have been limited. We apply Markov chain Monte Carlo methods to 1D burst models for the first time, and obtain system parameter estimations for the `Clocked Burster', GS 1826$-$238, by fitting multiple observed epochs simultaneously. We explore multiple parameters which are often held constant, including the neutron star mass, crustal heating rate, and hydrogen composition. To improve the computational efficiency, we precompute a grid of 3840 KEPLER models - the largest set of 1D burst simulations to date - and by interpolating over the model grid, we can rapidly sample burst predictions. We obtain estimates for a CNO metallicity of $Z_\mathrm{CNO} = 0.010^{+0.005}_{-0.004}$, a hydrogen fraction of $X_0 = 0.74^{+0.02}_{-0.03}$, a distance of $d \sqrt{ξ_\mathrm{b}} = 6.5^{+0.4}_{-0.6}\, \mathrm{kpc}$, and a system inclination of $i = {69^{+2}_{-3}}^{\circ}$.

preprint2020arXiv

On Monolithic Supermassive Stars

Supermassive stars have been proposed as the progenitors of the massive ($\sim 10^{9}\,\rm{M}_{\odot}$) quasars observed at $z\sim7$. Prospects for directly detecting supermassive stars with next-generation facilities depend critically on their intrinsic lifetimes, as well as their formation rates. We use the 1D stellar evolution code Kepler to explore the theoretical limiting case of zero-metallicity, non-rotating stars, formed monolithically with initial masses between $10\,\rm{kM}_{\odot}$ and $190\,\rm{kM}_{\odot}$. We find that stars born with masses between $\sim60\,\rm{kM}_{\odot}$ and $\sim150\,\rm{kM}_{\odot}$ collapse at the end of the main sequence, burning stably for $\sim1.5\,\rm{Myr}$. More massive stars collapse directly through the general relativistic instability after only a thermal timescale of $\sim3\,\rm{kyr}$--$4\,\rm{kyr}$. The expected difficulty in producing such massive, thermally-relaxed objects, together with recent results for currently preferred rapidly-accreting formation models, suggests that such ``truly direct'' or ``dark'' collapses may not be typical for supermassive objects in the early Universe. We close by discussing the evolution of supermassive stars in the broader context of massive primordial stellar evolution and the possibility of supermassive stellar explosions.

preprint2020arXiv

Production of Lithium in Primordial Supernovae

The first generation of stars is quite unique. The absence of metals likely affects their formation, with current models suggesting a much more top-heavy initial mass fraction than what we observe today, and some of their other properties, such as rotation rates and binarity, are largely unknown or constrained by direct observations. But even non-rotation single stars of a given mass will evolve quite differently due to the absence of the metals: the stars will mostly remain much more compact until their death, with the hydrogen-rich later reaching down ten teems deeper in radius then in modern stars. When they explode as supernovae, the exposure to the supernova neutrino flux is much enhanced, allowing for copious production of lithium. This production will not be constant for all stars but largely vary across the mass range. Such production even more challenges the presence of the Spite Plateau.

preprint2020arXiv

The Chemical Evolution of Iron-Peak Elements with Hypernovae

We calculate the mean evolution of the iron-peak abundance ratios [(Cr,Mn,Co,Zn)/Fe] in the Galaxy, using modern supernova and hypernova chemical yields and a Galactic Chemical Evolution code that assumes homogeneous chemical evolution. We investigate a range of hypernova occurrence rates and are able to produce a chemical composition that is a reasonable fit to the observed values in metal-poor stars. This requires a hypernova occurence rate that is large (50%) in the early Universe, decreasing throughout evolution to a value that is within present day observational constraints (>~ 1%). A large hypernova occurence rate is beneficial to matching the high [Zn/Fe] observed in the most metal-poor stars, although including hypernovae with progenitor mass >= 60 solar masses is detrimental to matching the observed [(Mn,Co)/Fe] evolution at low [Fe/H]. A significant contribution from HNe seems to be critical for producing supersolar [(Co,Zn)/Fe] at low metallicity, though more work will need to be done in order to match the most extreme values. We also emphasise the need to update models for the enrichment sources at higher metallicity, as the satisfactory recovery of the solar values of [(Cr,Mn,Co,Zn)/Fe] still presents a challenge.

preprint2020arXiv

The impact of fallback on the compact remnants and chemical yields of core-collapse supernovae

Fallback in core-collapse supernovae plays a crucial role in determining the properties of the compact remnants and of the ejecta composition. We perform three-dimensional simulations of mixing and fallback for selected non-rotating supernova models to study how explosion energy and asymmetries correlate with the remnant mass, remnant kick, and remnant spin. We find that the strongest kick and spin is imparted by partial fallback in an asymmetric explosion. Black hole (BH) kicks of several hundred $\mathrm{km}\,\mathrm{s}^{-1}$ and spin parameters of $\mathord{\sim}0.25$ can be obtained in this scenario. If the initial explosion energy barely exceeds the envelope binding energy, stronger fallback results, and the remnant kick and spin remain small. If the explosion energy is high with respect to the envelope binding energy, there is little fallback with a small effect on the remnant kick, but the spin-up by fallback can be substantial. For a non-rotating $12\,\mathrm{M}_\odot$ progenitor, we find that the neutron star (NS) is spun up to millisecond periods. The high specific angular momentum of the fallback material can also lead to disk formation around black holes. Fallback may thus be a pathway towards millisecond-magnetar or collapsar-type engines for hypernovae and gamma-ray bursts that does not require rapid progenitor rotation. Within our small set of simulations, none reproduced the peculiar layered fallback necessary to explain the metal-rich iron-poor composition of many carbon-enhanced metal-poor (CEMP) stars. Models with different explosion energy and different realisations of asymmetries may, however, be compatible with CEMP abundance patterns.

preprint2019arXiv

Large-Scale Mixing in a Violent Oxygen-Neon Shell Merger Prior to a Core-Collapse Supernova

We present a seven-minute long $4π$-3D simulation of a shell merger event in a non-rotating $18.88\, M_\odot$ supernova progenitor before the onset of gravitational collapse. The key motivation is to capture the large-scale mixing and asymmetries in the wake of the shell merger before collapse using a self-consistent approach. The $4π$ geometry is crucial as it allows us to follow the growth and evolution of convective modes on the largest possible scales. We find significant differences between the kinematic, thermodynamic and chemical evolution of the 3D and the 1D model. The 3D model shows vigorous convection leading to more efficient mixing of nuclear species. In the 3D case the entire oxygen shell attains convective Mach numbers of $\mathord{\approx}\, 0.1$, whereas in the 1D model, the convective velocities are much lower and there is negligible overshooting across convective boundaries. In the 3D case, the convective eddies entrain nuclear species from the neon (and carbon) layers into the deeper part of the oxygen burning shell, where they burn and power a violent convection phase with outflows. This is a prototypical model of a convective-reactive system. Due to the strong convection and the resulting efficient mixing, the interface between the neon layer and the silicon-enriched oxygen layer disappears during the evolution, and silicon is mixed far out into merged oxygen/neon shell. Neon entrained inwards by convective downdrafts burns, resulting in lower neon mass in the 3D model compared to the 1D model at time of collapse. In addition, the 3D model develops remarkable large-scale, large-amplitude asymmetries, which may have important implications for the impending gravitational collapse and the subsequent explosion.

preprint2016arXiv

Evidence from stable isotopes and Be-10 for solar system formation triggered by a low-mass supernova

About 4.6 billion years ago, some event disturbed a cloud of gas and dust, triggering the gravitational collapse that led to the formation of the solar system. A core-collapse supernova, whose shock wave is capable of compressing such a cloud, is an obvious candidate for the initiating event. This hypothesis can be tested because supernovae also produce telltale patterns of short-lived radionuclides, which would be preserved today as isotopic anomalies. Previous studies of the forensic evidence have been inconclusive, finding a pattern of isotopes differing from that produced in conventional supernova models. Here we argue that these difficulties either do not arise or are mitigated if the initiating supernova was a special type, low in mass and explosion energy. Key to our conclusion is the demonstration that short-lived Be-10 can be readily synthesized in such supernovae by neutrino interactions, while anomalies in stable isotopes are suppressed.

preprint2016arXiv

Fast evolving pair-instability supernova models: evolution, explosion, light curves

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

preprint2016arXiv

Mass and Metallicity Requirement in Stellar Models for Galactic Chemical Evolution Applications

We used a one-zone chemical evolution model to address the question of how many masses and metallicities are required in grids of massive stellar models in order to ensure reliable galactic chemical evolution predictions. We used a set of yields that includes seven masses between 13 and 30 Msun, 15 metallicities between 0 and 0.03 in mass fraction, and two different remnant mass prescriptions. We ran several simulations where we sampled subsets of stellar models to explore the impact of different grid resolutions. Stellar yields from low- and intermediate-mass stars and from Type Ia supernovae have been included in our simulations, but with a fixed grid resolution. We compared our results with the stellar abundances observed in the Milky Way for O, Na, Mg, Si, Ca, Ti, and Mn. Our results suggest that the range of metallicity considered is more important than the number of metallicities within that range, which only affects our numerical predictions by about 0.1 dex. We found that our predictions at [Fe/H] < -2 are very sensitive to the metallicity range and the mass sampling used for the lowest metallicity included in the set of yields. Variations between results can be as high as 0.8 dex, for any remnant mass prescription. At higher [Fe/H], we found that the required number of masses depends on the element of interest and on the remnant mass prescription. With a monotonic remnant mass prescription where every model explodes as a core-collapse supernova, the mass resolution induces variations of 0.2 dex on average. But with a remnant mass prescription that includes islands of non-explodability, the mass resolution can cause variations of about 0.2 to 0.7 dex depending on the choice of metallicity range. With such a prescription, explosive or non-explosive models can be missed if not enough masses are selected, resulting in over- or under-estimations of the mass ejected by massive stars.

preprint2016arXiv

Observational Constraints on First-Star Nucleosynthesis. II. Spectroscopy of an Ultra Metal-Poor CEMP-no Star

We report on the first high-resolution spectroscopic analysis of HE0020-1741, a bright (V=12.9), ultra metal-poor ([Fe/H] = -4.1), carbon-enhanced ([C/Fe] = +1.7) star selected from the Hamburg/ESO Survey. This star exhibits low abundances of neutron-capture elements ([Ba/Fe] = -1.1), and an absolute carbon abundance A(C) = 6.1; based on either criterion, HE0020-1741 is sub-classified as a CEMP-no star. We show that the light-element abundance pattern of HE0020-1741 is consistent with predicted yields from a massive (M = 21.5 Mo), primordial composition, supernova (SN) progenitor. We also compare the abundance patterns of other ultra metal-poor stars from the literature with available measures of C, N, Na, Mg, and Fe abundances with an extensive grid of SN models (covering the mass range 10 Mo - 100 Mo), in order to probe the nature of their likely stellar progenitors. Our results suggest that at least two classes of progenitors are required at [Fe/H] < -4.0, as the abundance patterns for more than half of the sample studied in this work (7 out of 12 stars) cannot be easily reproduced by the predicted yields.

preprint2016arXiv

The Last Minutes of Oxygen Shell Burning in a Massive Star

We present the first 3D simulation of the last minutes of oxygen shell burning in an 18 solar mass supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a 1D stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ~0.1 at collapse, and an l=2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 to 0.56 solar masses due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12--24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.

preprint2015arXiv

Are Models for Core-Collapse Supernova Progenitors Consistent with the Properties of Supernova Remnants?

The recent discovery that the Fe-K line luminosities and energy centroids observed in nearby SNRs are a strong discriminant of both progenitor type and circumstellar environment has implications for our understanding of supernova progenitor evolution. Using models for the chemical composition of core-collapse supernova ejecta, we model the dynamics and thermal X-ray emission from shocked ejecta and circumstellar material, modeled as an $r^{-2}$ wind, to ages of 3000 years. We compare the X-ray spectra expected from these models to observations made with the Suzaku satellite. We also model the dynamics and X-ray emission from Type Ia progenitor models. We find a clear distinction in Fe-K line energy centroid between core-collapse and Type Ia models. The core-collapse supernova models predict higher Fe-K line centroid energies than the Type Ia models, in agreement with observations. We argue that the higher line centroids are a consequence of the increased densities found in the circumstellar environment created by the expansion of the slow-moving wind from the massive progenitors.

preprint2015arXiv

Can direct collapse black holes launch gamma-ray bursts and grow to supermassive black holes?

The existence of black holes (BHs) of mass ~ 10^{9} M_sun at z > 6 is a big puzzle in astrophysics because even optimistic estimates of the accretion time are insufficient for stellar mass BHs of ~ 10 M_sun to grow into such supermassive BHs. A resolution of this puzzle might be the direct collapse of supermassive stars with mass M ~ 10^{5} M_sun into massive seed BHs. We find that if a jet is launched from the accretion disk around the central BH, the jet can break out the star because of the structure of the radiation pressure-dominated envelope. Such ultra-long gamma-ray bursts with duration of ~ 10^{4} - 10^{6} s and flux of 10^{-11} - 10^{-8} erg s^{-1} cm^{-2} could be detectable by Swift. We estimate an event rate of < 1yr^{-1}. The total explosion energy is > 10^{55} - 10^{56} erg. The resulting negative feedback delays the growth of the remnant BH by about 70 Myr or evacuates the host galaxy completely.

preprint2015arXiv

Cosmological Impact of Population III Binaries

We present the results of the stellar feedback from Pop III binaries by employing improved, more realistic Pop III evolutionary stellar models. To facilitate a meaningful comparison, we consider a fixed mass of 60 solar masses (Msun) incorporated in Pop III stars, either contained in a single star, or split up in binary stars of 30 Msun each or an asymmetric case of one 45 Msun and one 15 Msun star. Whereas the sizes of the resulting HII regions are comparable across all cases, the HeIII regions around binary stars are significantly smaller than that of the single star. Consequently, the He$^{+}$ 1640 angstrom recombination line is expected to become much weaker. Supernova feedback exhibits great variety due to the uncertainty in possible explosion pathways. If at least one of the component stars dies as a hypernova about ten times more energetic than conventional core-collapse supernovae, the gas inside the host minihalo is effectively blown out, chemically enriching the intergalactic medium (IGM) to an average metallicity of $10^{-4}-10^{-3}$ solar metallicity (Zsun), out to $\sim 2$ kpc. The single star, however, is more likely to collapse into a black hole, accompanied by at most very weak explosions. The effectiveness of early chemical enrichment would thus be significantly reduced, in difference from the lower mass binary stars, where at least one component is likely to contribute to heavy element production and dispersal. Important new feedback physics is also introduced if close binaries can form high-mass x-ray binaries, leading to the pre-heating and -ionization of the IGM beyond the extent of the stellar HII regions.

preprint2015arXiv

Metal-Poor Stars Observed with the Magellan Telescope. III. New Extremely and Ultra Metal-Poor Stars from SDSS/SEGUE and Insights on the Formation of Ultra Metal-Poor Stars

We report the discovery of one extremely metal-poor (EMP; [Fe/H]<-3) and one ultra metal-poor (UMP; [Fe/H]<-4) star selected from the SDSS/SEGUE survey. These stars were identified as EMP candidates based on their medium-resolution (R~2,000) spectra, and were followed-up with high-resolution (R~35,000) spectroscopy with the Magellan-Clay Telescope. Their derived chemical abundances exhibit good agreement with those of stars with similar metallicities. We also provide new insights on the formation of the UMP stars, based on comparison with a new set of theoretical models of supernovae nucleosynthesis. The models were matched with 20 UMP stars found in the literature, together with one of the program stars (SDSS J1204+1201), with [Fe/H]=-4.34. From fitting their abundances, we find that the supernovae progenitors, for stars where carbon and nitrogen are measured, had masses ranging from 20.5 M_sun to 28 M_sun and explosion energies from 0.3 to 0.9x10^51 erg. These results are highly sensitive to the carbon and nitrogen abundance determinations, which is one of the main drivers for future high-resolution follow-up of UMP candidates. In addition, we are able to reproduce the different CNO abundance patterns found in UMP stars with a single progenitor type, by varying its mass and explosion energy.

preprint2015arXiv

Neutrino-Induced Nucleosynthesis in Helium Shells of Early Core-Collapse Supernovae

We summarize our studies on neutrino-driven nucleosynthesis in He shells of early core-collapse supernovae with metallicities of $Z\lesssim 10^{-3}\,Z_\odot$. We find that for progenitors of $\sim 11$--$15\,{\rm M}_\odot$, the neutrons released by $^4{\rm He}(\barν_e,e^+n)^3{\rm H}$ in He shells can be captured to produce nuclei with mass numbers up to $A \sim 200$. This mechanism is sensitive to neutrino emission spectra and flavor oscillations. In addition, we find two new primary mechanisms for neutrino-induced production of $^{9}$Be in He shells. The first mechanism produces $^9$Be via $^7{\rm Li}(n,γ)^8{\rm Li}(n,γ)^9{\rm Li}(e^-\barν_e)^9{\rm Be}$ and relies on a low explosion energy for its survival. The second mechanism operates in progenitors of $\sim 8\,{\rm M}_\odot$, where $^9$Be can be produced directly via $^7{\rm Li}(^3{\rm H},n_0)^9{\rm Be}$ during the rapid expansion of the shocked He-shell material. The light nuclei $^7$Li and $^3$H involved in these mechanisms are produced by neutrino interactions with $^4$He. We discuss the implications of neutrino-induced nucleosynthesis in He shells for interpreting the elemental abundances in metal-poor stars.

preprint2015arXiv

Nucleosynthesis in a Primordial Supernova: Carbon and Oxygen Abundances in SMSS J031300.36-670839.31

SMSS J031300.36-670839.3 (hereafter SM0313-6708) is a sub-giant halo star, with no detectable Fe lines and large overabundances of C and Mg relative to Ca. We obtained VLT-UVES spectra extending to 3060 Angstroms showing strong OH A-X band lines enabling an oxygen abundance to be derived. The OH A-X band lines in SM0313-6708 are much stronger than the CH C-X band lines. Spectrum synthesis fits indicate an [O/C] ratio of 0.02 +- 0.175. Our high S/N UVES data also enabled us to lower the Fe abundance limit to [Fe/H]{3D},NLTE < -7.52 (3 sigma). These data support our previous suggestion that the star formed from the iron-poor ejecta of a single massive star Population III supernova.

preprint2015arXiv

Pair-Instability Supernovae of Non-Zero Metallicity Stars

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.

preprint2015arXiv

The influence of accretion rate and metallicity on thermonuclear bursts: predictions from KEPLER models

Using the KEPLER hydrodynamics code, 464 models of thermonuclear X-ray bursters were performed across a range of accretion rates and compositions. We present the library of simulated burst profiles from this sample, and examine variations in the simulated lightcurve for different model conditions. We find that the recurrence time varies as a power law against accretion rate, and measure its slope while mixed H/He burning is occurring for a range of metallicities, finding the power law gradient to vary from $η= 1.1$ to $1.24$. We also identify the accretion rates at which mixed H/He burning stops and a transition occurs to different burning regimes. We also explore how varying the accretion rate and metallicity affects burst morphology in both the rise and tail.

preprint2015arXiv

The interplay between chemistry and nucleation in the formation of carbonaceous dust in supernova ejecta

Core-collapse supernovae are considered to be important contributors to the primitive dust enrichment of the interstellar medium in the high-redshift universe. Theoretical models of dust formation in stellar explosions have so far provided controversial results and a generally poor fit to the observations of dust formation in local supernovae. We present a new methodology for the calculation of carbonaceous dust formation in young supernova remnants. Our new technique uses both the nucleation theory and a chemical reaction network to allow us to compute the dust growth beyond the molecular level as well as to consider chemical erosion of the forming grains. We find that carbonaceous dust forms efficiently in the core of the ejecta, but takes several years to condensate, longer than previously estimated. It forms unevenly and remains concentrated in the inner part of the remnant. These results support the role of core-collapse supernovae as dust factories and provide new insight on the observations of SN 1987A, in which large amounts of dust have been detected to form on a timescale of years after core collapse.

preprint2015arXiv

The Remarkable Deaths of 9 - 11 Solar Mass Stars

The post-helium burning evolution of stars from 7 to 11 solar masses is complicated by the lingering effects of degeneracy and off-center ignition. Here stars in this mass range are studied using a standard set of stellar physics. Two important aspects of the study are the direct coupling of a reaction network of roughly 220 nuclei to the structure calculation at all stages and the use of a sub grid model to describe the convective bounded flame that develops during neon and oxygen burning. Below 9.0 solar masses, degenerate oxygen-neon cores form that may become either white dwarfs or electron-capture supernovae. Above 10.3 solar masses the evolution proceeds "normally" to iron-core collapse, without composition inversions or degenerate flashes. Emphasis here is upon the stars in between which typically ignite oxygen burning off center. After oxygen burns in a convectively bounded flame, silicon burning ignites in a degenerate flash that commences closer to the stellar center and with increasing violence for stars of larger mass. In some cases the silicon flash is so violent that it could lead to the early ejection of the hydrogen envelope. This might have interesting observable consequences. For example, the death of a 10.0 solar mass star could produce two supernova-like displays, a faint low energy event due to the silicon flash, and an unusually bright supernova many months later as the low energy ejecta from core collapse collides with the previously ejected envelope. The potential relation to the Crab supernova is discussed.

preprint2014arXiv

Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.

preprint2014arXiv

Effective Helium Burning Rates and the Production of the Neutrino Nuclei

Effective values for the key helium burning reaction rates, triple-alpha and 12C(alpha,gamma)16O, are obtained by adjusting their strengths so as to obtain the best match with the solar abundance pattern of isotopes uniquely or predominately made in core collapse supernovae. These effective rates are then used to determine the production of the neutrino isotopes. The use of effective rates considerably reduces the uncertainties in the production factors arising from uncertainties in the helium burning rates, and improves our ability to use the production of 11B to constrain the neutrino emission from supernovae.

preprint2014arXiv

General Relativistic Instability Supernova of a Supermassive Population III Star

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.

preprint2014arXiv

Numerical approaches for multidimensional simulations of stellar explosions

We introduce numerical algorithms for initializing multidimensional simulations of stellar explosions with 1D stellar evolution models. The initial mapping from 1D profiles onto multidimensional grids can generate severe numerical artifacts, one of the most severe of which is the violation of conservation laws for physical quantities. We introduce a numerical scheme for mapping 1D spherically-symmetric data onto multidimensional meshes so that these physical quantities are conserved. We verify our scheme by porting a realistic 1D Lagrangian stellar profile to the new multidimensional Eulerian hydro code CASTRO. Our results show that all important features in the profiles are reproduced on the new grid and that conservation laws are enforced at all resolutions after mapping. We also introduce a numerical scheme for initializing multidimensional supernova simulations with realistic perturbations predicted by 1D stellar evolution models. Instead of seeding 3D stellar profiles with random perturbations, we imprint them with velocity perturbations that reproduce the Kolmogorov energy spectrum expected for highly turbulent convective regions in stars. Our models return Kolmogorov energy spectra and vortex structures like those in turbulent flows before the modes become nonlinear. Finally, we describe approaches to determining the resolution for simulations required to capture fluid instabilities and nuclear burning. Our algorithms are applicable to multidimensional simulations besides stellar explosions that range from astrophysics to cosmology.

preprint2014arXiv

Pair Instability Supernovae of Very Massive Population III Stars

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.

preprint2014arXiv

Pair-Instability Supernovae in the Local Universe

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}$.

preprint2014arXiv

Stellar origin of the 182Hf cosmochronometer and the presolar history of solar system matter

Among the short-lived radioactive nuclei inferred to be present in the early solar system via meteoritic analyses, there are several heavier than iron whose stellar origin has been poorly understood. In particular, the abundances inferred for 182Hf (half-life = 8.9 million years) and 129I (half-life = 15.7 million years) are in disagreement with each other if both nuclei are produced by the rapid neutron-capture process. Here, we demonstrate that contrary to previous assumption, the slow neutron-capture process in asymptotic giant branch stars produces 182Hf. This has allowed us to date the last rapid and slow neutron-capture events that contaminated the solar system material at roughly 100 million years and 30 million years, respectively, before the formation of the Sun.

preprint2014arXiv

The Deaths of Very Massive Stars

The theory underlying the evolution and death of stars heavier than 10 Msun on the main sequence is reviewed with an emphasis upon stars much heavier than 30 Msun. These are stars that, in the absence of substantial mass loss, are expected to either produce black holes when they die, or, for helium cores heavier than about 35 Msun, encounter the pair instability. A wide variety of outcomes is possible depending upon the initial composition of the star, its rotation rate, and the physics used to model its evolution. These heavier stars can produce some of the brightest supernovae in the universe, but also some of the faintest. They can make gamma-ray bursts or collapse without a whimper. Their nucleosynthesis can range from just CNO to a broad range of elements up to the iron group. Though rare nowadays, they probably played a disproportionate role in shaping the evolution of the universe following the formation of its first stars.

preprint2014arXiv

Two-Dimensional Simulations of Pulsational Pair-Instability Supernovae

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.

preprint2013arXiv

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 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.

preprint2013arXiv

Finding the First Cosmic Explosions II: Core-Collapse Supernovae

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.

preprint2013arXiv

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 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.

preprint2013arXiv

Metallicity-Dependent Galactic Isotopic Decomposition for Nucleosynthesis

All stellar evolution models for nucleosynthesis require an initial isotopic abundance set to use as a starting point. Generally, our knowledge of isotopic abundances of stars is fairly incomplete except for the Solar System. We present a first model for a complete average isotopic decomposition as a function of metallicity. Our model is based on the underlying nuclear astrophysics processes, and is fitted to observational data, rather than traditional forward galactic chemical evolution modeling which integrates stellar yields beginning from big bang nucleosynthesis. We first decompose the isotopic solar abundance pattern into contributions from astrophysical sources. Each contribution is then assumed to scale as a function of metallicity. The resulting total isotopic abundances are summed into elemental abundances and fitted to available halo and disk stellar data to constrain the model's free parameter values. This procedure allows us to use available elemental observational data to reconstruct and constrain both the much needed complete isotopic evolution that is not accessible to current observations, and the underlying astrophysical processes. Our model finds a best fit for Type Ia contributing \simeq 0.7 to the solar Fe abundance, and and Type Ia onset occurring at [Fe/H]\simeq -1.2.

preprint2013arXiv

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 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.

preprint2013arXiv

The Biggest Explosions in the Universe. II

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.

preprint2013arXiv

The Impact of Helium-Burning Reaction Rates on Massive Star Evolution and Nucleosynthesis

We study the sensitivity of presupernova evolution and supernova nucleosynthesis yields of massive stars to variations of the helium-burning reaction rates within the range of their uncertainties. We use the current solar abundances from Lodders (2009) for the initial stellar composition. We compute a grid of 12 initial stellar masses and 176 models per stellar mass to explore the effects of independently varying the 12^C(a,g)16^O and 3a reaction rates, denoted R_a12 and R_3a, respectively. The production factors of both the intermediate-mass elements (A=16-40) and the s-only isotopes along the weak s-process path (70Ge, 76Se, 80Kr, 82Kr, 86Sr, and 87Sr) were found to be in reasonable agreement with predictions for variations of R_3a and R_a12 of +/-25%; the s-only isotopes, however, tend to favor higher values of R_3a than the intermediate-mass isotopes. The experimental uncertainty (one standard deviation) in R_3a(R_a12) is approximately +/-10%(+/-25%). The results show that a more accurate measurement of one of these rates would decrease the uncertainty in the other as inferred from the present calculations. We also observe sharp changes in production factors and standard deviations for small changes in the reaction rates, due to differences in the convection structure of the star. The compactness parameter was used to assess which models would likely explode as successful supernovae, and hence contribute explosive nucleosynthesis yields. We also provide the approximate remnant masses for each model and the carbon mass fractions at the end of core-helium burning as a key parameter for later evolution stage

preprint2013arXiv

The Metallicity Dependence of the Minimum Mass for Core-Collapse Supernovae

Understanding the progenitors of core collapse supernovae and their population statistics is a key ingredient for many current studies in astronomy but as yet this remains elusive. Using the MESA stellar evolution code we study the dependence of the lower mass limit for making core collapse supernovae (SNe) as function of initial stellar metallicity. We find that this mass limit is smallest at approximately [Z] = -2 with a value of ~ 8.3 Msun. At [Z] = 0 the limit is ~ 9.5 Msun and continues to rise with higher metallicity. As a consequence, for a fixed initial mass function the supernova rate may be 20% to 25% higher at [Z] = -2. This affects the association of observed SN rates as a probe for the cosmological star formation rate, rate predictions for supernova surveys, and population synthesis studies.

preprint2013arXiv

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 ~ 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.

preprint2013arXiv

Very Massive Stars (VMS) in the Local Universe

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

preprint2012arXiv

Seeing the First Supernovae at the Edge of the Universe with JWST

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.

preprint2011arXiv

Conservative Initial Mapping For Multidimensional Simulations of Stellar Explosions

Mapping one-dimensional stellar profiles onto multidimensional grids as initial conditions for hydrodynamics calculations can lead to numerical artifacts, one of the most severe of which is the violation of conservation laws for physical quantities such as energy and mass. Here we introduce a numerical scheme for mapping one-dimensional spherically-symmetric data onto multidimensional meshes so that these physical quantities are conserved. We validate our scheme by porting a realistic 1D Lagrangian stellar profile to the new multidimensional Eulerian hydro code CASTRO. Our results show that all important features in the profiles are reproduced on the new grid and that conservation laws are enforced at all resolutions after mapping.

preprint2011arXiv

Long Gamma-Ray Transients from Collapsars

In the collapsar model for common gamma-ray bursts, the formation of a centrifugally supported disk occurs during the first $\sim$10 seconds following the collapse of the iron core in a massive star. This only occurs in a small fraction of massive stellar deaths, however, and requires unusual conditions. A much more frequent occurrence could be the death of a star that makes a black hole and a weak or absent outgoing shock, but in a progenitor that only has enough angular momentum in its outermost layers to make a disk. We consider several cases where this is likely to occur - blue supergiants with low mass loss rates, tidally-interacting binaries involving either helium stars or giant stars, and the collapse to a black hole of very massive pair-instability supernovae. These events have in common the accretion of a solar mass or so of material through a disk over a period much longer than the duration of a common gamma-ray burst. A broad range of powers is possible, $10^{47}$ to $10^{50}\,$erg s$^{-1}$, and this brightness could be enhanced by beaming. Such events were probably more frequent in the early universe where mass loss rates were lower. Indeed this could be one of the most common forms of gamma-ray transients in the universe and could be used to study first generation stars. Several events could be active in the sky at any one time. A recent example of this sort of event may have been the SWIFT transient Sw-1644+57.

preprint2011arXiv

Multidimensional Simulations of Thermonuclear Supernovae from the First Stars

Theoretical models suggest that the first stars in the universe could have been very massive, with typical masses $\gtrsim$ 100 \Msun. Many of them might have died as energetic thermonuclear explosions known as pair-instability supernovae (PSNe). We present multidimensional numerical simulations of PSNe with the new radiation-hydrodynamics code CASTRO. Our models capture all explosive burning and follow the explosion until the shock breaks out from the stellar surface. 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 sec after the explosion begins. Later, when the shock reaches the hydrogen envelope a strong reverse shock forms that rapidly develops additional Rayleigh-Taylor instabilities. In red supergiant progenitors, the amplitudes of these instabilities are sufficient to mix the supernova's ejecta and alter its observational signature. Our results provide useful predictions for the detection of PSNe by forthcoming telescopes.

preprint2011arXiv

Pair Instability Supernovae: Light Curves, Spectra, and Shock Breakout

For the initial mass range (140 < M < 260 Msun) stars die in a thermonuclear runaway triggered by the pair-production instability. The supernovae they make can be remarkably energetic (up to ~10^53 ergs) and synthesize considerable amounts of radioactive isotopes. Here we model the evolution, explosion, and observational signatures of representative pair-instability supernovae (PI SNe) spanning a range of initial masses and envelope structures. The predicted light curves last for hundreds of days and range in luminosity, from very dim to extremely bright, L ~ 10^44 ergs/s. The most massive events are bright enough to be seen at high redshift, but the extended light curve duration (~1 year) -- prolonged by cosmological time-dilation -- may make it difficult to detect them as transients. An alternative approach may be to search for the brief and luminous outbreak occurring when the explosion shock wave reaches the stellar surface. Using a multi-wavelength radiation-hydrodynamics code we calculate that, in the rest-frame, the shock breakout transients of PI SNe reach luminosities of 10^45-10^46 ergs/s, peak at wavelengths ~30-170 Angstroms, and last for several hours. We explore the detectability of PI SNe emission at high redshift, and discuss how observations of the light curves, spectra, and breakout emission can be used to constrain the mass, radius, and metallicity of the progenitor.

preprint2010arXiv

Multi-Dimensional Simulations of Pair-Instability Supernovae

We present preliminary results from multidimensional numerical studies of pair instability supernova (PSN), studying the fluid instabilities that occur in multiple spatial dimensions. We use the new radiation-hydrodynamics code, CASTRO, and introduce a new mapping procedure that defines the initial conditions for the multidimensional runs in such a way that conservation of physical quantities is guaranteed at any level of resolution.

preprint2010arXiv

Production of 26Al, 44Ti and 60Fe in Supernovae-sensitivity to the helium burning rates

We have studied the sensitivity of supernova production of the gamma emitting nuclei 26Al, 44Ti and 60Fe to variations of the rates of the triple alpha and 12C(alpha, gamma) reactions. Over a range of twice their experimental uncertainties we find variations in the production of 60Fe by more than a factor of five. Smaller variations, about a factor of two to three, were observed for 26Al and 44Ti. The yields of these isotopes change significantly when the abundances of Lodders (2003) are used instead of those of Anders and Grevesse (1989). These sensitivities will limit conclusions based on a comparison of observed gamma ray intensities and stellar models until the helium burning rates are better known. Prospects for improving the helium burning rates are discussed and a new version of the Boyes rate for 12C(alpha, gamma} is presented.

preprint2010arXiv

Production of 26Al, 44Ti, and 60Fe in Core-Collapse Supernovae: Sensitivity to the Rates of the Triple Alpha and 12C(a,g)16O Reactions

We have studied the sensitivity to variations in the the triple alpha and 12C(a,g)16O reaction rates of the production of 26Al, 44Ti, and 60Fe in core-collapse supernovae. We used the KEPLER code to model the evolution of 15, 20, and 25 solar mass stars to the onset of core collapse and simulated the ensuing supernova explosion using a piston model for the explosion. Calculations were performed for the Anders and Grevesse (1989) and Lodders (2003) abundances. Over a range of twice the experimental uncertainty, sigma, for each helium-burning rate, the production of 26Al, 60Fe, and their ratio vary by factors of five or more. For some species, similar variations were observed for much smaller rate changes, 0.5 sigma or less. The production of 44Ti was less sensitive to changes in the helium-burning rates. Production of all three isotopes depended on the solar abundance set used for the initial stellar composition.

preprint2010arXiv

The Nucleosynthetic Imprint of 15-40 Solar Mass Primordial Supernovae on Metal-Poor Stars

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.

preprint2010arXiv

Two Dimensional Simulations of Pair-Instability Supernovae

We present preliminary results from two dimensional numerical studies of pair instability supernova (PSN). We study nuclear burning, hydrodynamic instabilities and explosion of very massive stars. Use a new radiation-hydrodynamics code, CASTRO.

preprint2008arXiv

Detailed Abundances for 28 Metal-poor Stars: Stellar Relics in the Milky Way

We present the results of an abundance analysis for a sample of stars with $-4<$[Fe/H]$<-2$. The data were obtained with the HIRES spectrograph at Keck Observatory. The set includes 28 stars, with effective temperature ranging from 4800 to 6600 K. For 13 stars with [Fe/H]$<-2.6$, including nine with [Fe/H]$<-3.0$ and one with [Fe/H]$=-4.0$, these are the first reported detailed abundances. For the most metal-poor star in our sample, CS 30336-049, we measure an abundance pattern that is very similar to stars in the range [Fe/H]$\sim-3.5$, including a normal C+N abundance. We also find that it has very low but measurable Sr and Ba, indicating some neutron-capture activity even at this low of a metallicity. We explore this issue further by examining other very neutron-capture-deficient stars, and find that at the lowest levels, [Ba/Sr] exhibits the ratio of the main r-process. We also report on a new r-process-enhanced star, CS 31078-018. This star has [Fe/H]$=-2.85$, [Eu/Fe]$=1.23$, and [Ba/Eu]$=-0.51$. CS 31078-018 exhibits an ``actinide boost'', i.e. much higher [Th/Eu] than expected and at a similar level to CS 31082-001. Our spectra allow us to further constrain the abundance scatter at low metallicities, which we then use to fit to the zero-metallicity Type II supernova yields of Heger & Woosley (2008). We find that supernovae with progenitor masses between 10 and 20 M$_{\odot}$ provide the best matches to our abundances.

preprint2008arXiv

First Stars III Conference Summary

The understanding of the formation, life, and death of Population III stars, as well as the impact that these objects had on later generations of structure formation, is one of the foremost issues in modern cosmological research and has been an active area of research during the past several years. We summarize the results presented at "First Stars III," a conference sponsored by Los Alamos National Laboratory, the Kavli Institute for Particle Astrophysics and Cosmology, and the Joint Institute for Nuclear Astrophysics. This conference, the third in a series, took place in July 2007 at the La Fonda Hotel in Santa Fe, New Mexico, U.S.A.

preprint2008arXiv

Mixing in Zero and Solar Metallicity Supernovae

Two-dimensional simulations of mixing and fall back in non-rotating massive stars have been carried out using realistic initial models for the presupernova star and assuming standard spherically symmetric explosions of 1.2*10^51ergs. Stars of 15 and 25 solar masses with both primordial and solar composition were modeled. The zero metallicity supernova progenitors were compact blue stars and the amount of Rayleigh-Taylor induced mixing in them was greatly reduced compared with what was seen in the red supergiants with solar metallicity. The compact zero-metal stars also experienced more fallback than their solar metallicity counterparts. As a result, the ejected nucleosynthesis from the two populations was very different. For the simple explosion model assumed, low metallicity stars ejected ejected too little iron and intermediate mass elements even to explain the abundance patterns in the most iron-poor stars found to date, suggesting that some important ingredient is missing. Rotation is likely to alter these conclusions by producing a greater fraction of red supergiants among Population III stars. The velocities of the heavy elements in all models considered - both red and blue supergiants - were less than observed in SN 1987A, suggesting that at least occasionally, asymmetric aspects of the explosion mechanism and fallback play a major role in mixing.

preprint2008arXiv

Nucleosynthesis and Evolution of Massive Metal-Free Stars

The evolution and explosion of metal-free stars with masses 10--100 solar masses are followed, and their nucleosynthetic yields, light curves, and remnant masses determined. When the supernova yields are integrated over a Salpeter initial mass function, the resulting elemental abundance pattern is qualitatively solar, but with marked deficiencies of odd-Z elements with 7 <= Z <= 13. Neglecting the contribution of the neutrino wind from the neutron stars that they make, no appreciable abundances are made for elements heavier than germanium. The computed pattern compares favorably with what has been observed in metal-deficient stars with [Z] ~< -3. Most of the stars end their lives as blue supergiants and make supernovae with distinctive light curves resembling SN 1987A, but some produce primary nitrogen by dredge up and become red supergiants. A novel automated fitting algorithm is developed for determining optimal combinations of explosion energy, mixing, and initial mass function in the large model data base to agree with specified data sets. The model is applied to the low metallicity sample of Cayrel et al. (2004) and the two ultra-iron-poor stars HE0107-5240 and HE1327-2326. Best agreement with these low metallicity stars is achieved with very little mixing, and none of the metal-deficient data sets considered show the need for a high energy explosion component. To the contrary, explosion energies somewhat less than 1.2 B seem to be preferred in most cases. (abbreviated)

preprint2008arXiv

The Impact of Neutrino Magnetic Moments on the Evolution of Massive Stars

We explore the sensitivity of massive stars to neutrino magnetic moments. We find that the additional cooling due to the neutrino magnetic moments bring about qualitative changes to the structure and evolution of stars in the mass window 7 Msun < M < 18 Msun, rather than simply changing the time scales for the burning. We describe some of the consequences of this modified evolution: the shifts in the threshold masses for creating core-collapse supernovae and oxygen-neon-magnesium white dwarfs and the appearance of a new type of supernova in which a partial carbon-oxygen core explodes within a massive star. The resulting sensitivity to the magnetic moment is at the level of (2-4) * 10^{-11} μ_B.

preprint2006arXiv

The Lack of Gamma-Ray Bursts from Population III Binaries

We study the evolution of first star (Population III) binaries. Under specific conditions, these stars may produce high redshift gamma-ray bursts (GRBs). We demonstrate that the occurrence rate of GRBs does not depend sensitively on evolutionary parameters in the population synthesis models. We show that the first binaries may form a very small group (< 1%) of fast rotating stars through binary tidal interactions that make GRBs. This finding is contrary to the Bromm & Loeb assumption that all stars in close Population III binaries will be spun up by tides and produce a GRB. We find that there is simply not enough fast rotating stars in Population III binaries to expect detection with SWIFT. Predicted detection rates, even with very optimistic assumptions on binary fraction, evolutionary parameters and GRB detection, are very small: 0.1-0.01 per year.

preprint2006arXiv

The Supernova Gamma-Ray Burst Connection

The chief distinction between ordinary supernovae and long-soft gamma-ray bursts (GRBs) is the degree of differential rotation in the inner several solar masses when a massive star dies, and GRBs are rare mainly because of the difficulty achieving the necessary high rotation rate. Models that do provide the necessary angular momentum are discussed, with emphasis on a new single star model whose rapid rotation leads to complete mixing on the main sequence and avoids red giant formation. This channel of progenitor evolution also gives a broader range of masses than previous models, and allows the copious production of bursts outside of binaries and at high redshifts. However, even the production of a bare helium core rotating nearly at break up is not, by itself, a sufficient condition to make a gamma-ray burst. Wolf-Rayet mass loss must be low, and will be low in regions of low metallicity. This suggests that bursts at high redshift (low metallicity) will, on the average, be more energetic, have more time structure, and last longer than bursts nearby. Every burst consists of three components: a polar jet (~0.1 radian), high energy, subrelativistic mass ejection (~1 radian), and low velocity equatorial mass that can fall back after the initial explosion. The relative proportions of these three components can give a diverse assortment of supernovae and high energy transients whose properties may vary with redshift.

preprint2005arXiv

Millihertz Quasi-Periodic Oscillations from Marginally Stable Nuclear Burning on an Accreting Neutron Star

We investigate marginally stable nuclear burning on the surface of accreting neutron stars as an explanation for the mHz quasi-periodic oscillations (QPOs) observed from three low mass X-ray binaries. At the boundary between unstable and stable burning, the temperature dependence of the nuclear heating rate and cooling rate almost cancel. The result is an oscillatory mode of burning, with an oscillation period close to the geometric mean of the thermal and accretion timescales for the burning layer. We describe a simple one-zone model which illustrates this basic physics, and then present detailed multizone hydrodynamical calculations of nuclear burning close to the stability boundary using the KEPLER code. Our models naturally explain the characteristic 2 minute period of the mHz QPOs, and why they are seen only in a very narrow range of X-ray luminosities. The oscillation period is sensitive to the accreted hydrogen fraction and the surface gravity, suggesting a new way to probe these parameters. A major puzzle is that the accretion rate at which the oscillations appear in the theoretical models is an order of magnitude larger than the rate implied by the X-ray luminosity when the mHz QPOs are seen. We discuss the implications for our general understanding of nuclear burning on accreting neutron stars. One possibility is that the accreted material covers only part of the neutron star surface at luminosities Lx > ~1E37 erg/s.

preprint1997arXiv

Pulsations in red supergiants with high L/M ratio -- Implications for the stellar and circumstellar structure of supernova progenitors

We investigate the pulsational properties of RSG models --- which we evolve from ZAMS masses in the range 10 to $20 \Msun$ --- by means of linear and non-linear calculations. We find period and growth rate of the dominant fundamental mode to increase with increasing luminosity-to-mass ratio $L/M$. Our models obtain relatively large $L/M$ values due to the inclusion of rotation in the evolutionary calculations; however, the largest values are obtained at and beyond central He-exhaustion due to major internal rearrangements of the nuclear burning regions. Our non-linear calculations as well as the behavior of the linear period and growth rate of the pulsations for periods approaching the Kelvin-Helmholtz time scale of the H-rich stellar envelope point towards the possibility of large amplitude pulsations. Such properties are similar to that found in AGB stars and suggest the possibility of a ``superwind'' to occur before the RSGs explode as supernovae. We conclude that changes in global stellar properties during the last few $10^4 \yr$ before core collapse may lead to drastic changes in the pulsational and wind properties of pre-supernova stars, with marked consequences for the immediate pre-supernova structure of the star and the circumstellar medium. We compare our results with observations of long-period OH/IR variables and discuss observational evidence for our scenario from observed supernova light curves, spectra and remnants.

Alexander Heger

What is connected

Connect this record

See the researcher in context

Building this map preview

67 published item(s)

Spectacular nucleosynthesis from early massive stars

Long-term evolution of post-explosion Helium-star Companions of Type Iax Supernovae

The Impact of the New $^{65\!}$As(p,$γ$)$^{66\!}$Se Reaction Rate on the Two-Proton Sequential Capture of $^{64}\!$Ge, Weak GeAs Cycles, and Type-I X-Ray Bursts such as the Clocked Burster GS 1826$-$24

The Regulated NiCu Cycles with the new $^{57}$Cu(p,$γ$)$^{58}$Zn reaction rate and the Influence on Type-I X-Ray Bursts: GS 1826$-$24 Clocked Burster

The Role of the Hadron-Quark Phase Transition in Core-Collapse Supernovae

On the origin of nitrogen at low metallicity

A Minimum Dilution Scenario for Supernovae and Consequences for Extremely Metal-Poor Stars

Helium and Nitrogen Enrichment in Massive Main Sequence Stars: Mechanisms and Implications for the Origin of WNL Stars

Multi-epoch X-ray burst modelling: MCMC with large grids of 1D simulations

On Monolithic Supermassive Stars

Production of Lithium in Primordial Supernovae

The Chemical Evolution of Iron-Peak Elements with Hypernovae

The impact of fallback on the compact remnants and chemical yields of core-collapse supernovae

Large-Scale Mixing in a Violent Oxygen-Neon Shell Merger Prior to a Core-Collapse Supernova

Evidence from stable isotopes and Be-10 for solar system formation triggered by a low-mass supernova

Fast evolving pair-instability supernova models: evolution, explosion, light curves

Mass and Metallicity Requirement in Stellar Models for Galactic Chemical Evolution Applications

Observational Constraints on First-Star Nucleosynthesis. II. Spectroscopy of an Ultra Metal-Poor CEMP-no Star

The Last Minutes of Oxygen Shell Burning in a Massive Star

Are Models for Core-Collapse Supernova Progenitors Consistent with the Properties of Supernova Remnants?

Can direct collapse black holes launch gamma-ray bursts and grow to supermassive black holes?

Cosmological Impact of Population III Binaries

Metal-Poor Stars Observed with the Magellan Telescope. III. New Extremely and Ultra Metal-Poor Stars from SDSS/SEGUE and Insights on the Formation of Ultra Metal-Poor Stars

Neutrino-Induced Nucleosynthesis in Helium Shells of Early Core-Collapse Supernovae

Nucleosynthesis in a Primordial Supernova: Carbon and Oxygen Abundances in SMSS J031300.36-670839.31

Pair-Instability Supernovae of Non-Zero Metallicity Stars

The influence of accretion rate and metallicity on thermonuclear bursts: predictions from KEPLER models

The interplay between chemistry and nucleation in the formation of carbonaceous dust in supernova ejecta

The Remarkable Deaths of 9 - 11 Solar Mass Stars

Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning

Effective Helium Burning Rates and the Production of the Neutrino Nuclei

General Relativistic Instability Supernova of a Supermassive Population III Star

Numerical approaches for multidimensional simulations of stellar explosions

Pair Instability Supernovae of Very Massive Population III Stars

Pair-Instability Supernovae in the Local Universe

Stellar origin of the 182Hf cosmochronometer and the presolar history of solar system matter

The Deaths of Very Massive Stars

Two-Dimensional Simulations of Pulsational Pair-Instability Supernovae

Finding the First Cosmic Explosions I: Pair-Instability Supernovae

Finding the First Cosmic Explosions II: Core-Collapse Supernovae

Finding the First Cosmic Explosions. III. Pulsational Pair-Instability Supernovae

Metallicity-Dependent Galactic Isotopic Decomposition for Nucleosynthesis

Supermassive Population III Supernovae and the Birth of the First Quasars

The Biggest Explosions in the Universe. II

The Impact of Helium-Burning Reaction Rates on Massive Star Evolution and Nucleosynthesis

The Metallicity Dependence of the Minimum Mass for Core-Collapse Supernovae

The Supernova that Destroyed a Protogalaxy: Prompt Chemical Enrichment and Supermassive Black Hole Growth

Very Massive Stars (VMS) in the Local Universe

Seeing the First Supernovae at the Edge of the Universe with JWST

Conservative Initial Mapping For Multidimensional Simulations of Stellar Explosions

Long Gamma-Ray Transients from Collapsars

Multidimensional Simulations of Thermonuclear Supernovae from the First Stars

Pair Instability Supernovae: Light Curves, Spectra, and Shock Breakout

Multi-Dimensional Simulations of Pair-Instability Supernovae

Production of 26Al, 44Ti and 60Fe in Supernovae-sensitivity to the helium burning rates

Production of 26Al, 44Ti, and 60Fe in Core-Collapse Supernovae: Sensitivity to the Rates of the Triple Alpha and 12C(a,g)16O Reactions

The Nucleosynthetic Imprint of 15-40 Solar Mass Primordial Supernovae on Metal-Poor Stars

Two Dimensional Simulations of Pair-Instability Supernovae

Detailed Abundances for 28 Metal-poor Stars: Stellar Relics in the Milky Way

First Stars III Conference Summary

Mixing in Zero and Solar Metallicity Supernovae

Nucleosynthesis and Evolution of Massive Metal-Free Stars

The Impact of Neutrino Magnetic Moments on the Evolution of Massive Stars

The Lack of Gamma-Ray Bursts from Population III Binaries

The Supernova Gamma-Ray Burst Connection

Millihertz Quasi-Periodic Oscillations from Marginally Stable Nuclear Burning on an Accreting Neutron Star

Pulsations in red supergiants with high L/M ratio -- Implications for the stellar and circumstellar structure of supernova progenitors