Arxiv:1710.04254V5 [Astro-Ph.SR] 25 Jun 2019 Hr Xssawd Iest Nsei Se
Total Page:16
File Type:pdf, Size:1020Kb
Accepted for publication in Astrophysical Journal Supplement, submitted on 19 September 2017, revised on 16 April 2018 Typeset using LATEX twocolumn style in AASTeX62 Explosive Nucleosynthesis in Near-Chandrasekhar Mass White Dwarf Models for Type Ia supernovae: Dependence on Model Parameters Shing-Chi Leung, Ken’ichi Nomoto1 1Kavli Institute for the Physics and Mathematics of the Universe (WPI), The University of Tokyo Institutes for Advanced Study The University of Tokyo, Kashiwa, Chiba 277-8583, Japan ABSTRACT We present two-dimensional hydrodynamics simulations of near-Chandrasekhar mass white dwarf (WD) models for Type Ia supernovae (SNe Ia) using the turbulent deflagration model with deflagration- detonation transition (DDT). We perform a parameter survey for 41 models to study the effects of the initial central density (i.e., WD mass), metallicity, flame shape, DDT criteria, and turbulent flame formula for a much wider parameter space than earlier studies. The final isotopic abundances of 11C to 91Tc in these simulations are obtained by post-process nucleosynthesis calculations. The survey includes SNe Ia models with the central density from 5 108 g cm−3 to 5 109 g cm−3 (WD masses × × of 1.30 - 1.38 M⊙), metallicity from 0 to 5 Z⊙, C/O mass ratio from 0.3 - 1.0 and ignition kernels including centered and off-centered ignition kernels. We present the yield tables of stable isotopes from 12C to 70Zn as well as the major radioactive isotopes for 33 models. Observational abundances of 55Mn, 56Fe, 57Fe and 58Ni obtained from the solar composition, well-observed SNe Ia and SN Ia remnants are used to constrain the explosion models and the supernova progenitor. The connection between the pure turbulent deflagration model and the subluminous SNe Iax is discussed. We find that dependencies of the nucleosynthesis yields on the metallicity and the central density (WD mass) are large. To fit these observational abundances and also for the application of galactic chemical evolution modeling, these dependencies on the metallicity and WD mass should be taken into account. Keywords: (stars:)nuclear reactions, nucleosynthesis, abundances, hydrodynamics 1. INTRODUCTION ously known variations ranging from luminous SNe Ia Type Ia supernovae (SNe Ia) have been shown (super-Chandra and SN 1991T-like) to SN 1991bg-like to be the major source of many iron-peak elements faint SNe Ia, a new sub-class of SN 2002cx-like, or in the galaxies (e.g., Arnett 1996; Matteucci 2001, Type Iax SNe (see, e.g., a recent review by Jha 2017) 2012; Nomoto et al. 1984, 1997; Nomoto & Leung 2017, and SN 2002es-like SNe (e.g., Taubenberger 2017) have 2018). To understand how SNe Ia contribute to the been reported. Such brightness variations imply a large 56 variation of the Ni mass ( 0.1 M⊙ to 1.4 M⊙) metal enrichment process in the galaxies, and to explain ∼ ∼ the growing diversities of the observational results, sim- produced in the explosions. ulations of SN Ia models with much wider parameter Studies of nucleosynthesis of one-dimensional mod- els have shown some important dependencies on the arXiv:1710.04254v5 [astro-ph.SR] 25 Jun 2019 ranges need to be done. SNe Ia have been well-modeled by the thermonu- model parameters. For example, Fe-peak elements clear explosions of carbon-oxygen white dwarfs (WDs) synthesized in the Chandrasekhar mass model W7 (e.g., Hillebrandt & Niemeyer 2000). Both the Chan- (Nomoto et al. 1984) is shown to be consistent with drasekhar mass WD model and sub-Chandra mass WD the solar abundances, except for a significant over- 58 model can be consistent with the similarity of SN Ia production of Ni, where the rate of electron capture is light curves (e.g., Branch & Wheeler 2017). The empir- important (e.g., Thielemann et al. 1986; Iwamoto et al. ical similarity later leads to the discovery of accelerating 1999; Brachwitz et al. 2000; Mori et al. 2016). The so- cosmological expansion and the existence of dark energy lar abundance pattern of Fe-peak elements can also be (Riess et al. 1998; Perlmutter et al. 1999). reproduced by the sub-Chandra models with solar met- However, recent observations have suggested that allcity (e.g., Shigeyama et al. 1992; Nomoto et al. 1994), there exists a wide diversity in SNe Ia (see, e.g., a although the Ni/Fe ratio (which depends on metallicity) review by Taubenberger 2017). In addition to previ- tends to be under solar because of the low central den- 2 sity of exploding WDs. As shown in [Ni/Fe] mentioned Such results indicate that properties of these SNe Ia above, the central density of the WD is an important might have important metallicity effects. Nucleosynthe- parameter because of the level of electron capture. If sis in SNe Ia has been studied extensively but still only the central density is high enough, synthesis of certain a small parameter space has been explored. neutron-rich isotopes, such as 48Ca, 50Ti, and 54Cr can In view of this background, we make systematic mod- be significant (Woosley 1997). eling of SNe Ia for various explosion configuration and Another interesting example is [Mn/Fe], which in- setting to see how the model parameters of SNe Ia af- creases with increasing [Fe/H] in the Galactic Halo (e.g., fect the WD explosions and their chemical yields for Hinkel et al. 2014; Mishenina et al. 2015). [Mn/Fe] Chandrasekhar mass WD for much more wide range of in metal-poor stars in dwarf-spheroidal galaxies (e.g., parameters (i.e., WD masses (central densities), metal- Larsen et al. 2014; Sbordone et al. 2015) provides an- licity, flame structure). In the forthcoming paper, we other information on metallicity dependence. In order will present our sub-Chandrasekhar mass models. The to calculate the chemical evolution of dwarf-spheroidal chemical yields of SNe Ia, which depend on the model galaxies, metallicity-dependent SN Ia yields are neces- parameters and environmental conditions, can be con- sary (e.g., Kobayashi et al. 2015). strained by the observed abundance patterns of Fe-peak Recent observations of SNe Ia remnants in the nebu- elements in various stars and systems. lar phase have provided important insights to the models We use our own 2D hydrodynamics code (see Ap- of SNe Ia. They include Tycho (Yamaguchi et al. 2015), pendix A), because our 2D hydro code is suitable to Kepler (Park et al. 2013), and 3C 397 (Yamaguchi et al. calculate many models for a wide range of parameters 2015). The relative X-ray flux of iron-peak elements than 3D hydro code. The typical running time for one can give promising constraints on the explosion condi- 2D model is days on a single machine; while it takes ∼ tions. For example, these three remnants have been sug- weeks to months for a cluster to calculate the explosion gested to have progenitors with super-solar metallicity phase of a model in 3D. Certainly, the 2D simulations (Yamaguchi et al. 2015). have some qualitative differences from the 3D simula- These would suggest the importance to obtain the SN tions in two ways. First, the flame in 2D models tends Ia yields for a wide range of environmental conditions, to propagate faster than in 3D models because of the such as metallicity and the mass (and thus the central larger surface area for the same 2D-projection. Second, density) of the WDs. the imposed symmetry may enhance the growth of hy- Such study will be important for the future use of drodynamical instabilities owing to the imposed reflec- galactic chemical evolution for an accurate modeling of tive boundaries, which stimulate the growth of boundary isotopes as a function of metallicity. Nucleosynthesis of flows. multi-dimensional hydrodynamics simulations was made In Section 3 we construct the benchmark to be a typ- in Travaglio et al. (2004) with the use of tracer parti- ical SNe Ia model. In Section 4, we present the nucle- cle scheme (Seitenzahl et al. 2010). The effects of ini- osynthesis yields of our models and show how large the tial flame structure to the chemical yield was studied in effects of each model parameter are. We then discuss Maeda et al. (2010); Fink et al. (2014); Seitenzahl et al. how our results can be applied to observational data to (2013) for different explosion models. constrain the model parameter. In Section 5 we com- A few more recent works have studied these objects. pare our results with earlier calculations. In Appendix, In Shen et al. (2017) the sub-Chandrasekhar SN Ia mod- we summarize the numerical code which is specifically els are revisited and showed that the sub-Chandrasekhar developed for modeling SNe Ia (Leung et al. 2015a). We SNe Ia can be connected to the remnant 3C 397 when ap- describe the updates and changes in the hydrodynamics propriate amount of reverse shock-heating is considered. and nucleosynthesis. Finally, in order to apply for the Similar explorations were done in Dave et al. (2017), chemical evolution modeling and comparison with obser- where some representative models of pure turbulent vational data, we present tables of the nucleosynthesis deflagration, deflagration-detonation transition (DDT) yields of our 24 models. and gravitationally confined detonation are explored. It is shown that for a sub-set of central densities, C/O 2. INITIAL MODELS AND METHODS ratio and high offset in the initial flame allow models 2.1. Initial Models to produce super-solar [Mn/Fe] to match the observed We construct the initial C+O WD models at the cen- data. See, e.g. Nomoto & Leung (2017) for a general tral carbon ignition with the masses from M =1.30M⊙ review of nucleosynthesis pattern and its connection to 8 to 1.39 M⊙ (and thus the central densities from 5 10 g explosion mechanism.