Normal type Ia supernovae from violent mergers of white dwarf binaries

One of the most important questions regarding the progenitor systems of Type Ia supernovae (SNe Ia) is whether mergers of two white dwarfs can lead to explosions that reproduce observations of normal events. Here we present a fully three-dimensional simulation of a violent merger of two carbon-oxygen white dwarfs with masses of 0.9 M and 1.1 M combining very high resolution and exact initial conditions. A well-tested combination of codes is used to study the system. We start with the dynamical inspiral phase and follow the subsequent thermonuclear explosion under the plausible assumption that a detonation forms in the process of merging. We then perform detailed nucleosynthesis calculations and radiative transfer simulations to predict synthetic observables from the homologously expanding supernova ejecta. We find that synthetic color light curves of our merger, which produces about 0.62 M of Ni, show good agreement with those observed for normal SNe Ia in all wave bands from U to K. Line velocities in synthetic spectra around maximum light also agree well with observations. We conclude that violent mergers of massive white dwarfs can closely resemble normal SNe Ia. Therefore, depending on the number of such massive systems available these mergers may contribute at least a small fraction to the observed population of normal SNe Ia. © 2012 The American Astronomical Society. All rights reserved.

Monte Carlo radiation hydrodynamics:Methods, tests and application to Type Ia supernova ejecta

In astrophysical systems, radiation-matter interactions are important in transferring energy and momentum between the radiation field and the surrounding material. This coupling often makes it necessary to consider the role of radiation when modelling the dynamics of astrophysical fluids. During the last few years, there have been rapid developments in the use of Monte Carlo methods for numerical radiative transfer simulations. Here, we present an approach to radiation hydrodynamics that is based on coupling Monte Carlo radiative transfer techniques with finite-volume hydrodynamical methods in an operator-split manner. In particular, we adopt an indivisible packet formalism to discretize the radiation field into an ensemble of Monte Carlo packets and employ volume-based estimators to reconstruct the radiation field characteristics. In this paper the numerical tools of this method are presented and their accuracy is verified in a series of test calculations. Finally, as a practical example, we use our approach to study the influence of the radiation-matter coupling on the homologous expansion phase and the bolometric light curve of Type Ia supernova explosions. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

On the brightness distribution of type ia supernovae from violent white dwarf mergers

We investigate the brightness distribution expected for thermonuclear explosions that might result from the ignition of a detonation during the violent merger of white dwarf (WD) binaries. Violent WD mergers are a subclass of the canonical double degenerate scenario where two carbon-oxygen (CO) WDs merge when the larger WD fills its Roche lobe. Determining their brightness distribution is critical for evaluating whether such an explosion model could be responsible for a significant fraction of the observed population of Type Ia supernovae (SNe Ia). We argue that the brightness of an explosion realized via the violent merger model is mainly determined by the mass of Ni produced in the detonation of the primary COWD. To quantify this link, we use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass (m) and expected peak bolometric brightness (M). We use this m-M relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We also investigate the sensitivity of our results to assumptions about the conditions required to realize a detonation during violent mergers ofWDs. We find a striking similarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a M distribution that roughly covers the range and matches the shape of the one observed for SNe Ia. However, this agreement hinges on a particular phase of mass accretion during binary evolution: the primary WD gains ~0.15-0.35M? from a slightly evolved helium star companion. In our standard binary evolution model, such an accretion phase is predicted to occur for about 43 per cent of all binary systems that ultimately give rise to binary CO WD mergers. We also find that with high probability, violent WD mergers involving the most massive primaries (?1.3M?, which should produce bright SNe) have delay times ?500 Myr. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

Three-dimensional delayed-detonation models with nucleosynthesis for type ia supernovae

