Detailed Nucleosynthesis Research Articles

Chemical evolution of the Galactic bulge with different stellar populations

Context. The metallicity distribution function (MDF) of the Galactic bulge is characterized by a multi-peak shape, with a metal-poor peak centered at [Fe/H] ∼ −0.3 dex and a metal-rich peak centered at [Fe/H] ∼ +0.3 dex. The bimodality of the MDF is also reflected in the [α/Fe] versus [Fe/H] abundance ratios, suggesting the presence of different stellar populations in the bulge. Aims. In this work we aim to reproduce the observed MDF of the Galactic bulge by testing a scenario in which the metal-poor component of the bulge is formed by stars formed in situ, during a strong burst of star formation, while the metal-rich population is formed by stars created in situ during a second burst of star formation and/or stars accreted from the innermost part of the Galactic disk as an effect of a growing bar. Methods. We adopted a chemical evolution model that is able to follow the evolution of several chemical species with detailed nucleosynthesis prescriptions. In particular, because of the importance of the production of Fe in constraining the MDF, close attention is paid to the production of this element in both Type Ia supernovae and massive stars. In particular, we included yields from rotating massive stars with different rotational velocity prescriptions. Our model also takes the infall and outflow of gas into account, as well as the effect of stellar migration. Results are compared to ∼13 000 stars from the SDSS/APOGEE survey that belong to the region located at a Galactocentric distance RGC ≤ 3.5 kpc. Results. We successfully reproduce the observed double-peak shape of the bulge MDF as well as the abundance trends of the α elements relative to Fe by assuming both (i) a multi-burst star formation history with a quenching of the first burst of ∼102 Myr and (ii) migration of stars from the innermost part of the Milky Way disk, as an effect of a growing bar. According to our results, the fraction of the stellar mass of the bulge-bar that belongs to the inner disk is ∼40%. In terms of the nucleosynthesis, we conclude that models that assume either no rotation for massive stars or a distribution of rotational velocities that favors slow rotation at high metallicities best reproduce the observed MDF as well as the [α/Fe] and the [Ce/Fe] versus [Fe/H] abundance patterns.

On the core-collapse supernova explanation for LAMOST J1010 + 2358

ABSTRACT Low-metallicity very massive stars with an initial mass of ∼140–$260\, \mathrm{M}_\odot$ are expected to end their lives as pair-instability supernovae (PISNe). The abundance pattern resulting from a PISN differs drastically from regular core-collapse supernova (CCSN) models and is expected to be seen in very metal-poor (VMP) stars of [Fe/H] ≲ −2. Despite the routine discovery of many VMP stars, the unique abundance pattern expected from PISNe has not been unambiguously detected. The recently discovered VMP star LAMOST J1010 + 2358, however, shows a peculiar abundance pattern that is remarkably well fit by a PISN, indicating the potential first discovery of a bonafide star born from gas polluted by a PISN. In this paper, we study the detailed nucleosynthesis in a large set of models of CCSN of Pop III and Pop II star of metallicity [Fe/H] = −3 with masses ranging from 12 to $30\, \mathrm{M}_\odot$. We find that the observed abundance pattern in LAMOST J1010 + 2358 can be fit at least equally well by CCSN models of ∼12–$14\, \mathrm{M}_\odot$ that undergo negligible fallback following the explosion. The best-fitting CCSN models provide a fit that is even marginally better than the best-fitting PISN model. We conclude the measured abundance pattern in LAMOST J1010 + 2358 could have originated from a CCSN and therefore cannot be unambiguously identified with a PISN given the set of elements measured in it to date. We identify key elements that need to be measured in future detections in stars like LAMOST J1010 + 2358 that can differentiate between CCSN and PISN origin.

Chemical evolution of elliptical galaxies I: supernovae and AGN feedback

