Common envelope (CE) phases in binary systems where the primary star reaches the tip of the red giant branch are discussed as a formation scenario for hot subluminous B-type (sdB) stars. For some of these objects, observations point to very low-mass companions. In hydrodynamical CE simulations with the moving-mesh code AREPO, we test whether low-mass objects can successfully unbind the envelope. The success of envelope removal in our simulations critically depends on whether or not the ionization energy released by recombination processes in the expanding material is taken into account. If this energy is thermalized locally, envelope ejection eventually leading to the formation of an sdB star is possible with companion masses down to the brown dwarf range. For even lower companion masses approaching the regime of giant planets, however, envelope removal becomes increasingly difficult or impossible to achieve. Our results are consistent with current observational constraints on companion masses of sdB stars. Based on a semianalytic model, we suggest a new criterion for the lowest companion mass that is capable of triggering a dynamical response of the primary star thus potentially facilitating the ejection of a common envelope. This gives an estimate consistent with the findings of our hydrodynamical simulations.

About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we follow the long-term evolution of such a 3D merger product in a 1D stellar evolution code. During a thermal relaxation phase after the coalescence, the merger product reaches critical surface rotation, sheds mass and then spins down primarily because of internal mass readjustments. The spin of the merger product after thermal relaxation is mainly set by the co-evolution of the star-torus structure left after coalescence. This evolution is still uncertain, so we also consider magnetic braking and other angular-momentum-gain and -loss mechanisms that may influence the final spin of the merged star. Because of core compression and mixing of carbon and nitrogen in the merger, enhanced nuclear burning drives a transient convective core that greatly contributes to the rejuvenation of the star. Once the merger product relaxed back to the main sequence, it continues its evolution similar to that of a genuine single star of comparable mass. It is a slow rotator that matches the magnetic blue straggler $\tau$ Sco. Our results show that merging is a promising mechanism to explain some magnetic massive stars and it may also be key to understand the origin of the strong magnetic fields of highly-magnetic white dwarfs and magnetars.

Sub-Chandrasekhar mass white dwarfs accreting a helium shell on a carbon-oxygen core are potential progenitors of normal Type Ia supernovae. This work focuses on the details of the onset of the carbon detonation in the double detonation sub-Chandrasekhar model. In order to simulate the influence of core-shell mixing on the carbon ignition mechanism, the helium shell and its detonation are followed with an increased resolution compared to the rest of the star treating the propagation of the detonation wave more accurately. This significantly improves the predictions of the nucleosynthetic yields from the helium burning. The simulations were carried out with the AREPO code. A carbon-oxygen core with a helium shell was set up in one dimension and mapped to three dimensions. We ensured the stability of the white dwarf with a relaxation step before the hydrodynamic detonation simulation started. Synthetic observables were calculated with the radiative transfer code ARTIS. An ignition mechanism of the carbon detonation was observed, which received little attention before. In this "scissors mechanism", the impact the helium detonation wave has on unburnt material when converging opposite to its ignition spot is strong enough to ignite a carbon detonation. This is possible in a carbon enriched transition region between the core and shell. The detonation mechanism is found to be sensitive to details of the core-shell transition and our models illustrate the need to consider core-shell mixing taking place during the accretion process. Even though the detonation ignition mechanism differs form the converging shock mechanism, the differences in the synthetic observables are not significant. Though they do not fit observations better than previous simulations, they illustrate the need for multi-dimensional simulations.

A significant fraction of stars between 7-11 solar masses are thought to become supernovae, but the explosion mechanism is unclear. The answer depends critically on the rate of electron capture on $^{20}$Ne in the degenerate oxygen-neon stellar core. However, due to the unknown strength of the transition between the ground states of $^{20}$Ne and $^{20}$F, it has not previously been possible to fully constrain the rate. By measuring the transition, we have established that its strength is exceptionally large and enhances the capture rate by several orders of magnitude. This has a decisive impact on the evolution of the core, increasing the likelihood that the star is (partially) disrupted by a thermonuclear explosion rather than collapsing to form a neutron star. Importantly, our measurement resolves the last remaining nuclear physics uncertainty in the final evolution of degenerate oxygen-neon stellar cores, allowing future studies to address the critical role of convection, which at present is poorly understood.

Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars (<1.5 solar masses) are generated in-situ by a dynamo process in their turbulent, convective envelopes. Unlike such stars, intermediate-mass and high-mass stars (>1.5 solar masses; referred to as "massive" stars here) have relatively quiet, radiative envelopes where a solar-like dynamo cannot operate. However, about 10% of them, including 78 Vir, have strong, large-scale surface magnetic fields whose origin has remained a major mystery. The massive star $\tau$ Sco is a prominent member of this group and appears to be surprisingly young compared to other presumably coeval members of the Upper Scorpius association. Here, we present the first 3D magneto-hydrodynamical simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain $\tau$ Sco's magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. These stars are promising progenitors of those neutron stars that host the strongest magnetic fields in the Universe, so-called magnetars, and that may give rise to some of the enigmatic fast radio bursts. Strong magnetic fields affect the explosions of core-collapse supernovae and, moreover, those magnetic stars that have rapidly-rotating cores at the end of their lives might provide the right conditions to power long-duration gamma-ray bursts and super-luminous supernovae.

