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 similar to 1.1 M-aS (TM), hybrid WDs in a binary system may easily approach the Chandrasekhar mass (M-Ch) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near-M-Ch 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 M-aS (TM) 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-M-Ch WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near-M-Ch 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 M-Ch to be of the order of 1 per cent of the Galactic SN Ia rate.

The impact of accretion disc winds on the optical spectra of cataclysmic variables

Many high-state non-magnetic cataclysmic variables (CVs) exhibit blueshifted absorption or P-Cygni profiles associated with ultraviolet (UV) resonance lines. These features imply the existence of powerful accretion disc winds in CVs. Here, we use our Monte Carlo ionization and radiative transfer code to investigate whether disc wind models that produce realistic UV line profiles are also likely to generate observationally significant recombination line and continuum emission in the optical waveband. We also test whether outflows may be responsible for the single-peaked emission line profiles often seen in high-state CVs and for the weakness of the Balmer absorption edge (relative to simple models of optically thick accretion discs). We find that a standard disc wind model that is successful in reproducing the UV spectra of CVs also leaves a noticeable imprint on the optical spectrum, particularly for systems viewed at high inclination. The strongest optical wind-formed recombination lines are H alpha and He ii lambda 4686. We demonstrate that a higher density outflow model produces all the expected H and He lines and produces a recombination continuum that can fill in the Balmer jump at high inclinations. This model displays reasonable verisimilitude with the optical spectrum of RW Trianguli. No single-peaked emission is seen, although we observe a narrowing of the double-peaked emission lines from the base of the wind. Finally, we show that even denser models can produce a single-peaked H alpha line. On the basis of our results, we suggest that winds can modify, and perhaps even dominate, the line and continuum emission from CVs.