ABSTRACT We study the formation and evolution of elliptical galaxies and how they suppress star formation and maintain it quenched. A one-zone chemical model which follows in detail the time evolution of gas mass and its chemical abundances during the active and passive evolution is adopted. The model includes both gas infall and outflow as well as detailed stellar nucleosynthesis. Elliptical galaxies with different infall masses, following a down-sizing in star formation scenario, are considered. In the chemical evolution simulation, we include a novel calculation of the feedback processes. We include heating by stellar wind, core-collapse supernovae (SNe), Type Ia SNe (usually not highlighted in galaxy formation simulations), and active galactic nucleus (AGN) feedback. The AGN feedback is a novelty in this kind of models and is computed by considering a Bondi-Eddington limited accretion onto the central supermassive black hole. We successfully reproduce several observational features, such as the [α/Fe] ratios increasing with galaxy mass, mass-metallicity, MBH–σ and MBH–M* relations. Moreover, we show that stellar feedback and in particular Type Ia SNe, has a main role in maintaining quenched the star formation after the occurrence of the main galactic wind, especially in low-mass ellipticals. For larger systems, the contribution from AGN to thermal energy of gas appears to be necessary. However, the effect of the AGN on the development of the main galactic wind is negligible, unless an unreasonable high-AGN efficiency or an extremely low-stellar feedback are assumed. We emphasize the important role played by Type Ia SNe in the energy budget of early-type galaxies.

Chemical evolution of 26Al and 60Fe in the Milky Way

ABSTRACT We present theoretical mass estimates of 26Al and 60Fe throughout the Galaxy, performed with a numerical chemical evolution model including detailed nucleosynthesis prescriptions for stable and radioactive nuclides. We compared the results for several sets of stellar yields taken from the literature, for massive, low and intermediate mass stars, nova systems (only for 26Al) and supernovae Type Ia. We then computed the total masses of 26Al and 60Fe in the Galaxy. We studied the bulge and the disc of the Galaxy in a Galactocentric radius range of 0–22 kpc. We assumed that the bulge region (within 2 kpc) evolved quickly suffering a strong star formation burst, while the disc formed more slowly and inside–out. We compared our results with the 26Al mass observed by the gamma-ray surveys COMPTEL and INTEGRAL to select the best model. Concerning 60Fe, we do not have any observed mass value so we just performed a theoretical prediction for future observations. In conclusion, low, intermediate mass stars, and Type Ia supernovae contribute negligibly to the two isotopes, while massive stars are the dominant source. The nova contribution is, however, necessary to reproduce the observations of 26Al. Our best model predicts 2.12 M⊙ of 26Al, in agreement with observations, while for 60Fe our best mass estimate is ∼1.05 M⊙. We also predicted the present injection rate of 26Al and 60Fe in the Galaxy and compared it with previous results, and we found a larger present time injection rate along the disc.

PUSHing Core-collapse Supernovae to Explosions in Spherical Symmetry. V. Equation of State Dependency of Explosion Properties, Nucleosynthesis Yields, and Compact Remnants

Abstract In this fifth paper of the series, we use the parameterized, spherically symmetric explosion method PUSH to investigate the impact of eight different nuclear equations of state (EOS). We present and discuss the explosion properties and the detailed nucleosynthesis yields, and predict the remnant (neutron star or black hole) for all our simulations. For this, we perform two sets of simulations. First, a complete study of nonrotating stars from 11 to 40 M ⊙ at three different metallicities using the SFHo EOS; and, second, a suite of simulations for four progenitors (16 M ⊙ at three metallicities and 25 M ⊙ at solar metallicity) for eight different nuclear EOS. We compare our predicted explosion energies and yields to observed supernovae and to the metal-poor star HD 84937. We find EOS-dependent differences in the explosion properties and the nucleosynthesis yields. However, when comparing to observations, these differences are not large enough to rule out any EOS considered in this work.

Type Iax supernovae from deflagrations in Chandrasekhar mass white dwarfs

