Binary stellar evolution calculations predict that Chandrasekhar-mass carbon/oxygen white dwarfs (WDs) show a radially varying profile for the composition with a carbon depleted core. Many recent ...multi-dimensional simulations of Type Ia supernovae (SNe Ia), however, assume the progenitor WD has a homogeneous chemical composition. In this work, we explore the impact of different initial carbon profiles of the progenitor WD on the explosion phase and on synthetic observables in the Chandrasekhar-mass delayed detonation model. Spectra and light curves are compared to observations to judge the validity of the model. The explosion phase is simulated using the finite volume supernova code LEAFS, which is extended to treat different compositions of the progenitor WD. The synthetic observables are computed with the Monte Carlo radiative transfer code ARTIS. Differences in binding energies of carbon and oxygen lead to a lower nuclear energy release for carbon depleted material; thus, the burning fronts that develop are weaker and the total nuclear energy release is smaller. For otherwise identical conditions, carbon depleted models produce less Ni-56. Comparing different models with similar Ni-56 yields shows lower kinetic energies in the ejecta for carbon depleted models, but only small differences in velocity distributions and line velocities in spectra. The light curve width-luminosity relation (WLR) obtained for models with differing carbon depletion is roughly perpendicular to the observed WLR, hence the carbon mass fraction is probably only a secondary parameter in the family of SNe Ia.
Context: Manganese is predominantly synthesised in Type Ia supernova (SN Ia) explosions. Owing to the entropy dependence of the Mn yield in explosive thermonuclear burning, SNe Ia involving near ...Chandrasekhar-mass white dwarfs (WDs) are predicted to produce Mn to Fe ratios significantly exceeding those of SN Ia explosions involving sub-Chandrasekhar mass primary WDs. Of all current supernova explosion models, only SN Ia models involving near-Chandrasekhar mass WDs produce Mn/Fe > 0.0. Aims: Using the specific yields for competing SN Ia scenarios, we aim to constrain the relative fractions of exploding near-Chandrasekhar mass to sub-Chandrasekhar mass primary WDs in the Galaxy. Methods: We extract the Mn yields from three-dimensional thermonuclear supernova simulations referring to different initial setups and progenitor channels. We then compute the chemical evolution of Mn in the Solar neighborhood, assuming SNe Ia are made up of different relative fractions of the considered explosion models. Results: We find that due to the entropy dependence of freeze-out yields from nuclear statistical equilibrium, Mn/Fe strongly depends on the mass of the exploding WD, with near-Chandraskher mass WDs producing substantially higher Mn/Fe than sub-Chandrasekhar mass WDs. Of all nucleosynthetic sources potentially influencing the chemical evolution of Mn, only explosion models involving the thermonuclear incineration of near-Chandrasekhar mass WDs predict solar or super-solar Mn/Fe. Consequently, we find in our chemical evolution calculations that the observed Mn/Fe in the Solar neighborhood at Fe/H > 0.0 cannot be reproduced without near-Chandrasekhar mass SN Ia primaries. Assuming that 50 per cent of all SNe Ia stem from explosive thermonuclear burning in near-Chandrasekhar mass WDs results in a good match to data.
We use a high resolution cosmological zoom simulation of a Milky Way-sized halo to study the observable features in velocity and metallicity space associated with the dynamical influence of spiral ...arms. For the first time, we demonstrate that spiral arms, that form in a disc in a fully cosmological environment with realistic galaxy formation physics, drive large-scale systematic streaming motions. In particular, on the trailing edge of the spiral arms the peculiar galacto-centric radial and azimuthal velocity field is directed radially outward and azimuthally backward, whereas it is radially inward and azimuthally forward on the leading edge. Owing to the negative radial metallicity gradient, this systematic motion drives, at a given radius, an azimuthal variation in the residual metallicity that is characterised by a metal rich trailing edge and a metal poor leading edge. We show that these signatures are theoretically observable in external galaxies with Integral Field Unit instruments such as VLT/MUSE, and if detected, would provide evidence for large-scale systematic radial migration driven by spiral arms.
The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically ...important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly \(100\,\rm{Myr}\) in the center of halos until saturation occurs around \(z=2-3\), when the magnetic energy density reaches about \(10\%\) of the turbulent energy density with a typical strength of \(10-50\,\rm{\mu G}\). In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until \(z=0\). At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around \(z=0.5\) to \(z=0\). The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.
The progenitor problem of Type Ia supernovae (SNe Ia) is still unsolved. Most of these events are thought to be explosions of carbon-oxygen (CO) white dwarfs (WDs), but for many of the explosion ...scenarios, particularly those involving the externally triggered detonation of a sub-Chandrasekhar mass WD (sub-M Ch WD), there is also a possibility of having an oxygen-neon (ONe) WD as progenitor. We simulate detonations of ONe WDs and calculate synthetic observables from these models. The results are compared with detonations in CO WDs of similar mass and observational data of SNe Ia. We perform hydrodynamic explosion simulations of detonations in initially hydrostatic ONe WDs for a range of masses below the Chandrasekhar mass (M Ch), followed by detailed nucleosynthetic postprocessing with a 384-isotope nuclear reaction network. The results are used to calculate synthetic spectra and light curves, which are then compared with observations of SNe Ia. We also perform binary evolution calculations to determine the number of SNe Ia involving ONe WDs relative to the number of other promising progenitor channels. The ejecta structures of our simulated detonations in sub-M Ch ONe WDs are similar to those from CO WDs. There are, however, small systematic deviations in the mass fractions and the ejecta velocities. These lead to spectral features that are systematically less blueshifted. Nevertheless, the synthetic observables of our ONe WD explosions are similar to those obtained from CO models. Our binary evolution calculations show that a significant fraction (3-10%) of potential progenitor systems should contain an ONe WD. The comparison of our ONe models with our CO models of comparable mass (1.2 Msun) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.
The gravitationally confined detonation (GCD) model has been proposed as a possible explosion mechanism for Type Ia supernovae in the single-degenerate evolution channel. Driven by buoyancy, a ...deflagration flame rises in a narrow cone towards the surface. For the most part, the flow of the expanding ashes remains radial, but upon reaching the outer, low-pressure layers of the white dwarf, an additional lateral component develops. This makes the deflagration ashes converge again at the opposite side, where the compression heats fuel and a detonation may be launched. To test the GCD explosion model, we perform a 3D simulation for a model with an ignition spot offset near the upper limit of what is still justifiable, 200 km. This simulation meets our deliberately optimistic detonation criteria and we initiate a detonation. The detonation burns through the white dwarf and leads to its complete disruption. We determine nucleosynthetic yields by post-processing 10^6 tracer particles with a 384 nuclide reaction network and we present multi-band light curves and time-dependent optical spectra. We find that our synthetic observables show a prominent viewing-angle sensitivity in UV and blue bands, which is in tension with observed SNe Ia. The strong dependence on viewing-angle is caused by the asymmetric distribution of the deflagration ashes in the outer ejecta layers. Finally, we perform a comparison of our model to SN 1991T. The overall flux-level of the model is slightly too low and the model predicts pre-maximum light spectral features due to Ca, S, and Si that are too strong. Furthermore, the model chemical abundance stratification qualitatively disagrees with recent abundance tomography results in two key areas: our model lacks low velocity stable Fe and instead has copious amounts of high-velocity 56Ni and stable Fe. We therefore do not find good agreement of the model with SN 1991T.
Extremely luminous, super-Chandrasekhar (SC) Type Ia Supernovae (SNe Ia) are as yet an unexplained phenomenon. We analyse a well-observed SN of this class, SN 2009dc, by modelling its photospheric ...spectra with a spectral synthesis code, using the technique of 'Abundance Tomography'. We present spectral models based on different density profiles, corresponding to different explosion scenarios, and discuss their consistency. First, we use a density structure of a simulated explosion of a 2 M_sun rotating C-O white dwarf (WD), which is often proposed as a possibility to explain SC SNe Ia. Then, we test a density profile empirically inferred from the evolution of line velocities (blueshifts). This model may be interpreted as a core-collapse SN with an ejecta mass ~ 3 M_sun. Finally, we calculate spectra assuming an interaction scenario. In such a scenario, SN 2009dc would be a standard WD explosion with a normal intrinsic luminosity, and this luminosity would be augmented by interaction of the ejecta with a H-/He-poor circumstellar medium. We find that no model tested easily explains SN 2009dc. With the 2 M_sun WD model, our abundance analysis predicts small amounts of burning products in the intermediate-/high-velocity ejecta (v > 9000 km/s). However, in the original explosion simulations, where the nuclear energy release per unit mass is large, burned material is present at high v. This contradiction can only be resolved if asymmetries strongly affect the radiative transfer or if C-O WDs with masses significantly above 2 M_sun exist. In a core-collapse scenario, low velocities of Fe-group elements are expected, but the abundance stratification in SN 2009dc seems 'SN Ia-like'. The interaction-based model looks promising, and we have some speculations on possible progenitor configurations. However, radiation-hydro simulations will be needed to judge whether this scenario is realistic at all.
We present results for a suite of fourteen three-dimensional, high resolution hydrodynamical simulations of delayed-detonation modelsof Type Ia supernova (SN Ia) explosions. This model suite ...comprises the first set of three-dimensional SN Ia simulations with detailed isotopic yield information. As such, it may serve as a database for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ a physically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration to detonation transition (DDT) probabilities. To obtain different strengths of the deflagration phase and thereby different degrees of pre-expansion, we have chosen a sequence of initial models with 1, 3, 5, 10, 20, 40, 100, 150, 200, 300, and 1600 (two different realizations) ignition kernels in a hydrostatic white dwarf with central density of 2.9 x 10^9 gcc, plus in addition one high central density (5.5 x 10^9 gcc), and one low central density (1.0 x 10^9 gcc) rendition of the 100 ignition kernel configuration. For each simulation we determined detailed nucleosynthetic yields by post-processing 10^6 tracer particles with a 384 nuclide reaction network. All delayed detonation models result in explosions unbinding the white dwarf, producing a range of 56Ni masses from 0.32 to 1.11 solar masses. As a general trend, the models predict that the stable neutron-rich iron group isotopes are not found at the lowest velocities, but rather at intermediate velocities (~3,000 - 10,000 km/s) in a shell surrounding a 56Ni-rich core. The models further predict relatively low velocity oxygen and carbon, with typical minimum velocities around 4,000 and 10,000 km/s, respectively.