When a star approaches a black hole closely, it may be pulled apart by gravitational forces in a tidal disruption event (TDE). The flares produced by TDEs are unique tracers of otherwise quiescent supermassive black holes (SMBHs) located at the centre of most galaxies. In particular, the appearance of such flares and the subsequent decay of the light curve are both sensitive to whether the star is partially or totally destroyed by the tidal field. However, the physics of the disruption and the fall-back of the debris are still poorly understood. We are here modelling the hydrodynamical evolution of realistic stars as they approach a SMBH on parabolic orbits, using for the first time the moving-mesh code AREPO, which is particularly well adapted to the problem through its combination of quasi-Lagrangian behaviour, low advection errors, and high accuracy typical of mesh-based techniques. We examine a suite of simulations with different impact parameters, allowing us to determine the critical distance at which the star is totally disrupted, the energy distribution and the fallback rate of the debris, as well as the hydrodynamical evolution of the stellar remnant in the case of a partial disruption. Interestingly, we find that the internal evolution of the remnant's core is strongly influenced by persistent vortices excited in the tidal interaction. These should be sites of strong magnetic field amplification, and the associated mixing may profoundly alter the subsequent evolution of the tidally pruned star.

(Abridged) The explosion mechanism of electron-capture supernovae (ECSNe) remains equivocal. We attempt to constrain the explosion mechanism (neutron-star-forming implosion or thermonuclear explosion) and the frequency of occurrence of ECSNe using nucleosynthesis simulations of the latter scenario, population synthesis, the solar abundance distribution, pre-solar meteoritic oxide grain isotopic ratio measurements and the white dwarf mass-radius relation. Tracer particles from 3d hydrodynamic simulations were post-processed with a large nuclear reaction network in order to determine the complete compositional state of the bound ONeFe remnant and the ejecta, and population synthesis simulations were performed in order to estimate the ECSN rate with respect to the CCSN rate. The 3d deflagration simulations drastically overproduce the neutron-rich isotopes $^{48}$Ca, $^{50}$Ti, $^{54}$Cr, $^{60}$Fe and several of the Zn isotopes relative to their solar abundances. Using the solar abundance distribution as our constraint, we place an upper limit on the frequency of thermonuclear ECSNe as 1$-$3~\% the frequency at which core-collapse supernovae (FeCCSNe) occur. This is on par with or 1~dex lower than the estimates for ECSNe from single stars. The upper limit from the yields is also in relatively good agreement with the predictions from our population synthesis simulations. The $^{54}$Cr/$^{52}$Cr and $^{50}$Ti/$^{48}$Ti isotopic ratios in the ejecta are a near-perfect match with recent measurements of extreme pre-solar meteoritc oxide grains, and $^{53}$Cr/$^{52}$Cr can also be matched if the ejecta condenses before mixing with the interstellar medium. Theoretical mass-radius relations for the bound ONeFe WD remnants of these explosions are apparently consistent with several observational WD candidates.

Despite the importance of Type Ia supernovae (SNe Ia) for modern astrophysics, their detailed mechanism is still not fully understood. In this contribution, we present recent findings from numerical explosion models in the context of the observed diversity of SNe Ia and we discuss how these models can help to shed light on the explosion mechanism and the progenitor stars of SNe Ia. In addition, we introduce the Heidelberg Supernova Model Archive (HESMA), a new online data base where we provide integrated isotopic abundances and radially averaged ejecta profiles and synthetic observables for a wide range of state-of-the-art explosion models.