We present results for a suite of 14 three-dimensional, high-resolution hydrodynamical simulations of delayed-detonation models of 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 data base for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ aphysically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration-to-detonation transition 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 a central density of 2.9 × 10 g cm, as well as one high central density (5.5 × 10 g cm) and one low central density (1.0 × 10 g cm) rendition of the 100 ignition kernel configuration. For each simulation, we determined detailed nucleosynthetic yields by postprocessing10 tracer particles with a 384 nuclide reaction network. All delayed-detonation models result in explosions unbinding thewhite dwarf, producing a range of 56Ni masses from 0.32 to 1.11M. As a general trend, the models predict that the stableneutron-rich iron-group isotopes are not found at the lowest velocities, but rather at intermediate velocities (~3000×10 000 km s) in a shell surrounding a Ni-rich core. The models further predict relatively low-velocity oxygen and carbon, with typical minimum velocities around 4000 and 10 000 km s, respectively. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

Sub-luminous type Ia supernovae from the mergers of equal-mass white dwarfs with mass 0.9M

Type Ia supernovae are thought to result from thermonuclear explosions of carbong'oxygen white dwarf stars. Existing models generally explain the observed properties, with the exception of the sub-luminous 1991bg-like supernovae. It has long been suspected that the merger of two white dwarfs could give rise to a type Ia event, but hitherto simulations have failed to produce an explosion. Here we report a simulation of the merger of two equal-mass white dwarfs that leads to a sub-luminous explosion, although at the expense of requiring a single common-envelope phase, and component masses of 0.9M. The light curve is too broad, but the synthesized spectra, red colour and low expansion velocities are all close to what is observed for sub-luminous 1991bg-like events. Although the mass ratios can be slightly less than one and still produce a sub-luminous event, the masses have to be in the range 0.83M to 0.9M. © 2010 Macmillan Publishers Limited. All rights reserved.

Galactic Bulge PNe:Carbon molecules in oxygen-rich environments

Galactic bulge planetary nebulae show evidence of mixed chemistry with emission from both silicate dust and PAHs. This mixed chemistry is unlikely to be related to carbon dredge up, as third dredge-up is not expected to occur in the low mass bulge stars. We show that the phenomenon is widespread, and is seen in 30 nebulae out of our sample of 40. A strong correlation is found between strength of the PAH bands and morphology, in particular, the presence of a dense torus. A chemical model is presented which shows that hydrocarbon chains can form within oxygen-rich gas through gas-phase chemical reactions. We conclude that the mixed chemistry phenomenon occurring in the galactic bulge planetary nebulae is best explained through hydrocarbon chemistry in an UV-irradiated, dense torus. © 2012 International Astronomical Union.

Infrared non-detection of fomalhaut b: Implications for the planet interpretation

The nearby A4-type star Fomalhaut hosts a debris belt in the form of an eccentric ring, which is thought to be caused by dynamical influence from a giant planet companion. In 2008, a detection of a point source inside the inner edge of the ring was reported and was interpreted as a direct image of the planet, named Fomalhaut b. The detection was made at 600-800nm, but no corresponding signatures were found in the near-infrared range, where the bulk emission of such a planet should be expected. Here, we present deep observations of Fomalhaut with Spitzer/IRAC at 4.5 µm, using a novel point-spread function subtraction technique based on angular differential imaging and Locally Optimized Combination of Images, in order to substantially improve the Spitzer contrast at small separations. The results provide more than an order ofmagnitude improvement in the upper flux limit of Fomalhaut b and exclude the possibility that any flux from a giant planet surface contributes to the observed flux at visible wavelengths. This renders any direct connection between the observed light source and the dynamically inferred giant planet highly unlikely. We discuss several possible interpretations of the total body of observations of the Fomalhaut system and find that the interpretation that best matches the available data for the observed source is scattered light from a transient or semi-transient dust cloud. © 2012 The American Astronomical Society. All rights reserved.

The host galaxy of the super-luminous SN 2010gx and limits on explosive nickel-56 production

