BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. A. Sim

S. A. Sim appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRastro-ph.COastro-ph.IMastro-ph.GA

Catalog footprint

What is connected

60works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

60 published item(s)

preprint2025arXiv

NLTE spectral modelling for a carbon-oxygen and helium white-dwarf merger as a Ca-rich transient candidate

We carry out NLTE (non local thermodynamic equilibrium) radiative transfer simulations to determine whether explosion during the merger of a carbon-oxygen (CO) white dwarf (WD) with a helium (He) WD can reproduce the characteristic Ca II/[Ca II] and He I lines observed in Ca-rich transients. Our study is based on a 1D representation of a hydrodynamic simulation of a 0.6 $M_{\odot}$ CO + 0.4 $M_{\odot}$ He WD merger. We calculate both photospheric and nebular-phase spectra including treatment for non-thermal electrons, as is required for accurate modelling of He I and [Ca II]. Consistent with Ca-rich transients, our simulation predicts a nebular spectrum dominated by emission from [Ca II] 7291, 7324 angstrom and the Ca II near-infrared (NIR) triplet. The photospheric-phase synthetic spectrum also exhibits a strong Ca II NIR triplet, prominent optical absorption due to He I 5876 angstrom and He I 10830 angstrom in the NIR, as is commonly observed for Ca-rich transients. Overall, our results therefore suggest that CO+He WD mergers are a promising channel for Ca-rich transients. However, the current simulation overpredicts some He I features, in particular both He I 6678 and 7065 angstrom and shows a significant contribution from Ti II, which results in a spectral energy distribution that is substantially redder than most Ca-rich transients at peak. Additionally the Ca II nebular emission features are too broad. Future work should investigate if these discrepancies can be resolved by considering full 3D models and exploring a range of CO+He WD binary configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A New Emulated Monte Carlo Radiative Transfer Disk-Wind Model: X-Ray Accretion Disk-wind Emulator -- XRADE

We present a new X-Ray Accretion Disk-wind Emulator (\textsc{xrade}) based on the 2.5D Monte Carlo radiative transfer code which provides a physically-motivated, self-consistent treatment of both absorption and emission from a disk-wind by computing the local ionization state and velocity field within the flow. \textsc{xrade} is then implemented through a process that combines X-ray tracing with supervised machine learning. We develop a novel emulation method consisting in training, validating, and testing the simulated disk-wind spectra into a purposely built artificial neural network. The trained emulator can generate a single synthetic spectrum for a particular parameter set in a fraction of a second, in contrast to the few hours required by a standard Monte Carlo radiative transfer pipeline. The emulator does not suffer from interpolation issues with multi-dimensional spaces that are typically faced by traditional X-ray fitting packages such as \textsc{xspec}. \textsc{xrade} will be suitable to a wide number of sources across the black-hole mass, ionizing luminosity, and accretion rate scales. As an example, we demonstrate the applicability of \textsc{xrade} to the physical interpretation of the X-ray spectra of the bright quasar PDS 456, which hosts the best-established accretion-disk wind observed to date. We anticipate that our emulation method will be an indispensable tool for the development of high-resolution theoretical models, with the necessary flexibility to be optimized for the next generation micro-calorimeters on board future missions, like \textit{XRISM/resolve} and \textit{Athena/X-IFU}. This tool can also be implemented across a wide variety of X-ray spectral models and beyond.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modelling the spectra of the kilonova AT2017gfo -- I: The photospheric epochs