Context. Due to the ever increasing number of observations during the past decades, Type Ia supernovae are nowadays regarded as a heterogeneous class of optical transients consisting of several subtypes. One of the largest of these subclasses is the class of Type Iax supernovae. They have been suggested to originate from pure deflagrations in carbon-oxygen Chandrasekhar mass white dwarfs because the outcome of this explosion scenario is in general agreement with their subluminous nature. Aims. Although a few deflagration studies have already been carried out, the full diversity of the class has not been captured yet. This, in particular, holds for the faint end of the subclass. We therefore present a parameter study of single-spot ignited deflagrations in Chandrasekhar mass white dwarfs varying the location of the ignition spark, the central density, the metallicity, and the composition of the white dwarf. We also explore a rigidly rotating progenitor to investigate whether the effect of rotation can spawn additional trends. Methods. We carried out three-dimensional hydrodynamic simulations employing the LEAFS code. Subsequently, detailed nucleosynthesis results were obtained with the nuclear network code YANN. In order to compare our results to observations, we calculated synthetic spectra and light curves with the ARTIS code. Results. The new set of models extends the range in brightness covered by previous studies to the lower end. Our single-spot ignited explosions produce 56Ni masses from 5.8 × 10−3 to 9.2 × 10−2 M⊙. In spite of the wide exploration of the parameter space, the main characteristics of the models are primarily driven by the mass of 56Ni and form a one-dimensional sequence. Secondary parameters seem to have too little impact to explain the observed trend in the faint part of the Type Iax supernova class. We report kick velocities of the gravitationally bound explosion remnants from 6.9 to 369.8 km s−1. The magnitude as well as the direction of the natal kick is found to depend on the strength of the deflagration. Conclusions. This work corroborates the results of previous studies of deflagrations in Chandrasekhar mass white dwarfs. The wide exploration of the parameter space in initial conditions and viewing angle effects in the radiative transfer lead to a significant spread in the synthetic observables. The trends in observational properties toward the faint end of the class are, however, not reproduced. This motivates a quantification of the systematic uncertainties in the modeling procedure and the influence of the 56Ni-rich bound remnant to get to the bottom of these discrepancies. Moreover, while the pure deflagration scenario remains a favorable explanation for bright and intermediate luminosity Type Iax supernovae, our results suggest that other mechanisms also contribute to this class of events.

Production of Very Light Elements and Strontium in the Early Ejecta of Neutron Star Mergers