Super-luminous supernovae have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. While these extremely luminous explosions have been observed from z=0.1 to 1.55, the closest explosions allow more detailed investigations of their host galaxies. We present a detailed analysis of the host galaxy of SN 2010gx (z=0.23), one of the best studied super-luminous type Ic supernovae. The host is a dwarf galaxy (M_g=-17.42+/-0.17) with a high specific star formation rate. It has a remarkably low metallicity of 12+log(O/H)=7.5+/-0.1 dex as determined from the detection of the [OIII] 4363 Angs line. This is the first reliable metallicity determination of a super-luminous stripped-envelope supernova host. We collected deep multi-epoch imaging with Gemini + GMOS between 240-560 days after explosion to search for any sign of radioactive nickel-56, which might provide further insights on the explosion mechanism and the progenitor's nature. We reach griz magnitudes of m_AB~26, but do not detect SN 2010gx at these epochs. The limit implies that any nickel-56 production was similar to or below that of SN 1998bw (a luminous type Ic SN that produced around 0.4 M_sun of nickel-56). The low volumetric rates of these supernovae (~10^-4 of the core-collapse population) could be qualitatively matched if the explosion mechanism requires a combination of low-metallicity (below 0.2 Z_sun), high progenitor mass (>60 M_sun) and high rotation rate (fastest 10% of rotators).

Red and dead:The progenitor of SN 2012aw in M95

Core-collapse supernovae (SNe) are the spectacular finale to massive stellar evolution. In this Letter, we identify a progenitor for the nearby core-collapse SN 2012aw in both ground based near-infrared, and space based optical pre-explosion imaging. The SN itself appears to be a normal Type II Plateau event, reaching a bolometric luminosity of 10$^{42}$ erg s$^{-1}$ and photospheric velocities of $\sim$11,000 \kms\ from the position of the H$\beta$ P-Cygni minimum in the early SN spectra. We use an adaptive optics image to show that the SN is coincident to within 27 mas with a faint, red source in pre-explosion HST+WFPC2, VLT+ISAAC and NTT+SOFI images. The source has magnitudes $F555W$=26.70$\pm$0.06, $F814W$=23.39$\pm$0.02, $J$=21.1$\pm$0.2, $K$=19.1$\pm$0.4, which when compared to a grid of stellar models best matches a red supergiant. Interestingly, the spectral energy distribution of the progenitor also implies an extinction of $A_V>$1.2 mag, whereas the SN itself does not appear to be significantly extinguished. We interpret this as evidence for the destruction of dust in the SN explosion. The progenitor candidate has a luminosity between 5.0 and 5.6 log L/\lsun, corresponding to a ZAMS mass between 14 and 26 \msun\ (depending on $A_V$), which would make this one of the most massive progenitors found for a core-collapse SN to date.

The superluminous supernova PS1-11ap: bridging the gap between low and high redshift

We present optical photometric and spectroscopic coverage of the superluminous supernova (SLSN) PS1-11ap, discovered with the Pan-STARRS1 Medium Deep Survey at z = 0.524. This intrinsically blue transient rose slowly to reach a peak magnitude of Mu = −21.4 mag and bolometric luminosity of 8 × 1043 erg s−1 before settling on to a relatively shallow gradient of decline. The observed decline is significantly slower than those of the SLSNe-Ic which have been the focus of much recent attention. Spectroscopic similarities with the lower redshift SN2007bi and a decline rate similar to 56Co decay time-scale initially indicated that this transient could be a candidate for a pair instability supernova (PISN) explosion. Overall the transient appears quite similar to SN2007bi and the lower redshift object PTF12dam. The extensive data set, from 30 d before peak to 230 d after, allows a detailed and quantitative comparison with published models of PISN explosions. We find that the PS1-11ap data do not match these model explosion parameters well, supporting the recent claim that these SNe are not pair instability explosions. We show that PS1-11ap has many features in common with the faster declining SLSNe-Ic, and the light-curve evolution can also be quantitatively explained by the magnetar spin-down model. At a redshift of z = 0.524, the observer-frame optical coverage provides comprehensive rest-frame UV data and allows us to compare it with the SLSNe recently found at high redshifts between z = 2 and 4. While these high-z explosions are still plausible PISN candidates, they match the photometric evolution of PS1-11ap and hence could be counterparts to this lower redshift transient.