Polarization spectral synthesis for Type Ia supernova explosion models

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multidimensional supernova explosion models. The approach utilizes 'virtual-packets' that are generated during the propagation of the Monte Carlo quanta and used to compute synthetic observables for specific observer orientations. Compared to extracting synthetic observables by direct binning of emergent Monte Carlo quanta, this virtual-packet approach leads to a substantial reduction in the Monte Carlo noise. This is not only vital for calculating synthetic spectropolarimetry (since the degree of polarization is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealized test code to simple geometries. We then describe its implementation in the Monte Carlo radiative transfer code ARTIS and present test calculations for simple models for Type Ia supernovae. Specifically, we use the well-known one-dimensional W7 model to verify that our scheme can accurately recover zero polarization from a spherical model, and to demonstrate the reduction in Monte Carlo noise compared to a simple packet-binning approach. To investigate the impact of aspherical ejecta on the polarization spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.

Late-time spectral line formation in Type IIb supernovae, with application to SN 1993J, SN 2008ax, and SN 2011dh

We investigate line formation processes in Type IIb supernovae (SNe) from 100 to 500 days post-explosion using spectral synthesis calculations. The modelling identifies the nuclear burning layers and physical mechanisms that produce the major emission lines, and the diagnostic potential of these. We compare the model calculations with data on the three best observed Type IIb SNe to-date - SN 1993J, SN 2008ax, and SN 2011dh. Oxygen nucleosynthesis depends sensitively on the main-sequence mass of the star and modelling of the [O I] lambda lambda 6300, 6364 lines constrains the progenitors of these three SNe to the M-ZAMS = 12-16 M-circle dot range (ejected oxygen masses 0.3-0.9 M-circle dot), with SN 2011dh towards the lower end and SN 1993J towards the upper end of the range. The high ejecta masses from M-ZAMS greater than or similar to 17 M-circle dot progenitors give rise to brighter nebular phase emission lines than observed. Nucleosynthesis analysis thus supports a scenario of low-to-moderate mass progenitors for Type IIb SNe, and by implication an origin in binary systems. We demonstrate how oxygen and magnesium recombination lines may be combined to diagnose the magnesium mass in the SN ejecta. For SN 2011dh, a magnesium mass of 0.02-0.14 M-circle dot is derived, which gives a Mg/O production ratio consistent with the solar value. Nitrogen left in the He envelope from CNO burning gives strong [N II] lambda lambda 6548, 6583 emission lines that dominate over Ha emission in our models. The hydrogen envelopes of Type IIb SNe are too small and dilute to produce any noticeable H alpha emission or absorption after similar to 150 days, and nebular phase emission seen around 6550 angstrom is in many cases likely caused by [N II] lambda lambda 6548, 6583. Finally, the influence of radiative transport on the emergent line profiles is investigated. Significant line blocking in the metal core remains for several hundred days, which affects the emergent spectrum. These radiative transfer effects lead to early-time blueshifts of the emission line peaks, which gradually disappear as the optical depths decrease with time. The modelled evolution of this effect matches the observed evolution in SN 2011dh.

Developed turbulence and nonlinear amplification of magnetic fields in laboratory and astrophysical plasmas

The visible matter in the universe is turbulent and magnetized. Turbulence in galaxy clusters is produced by mergers and by jets of the central galaxies and believed responsible for the amplification of magnetic fields. We report on experiments looking at the collision of two laser-produced plasma clouds, mimicking, in the laboratory, a cluster merger event. By measuring the spectrum of the density fluctuations, we infer developed, Kolmogorov-like turbulence. From spectral line broadening, we estimate a level of turbulence consistent with turbulent heating balancing radiative cooling, as it likely does in galaxy clusters. We show that the magnetic field is amplified by turbulent motions, reaching a nonlinear regime that is a precursor to turbulent dynamo. Thus, our experiment provides a promising platform for understanding the structure of turbulence and the amplification of magnetic fields in the universe.

The nebular spectra of SN 2012aw and constraints on stellar nucleosynthesis from oxygen emission lines

We present nebular-phase optical and near-infrared spectroscopy of the Type IIP supernova SN 2012aw combined with non-local thermodynamic equilibrium radiative transfer calculations applied to ejecta from stellar evolution/explosion models. Our spectral synthesis models generally show good agreement with the ejecta from a M_ZAMS = 15 M_⊙progenitor star. The emission lines of oxygen, sodium, and magnesium are all consistent with the nucleosynthesis in a progenitor in the 14-18 M_⊙ range.We also demonstrate how the evolution of the oxygen cooling lines of [O I] λ5577, [O I] λ6300, and [O I] λ6364 can be used to constrain the mass of oxygen in the non-molecularly cooled ashes to < 1 M_⊙, independent of the mixing in the ejecta. This constraint implies that any progenitor model of initial mass greater than 20 M_⊙ would be difficult to reconcile with the observed line strengths. A stellar progenitor of around M_ZAMS = 15 M_⊙ can consistently explain the directly measured luminosity of the progenitor star, the observed nebular spectra, and the inferred pre-supernova mass-loss rate.We conclude that there is still no convincing example of a Type IIP supernova showing the nucleosynthesis products expected from an M_ZAMS > 20 M_⊙ progenitor. © 2014 The Author. Published by Oxford University Press on behalf of the Royal Astronomical Society.

Constraining the physical properties of Type II-Plateau supernovae using nebular phase spectra

We present a study of the nebular phase spectra of a sample of Type II-Plateau supernovae with identified progenitors or restrictive limits. The evolution of line fluxes, shapes and velocities is compared within the sample, and interpreted by the use of a spectral synthesis code. The small diversity within the data set can be explained by strong mixing occurring during the explosion, and by recognizing that most lines have significant contributions from primordial metals in the H envelope, which dominates the total ejecta mass in these types of objects. In particular, when using the [O I] 6300, 6364 Å doublet for estimating the core mass of the star, care has to be taken to account for emission from primordial O in the envelope. Finally, a correlation between the Hα line width and the mass of 56Ni is presented, suggesting that higher energy explosions are associated with higher 56Ni production.

On the diversity of superluminous supernovae: ejected mass as the dominant factor

We assemble a sample of 24 hydrogen-poor superluminous supernovae(SLSNe). Parameterizing the light-curve shape through rise and declinetime-scales shows that the two are highly correlated. Magnetar-poweredmodels can reproduce the correlation, with the diversity in rise anddecline rates driven by the diffusion time-scale. Circumstellarinteraction models can exhibit a similar rise-decline relation, but onlyfor a narrow range of densities, which may be problematic for thesemodels. We find that SLSNe are approximately 3.5 mag brighter and havelight curves three times broader than SNe Ibc, but that the intrinsicshapes are similar. There are a number of SLSNe with particularly broadlight curves, possibly indicating two progenitor channels, butstatistical tests do not cleanly separate two populations. The generalspectral evolution is also presented. Velocities measured from Fe II aresimilar for SLSNe and SNe Ibc, suggesting that diffusion timedifferences are dominated by mass or opacity. Flat velocity evolution inmost SLSNe suggests a dense shell of ejecta. If opacities in SLSNe aresimilar to other SNe Ibc, the average ejected mass is higher by a factor2-3. Assuming κ = 0.1 cm2 g-1, we estimate amean (median) SLSN ejecta mass of 10 M⊙ (6M⊙), with a range of 3-30 M⊙. Doubling theassumed opacity brings the masses closer to normal SNe Ibc, but with ahigh-mass tail. The most probable mechanism for generating SLSNe seemsto be the core collapse of a very massive hydrogen-poor star, forming amillisecond magnetar.

Type Ia supernovae from exploding oxygen-neon white dwarfs

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

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

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

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

Conclusions. 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 (similar to 1.2 M-circle dot) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.

Self-consistent modelling of line-driven hot-star winds with Monte Carlo radiation hydrodynamics

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multiline effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach are demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally, we discuss some of the current limitations of our method, some possible extensions and potential future applications.

5.9-keV Mn K-shell X-ray luminosity from the decay of 55Fe in Type Ia supernova models

We show that the X-ray line flux of the Mn Kα line at 5.9 keV from the decay of 55Fe is a promising diagnostic to distinguish between Type Ia supernova (SN Ia) explosion models. Using radiation transport calculations, we compute the line flux for two three-dimensional explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two (1.1 and 0.9 M⊙) white dwarfs. Both models are based on solar metallicity zero-age main-sequence progenitors. Due to explosive nuclear burning at higher density, the delayed-detonation model synthesizes ˜3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kα line flux of the delayed-detonation model exceeds that of the merger model by a factor of ˜4.5. Since in both models the 5.9-keV X-ray flux peaks five to six years after the explosion, a single measurement of the X-ray line emission at this time can place a constraint on the explosion physics that is complementary to those derived from earlier phase optical spectra or light curves. We perform detector simulations of current and future X-ray telescopes to investigate the possibilities of detecting the X-ray line at 5.9 keV. Of the currently existing telescopes, XMM-Newton/pn is the best instrument for close (≲1-2 Mpc), non-background limited SNe Ia because of its large effective area. Due to its low instrumental background, Chandra/ACIS is currently the best choice for SNe Ia at distances above ˜2 Mpc. For the delayed-detonation scenario, a line detection is feasible with Chandra up to ˜3 Mpc for an exposure time of 106 s. We find that it should be possible with currently existing X-ray instruments (with exposure times ≲5 × 105 s) to detect both of our models at sufficiently high S/N to distinguish between them for hypothetical events within the Local Group. The prospects for detection will be better with future missions. For example, the proposed Athena/X-IFU instrument could detect our delayed-detonation model out to a distance of ˜5 Mpc. This would make it possible to study future events occurring during its operational life at distances comparable to those of the recent supernovae SN 2011fe (˜6.4 Mpc) and SN 2014J (˜3.5 Mpc).

The white dwarf's carbon fraction as a secondary parameter of Type Ia supernovae

Context. Binary stellar evolution calculations predict thatChandrasekhar-mass carbon/oxygen white dwarfs (WDs) show a radiallyvarying profile for the composition with a carbon depleted core. Manyrecent multi-dimensional simulations of Type Ia supernovae (SNe Ia),however, assume the progenitor WD has a homogeneous chemicalcomposition.
Aims: In this work, we explore the impact ofdifferent initial carbon profiles of the progenitor WD on the explosionphase and on synthetic observables in the Chandrasekhar-mass delayeddetonation model. Spectra and light curves are compared to observationsto judge the validity of the model.
Methods: The explosion phaseis simulated using the finite volume supernova code Leafs, which isextended to treat different compositions of the progenitor WD. Thesynthetic observables are computed with the Monte Carlo radiativetransfer code Artis. Results: Differences in binding energies ofcarbon and oxygen lead to a lower nuclear energy release for carbondepleted material; thus, the burning fronts that develop are weaker andthe total nuclear energy release is smaller. For otherwise identicalconditions, carbon depleted models produce less 56Ni.Comparing different models with similar 56Ni yields showslower kinetic energies in the ejecta for carbon depleted models, butonly small differences in velocity distributions and line velocities inspectra. The light curve width-luminosity relation (WLR) obtained formodels with differing carbon depletion is roughly perpendicular to theobserved WLR, hence the carbon mass fraction is probably only asecondary parameter in the family of SNe Ia.
Tables 3 and 4 are available in electronic form at http://www.aanda.org

Gamma-ray diagnostics of Type Ia supernovae. Predictions of observables from three-dimensional modeling

Context. Although the question of progenitor systems and detailed explosion mechanisms still remains a matter of discussion, it is commonly believed that Type Ia supernovae (SNe Ia) are production sites of large amounts of radioactive nuclei. Even though the gamma-ray emission due to radioactive decays is responsible for powering the light curves of SNe Ia, gamma rays themselves are of particular interest as a diagnostic tool because they directly lead to deeper insight into the nucleosynthesis and the kinematics of these explosion events. Aims: We study the evolution of gamma-ray line and continuum emission of SNe Ia with the objective of analyzing the relevance of observations in this energy range. We seek to investigate the chances for the success of future MeV missions regarding their capabilities for constraining the intrinsic properties and the physical processes of SNe Ia. Methods: Focusing on two of the most broadly discussed SN Ia progenitor scenarios - a delayed detonation in a Chandrasekhar-mass white dwarf (WD) and a violent merger of two WDs - we used three-dimensional explosion models and performed radiative transfer simulations to obtain synthetic gamma-ray spectra. Both chosen models produce the same mass of 56Ni and have similar optical properties that are in reasonable agreement with the recently observed supernova SN 2011fe. We examine the gamma-ray spectra with respect to their distinct features and draw connections to certain characteristics of the explosion models. Applying diagnostics, such as line and hardness ratios, the detection prospects for future gamma-ray missions with higher sensitivities in the MeV energy range are discussed. Results: In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be quite different. The almost direct connection of the emission of gamma rays to fundamental physical processes occurring in SNe Ia permits additional constraints concerning several explosion model properties that are not easily accessible within other wavelength ranges. Proposed future MeV missions such as GRIPS will resolve all spectral details only for nearby SNe Ia, but hardness ratio and light curve measurements still allow for a distinction of the two different models at 10 Mpc and 16 Mpc for an exposure time of 106 s. The possibility of detecting the strongest line features up to the Virgo distance will offer the opportunity to build up a first sample of SN Ia detections in the gamma-ray energy range and underlines the importance of future space observatories for MeV gamma rays.

Self-regulated Shocks in Massive Star Binary Systems

In an early-type, massive star binary system, X-ray bright shocks result from the powerful collision of stellar winds driven by radiation pressure on spectral line transitions. We examine the influence of the X-rays from the wind-wind collision shocks on the radiative driving of the stellar winds using steady-state models that include a parameterized line force with X-ray ionization dependence. Our primary result is that X-ray radiation from the shocks inhibits wind acceleration and can lead to a lower pre-shock velocity, and a correspondingly lower shocked plasma temperature, yet the intrinsic X-ray luminosity of the shocks, L X, remains largely unaltered, with the exception of a modest increase at small binary separations. Due to the feedback loop between the ionizing X-rays from the shocks and the wind driving, we term this scenario as self-regulated shocks. This effect is found to greatly increase the range of binary separations at which a wind-photosphere collision is likely to occur in systems where the momenta of the two winds are significantly different. Furthermore, the excessive levels of X-ray ionization close to the shocks completely suppress the line force, and we suggest that this may render radiative braking less effective. Comparisons of model results against observations reveal reasonable agreement in terms of log (L X/L bol). The inclusion of self-regulated shocks improves the match for kT values in roughly equal wind momenta systems, but there is a systematic offset for systems with unequal wind momenta (if considered to be a wind-photosphere collision).

Generation of cluster states in optomechanical quantum systems

We consider an optomechanical quantum system composed of a single cavity mode interacting with N mechanical resonators. We propose a scheme for generating continuous-variable graph states of arbitrary size and shape, including the so-called cluster states for universal quantum computation. The main feature of this scheme is that, differently from previous approaches, the graph states are hosted in the mechanical degrees of freedom rather than in the radiative ones. Specifically, via a 2N-tone drive, we engineer a linear Hamiltonian which is instrumental to dissipatively drive the system to the desired target state. The robustness of this scheme is assessed against finite interaction times and mechanical noise, confirming it as a valuable approach towards quantum state engineering for continuous-variable computation in a solid-state platform.