Despite their use as cosmological distance indicators and their importance in the chemical evolution of galaxies, the unequivocal identification of the progenitor systems and explosion mechanism of normal type Ia supernova (SN Ia) remains elusive. The leading hypothesis is that such a supernova is a thermonuclear explosion of a carbon-oxygen white dwarf, but the exact explosion mechanism is still a matter of debate. Observation of a galactic SN Ia would be of immense value in answering the many open questions related to these events. One potentially useful source of information about the explosion mechanism and progenitor is the neutrino signal. In this paper we compute the expected neutrino signal from a gravitationally confined detonation (GCD) explosion scenario for a SN~Ia and show how the flux at Earth contains features in time and energy unique to this scenario. We then calculate the expected event rates in the Super-K, Hyper-K, JUNO, DUNE, and IceCube detectors and find both Hyper-K and IceCube would see a few events for a GCD supernova at 1 kpc or closer, while Super-K, JUNO, and DUNE would see a events if the supernova were closer than ${\sim}0.3$ kpc. The distance and detector criteria needed to resolve the time and spectral features arising from the explosion mechanism, neutrino production, and neutrino oscillation processes are also discussed. The neutrino signal from the GCD is then compared with the signal from a deflagration-to-detonation transition (DDT) explosion model computed previously. We find the overall event rate is the most discriminating feature between the two scenarios followed by the event rate time structure. Using the event rate in the Hyper-K detector alone, the DDT can be distinguished from the GCD at 2$\sigma$ if the distance to the supernova is less than $2.3\;{\rm kpc}$ for a normal mass ordering and $3.6\;{\rm kpc}$ for an inverted ordering.

Hydrodynamical simulations of stellar interactions require stable models of stars as initial conditions. Such initial models, however, are difficult to construct for giant stars because of the wide range in spatial scales of the hydrostatic equilibrium and in dynamical timescales between the core and the envelope of the giant. They are needed for, e.g., modeling the common envelope phase where a giant envelope encompasses both the giant core and a companion star. Here, we present a new method of approximating and reconstructing giant profiles from a stellar evolution code to produce stable models for multi-dimensional hydrodynamical simulations. We determine typical stellar stratification profiles with the 1D stellar evolution code MESA. After an appropriate mapping, hydrodynamical simulations are conducted using the moving-mesh code AREPO. The giant profiles are approximated by replacing the core of the giant with a point mass and by constructing a suitable continuation of the profile to the center. Different reconstruction methods are tested that can specifically control the convective behaviour of the model. After mapping to a grid, a relaxation procedure that includes damping of spurious velocities yields stable models in three-dimensional hydrodynamical simulations. Initially convectively stable configurations lead to stable hydrodynamical models while for stratifications that are convectively unstable in the stellar evolution code, simulations recover the convective behaviour of the initial model and show large convective plumes with Mach numbers up to 0.8. Examples are shown for a $2M_\odot$ red giant and a $0.67M_\odot$ asymptotic giant branch star. A detailed analysis shows that the improved method reliably provides stable models of giant envelopes that can be used as initial conditions for subsequent hydrodynamical simulations of stellar interactions involving giant stars.

During the common envelope (CE) phase, a giant star in a binary system overflows its Roche lobe and unstable mass transfer leads to a spiral-in of the companion, resulting in a close binary system or in a merger of the stellar cores. Dynamo processes during the CE phase have been proposed as a mechanism to generate magnetic fields that are important for forming magnetic white dwarfs (MWDs) and for shaping planetary nebulae. Here, we present the first magnetohydrodynamics simulations of the dynamical spiral-in during a CE phase. We find that magnetic fields are strongly amplified in the accretion stream around the $1M_\odot$ companion as it spirals into the envelope of a $2M_\odot$ RG. This leads to field strengths of 10 to 100 kG throughout the envelope after 120 d. The magnetic field amplification is consistent with being driven by the magnetorotational instability. The field strengths reached in our simulation make the magnetic field interesting for diagnostic purposes, but they are dynamically irrelevant. They are also too small to explain the formation of the highest fields found in MWDs, but may be relevant for luminous red novae, and detecting magnetic fields in these events would support the scenario as proposed here.

In the classical picture, electron-capture supernovae and the accretion-induced collapse of oxygen-neon white dwarfs (ONeWDs) undergo an oxygen deflagration phase before gravitational collapse produces a neutron star (NS). These types of core collapse events are postulated to explain several astronomical phenomena. In this work, the deflagration phase is simulated for the first time using multidimensional hydrodynamics, with the aim of gaining new insight into the explosive deaths of $8-10~M_\odot$ stars and ONeWDs that accrete material from a binary companion star. The main aim is to determine whether these events are thermonuclear or core-collapse supernova explosions, and hence whether NSs are formed by such phenomena. The deflagration is simulated in ONe cores with three different central ignition densities. The intermediate density case is perhaps the most realistic, being based on recent nuclear physics calculations and 1D stellar models. The 3D hydrodynamic simulations presented in this work begin from a centrally confined flame structure using a level-set-based flame approach and are performed in $256^3$ and $512^3$ numerical resolutions. In the simulations with intermediate and low ignition density, the cores do not appear to collapse into NSs. Instead, almost a solar mass of material becomes unbound from the cores, leaving bound remnants. These simulations represent the case in which semiconvective mixing during the electron-capture phase preceding the deflagration is inefficient. The masses of the bound remnants double when Coulomb corrections are included in the EoS, however they still do not exceed the effective Chandrasekhar mass and, hence, would not collapse into NSs. The simulations with the highest ignition density ($\log_{10}\rho_{\rm c}=10.3$), representing the case where semiconvective mixing is very efficient, show clear signs that the core will collapse into a NS.

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.