The kilonova (KN) associated with the binary neutron star (BNS) merger GW170817 is the only known electromagnetic counterpart to a gravitational wave source. Here we produce a sequence of radiative transfer models (using $\textsc{tardis}$) with updated atomic data, and compare them to accurately calibrated spectra. We use element compositions from nuclear network calculations based on a realistic hydrodynamical simulation of a BNS merger. We show that the blue spectrum at +1.4 days after merger requires a nucleosynthetic trajectory with a high electron fraction. Our best-fitting model is composed entirely of first $r$-process peak elements (Sr & Zr) and the strong absorption feature is reproduced well by Sr$\,\textsc{ii}$ absorption. At this epoch, we set an upper limit on the lanthanide mass fraction of $X_{\textsc{ln}} \lesssim 5 \times 10^{-3}$. In contrast, all subsequent spectra from $+2.4 - 6.4$ days require the presence of a modest amount of lanthanide material ($X_{\textsc{ln}} \simeq 0.05^{+0.05}_{-0.02}$), produced by a trajectory with $Y_{\rm e} = 0.29$. This produces lanthanide-induced line blanketing below 6000$\,$Å, and sufficient light $r$-process elements to explain the persistent strong feature at $\sim 0.7 - 1.0 \, \micron$ (Sr$\,\textsc{ii}$). The composition gives good matches to the observed data, indicating that the strong blue flux deficit results in the near-infrared (NIR) excess. The disjoint in composition between the first epoch and all others indicates either ejecta stratification, or the presence of two distinct components of material. This further supports the `two-component' kilonova model, and constrains the element composition from nucleosynthetic trajectories. The major uncertainties lie in availability of atomic data and the ionisation state of the expanding material.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Models of pulsationally assisted gravitationally confined detonations with different ignition conditions

Over the past decades, many explosion scenarios for Type Ia supernovae have been proposed and investigated including various combinations of deflagrations and detonations in white dwarfs of different masses up to the Chandrasekhar mass. One of these is the gravitationally confined detonation model. In this case a weak deflagration burns to the surface, wraps around the bound core, and collides at the antipode. A subsequent detonation is then initiated in the collision area. Since the parameter space for this scenario, that is, varying central densities and ignition geometries, has not been studied in detail, we used pure deflagration models of a previous parameter study dedicated to Type Iax supernovae as initial models to investigate the gravitationally confined detonation scenario. We aim to judge whether this channel can account for one of the many subgroups of Type Ia supernovae, or even normal events. To this end, we employed a comprehensive pipeline for three-dimensional Type Ia supernova modeling that consists of hydrodynamic explosion simulations, nuclear network calculations, and radiative transfer. The observables extracted from the radiative transfer are then compared to observed light curves and spectra. The study produces a wide range in masses of synthesized 56 Ni ranging from 0.257 to 1.057 $M_\odot$ , and, thus, can potentially account for subluminous as well as overluminous Type Ia supernovae in terms of brightness. However, a rough agreement with observed light curves and spectra can only be found for 91T-like objects. Although several discrepancies remain, we conclude that the gravitationally confined detonation model cannot be ruled out as a mechanism to produce 91T-like objects. However, the models do not provide a good explanation for either normal Type Ia supernovae or Type Iax supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

SN 2020kyg and the rates of faint Iax Supernovae from ATLAS

We present multi-wavelength follow-up observations of the ATLAS discovered faint Iax supernova SN 2020kyg that peaked at an absolute magnitude of $M_g \approx -14.9 \pm 0.2$, making it another member of the faint Iax supernova population. The bolometric light curve requires only $\approx 7 \times 10^{-3}$ M$_{\odot}$ of radioactive $^{56}$Ni, with an ejected mass of $M_{\rm ej} \sim 0.4$ M$_{\odot}$ and a low kinetic energy of $E \approx 0.05 \pm 0.02 \times 10^{51}$ erg. We construct a homogeneous volume-limited sample of 902 transients observed by ATLAS within 100 Mpc during a 3.5 year span. Using this sample, we constrain the rates of faint Iax ($M_r \gtrsim -16$) events within 60 Mpc at $12^{+14}_{-8}\%$ of the SN Ia rate. The overall Iax rate, at $15^{+17}_{-9}\%$ of the Ia rate, is dominated by the low-luminosity events, with luminous SNe Iax ($M_r \lesssim -17.5$) like 2002cx and 2005hk accounting for only $0.9^{+1.1}_{-0.5}\%$ of the Ia rate (a 2$σ$ upper limit of approximately 3\%). We favour the hybrid CONe WD + He star progenitor channel involving a failed deflagration of a near Chandrasekhar mass white dwarf, expected to leave a bound remnant and a surviving secondary companion, as a candidate explanation for faint Iax explosions. This scenario requires short delay times, consistent with the observed environments of SNe Iax. Furthermore, binary population synthesis calculations have suggested rates of $1-18\%$ of the SN Ia rate for this channel, consistent with our rate estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Exploring the diversity of double detonation explosions for type Ia supernovae: Effects of the post-explosion helium shell composition

The detonation of a helium shell on top of a carbon-oxygen white dwarf has been argued as a potential explosion mechanism for type Ia supernovae (SNe~Ia). The ash produced during helium shell burning can lead to light curves and spectra that are inconsistent with normal SNe~Ia, but may be viable for some objects showing a light curve bump within the days following explosion. We present a series of radiative transfer models designed to mimic predictions from double detonation explosion models. We consider a range of core and shell masses, and systematically explore multiple post-explosion compositions for the helium shell. We find that a variety of luminosities and timescales for early light curve bumps result from those models with shells containing $^{56}$Ni, $^{52}$Fe, or $^{48}$Cr. Comparing our models to SNe~Ia with light curve bumps, we find that these models can reproduce the shapes of almost all of the bumps observed, but only those objects with red colours around maximum light ($B-V \gtrsim 1$) are well matched throughout their evolution. Consistent with previous works, we also show that those models in which the shell does not contain iron-group elements provide good agreement with normal SNe~Ia of different luminosities from shortly after explosion up to maximum light. While our models do not amount to positive evidence in favour of the double detonation scenario, we show that provided the helium shell ash does not contain iron-group elements, it may be viable for a wide range of normal SNe~Ia.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

HOLISMOKES -- III. Achromatic Phase of Strongly Lensed Type Ia Supernovae

To use strongly lensed Type Ia supernovae (LSNe Ia) for cosmology, a time-delay measurement between the multiple supernova (SN) images is necessary. The sharp rise and decline of SN Ia light curves make them promising for measuring time delays, but microlensing can distort these light curves and therefore add large uncertainties to the measurements. An alternative approach is to use color curves where uncertainties due to microlensing are significantly reduced for a certain period of time known as the achromatic phase. In this work, we investigate in detail the achromatic phase, testing four different SN Ia models with various microlensing configurations. We find on average an achromatic phase of around three rest-frame weeks or longer for most color curves but the spread in the duration of the achromatic phase (due to different microlensing maps and filter combinations) is quite large and an achromatic phase of just a few days is also possible. Furthermore, the achromatic phase is longer for smoother microlensing maps, lower macro-magnifications and larger mean Einstein radii of microlenses. From our investigations, we do not find a strong dependency on the model or on asymmetries in the SN ejecta. Further, we find that three independent LSST color curves exhibit features such as extreme points or turning points within the achromatic phase, which make them promising for time-delay measurements. These curves contain combinations of rest-frame bands $u$, $g$, $r$, and $i$ and to observe them for typical LSN Ia redshifts, it would be ideal to cover (observer-frame) filters $r$, $i$, $z$, $y$, $J$, and $H$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

HOLISMOKES -- VII. Time-delay measurement of strongly lensed Type Ia supernovae using machine learning

The Hubble constant ($H_0$) is one of the fundamental parameters in cosmology, but there is a heated debate around the $>$4$σ$ tension between the local Cepheid distance ladder and the early Universe measurements. Strongly lensed Type Ia supernovae (LSNe Ia) are an independent and direct way to measure $H_0$, where a time-delay measurement between the multiple supernova (SN) images is required. In this work, we present two machine learning approaches for measuring time delays in LSNe Ia, namely, a fully connected neural network (FCNN) and a random forest (RF). For the training of the FCNN and the RF, we simulate mock LSNe Ia from theoretical SN Ia models that include observational noise and microlensing. We test the generalizability of the machine learning models by using a final test set based on empirical LSN Ia light curves not used in the training process, and we find that only the RF provides a low enough bias to achieve precision cosmology; as such, RF is therefore preferred over our FCNN approach for applications to real systems. For the RF with single-band photometry in the $i$ band, we obtain an accuracy better than 1\% in all investigated cases for time delays longer than 15 days, assuming follow-up observations with a 5$σ$ point-source depth of 24.7, a two day cadence with a few random gaps, and a detection of the LSNe Ia 8 to 10 days before peak in the observer frame. In terms of precision, we can achieve an approximately 1.5-day uncertainty for a typical source redshift of $\sim$0.8 on the $i$ band under the same assumptions. To improve the measurement, we find that using three bands, where we train a RF for each band separately and combine them afterward, helps to reduce the uncertainty to $\sim$1.0 day. We have publicly released the microlensed spectra and light curves used in this work.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Type Iax supernovae from deflagrations in Chandrasekhar mass white dwarfs

Due to the increasing number of observations Type Ia supernovae are nowadays regarded as a heterogeneous class of objects consisting of several subclasses. One of the largest of these is the class of Type Iax supernovae (SNe Iax) which have been suggested to originate from pure deflagrations in CO Chandrasekhar-mass white dwarfs (WDs). Although a few deflagration studies have been carried out, the full diversity of the class is not captured yet. We therefore present a parameter study of single-spot ignited deflagrations with varying ignition locations, central densities, metallicities and compositions. We also explore a rigidly rotating progenitor and carry out 3D hydrodynamic simulations, nuclear network calculations and radiative transfer. The new models extend the range in brightness covered by previous studies to the lower end. Our explosions produce $^{56}$Ni masses from $5.8 \times 10^{-3}$ to $9.2 \times 10^{-2}\,M_\odot$. In spite of the wide exploration of the parameter space the main characteristics of the models are primarily driven by the mass of $^{56}$Ni. Secondary parameters have too little impact to explain the observed trend among faint SNe~Iax. We report kick velocities of the bound explosion remnants from $6.9$ to $369.8\,$km$\,s^{-1}$. The wide exploration of the parameter space and viewing-angle effects in the radiative transfer lead to a significant spread in the synthetic observables. The trends towards the faint end of the class are, however, not reproduced. This motivates a quantification of the systematic uncertainties in the modeling procedure and the influence of the $^{56}$Ni-rich bound remnant. While the pure deflagration scenario remains a favorable explanation for bright and intermediate luminosity SNe~Iax, the possibility that SNe~Iax do not consist of a single explosion scenario needs to be considered.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HOLISMOKES -- I. Highly Optimised Lensing Investigations of Supernovae, Microlensing Objects, and Kinematics of Ellipticals and Spirals

We present the HOLISMOKES programme on strong gravitational lensing of supernovae as a probe of supernova (SN) physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra using four different SN explosion models, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible ($<$1% distortion within the 1$σ$ spread, and $\lesssim$10% distortion within the 2$σ$ spread). This shows great prospects of using lensed SNe Ia to obtain intrinsic early-phase SN spectra for deciphering SN Ia progenitors. As a demonstration of the usefulness of lensed SNe Ia for cosmology, we simulate a sample of mock lensed SN Ia systems that are expected to have accurate and precise time-delay measurements in the era of the Rubin Observatory Legacy Survey of Space and Time (LSST). Adopting realistic yet conservative uncertainties on their time-delay distances and lens angular diameter distances (of 6.6% and 5%, respectively), we find that a sample of 20 lensed SNe Ia would allow a constraint on the Hubble constant ($H_0$) with 1.3% uncertainty in the flat $Λ$CDM cosmology. We find a similar constraint on $H_0$ in an open $Λ$CDM cosmology, while the constraint degrades to $3\%$ in a flat $w$CDM cosmology. We anticipate lensed SNe to be an independent and powerful probe of SN physics and cosmology in the upcoming LSST era.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

White dwarf deflagrations for Type Iax supernovae: Polarisation signatures from the explosion and companion interaction

Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths ($\lesssim$ 5000 $\unicode{x212B}$). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal $P\sim0.2$ per cent is found at blue wavelengths for orientations 30$^\circ$ and 45$^\circ$ away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Determining the $^{56}$Ni distribution of type Ia supernovae from observations within days of explosion

Recent studies have shown how the distribution of $^{56}$Ni within the ejecta of type Ia supernovae can have profound consequences on the observed light curves. Observations at early times can therefore provide important details on the explosion physics in thermonuclear supernovae. We present a series of radiative transfer calculations that explore variations in the $^{56}$Ni distribution. Our models also show the importance of the density profile in shaping the light curve, which is often neglected in the literature. Using our model set, we investigate the observations that are necessary to determine the $^{56}$Ni distribution as robustly as possible within the current model set. We find that this includes observations beginning at least $\sim$14 days before $B$-band maximum, extending to approximately maximum light with a high ($\lesssim$3 day) cadence, and in at least one blue and one red band are required (such as $B$ and $R$, or $g$ and $r$). We compare a number of well-observed type Ia supernovae that meet these criteria to our models and find that the light curves of $\sim$70-80\% of objects in our sample are consistent with being produced solely by variations in the $^{56}$Ni distributions. The remaining supernovae show an excess of flux at early times, indicating missing physics that is not accounted for within our model set, such as an interaction or the presence of short-lived radioactive isotopes. Comparing our model light curves and spectra to observations and delayed detonation models demonstrates that while a somewhat extended $^{56}$Ni distribution is necessary to reproduce the observed light curve shape, this does not negatively affect the spectra at maximum light. Investigating current explosion models shows that observations typically require a shallower decrease in the $^{56}$Ni mass towards the outer ejecta than is produced for models of a given $^{56}$Ni mass.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Monte Carlo radiative transfer for the nebular phase of Type Ia supernovae

We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Spectral modeling of type II supernovae II. A machine learning approach to quantitative spectroscopic analysis

There are now hundreds of publicly available supernova spectral time series. Radiative transfer modeling of this data gives insights into the physical properties of these explosions such as the composition, the density structure, or the intrinsic luminosity---this is invaluable for understanding the supernova progenitors, the explosion mechanism, or for constraining the supernova distance. However, a detailed parameter study of the available data has been out of reach due to the high dimensionality of the problem coupled with the still significant computational expense. We tackle this issue through the use of machine-learning emulators, which are algorithms for high-dimensional interpolation. These use a pre-calculated training dataset to mimic the output of a complex code but with run times orders of magnitude shorter. We present the application of such an emulator to synthetic type II supernova spectra generated with the TARDIS radiative transfer code. The results show that with a relatively small training set of 780 spectra we can generate emulated spectra with interpolation uncertainties of less than one percent. We demonstrate the utility of this method by automatic spectral fitting of two well-known type IIP supernovae; as an exemplary application, we determine the supernova distances from the spectral fits using the tailored-expanding-photosphere method. We compare our results to previous studies and find good agreement. This suggests that emulation of TARDIS spectra can likely be used to perform automatic and detailed analysis of many transient classes putting the analysis of large data repositories within reach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

LSQ13fn: A type II-Plateau supernova with a possibly low metallicity progenitor that breaks the standardised candle relation

We present optical imaging and spectroscopy of supernova (SN) LSQ13fn, a type II supernova with several hitherto-unseen properties. Although it initially showed strong symmetric spectral emission features attributable to \ion{He}{ii}, \ion{N}{iii}, and \ion{C}{iii}, reminiscent of some interacting SNe, it transitioned into an object that would fall more naturally under a type II-Plateau (IIP) classification. However, its spectral evolution revealed several unusual properties: metal lines appeared later than expected, were weak, and some species were conspicuous by their absence. Furthermore, the line velocities were found to be lower than expected given the plateau brightness, breaking the SNe~IIP standardised candle method for distance estimates. We found that, in combination with a short phase of early-time ejecta-circumstellar material interaction, metal-poor ejecta, and a large progenitor radius could reasonably account for the observed behaviour. Comparisons with synthetic model spectra of SNe~IIP of a given progenitor mass would imply a progenitor star metallicity as low as 0.1\,Z$_{\odot}$. LSQ13fn highlights the diversity of SNe~II and the many competing physical effects that come into play towards the final stages of massive star evolution immediately preceding core-collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the diversity of super-luminous supernovae: ejected mass as the dominant factor

We assemble a sample of 24 hydrogen-poor super-luminous supernovae (SLSNe). Parameterizing the light curve shape through rise and decline timescales shows that the two are highly correlated. Magnetar-powered models can reproduce the correlation, with the diversity in rise and decline rates driven by the diffusion timescale. Circumstellar interaction models can exhibit a similar rise-decline relation, but only for a narrow range of densities, which may be problematic for these models. We find that SLSNe are approximately 3.5 magnitudes brighter and have light curves 3 times broader than SNe Ibc, but that the intrinsic shapes are similar. There are a number of SLSNe with particularly broad light curves, possibly indicating two progenitor channels, but statistical tests do not cleanly separate two populations. The general spectral evolution is also presented. Velocities measured from Fe II are similar for SLSNe and SNe Ibc, suggesting that diffusion time differences are dominated by mass or opacity. Flat velocity evolution in most SLSNe suggests a dense shell of ejecta. If opacities in SLSNe are similar to other SNe Ibc, the average ejected mass is higher by a factor 2-3. Assuming $κ=0.1\,$cm$^2\,$g$^{-1}$, we estimate a mean (median) SLSN ejecta mass of 10$\,$M$_\odot$ (6$\,$M$_\odot$), with a range of 3-30$\,$M$_\odot$. Doubling the assumed opacity brings the masses closer to normal SNe Ibc, but with a high-mass tail. The most probable mechanism for generating SLSNe seems to be the core-collapse of a very massive hydrogen-poor star, forming a millisecond magnetar.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Predicting polarization signatures for double-detonation and delayed-detonation models of Type Ia supernovae

Calculations of synthetic spectropolarimetry are one means to test multi-dimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M$_{\odot}$ white dwarf binary system is too asymmetric to explain the low polarization levels commonly observed in normal Type Ia supernovae. Here, we present polarization simulations for two alternative scenarios: the sub-Chandrasekhar mass double-detonation and the Chandrasekhar mass delayed-detonation model. Specifically, we study a two-dimensional double-detonation model and a three-dimensional delayed-detonation model, and calculate polarization spectra for multiple observer orientations in both cases. We find modest polarization levels ($<$ 1 per cent) for both explosion models. Polarization in the continuum peaks at $\sim$ 0.1$-$0.3 per cent and decreases after maximum light, in excellent agreement with spectropolarimetric data of normal Type Ia supernovae. Higher degrees of polarization are found across individual spectral lines. In particular, the synthetic Si ii λ6355 profiles are polarized at levels that match remarkably well the values observed in normal Type Ia supernovae, while the low degrees of polarization predicted across the O i λ7774 region are consistent with the non-detection of this feature in current data. We conclude that our models can reproduce many of the characteristics of both flux and polarization spectra for well-studied Type Ia supernovae, such as SN 2001el and SN 2012fr. However, the two models considered here cannot account for the unusually high level of polarization observed in extreme cases such as SN 2004dt.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Slowly fading super-luminous supernovae that are not pair-instability explosions

Super-luminous supernovae that radiate more than 10^44 ergs per second at their peak luminosity have recently been discovered in faint galaxies at redshifts of 0.1-4. Some evolve slowly, resembling models of 'pair-instability' supernovae. Such models involve stars with original masses 140-260 times that of the Sun that now have carbon-oxygen cores of 65-30 solar masses. In these stars, the photons that prevent gravitational collapse are converted to electron-positron pairs, causing rapid contraction and thermonuclear explosions. Many solar masses of 56Ni are synthesized; this isotope decays to 56Fe via 56Co, powering bright light curves. Such massive progenitors are expected to have formed from metal-poor gas in the early Universe. Recently, supernova 2007bi in a galaxy at redshift 0.127 (about 12 billion years after the Big Bang) with a metallicity one-third that of the Sun was observed to look like a fading pair-instability supernova. Here we report observations of two slow-to-fade super-luminous supernovae that show relatively fast rise times and blue colours, which are incompatible with pair-instability models. Their late-time light-curve and spectral similarities to supernova 2007bi call the nature of that event into question. Our early spectra closely resemble typical fast-declining super-luminous supernovae, which are not powered by radioactivity. Modelling our observations with 10-16 solar masses of magnetar-energized ejecta demonstrates the possibility of a common explosion mechanism. The lack of unambiguous nearby pair-instability events suggests that their local rate of occurrence is less than 6x10^-6 times that of the core-collapse rate.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spectropolarimetry of superluminous supernovae: insight into their geometry

We present the first spectropolarimetric observations of a hydrogen-free superluminous supernova at z=0.1136, namely SN 2015bn. The transient shows significant polarization at both the observed epochs: one 24 days before maximum light in the rest-frame, and the subsequent at 27 days after peak luminosity. Analysis of the Q-U plane suggests the presence of a dominant axis and no physical departure from the main axis at either epoch. The polarization spectrum along the dominant axis is characterized by a strong wavelength dependence and an increase in the signal from the first to the second epoch. We use a Monte Carlo code to demonstrate that these properties are consistent with a simple toy model that adopts an axi-symmetric ellipsoidal configuration for the ejecta. We find that the wavelength dependence of the polarisation is possibly due to a strong wavelength dependence in the line opacity, while the higher level of polarisation at the second epoch is a consequence of the increase in the asphericity of the inner layers of the ejecta or the fact that the photosphere recedes into less spherical layers. The geometry of the superluminous supernova results similar to those of stripped-envelope core collapse SNe connected with GRB, while the overall evolution of the ejecta shape could be consistent with a central engine.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The peculiar Type Ia supernova iPTF14atg: Chandrasekhar-mass explosion or violent merger?

iPTF14atg, a subluminous peculiar Type Ia supernova (SN Ia) similar to SN 2002es, is the first SN Ia for which a strong UV flash was observed in the early-time light curves. This has been interpreted as evidence for a single-degenerate (SD) progenitor system where such a signal is expected from interactions between the SN ejecta and the non-degenerate companion star. Here, we compare synthetic observables of multi-dimensional state-of-the-art explosion models for different progenitor scenarios to the light curves and spectra of iPTF14atg. From our models, we have difficulties explaining the spectral evolution of iPTF14atg within the SD progenitor channel. In contrast, we find that a violent merger of two carbon-oxygen white dwarfs with 0.9 and 0.76 solar masses, respectively, provides an excellent match to the spectral evolution of iPTF14atg from 10d before to several weeks after maximum light. Our merger model does not naturally explain the initial UV flash of iPTF14atg. We discuss several possibilities like interactions of the SN ejecta with the circum-stellar medium and surface radioactivity from a He ignited merger that may be able to account for the early UV emission in violent merger models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The type Iax supernova, SN 2015H: a white dwarf deflagration candidate

We present results based on observations of SN 2015H which belongs to the small group of objects similar to SN 2002cx, otherwise known as type Iax supernovae. The availability of deep pre-explosion imaging allowed us to place tight constraints on the explosion epoch. Our observational campaign began approximately one day post-explosion, and extended over a period of about 150 days post maximum light, making it one of the best observed objects of this class to date. We find a peak magnitude of M$_r$ = -17.27 $\pm$ 0.07, and a ($Δm_{15})_r$ = 0.69 $\pm$ 0.04. Comparing our observations to synthetic spectra generated from simulations of deflagrations of Chandrasekhar mass carbon-oxygen white dwarfs, we find reasonable agreement with models of weak deflagrations that result in the ejection of ~0.2 M$_{\odot}$ of material containing ~0.07 M$_{\odot}$ of 56Ni. The model light curve however, evolves more rapidly than observations, suggesting that a higher ejecta mass is to be favoured. Nevertheless, empirical modelling of the pseudo-bolometric light curve suggests that $\lesssim$0.6 M_sun of material was ejected, implying that the white dwarf is not completely disrupted, and that a bound remnant is a likely outcome.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraints on explosive silicon burning in core-collapse supernovae from measured Ni/Fe ratios

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54Fe + 56Ni), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova that produced a Ni/Fe ratio of $3.4\pm1.2$ times solar, we find that burning of a fuel with neutron excess $η\approx 6\times 10^{-3}$ is required. Unless the progenitor metallicity is over 5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Supernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of $M_{\rm ZAMS} \lesssim 13$ Msun stars exploding with a delay time of less than one second ($M_{\rm cut} < 1.5$ Msun) are able to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and neutrino-neutronization obtained in electron-capture models remains the only viable explanation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

LSQ14bdq: A Type Ic super-luminous supernova with a double-peaked light curve

We present data for LSQ14bdq, a hydrogen-poor super-luminous supernova (SLSN) discovered by the La Silla QUEST survey and classified by the Public ESO Spectroscopic Survey of Transient Objects. The spectrum and light curve are very similar to slow-declining SLSNe such as PTF12dam. However, detections within ~1 day after explosion show a bright and relatively fast initial peak, lasting for ~15 days, prior to the usual slow rise to maximum light. The broader, main peak can be fit with either central engine or circumstellar interaction models. We discuss the implications of the precursor peak in the context of these models. It is too bright and narrow to be explained as a normal 56Ni-powered SN, and we suggest that interaction models may struggle to fit the two peaks simultaneously. We propose that the initial peak may arise from the post-shock cooling of extended stellar material, and reheating by a central engine drives the second peak. In this picture, we show that an explosion energy of ~2x10^{52} erg and a progenitor radius of a few hundred solar radii would be required to power the early emission. The competing engine models involve rapidly spinning magnetars (neutron stars) or fall-back accretion onto a central black hole. The prompt energy required may favour the black hole scenario. The bright initial peak effectively rules out a compact Wolf-Rayet star as a progenitor, since the inferred energies and ejected masses become unphysical.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Polarisation spectral synthesis for Type Ia supernova explosion models

We present a Monte Carlo radiative transfer technique for calculating synthetic spectropolarimetry for multi-dimensional supernova explosion models. The approach utilises "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 vital for calculating synthetic spectropolarimetry (since the degree of polarisation is typically very small) but also useful for calculations of light curves and spectra. We first validate our approach via application of an idealised 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 polarisation 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 polarisation spectra, we then use ARTIS to calculate synthetic observables for prolate and oblate ellipsoidal models with Type Ia supernova compositions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 multi-line 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 is 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Type Ia supernovae from violent mergers of carbon-oxygen white dwarfs: polarisation signatures

The violent merger of two carbon-oxygen white dwarfs has been proposed as a viable progenitor for some Type Ia supernovae. However, it has been argued that the strong ejecta asymmetries produced by this model might be inconsistent with the low degree of polarisation typically observed in Type Ia supernova explosions. Here, we test this claim by carrying out a spectropolarimetric analysis for the model proposed by Pakmor et al. (2012) for an explosion triggered during the merger of a 1.1 M$_{\odot}$ and 0.9 M$_{\odot}$ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarisation signal varies significantly with viewing angle. We find that polarisation levels for observers in the equatorial plane are modest ($\lesssim$ 1 per cent) and show clear evidence for a dominant axis, as a consequence of the ejecta symmetry about the orbital plane. In contrast, orientations out of the plane are associated with higher degrees of polarisation and departures from a dominant axis. While the particular model studied here gives a good match to highly-polarised events such as SN 2004dt, it has difficulties in reproducing the low polarisation levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarisation features that are not observed in the majority of currently available spectropolarimetric data of Type Ia supernovae. Future studies will map out the parameter space of the merger scenario to investigate if alternative models can provide better agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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 Kalpha 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 3D explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two 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 synthesises 3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kalpha 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. For the delayed-detonation scenario, a line detection is feasible with Chandra up to 3 Mpc for an exposure time of 10^6 s. We find that it should be possible with currently existing X-ray instruments (with exposure times 5x10^5 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).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

OGLE-2013-SN-079: a lonely supernova consistent with a helium shell detonation

We present observational data for a peculiar supernova discovered by the OGLE-IV survey and followed by the Public ESO Spectroscopic Survey for Transient Objects. The inferred redshift of $z=0.07$ implies an absolute magnitude in the rest-frame $I$-band of M$_{I}\sim-17.6$ mag. This places it in the luminosity range between normal Type Ia SNe and novae. Optical and near infrared spectroscopy reveal mostly Ti and Ca lines, and an unusually red color arising from strong depression of flux at rest wavelengths $<5000$ Å. To date, this is the only reported SN showing Ti-dominated spectra. The data are broadly consistent with existing models for the pure detonation of a helium shell around a low-mass CO white dwarf and "double-detonation" models that include a secondary detonation of a CO core following a primary detonation in an overlying helium shell.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The ejected mass distribution of type Ia supernovae: A significant rate of non-Chandrasekhar-mass progenitors

The ejected mass distribution of type Ia supernovae directly probes progenitor evolutionary history and explosion mechanisms, with implications for their use as cosmological probes. Although the Chandrasekhar mass is a natural mass scale for the explosion of white dwarfs as type Ia supernovae, models allowing type Ia supernovae to explode at other masses have attracted much recent attention. Using an empirical relation between the ejected mass and the light curve width, we derive ejected masses $M_\mathrm{ej}$ and $^{56}$Ni masses $M_\mathrm{Ni}$ for a sample of 337 type Ia supernovae with redshifts $z < 0.7$ used in recent cosmological analyses. We use hierarchical Bayesian inference to reconstruct the joint $M_\mathrm{ej}$-$M_\mathrm{Ni}$ distribution, accounting for measurement errors. The inferred marginal distribution of $M_\mathrm{ej}$ has a long tail towards sub-Chandrasekhar masses, but cuts off sharply above 1.4 $M_\odot$. Our results imply that 25\%-50\% of normal type Ia supernovae are inconsistent with Chandrasekhar-mass explosions, with almost all of these being sub-Chandrasekhar-mass; super-Chandrasekhar-mass explosions make up no more than 1\% of all spectroscopically normal type Ia supernovae. We interpret the type Ia supernova width-luminosity relation as an underlying relation between $M_\mathrm{ej}$ and $M_\mathrm{Ni}$, and show that the inferred relation is not naturally explained by the predictions of any single known explosion mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Three-dimensional pure deflagration models with nucleosynthesis and synthetic observables for Type Ia supernovae

We investigate whether pure deflagration models of Chandrasekhar-mass carbon-oxygen white dwarf stars can account for one or more subclass of the observed population of Type Ia supernova (SN Ia) explosions. We compute a set of 3D full-star hydrodynamic explosion models, in which the deflagration strength is parametrized using the multispot ignition approach. For each model, we calculate detailed nucleosynthesis yields in a post-processing step with a 384 nuclide nuclear network. We also compute synthetic observables with our 3D Monte Carlo radiative transfer code for comparison with observations. For weak and intermediate deflagration strengths (energy release E_nuc <~ 1.1 x 10^51 erg), we find that the explosion leaves behind a bound remnant enriched with 3 to 10 per cent (by mass) of deflagration ashes. However, we do not obtain the large kick velocities recently reported in the literature. We find that weak deflagrations with E_nuc ~ 0.5 x 10^51 erg fit well both the light curves and spectra of 2002cx-like SNe Ia, and models with even lower explosion energies could explain some of the fainter members of this subclass. By comparing our synthetic observables with the properties of SNe Ia, we can exclude the brightest, most vigorously ignited models as candidates for any observed class of SN Ia: their B - V colours deviate significantly from both normal and 2002cx-like SNe Ia and they are too bright to be candidates for other subclasses.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Type Ia supernova bolometric light curves and ejected mass estimates from the Nearby Supernova Factory

We present a sample of normal type Ia supernovae from the Nearby Supernova Factory dataset with spectrophotometry at sufficiently late phases to estimate the ejected mass using the bolometric light curve. We measure $^{56}$Ni masses from the peak bolometric luminosity, then compare the luminosity in the $^{56}$Co-decay tail to the expected rate of radioactive energy re- lease from ejecta of a given mass. We infer the ejected mass in a Bayesian context using a semi-analytic model of the ejecta, incorporating constraints from contemporary numerical models as priors on the density structure and distribution of $^{56}$Ni throughout the ejecta. We find a strong correlation between ejected mass and light curve decline rate, and consequently $^{56}$Ni mass, with ejected masses in our data ranging from 0.9-1.4 $M_\odot$. Most fast-declining (SALT2 $x_1 < -1$) normal SNe Ia have significantly sub-Chandrasekhar ejected masses in our fiducial analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Simple Disk Wind Model for Broad Absorption Line Quasars

Approximately 20% of quasi-stellar objects (QSOs) exhibit broad, blue-shifted absorption lines in their ultraviolet spectra. Such features provide clear evidence for significant outflows from these systems, most likely in the form of accretion disk winds. These winds may represent the "quasar" mode of feedback that is often invoked in galaxy formation/evolution models, and they are also key to unification scenarios for active galactic nuclei (AGN) and QSOs. To test these ideas, we construct a simple benchmark model of an equatorial, biconical accretion disk wind in a QSO and use a Monte Carlo ionization/radiative transfer code to calculate the ultraviolet spectra as a function of viewing angle. We find that for plausible outflow parameters, sightlines looking directly into the wind cone do produce broad, blue-shifted absorption features in the transitions typically seen in broad absorption line QSOs. However, our benchmark model is intrinsically X-ray weak in order to prevent overionization of the outflow, and the wind does not yet produce collisionally excited line emission at the level observed in non-BAL QSOs. As a first step towards addressing these shortcomings, we discuss the sensitivity of our results to changes in the assumed X-ray luminosity and mass-loss rate, Mdot(wind). In the context of our adopted geometry, Mdot(wind) ~ Mdot(acc) is required in order to produce significant BAL features. The kinetic luminosity and momentum carried by such outflows would be sufficient to provide significant feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

Besides the fact that the gamma-ray emission due to radioactive decays is responsible for powering the light curves of Type Ia supernovae (SNe Ia), gamma rays themselves are of particular interest as a diagnostic tool because they provide a direct way to obtain deeper insights into the nucleosynthesis and the kinematics of these explosion events. 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 use three-dimensional explosion models and perform 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. In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be rather different. The almost direct connection of the emission of gamma rays to fundamental physical processes occuring 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 and 16 Mpc for an exposure time of 10^6 s, respectively. The possibility to detect 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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, LX 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 suppresses the line force, and we suggest that this may render radiative braking less effective. Comparisons of model results against observations reveals reasonable agreement in terms of log(LX/Lbol). 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).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

SN 2010lp---a Type Ia supernova from a violent merger of two carbon--oxygen White Dwarfs

SN 2010lp is a subluminous Type Ia supernova (SN Ia) with slowly-evolving lightcurves. Moreover, it is the only subluminous SN Ia observed so far that shows narrow emission lines of OI in late-time spectra, indicating unburned oxygen close to the centre of the ejecta. Most explosion models for SNe Ia cannot explain the narrow OI emission. Here, we present hydrodynamic explosion and radiative transfer calculations showing that the violent merger of two carbon-oxygen white dwarfs of 0.9 and 0.76 solar masses, respectively, adequately reproduces the early-time observables of SN 2010lp. Moreover, our model predicts oxygen close to the centre of the explosion ejecta, a pre-requisite for narrow OI emission in nebular spectra as observed in SN 2010lp.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spectroscopic Observations of SN 2012fr: A Luminous Normal Type Ia Supernova with Early High Velocity Features and Late Velocity Plateau

We present 65 optical spectra of the Type Ia supernova SN 2012fr, of which 33 were obtained before maximum light. At early times SN 2012fr shows clear evidence of a high-velocity feature (HVF) in the Si II 6355 line which can be cleanly decoupled from the lower velocity "photospheric" component. This Si II 6355 HVF fades by phase -5; subsequently, the photospheric component exhibits a very narrow velocity width and remains at a nearly constant velocity of v~12,000 km/s until at least 5 weeks after maximum brightness. The Ca II infrared (IR) triplet exhibits similar evidence for both a photospheric component at v~12,000 km/s with narrow line width and long velocity plateau, as well as a high-velocity component beginning at v~31,000 km/s two weeks before maximum. SN 2012fr resides on the border between the "shallow silicon" and "core-normal" subclasses in the Branch et al. (2009) classification scheme, and on the border between normal and "high-velocity" SNe Ia in the Wang et al. (2009a) system. Though it is a clear member of the "low velocity gradient" (LVG; Benetii et al., 2005) group of SNe Ia and exhibits a very slow light-curve decline, it shows key dissimilarities with the overluminous SN 1991T or SN 1999aa subclasses of SNe Ia. SN 2012fr represents a well-observed SN Ia at the luminous end of the normal SN Ia distribution, and a key transitional event between nominal spectroscopic subclasses of SNe Ia.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Synthetic light curves and spectra for three-dimensional delayed-detonation models of Type Ia supernovae

In a companion paper, Seitenzahl et al. (2013) presented a set of three-dimensional delayed detonation models for thermonuclear explosions of near-Chandrasekhar mass white dwarfs (WDs). Here, we present multi-dimensional radiative transfer simulations that provide synthetic light curves and spectra for those models. The model sequence explores both changes in the strength of the deflagration phase (controlled by the ignition configuration) and the WD central density. In agreement with previous studies, we find that the strength of the deflagration significantly affects the explosion and the observables. Variations in the central density also have an influence on both brightness and colour, but overall it is a secondary parameter in our set of models. In many respects, the models yield a good match to normal Type Ia supernovae (SNe Ia): peak brightness, rise/decline time scales and synthetic spectra are all in reasonable agreement. There are, however, several differences. In particular, the models are too red around maximum light, manifest spectral line velocities that are a little too high and yield I-band light curves that do not match observations. Although some of these discrepancies may simply relate to approximations made in the modelling, some pose real challenges to the models. If viewed as a complete sequence, our models do not reproduce the observed light-curve width-luminosity relation (WLR) of SNe Ia: all our models show similar B-band decline rates, irrespective of peak brightness. This suggests that simple variations in the strength of the deflagration phase in Chandrasekhar-mass deflagration-to-detonation models do not readily explain the observed diversity of normal SNe Ia. This may imply that some other parameter within the Chandrasekhar-mass paradigm is key to the WLR, or that a substantial fraction of normal SNe Ia arise from an alternative explosion scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variability of the High Velocity Outflow in the Quasar PDS 456

We present a comparison of two Suzaku X-ray observations of the nearby (z=0.184), luminous ($L_{bol} \sim 10^{47}$ erg s$^{-1}$) type I quasar, PDS456. A new 125ks Suzaku observation in 2011 caught the quasar during a period of low X-ray flux and with a hard X-ray spectrum, in contrast to a previous 190ks Suzaku observation in 2007 when the quasar appeared brighter and had a steep ($Γ>2$) X-ray spectrum. The 2011 X-ray spectrum contains a pronounced trough near 9\,keV in the quasar rest frame, which can be modeled with blue-shifted iron K-shell absorption, most likely from the He and H-like transitions of iron. The absorption trough is observed at a similar rest-frame energy as in the earlier 2007 observation, which appears to confirm the existence of a persistent high velocity wind in PDS 456, at an outflow velocity of $0.25-0.30$c. The spectral variability between 2007 and 2011 can be accounted for by variations in a partial covering absorber, increasing in covering fraction from the brighter 2007 observation to the hard and faint 2011 observation. Overall the low flux 2011 observation can be explained if PDS 456 is observed at relatively low inclination angles through a Compton thick wind, originating from the accretion disk, which significantly attenuates the X-ray flux from the quasar.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

3D deflagration simulations leaving bound remnants: a model for 2002cx-like Type Ia supernovae

2002cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae. Their light curves and spectra are characterized by distinct features that indicate strong mixing of the explosion ejecta. Pure turbulent deflagrations have been shown to produce such mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative transfer calculations for a 3D full-star deflagration of a Chandrasekhar-mass white dwarf. Our model is able to reproduce the characteristic observational features of SN 2005hk (a proto-typical 2002cx-like supernova), not only in the optical, but also in the near-infrared. For that purpose we present, for the first time, five near-infrared spectra of SN 2005hk from -0.2 to 26.6 days with respect to B-band maximum. Since our model burns only small parts of the initial white dwarf, it fails to completely unbind the white dwarf and leaves behind a bound remnant of ~1.03 solar masses -- consisting mainly of unburned carbon and oxygen, but also enriched by some amount of intermediate-mass and iron-group elements from the explosion products that fall back on the remnant. We discuss possibilities for detecting this bound remnant and how it might influence the late-time observables of 2002cx-like SNe.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining Type Ia supernova models: SN 2011fe as a test case

The nearby supernova SN 2011fe can be observed in unprecedented detail. Therefore, it is an important test case for Type Ia supernova (SN Ia) models, which may bring us closer to understanding the physical nature of these objects. Here, we explore how available and expected future observations of SN 2011fe can be used to constrain SN Ia explosion scenarios. We base our discussion on three-dimensional simulations of a delayed detonation in a Chandrasekhar-mass white dwarf and of a violent merger of two white dwarfs-realizations of explosion models appropriate for two of the most widely-discussed progenitor channels that may give rise to SNe Ia. Although both models have their shortcomings in reproducing details of the early and near-maximum spectra of SN 2011fe obtained by the Nearby Supernova Factory (SNfactory), the overall match with the observations is reasonable. The level of agreement is slightly better for the merger, in particular around maximum, but a clear preference for one model over the other is still not justified. Observations at late epochs, however, hold promise for discriminating the explosion scenarios in a straightforward way, as a nucleosynthesis effect leads to differences in the 55Co production. SN 2011fe is close enough to be followed sufficiently long to study this effect.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Interpreting the near-infrared spectra of the 'golden standard' Type Ia supernova 2005cf

We present nine near-infrared (NIR) spectra of supernova (SN) 2005cf at epochs from -10 d to +42 d with respect to B-band maximum, complementing the existing excellent data sets available for this prototypical Type Ia SN at other wavelengths. The spectra show a time evolution and spectral features characteristic of normal Type Ia SNe, as illustrated by a comparison with SNe 1999ee, 2002bo and 2003du. The broad-band spectral energy distribution (SED) of SN 2005cf is studied in combined ultraviolet (UV), optical and NIR spectra at five epochs between ~ 8 d before and ~ 10 d after maximum light. We also present synthetic spectra of the hydrodynamic explosion model W7, which reproduce the key properties of SN 2005cf not only at UV-optical as previously reported, but also at NIR wavelengths. From the radiative-transfer calculations we infer that fluorescence is the driving mechanism that shapes the SED of SNe Ia. In particular, the NIR part of the spectrum is almost devoid of absorption features, and instead dominated by fluorescent emission of both iron-group material and intermediate-mass elements at pre-maximum epochs, and pure iron-group material after maximum light. A single P-Cygni feature of Mg II at early epochs and a series of relatively unblended Co II lines at late phases allow us to constrain the regions of the ejecta in which the respective elements are abundant.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling the Fe K Line Profiles in Type I AGN with a Compton-Thick Disk Wind

We have modeled a small sample of Seyfert galaxies that were previously identified as having simple X-ray spectra with little intrinsic absorption. The sources in this sample all contain moderately broad components of Fe K-shell emission and are ideal candidates for testing the applicability of a Compton-thick accretion-disk wind model to AGN emission components. Viewing angles through the wind allow the observer to see the absorption signature of the gas, whereas face-on viewing angles allow the observer to see the scattered light from the wind. We find that the Fe K emission line profiles are well described with a model of a Compton-thick accretion-disk wind of solar abundances, arising tens to hundred of gravitational radii from the central black hole. Further, the fits require a neutral component of Fe K alpha emission that is too narrow to arise from the inner part of the wind, and likely comes from a more distant reprocessing region. Our study demonstrates that a Compton-thick wind can have a profound effect on the observed X-ray spectrum of an AGN, even when the system is not viewed through the flow.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Monte Carlo radiation hydrodynamics: methods, tests and application to supernova Type Ia 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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 \mathrm{M_\odot}$ and $1.1 \mathrm{M_\odot}$ 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 lightcurves of our merger, which produces about $0.62 \mathrm{M_\odot}$ of $^{56}\mathrm{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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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. 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 SNe Ia. We argue that the brightness of an explosion realized via the violent merger model is determined by the mass of nickel 56 produced in the detonation of the primary WD. We use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass and expected peak bolometric brightness. We use this relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We find a striking similarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a brightness 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.35 Msun from a slightly-evolved helium star companion. We also find that with high probability, violent WD mergers involving the most massive primaries (which should produce bright SNe) have delay times <500 Myr (Abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. Some of these outflows may be driven off gas infalling towards the central black hole. We perform radiative transfer calculations to compute the gas ionization state and X-ray spectra for two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of outflowing gas. This outflow causes strong absorption for observer orientation angles of < 35 degrees. Particularly in 3D, the strength of this absorption varies significantly for different lines-of-sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features since the ionization state is very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow is predicted to give rise to several emission features for all observer orientations. For the hydrodynamical simulations considered we find both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. (Abridged.)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

2D simulations of the double-detonation model for thermonuclear transients from low-mass carbon-oxygen white dwarfs

Thermonuclear explosions may arise in binaries in which a CO white dwarf (WD) accretes He from a companion. If the accretion rate allows a sufficiently large mass of He to accumulate prior to ignition of nuclear burning, the He surface layer may detonate, giving rise to an astrophysical transient. Detonation of the He layer generates shock waves that propagate into the underlying CO WD. This might directly ignite a detonation at the edge of the CO WD or compress the core of the WD sufficiently to trigger a CO detonation near the centre. If either ignition mechanism works, the two detonations can release sufficient energy to completely unbind the WD. Here we extend our 2D studies of this double-detonation model to low-mass CO WDs. We investigate the feasibility of triggering a secondary core detonation by shock convergence in low-mass CO WDs and the observable consequences of such a detonation. Our results suggest that core detonation is probable, even for the lowest CO core masses realized in nature. We compute spectra and light curves for models in which either an edge-lit or compression-triggered CO detonation is assumed to occur and compare these to models in which no CO detonation was allowed to occur. If significant shock compression of the CO WD occurs prior to detonation, explosion of the CO WD can produce a sufficiently large mass of radioactive iron-group nuclei to affect the light curves. In particular, this can lead to relatively slow post-maximum decline. If the secondary detonation is edge-lit, however, the CO WD explosion primarily yields intermediate-mass elements that affect the observables more subtly. In this case, NIR observations and detailed spectroscopic analysis would be needed to determine whether core detonation occurred. We comment on the implications of our results for understanding peculiar astrophysical transients including SN 2002bj, SN 2010X and SN 2005E.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analytical approximations to numerical solutions of theoretical emission measure distributions

Emission line fluxes from cool stars are widely used to establish an apparent emission measure distribution, Emd-app(Te), between temperatures characteristic of the low transition region and the low corona. The true emission measure distribution, Emd-t(Te), is determined by the energy balance and geometry adopted and, with a numerical model, can be used to predict Emd-app(Te), to guide further modelling. The scaling laws that exist between coronal parameters arise from the dimensions of the terms in the energy balance equation. Here, analytical approximations to numerical solutions for Emd-t(Te) are presented, which show how the constants in the coronal scaling laws are determined. The apparent emission measure distributions show a minimum value at some temperature (T0) and a maximum at the mean coronal temperature Tc (although in some stars, emission from active regions can contribute). It is shown that, for the energy balance and geometry adopted, the analytical values of the emission measure and electron pressure at T0 and Tc, depend on only three parameters: the stellar surface gravity and the values of T0 and Tc. The results are tested against full numerical solutions for Epsilon Eri (K2 V) and are applied to Procyon (alpha CMi; F5 IV/V). The analytical approximations can be used to restrict the required range of full numerical solutions, to check the assumed geometry and to show where the adopted energy balance may not be appropriate.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Delay Times and Rates for Type Ia Supernovae and Thermonuclear Explosions from Double-detonation Sub-Chandrasekhar Mass Models

We present theoretical delay times and rates of thermonuclear explosions that are thought to produce Type Ia supernovae, including the double-detonation sub-Chandrasekhar mass model, using the population synthesis binary evolution code StarTrack. If detonations of sub-Chandrasekhar mass carbon-oxygen white dwarfs following a detonation in an accumulated layer of helium on the white dwarf's surface ("double-detonation" models) are able to produce thermonuclear explosions which are characteristically similar to those of SNe Ia, then these sub-Chandrasekhar mass explosions may account for at least some substantial fraction of the observed SN Ia rate. Regardless of whether all double-detonations look like 'normal' SNe Ia, in any case the explosions are expected to be bright and thus potentially detectable. Additionally, we find that the delay time distribution of double-detonation sub-Chandrasekhar mass SNe Ia can be divided into two distinct formation channels: the 'prompt' helium-star channel with delay times <500 Myr (~10% of all sub-Chandras), and the 'delayed' double white dwarf channel, with delay times >800 Myr spanning up to a Hubble time (~90%). These findings coincide with recent observationally-derived delay time distributions which have revealed that a large number of SNe Ia are prompt with delay times <500 Myr, while a significant fraction also have delay times spanning ~1 Gyr to a Hubble time.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Detonations in Sub-Chandrasekhar Mass C+O White Dwarfs

Explosions of sub-Chandrasekhar-mass white dwarfs are one alternative to the standard Chandrasekhar-mass model of Type Ia supernovae. They are interesting since binary systems with sub-Chandrasekhar-mass primary white dwarfs should be common and this scenario would suggest a simple physical parameter which determines the explosion brightness, namely the mass of the exploding white dwarf. Here we perform one-dimensional hydrodynamical simulations, associated post-processing nucleosynthesis and multi-wavelength radiation transport calculations for pure detonations of carbon-oxygen white dwarfs. The light curves and spectra we obtain from these simulations are in good agreement with observed properties of Type Ia supernovae. In particular, for white dwarf masses from 0.97 - 1.15 Msun we obtain 56Ni masses between 0.3 and 0.8 Msun, sufficient to capture almost the complete range of Type Ia supernova brightnesses. Our optical light curve rise times, peak colours and decline timescales display trends which are generally consistent with observed characteristics although the range of B-band decline timescales displayed by our current set of models is somewhat too narrow. In agreement with observations, the maximum light spectra of the models show clear features associated with intermediate mass elements and reproduce the sense of the observed correlation between explosion luminosity and the ratio of the Si II lines at 6355 and 5972 Angstroms. We therefore suggest that sub-Chandrasekhar mass explosions are a viable model for Type Ia supernovae for any binary evolution scenario leading to explosions in which the optical display is dominated by the material produced in a detonation of the primary white dwarf.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Double-detonation sub-Chandrasekhar supernovae: can minimum helium shell masses detonate the core?

The explosion of sub-Chandrasekhar mass white dwarfs via the double detonation scenario is a potential explanation for type Ia supernovae. In this scenario, a surface detonation in a helium layer initiates a detonation in the underlying carbon/oxygen core leading to an explosion. For a given core mass, a lower bound has been determined on the mass of the helium shell required for dynamical burning during a helium flash, which is a necessary prerequisite for detonation. For a range of core and corresponding minimum helium shell masses, we investigate whether an assumed surface helium detonation is capable of triggering a subsequent detonation in the core even for this limiting case. We carried out hydrodynamic simulations on a co-expanding Eulerian grid in two dimensions assuming rotational symmetry. The detonations are propagated using the level-set approach and a simplified scheme for nuclear reactions that has been calibrated with a large nuclear network. The same network is used to determine detailed nucleosynthetic abundances in a post-processing step. Based on approximate detonation initiation criteria in the literature, we find that secondary core detonations are triggered for all of the simulated models, ranging in core mass from 0.810 up to 1.385 M_solar with corresponding shell masses from 0.126 down to 0.0035 M_solar. This implies that, as soon as a detonation triggers in a helium shell covering a carbon/oxygen white dwarf, a subsequent core detonation is virtually inevitable.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Double-detonation sub-Chandrasekhar supernovae: synthetic observables for minimum helium shell mass models

Abridged. In the double detonation scenario for Type Ia supernovae (SNe Ia) a detonation initiates in a shell of He-rich material accreted from a companion star by a sub-Chandrasekhar-mass White Dwarf (WD). This shell detonation drives a shock front into the carbon-oxygen (C/O) WD that triggers a secondary detonation in the core. The core detonation results in a complete disruption of the WD. Earlier studies concluded that this scenario has difficulties in accounting for the observed properties of SNe Ia since the explosion ejecta are surrounded by the products of explosive He burning in the shell. Recently, it was proposed that detonations might be possible for much less massive He shells than previously assumed. Moreover, it was shown that even detonations of these minimum He shell masses robustly trigger detonations of the C/O core. Here we present time-dependent multi-wavelength radiative transfer calculations for models with minimum He shell mass and derive synthetic observables for both the optical and γ-ray spectral regions. These differ strongly from those found in earlier simulations of sub-Chandrasekhar-mass explosions in which more massive He shells were considered. Our models predict light curves which cover both the range of brightnesses and the rise and decline times of observed SNe Ia. However, their colours and spectra do not match the observations. In particular, their B-V colours are generally too red. We show that this discrepancy is mainly due to the composition of the burning products of the He shell of our models which contain significant amounts of Ti and Cr. Using a toy model, we also show that the burning products of the He shell depend crucially on its initial composition. This leads us to conclude that good agreement between sub-Chandrasekhar-mass explosions and observed SNe Ia may still be feasible but further study of the shell properties is required.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-dimensional modelling of X-ray spectra for AGN accretion-disk outflows II

Highly-ionized fast accretion-disk winds have been suggested as an explanation for a variety of observed absorption and emission features in the X-ray spectra of Active Galactic Nuclei. Simple estimates have suggested that these flows may be massive enough to carry away a significant fraction of the accretion energy and could be involved in creating the link between supermassive black holes and their host galaxies. However, testing these hypotheses, and quantifying the outflow signatures, requires high-quality theoretical spectra for comparison with observations. Here we describe extensions of our Monte Carlo radiative transfer code that allow us to generate realistic theoretical spectra for a much wider variety of disk wind models than possible in our previous work. In particular, we have expanded the range of atomic physics simulated by the code so that L- and M-shell ions can now be included. We have also substantially improved our treatment of both ionization and radiative heating such that we are now able to compute spectra for outflows containing far more diverse plasma conditions. We present example calculations that illustrate the variety of spectral features predicted by parametrized outflow models and demonstrate their applicability to the interpretation of data by comparison with observations of the bright quasar PG1211+143. We find that the major features in the observed 2 - 10 keV spectrum of this object can be well-reproduced by our spectra, confirming that it likely hosts a massive outflow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-dimensional modelling of X-ray spectra for AGN accretion-disk outflows III: application to a hydrodynamical simulation

We perform multi-dimensional radiative transfer simulations to compute spectra for a hydrodynamical simulation of a line-driven accretion disk wind from an active galactic nucleus. The synthetic spectra confirm expectations from parameterized models that a disk wind can imprint a wide variety of spectroscopic signatures including narrow absorption lines, broad emission lines and a Compton hump. The formation of these features is complex with contributions originating from many of the different structures present in the hydrodynamical simulation. In particular, spectral features are shaped both by gas in a successfully launched outflow and in complex flows where material is lifted out of the disk plane but ultimately falls back. We also confirm that the strong Fe Kalpha line can develop a weak, red-skewed line wing as a result of Compton scattering in the outflow. In addition, we demonstrate that X-ray radiation scattered and reprocessed in the flow has a pivotal part in both the spectrum formation and determining the ionization conditions in the wind. We find that scattered radiation is rather effective in ionizing gas which is shielded from direct irradiation from the central source. This effect likely makes the successful launching of a massive disk wind somewhat more challenging and should be considered in future wind simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Type Ia Supernovae and Accretion Induced Collapse

Using the population synthesis binary evolution code StarTrack, we present theoretical rates and delay times of Type Ia supernovae arising from various formation channels. These channels include binaries in which the exploding white dwarf reaches the Chandrasekhar mass limit (DDS, SDS, and helium-rich donor scenario) as well as the sub-Chandrasekhar mass scenario, in which a white dwarf accretes from a helium-rich companion and explodes as a SN Ia before reaching the Chandrasekhar mass limit. We find that using a common envelope parameterization employing energy balance with alpha=1 and lambda=1, the supernova rates per unit mass (born in stars) of sub-Chandrasekhar mass SNe Ia exceed those of all other progenitor channels at epochs t=0.7 - 4 Gyr for a burst of star formation at t=0. Additionally, the delay time distribution of the sub-Chandrasekhar model can be divided in to two distinct evolutionary channels: the `prompt' helium-star channel with delay times < 500 Myr, and the `delayed' double white dwarf channel with delay times > 800 Myr spanning up to a Hubble time. These findings are in agreement with recent observationally-derived delay time distributions which predict that a large number of SNe Ia have delay times < 1 Gyr, with a significant fraction having delay times < 500 Myr. We find that the DDS channel is also able to account for the observed rates of SNe Ia. However, detailed simulations of white dwarf mergers have shown that most of these mergers will not lead to SNe Ia but rather to the formation of a neutron star via accretion-induced collapse. If this is true, our standard population synthesis model predicts that the only progenitor channel which can account for the rates of SNe Ia is the sub-Chandrasekhar mass scenario, and none of the other progenitors considered can fully account for the observed rates.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