Abstract We study the production of very light elements (Z < 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10−6 M ⊙ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10−6 M ⊙ of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r-process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817.

Nebular phase properties of supernova Ibc from He-star explosions

Following our recent work on Type II supernovae (SNe), we present a set of 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc SNe starting from state-of-the-art explosion models with detailed nucleosynthesis. Our grid of progenitor models is derived from He stars that were subsequently evolved under the influence of wind mass loss. These He stars, which most likely form through binary mass exchange, synthesize less oxygen than their single-star counterparts with the same zero-age main sequence (ZAMS) mass. This reduction is greater in He-star models evolved with an enhanced mass loss rate. We obtain a wide range of spectral properties at 200 d. In models from He stars with an initial mass > 6 M⊙, the [O I] λλ 6300, 6364 is of a comparable or greater strength than [Ca II] λλ 7291, 7323 – the strength of [O I] λλ 6300, 6364 increases with the He-star initial mass. In contrast, models from lower mass He stars exhibit a weak [O I] λλ 6300, 6364, strong [Ca II] λλ 7291, 7323, and also strong N II lines and Fe II emission below 5500 Å. The ejecta density, which is modulated by the ejecta mass, the explosion energy, and clumping, has a critical impact on gas ionization, line cooling, and spectral properties. We note that Fe II dominates the emission below 5500 Å and is stronger at earlier nebular epochs. It ebbs as the SN ages, while the fractional flux in [O I] λλ 6300, 6364 and [Ca II] λλ 7291, 7323 increases with a similar rate as the ejecta recombine. Although the results depend on the adopted wind mass loss rate and pre-SN mass, we find that He-stars of 6–8 M⊙ initially (ZAMS mass of 23–28 M⊙) match the properties of standard SNe Ibc adequately. This finding agrees with the offset in progenitor masses inferred from the environments of SNe Ibc relative to SNe II. Our results for less massive He stars are more perplexing since the predicted spectra are not seen in nature. They may be missed by current surveys or associated with Type Ibn SNe in which interaction power dominates over decay power.

Gamma-Ray Emission Produced by r-process Elements from Neutron Star Mergers

Abstract The observation of a radioactively powered kilonova AT 2017gfo associated with the gravitational wave event GW170817 from a binary neutron star merger proves that these events are ideal sites for the production of heavy r-process elements. The gamma-ray photons produced by the radioactive decay of heavy elements are unique probes for the detailed nuclide compositions. Based on the detailed r-process nucleosynthesis calculations and considering radiative transport calculations for the gamma rays in different shells, we study the gamma-ray emission in a merger ejecta on a timescale of a few days. It is found that the total gamma-ray energy generation rate evolution is roughly depicted as E ̇ ∝ t − 1.3 . For the dynamical ejecta with a low electron fraction (Y e ≲ 0.20), the dominant contributors of gamma-ray energy are the nuclides around the second r-process peak (A ∼ 130) and the decay chain of 132Te (t 1/2 = 3.21 days) → 132I (t 1/2 = 0.10 days) → 132Xe produces gamma-ray lines at 228, 668, and 773 keV. For the case of a wind ejecta with Y e ≳ 0.30, the dominant contributors of gamma-ray energy are the nuclides around the first r-process peak (A ∼ 80) and the decay chain of 72Zn (t 1/2 = 1.93 days) → 72Ga (t 1/2 = 0.59 days) → 72Ge produces gamma-ray lines at 145, 834, 2202, and 2508 keV. The peak fluxes of these lines are 10−9 ∼ 10−7 ph cm−2 s−1, which are marginally detectable with the next-generation MeV gamma-ray detector ETCC if the source is at a distance of 40 Mpc.

The explosion of 9–29M⊙stars as Type II supernovae: Results from radiative-transfer modeling at one year after explosion

We present a set of nonlocal thermodynamic equilibrium steady-state calculations of radiative transfer for one-year-old Type II supernovae (SNe) starting from state-of-the-art explosion models computed with detailed nucleosynthesis. This grid covers single-star progenitors with initial masses between 9 and 29M⊙, all evolved with the codeKEPLERat solar metallicity and ignoring rotation. The [O I]λλ6300, 6364 line flux generally grows with progenitor mass, and Hαexhibits an equally strong and opposite trend. The [Ca II]λλ7291, 7323 strength increases at low56Ni mass, at low explosion energy, or with clumping. This Ca IIdoublet, which forms primarily in the explosively produced Si/S zones, depends little on the progenitor mass but may strengthen if Ca+dominates in the H-rich emitting zones or if Ca is abundant in the O-rich zones. Indeed, Si–O shell merging prior to core collapse may boost the Ca IIdoublet at the expense of the O Idoublet, and may thus mimic the metal line strengths of a lower-mass progenitor. We find that the56Ni bubble effect has a weak impact, probably because it is too weak to induce much of an ionization shift in the various emitting zones. Our simulations compare favorably to observed SNe II, including SN 2008bk (e.g., the 9M⊙model), SN 2012aw (12M⊙model), SN 1987A (15M⊙model), or SN 2015bs (25M⊙model with no Si–O shell merging). SNe II with narrow lines and a low56Ni mass are well matched by the weak explosion of 9–11M⊙progenitors. The nebular-phase spectra of standard SNe II can be explained with progenitors in the mass range 12–15M⊙, with one notable exception for SN 2015bs. In the intermediate mass range, these mass estimates may increase by a fewM⊙, with allowance for clumping of the O-rich material or CO molecular cooling.

Chemical evolution of ultrafaint dwarf galaxies: testing the IGIMF

ABSTRACT We test the integrated galactic initial mass function (IGIMF) on the chemical evolution of 16 ultrafaint dwarf (UFD) galaxies discussing in detail the results obtained for three of them: Boötes I, Boötes II, and Canes Venatici I, taken as prototypes of the smallest and the largest UFDs. These objects have very small stellar masses (∼103–104 M⊙) and quite low metallicities ([Fe/H] < −1.0 dex). We consider four observational constraints: the present-day stellar mass, the [α/Fe] versus [Fe/H] relation, the stellar metallicity distribution function and the cumulative star formation history. Our model follows in detail the evolution of several chemical species (H, He, α-elements, and Fe). We take into account detailed nucleosynthesis and gas flows (in and out). Our results show that the IGIMF, coupled with the very low star formation rate predicted by the model for these galaxies (∼10−4–10−6 M⊙yr−1), cannot reproduce the main chemical properties, because it implies a negligible number of core-collapse SNe and even Type Ia SNe, the most important polluters of galaxies. On the other hand, a constant classical Salpeter IMF gives the best agreement with data, but we cannot exclude that other formulations of the IGIMF could reproduce the properties of these galaxies. Comparing with Galaxy data, we suggest that UFDs could not be the building blocks of the entire Galactic halo, although more data are necessary to draw firmer conclusions.

Rapid Transients Originating from Thermonuclear Explosions in Helium White Dwarf Tidal Disruption Events

Abstract We study the emission properties of thermonuclear explosions in a helium white dwarf (WD) tidal disruption event (TDE). We consider a TDE where a 0.2 M ⊙ helium WD is disrupted by a 102.5 M ⊙ intermediate-mass black hole (IMBH). The helium WD is not only tidally disrupted but is also detonated by the tidal compression and by succeeding shocks. We focus on the emission powered by radioactive nuclei in the unbound TDE ejecta. We perform hydrodynamic simulations coupled with nuclear reactions, post-process detailed nucleosynthesis calculations, and radiative transfer simulations. We thus derive multi-band light curves and spectra. The helium WD TDE shows rapid (Δt 1mag ≃ 5–10 days) and relatively faint ( ) light curves, because the ejecta mass and 56Ni mass are low (0.12 M ⊙ and 0.03 M ⊙, respectively). The spectra show strong calcium and Fe-peak features and very weak silicon features, reflecting the peculiar elemental abundance. The key feature is the Doppler shift of the spectral lines up to ≃ ±12,000 km s−1, depending on the viewing angle, due to the bulk motion of the ejecta. Our model matches well with two rapid and faint transients reported in Pursiainen et al. The particular model presented here does not match with observed SNe Iax, calcium-rich transients, or .Ia explosion candidates, either in the spectra or light curves. However, we expect a large variety of observational signatures once a wide range of the WD/black hole masses and orbital parameters are considered. This study contributes to the search for WD TDEs with current and upcoming surveys, and to the identification of IMBHs as disrupters in TDEs.

Nucleosynthesis for SN 1987A from single-star and binary-merger progenitors

We apply the parametrized, spherically symmetric explosion method PUSH to two sets of pre-explosion models suitable for SN 1987A: blue supergiants (BSGs) resulting from the merger of a main sequence star with a giant and red supergiants (RSGs) representing the end point of single-star stellar evolution. For each model, we perform a calibration of the PUSH method to the observational properties of SN 1987A and calculate the detailed explosive nucleosynthesis yields. We find that such a calibration to SN 1987A is only possible for one of the BSG models. We compare the yields from this model with the yields from the best-fit RSG model. The largest differences are found for nuclei in the mass range of which are mostly synthesized pre-explosion. We predict a neutron star with a gravitational mass of 1.48 M⊙ from the BSG model and a neutron star of 1.41 M⊙ from the RSG model.

Nucleosynthesis in 2D core-collapse supernovae of 11.2 and 17.0 M⊙ progenitors: implications for Mo and Ru production

Core-collapse supernovae (CCSNe) are the first polluters of heavy elements in the galactic history. As such, it is important to study the nuclear compositions of their ejecta, and understand their dependence on the progenitor structure (e.g., mass, compactness, metallicity). Here, we present a detailed nucleosynthesis study based on two long-term, two-dimensional CCSN simulations of a 11.2 M and a 17.0 M star. We find that in both models nuclei well beyond the iron group (up to ) can be produced, and discuss in detail also the nucleosynthesis of the p-nuclei 92,94Mo and 96,98Ru. While we observe the production of 92Mo and 94Mo in slightly neutron-rich conditions in both simulations, 96,98Ru can only be produced efficiently via the νp-process. Furthermore, the production of Ru in the νp-process heavily depends on the presence of very proton-rich material in the ejecta. This disentanglement of production mechanisms has interesting consequences when comparing to the abundance ratios between these isotopes in the solar system and in presolar grains.

Neutron Star Mergers and Nucleosynthesis of Heavy Elements

The existence of neutron star mergers has been supported since the discovery of the binary pulsar and the observation of its orbital energy loss, consistent with General Relativity. They are considered nucleosynthesis sites of the rapid neutron-capture process (r-process), which is responsible for creating approximately half of all heavy elements beyond Fe and is the only source of elements beyond Pb and Bi. Detailed nucleosynthesis calculations based on the decompression of neutron star matter are consistent with solar r-process abundances of heavy nuclei. Neutron star mergers have also been identified with short-duration [Formula: see text]-ray bursts via their IR afterglow. The high neutron densities in ejected matter permit a violent r-process, leading to fission cycling of the heaviest nuclei in regions far from (nuclear) stability. Uncertainties in several nuclear properties affect the abundance distributions. The modeling of astrophysical events also depends on the hydrodynamic treatment, the occurrence of a neutrino wind after the merger and before the possible emergence of a black hole, and the properties of black hole accretion disks. We discuss the effect of nuclear and modeling uncertainties and conclude that binary compact mergers are probably a (or the) dominant site of the production of r-process nuclei in our Galaxy.

Dark matter relic density in scalar-tensor gravity revisited

We revisit the calculation of dark matter relic abundances in scalar-tensor gravity using a generic form A(φ*) = eβφ*2/2 for the coupling between the scalar field φ* and the metric, for which detailed Big Bang Nucleosynthesis constraints are available. We find that BBN constraints restrict the modified expansion rate in these models to be almost degenerate with the standard expansion history at the time of dark matter decoupling. In this case the maximum level of enhancement of the dark matter relic density was found to be a factor of ∼ 3, several orders of magnitude below that found in previous investigations.

Deflagrations in hybrid CONe white dwarfs: a route to explain the faint Type Iax supernova 2008ha

Stellar evolution models predict the existence of hybrid white dwarfs (WDs) with a carbon-oxygen core surrounded by an oxygen-neon mantle. Being born with masses ~1.1 Msun, hybrid WDs in a binary system may easily approach the Chandrasekhar mass (MCh) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-MCh hybrid WD under the assumption that nuclear burning only occurs in carbon-rich material. Performing hydrodynamics simulations of the explosion and detailed nucleosynthesis post-processing calculations, we find that only 0.014 Msun of material is ejected while the remainder of the mass stays bound. The ejecta consist predominantly of iron-group elements, O, C, Si and S. We also calculate synthetic observables for our model and find reasonable agreement with the faint Type Iax SN 2008ha. This shows for the first time that deflagrations in near-MCh WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-MCh bound remnant opens the possibility for recurrent explosions or a subsequent accretion-induced collapse in faint Type Iax SNe, if further accretion episodes occur. From binary population synthesis calculations, we find the rate of hybrid WDs approaching MCh to be on the order of 1 percent of the Galactic SN Ia rate.

CARBON-ENHANCED METAL-POOR STARS: RELICS FROM THE DARK AGES

We use detailed nucleosynthesis calculations and a realistic prescription for the environment of the first stars to explore the first episodes of chemical enrichment that occurred during the dark ages. Based on these calculations, we propose a novel explanation for the increased prevalence of carbon-enhanced metal-poor (CEMP) stars with decreasing Fe abundance: The observed chemistry for the most metal-poor Galactic halo stars is the result of an intimate link between the explosions of the first stars and their host minihalo's ability to retain its gas. Specifically, high-energy supernovae produce a near solar ratio of C/Fe, but are effective in evacuating the gas from their host minihalo, thereby suppressing the formation of a second generation of stars. On the other hand, minihalos that host low-energy supernovae are able to retain their gas and form a second stellar generation but, as a result, the second stars are born with a supersolar ratio of C/Fe. Our models are able to accurately reproduce the observed distributions of [C/Fe] and [Fe/H], as well as the fraction of CEMP stars relative to non-CEMP stars as a function of [Fe/H] without any free parameters. We propose that the present lack of chemical evidence for very massive stars (>140 Msun), that ended their lives as a highly energetic pair-instability supernova, does not imply that such stars were rare or did not exist; the chemical products of these very massive first stars may have escaped from their host minihalo, and were never incorporated into subsequent generations of stars. Finally, our models suggest that the most Fe-poor stars currently known may have seen the enrichment from a small multiple of metal-free stars, and need not have been exclusively enriched by a solitary first star. These calculations support the idea that some of the surviving dwarf satellite galaxies of the Milky Way are relics of the first galaxies.

PREPARING FOR AN EXPLOSION: HYDRODYNAMIC INSTABILITIES AND TURBULENCE IN PRESUPERNOVAE

Both observations and direct numerical simulations are discordant with predictions of conventional stellar evolution codes for the latest stages of a massive star's life prior to core collapse. We suggest that the problem lies in the treatment of turbulent convection in these codes, which ignores finite amplitude fluctuations in velocity and temperature, and their nonlinear interaction with nuclear burning. The hydrodynamic instabilities that may arise prompt us to discuss a number of far-reaching implications for the fates of massive stars. In particular, we explore connections to enhanced presupernova mass loss, unsteady nuclear burning and consequent eruptions, swelling of the stellar radius that may trigger violent interactions with a companion star, and potential modifications to the core structure that could dramatically impact calculations of the core-collapse mechanism itself. These modifications may be of fundamental importance to the interpretation of measured isotopic anomalies in meteorites, changing the predictions of both mixing and detailed nucleosynthesis, and of young core-collapse supernova remnants. They may also make possible the development of an early warning system for the detection of impending core collapse.

Three-dimensional pure deflagration models with nucleosynthesis and synthetic observables for Type Ia supernovae

We investigate whether pure deflagration models of Chandrasekhar-mass carbon-oxygen white dwarf stars can account for one or more subclass of the observed population of Type Ia supernova (SN Ia) explosions. We compute a set of 3D full-star hydrodynamic explosion models, in which the deflagration strength is parametrized using the multispot ignition approach. For each model, we calculate detailed nucleosynthesis yields in a post-processing step with a 384 nuclide nuclear network. We also compute synthetic observables with our 3D Monte Carlo radiative transfer code for comparison with observations. For weak and intermediate deflagration strengths (energy release E_nuc <~ 1.1 x 10^51 erg), we find that the explosion leaves behind a bound remnant enriched with 3 to 10 per cent (by mass) of deflagration ashes. However, we do not obtain the large kick velocities recently reported in the literature. We find that weak deflagrations with E_nuc ~ 0.5 x 10^51 erg fit well both the light curves and spectra of 2002cx-like SNe Ia, and models with even lower explosion energies could explain some of the fainter members of this subclass. By comparing our synthetic observables with the properties of SNe Ia, we can exclude the brightest, most vigorously ignited models as candidates for any observed class of SN Ia: their B - V colours deviate significantly from both normal and 2002cx-like SNe Ia and they are too bright to be candidates for other subclasses.