There has been much debate about the origin of the diffuse $\gamma$--ray background in the MeV range. At lower energies, AGNs and Seyfert galaxies can explain the background, but not above $\simeq$0.3 MeV. Beyond $\sim$10 MeV blazars appear to account for the flux observed. That leaves an unexplained gap for which different candidates have been proposed, including annihilations of WIMPS. One candidate are Type Ia supernovae (SNe Ia). Early studies concluded that they were able to account for the $\gamma$--ray background in the gap, while later work attributed a significantly lower contribution to them. All those estimates were based on SN Ia explosion models which did not reflect the full 3D hydrodynamics of SNe Ia explosions. In addition, new measurements obtained since 2010 have provided new, direct estimates of high-z SNe Ia rates beyond $z\sim$2. We take into account these new advances to see the predicted contribution to the gamma--ray background. We use here a wide variety of explosion models and a plethora of new measurements of SNe Ia rates. SNe Ia still fall short of the observed background. Only for a fit, which would imply $\sim$150\% systematic error in detecting SNe Ia events, do the theoretical predictions approach the observed fluxes. This fit is, however, at odds at the highest redshifts with recent SN Ia rates estimates. Other astrophysical sources such as FSRQs do match the observed flux levels in the MeV regime, while SNe Ia make up to 30--50\% of the observed flux.

The progenitor system(s) and the explosion mechanism(s) of Type Ia supernovae (SNe Ia) are still under debate. Non-electromagnetic observables, in particular gravitational waves and neutrino emission, of thermonuclear supernovae are a complementary window to light curves and spectra for studying these enigmatic objects. A leading model for SNe Ia is the thermonuclear incineration of a near-Chandrasekhar mass carbon-oxygen white dwarf star in a "delayed-detonation". We calculate a three-dimensional hydrodynamic explosion for the N100 delayed-detonation model extensively discussed in the literature, taking the dynamical effects of neutrino emission from all important contributing source terms into account. Although neutrinos carry away $2 \times 10^{49}$ erg of energy, we confirm the common view that neutrino energy losses are dynamically not very important, resulting in only a modest reduction of the final kinetic energy by two per cent. We then calculate the gravitational wave signal from the time evolution of the quadrupole moment. Our model radiates $7 \times 10^{39}$ erg in gravitational waves and the spectrum has a pronounced peak around 0.4 Hz. Depending on viewing angle and polarization, we find that the future space-based gravitational wave missions DECIGO and BBO would be able to detect our source to a distance of 1.3 Mpc. We predict a clear signature of the deflagration-to-detonation transition in the neutrino and the gravitational wave signals. If observed, such a feature would be a strong indicator of the realization of delayed-detonations in near-Chandrasekhar mass white dwarfs.

Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.

The common envelope (CE) phase is an important stage in binary stellar evolution. It is needed to explain many close binary stellar systems, such as cataclysmic variables, Type Ia supernova progenitors, or X-ray binaries. To form the resulting close binary, the initial orbit has to shrink, thereby transferring energy to the primary giant's envelope that is hence ejected. The details of this interaction, however, are still not understood. Here, we present new hydrodynamic simulations of the dynamical spiral-in forming a CE system. We apply the moving-mesh code AREPO to follow the interaction of a $1M_\odot$ compact star with a $2M_\odot$ red giant possessing a $0.4M_\odot$ core. The nearly Lagrangian scheme combines advantages of smoothed particle hydrodynamics and traditional grid-based hydrodynamic codes and allows us to capture also small flow features at high spatial resolution. Our simulations reproduce the initial transfer of energy and angular momentum from the binary core to the envelope by spiral shocks seen in previous studies, but after about 20 orbits a new phenomenon is observed. Large-scale flow instabilities are triggered by shear flows between adjacent shock layers. These indicate the onset of turbulent convection in the common envelope, thus altering the transport of energy on longer time scales. At the end of our simulation, only 8% of the envelope mass is ejected. The failure to unbind the envelope completely may be caused by processes on thermal time scales or unresolved microphysics.

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.

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.

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.