3D non-LTE spectrum synthesis for Type Ia supernovae

Despite the importance of Type Ia supernovae as standard candles for cosmology and to the chemical evolution of the Universe, we still have no consistent picture of the nature of these events. Much progress has been made in the hydrodynamical explosion modelling of supernovae Ia in the last few years and fully 3-D explosion models are now available. However those simulations are not directly comparable to observations: to constrain explosion models, radiative transfer calculations must be carried out. We present a new 3-D Monte Carlo radiative transfer code which allows forward modelling of the spectral evolution of Type Ia supernovae from first principles, using hydrodynamical explosion models as input. Here, as a first application, we calculate line-of-sight dependent colour light curves for a toy model of an off-centre explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Late-time supernova light curves: The effect of internal conversion and Auger electrons

Energy release from radioactive decays contributes significantly to supernova light curves. Previous works, which considered the energy deposited by gamma-rays and positrons produced by 56Ni, 56Co, 57Ni, 57Co, 44Ti and 44Sc, have been quite successful in explaining the light curves of both core collapse and thermonuclear supernovae. We point out that Auger and internal conversion electrons together with the associated X-ray cascade, constitute an additional heat source. When a supernova is transparent to gamma-rays, these electrons can contribute significantly to light curves for reasonable nucleosynthetic yields. In particular, the electrons emitted in the decay of 57Co, which are largely due to internal conversion from a fortuitously low-lying 3/2- state in the daughter 57Fe, constitute an additional significant energy deposition channel. We show that when the heating by these electrons is accounted for, a slow-down in the lightcurve of SN 1998bw is naturally obtained for typical hypernova nucleosynthetic yields. Additionally, we show that for generic Type Ia supernova yields, the Auger electrons emitted in the ground-state to ground-state electron capture decay of 55Fe exceed the energy released by the 44Ti decay chain for many years after the explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nebular emission-line profiles of Type Ib/c Supernovae - probing the ejecta asphericity

In order to assess qualitatively the ejecta geometry of stripped-envelope core-collapse supernovae, we investigate 98 late-time spectra of 39 objects, many of them previously unpublished. We perform a Gauss-fitting of the [O I] 6300, 6364 feature in all spectra, with the position, full width at half maximum (FWHM) and intensity of the 6300 Gaussian as free parameters, and the 6364 Gaussian added appropriately to account for the doublet nature of the [O I] feature. On the basis of the best-fit parameters, the objects are organised into morphological classes, and we conclude that at least half of all Type Ib/c supernovae must be aspherical. Bipolar jet-models do not seem to be universally applicable, as we find too few symmetric double-peaked [O I] profiles. In some objects the [O I] line exhibits a variety of shifted secondary peaks or shoulders, interpreted as blobs of matter ejected at high velocity and possibly accompanied by neutron-star kicks to assure momentum conservation. At phases earlier than ~200d, a systematic blueshift of the [O I] 6300, 6364 line centroids can be discerned. Residual opacity provides the most convincing explanation of this phenomenon, photons emitted on the rear side of the SN being scattered or absorbed on their way through the ejecta. Once modified to account for the doublet nature of the oxygen feature, the profile of Mg I] 4571 at sufficiently late phases generally resembles that of [O I] 6300, 6364, suggesting negligible contamination from other lines and confirming that O and Mg are similarly distributed within the ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Time-dependent 3D spectrum synthesis for type Ia supernovae

A Monte Carlo code (ARTIS) for modelling time-dependent three-dimensional spectral synthesis in chemically inhomogeneous models of Type Ia supernova ejecta is presented. Following the propagation of gamma-ray photons, emitted by the radioactive decay of the nucleosynthesis products, energy is deposited in the supernova ejecta and the radiative transfer problem is solved self-consistently, enforcing the constraint of energy conservation in the co-moving frame. Assuming a photoionisation dominated plasma, the equations of ionisation equilibrium are solved together with the thermal balance equation adopting an approximate treatment of excitation. Since we implement a fully general treatment of line formation, there are no free parameters to adjust. Thus a direct comparison between synthetic spectra and light curves, calculated from hydrodynamic explosion models, and observations is feasible. The code is applied to the well known W7 explosion model and the results tested against other studies. Finally the effect of asymmetric ejecta on broad band light curves and spectra is illustrated using an elliptical toy model.

0 citations0 reviewsNo code linkedOpen access