Two-dimensional Monte Carlo simulations of H i line formation in massive young stellar object disc winds

Massive young stellar objects (YSOs) are powerful infrared Hi line emitters. It has been suggested that these lines form in an outflow from a disc surrounding the YSO. Here, new two-dimensional Monte Carlo radiative transfer calculations are described which test this hypothesis. Infrared spectra are synthesized for a YSO disc wind model based on earlier hydrodynamical calculations. The model spectra are in qualitative agreement with the observed spectra from massive YSOs, and therefore provide support for a disc wind explanation for the Hi lines. However, there are some significant differences: the models tend to overpredict the Bra/Br? ratio of equivalent widths and produce line profiles which are slightly too broad and, in contrast to typical observations, are double-peaked. The interpretation of these differences within the context of the disc wind picture and suggestions for their resolution via modifications to the assumed disc and outflow structure are discussed. © 2005 RAS.

Evidence for relativistic features in the X-ray spectrum of Mrk 335

We present an analysis of hard X-ray features in the spectrum of the bright Sy 1 galaxy Mrk 335 observed by the XMM-Newton satellite. Our analysis confirms the presence of a broad, ionized Fe Ka emission line in the spectrum, first found by Gondoin et al. The broad line can be modelled successfully by relativistic accretion disc reflection models. This interpretation is unusually robust in the case of Mrk 335 because of the lack of any ionized ('warm') absorber and the absence a clear narrow core to the line. Partial covering by neutral gas cannot, however, be ruled out statistically as the origin of the broad residuals. Regardless of the underlying continuum we report, for the first time in this source, the detection of a narrow absorption feature at the rest frame energy of ~5.9 keV. If the feature is identified with a resonance absorption line of iron in a highly ionized medium, the redshift of the line corresponds to an inflow velocity of ~0.11-0.15c. We present a simple model for the inflow, accounting approximately for relativistic and radiation pressure effects, and use Monte Carlo methods to compute synthetic spectra for qualitative comparison with the data. This modelling shows that the absorption feature can plausibly be reproduced by infalling gas providing that the feature is identified with Fe xxvi. We require the inflowing gas to extend over a limited range of radii at a few tens of r to match the observed feature. The mass accretion rate in the flow corresponds to 60 per cent of the Eddington limit, in remarkable agreement with the observed rate. The narrowness of the absorption line tends to argue against a purely gravitational origin for the redshift of the line, but given the current data quality we stress that such an interpretation cannot be ruled out. © 2006 The Authors.

Light curves for off-centre ignition models of Type Ia supernovae

Motivated by recent models involving off-centre ignition of Type Ia supernova explosions, we undertake three-dimensional time-dependent radiation transport simulations to investigate the range of bolometric light-curve properties that could be observed from supernovae in which there is a lop-sided distribution of the products from nuclear burning. We consider both a grid of artificial toy models which illustrate the conceivable range of effects and a recent three-dimensional hydrodynamical explosion model. We find that observationally significant viewing angle effects are likely to arise in such supernovae and that these may have important ramifications for the interpretation of the observed diversity of Type Ia supernova and the systematic uncertainties which relate to their use as standard candles in contemporary cosmology. © 2007 RAS.

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

Thermonuclear explosions may arise in binary star systems in which a carbon-oxygen (CO) white dwarf (WD) accretes helium-rich material from a companion star. If the accretion rate allows a sufficiently large mass of helium to accumulate prior to ignition of nuclear burning, the helium surface layer may detonate, giving rise to an astrophysical transient. Detonation of the accreted helium layer generates shock waves that propagate into the underlying CO WD. This might directly ignite a detonation of the CO WD at its surface (an edge-lit secondary detonation) or compress the core of the WD sufficiently to trigger a CO detonation near the centre. If either of these ignition mechanisms works, the two detonations (helium and CO) can then release sufficient energy to completely unbind the WD. These 'double-detonation' scenarios for thermonuclear explosion of WDs have previously been investigated as a potential channel for the production of Type Ia supernovae from WDs of ~ 1 M . Here we extend our 2D studies of the double-detonation model to significantly less massive CO WDs, the explosion of which could produce fainter, more rapidly evolving transients. 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 that are likely to be realized in nature. To quantify the observable signatures of core detonation, we compute spectra and light curves for models in which either an edge-lit or compression-triggered CO detonation is assumed to occur. We compare these to synthetic observables for 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 significantly 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, near-infrared observations and detailed spectroscopic analysis would be needed to determine whether a core detonation occurred. We comment on the implications of our results for understanding peculiar astrophysical transients including SN 2002bj, SN 2010X and SN 2005E. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

A theoretical color-velocity correlation for supernovae associated with gamma-ray bursts

We carry out the first multi-dimensional radiative transfer calculations to simultaneously compute synthetic spectra and light curves for models of supernovae driven by fast bipolar outflows. These allow us to make self-consistent predictions for the orientation dependence of both color evolution and spectral features. We compare models with different degrees of asphericity and metallicity and find significant observable consequences of both. In aspherical models, we find spectral and light curve features that vary systematically with observer orientation. In particular, we find that the early-phase light curves are brighter and bluer when viewed close to the polar axis but that the peak flux is highest for equatorial (off-axis) inclinations. Spectral line features also depend systematically on observer orientation, including the velocity of the Si II 6355 Å line. Consequently, our models predict a correlation between line velocity and color that could assist the identification of supernovae associated with off-axis jet-driven explosions. The amplitude and range of this correlation depends on the degree of asphericity, the metallicity, and the epoch of observation but we find that it is always present and acts in the same direction. © 2012. The American Astronomical Society. All rights reserved..

Normal type Ia supernovae from violent mergers of white dwarf binaries

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

On the γ-ray emission of Type Ia supernovae

A multidimension, time-dependent Monte Carlo code is used to compute sample ?-ray spectra to explore whether unambiguous constraints could be obtained from ?-ray observations of Type Ia supernovae. Both spherical and aspherical geometries are considered and it is shown that moderate departures from sphericity can produce viewing-angle effects that are at least as significant as those caused by the variation of key parameters in 1D models. Thus, ?-ray data could, in principle, carry some geometrical information, and caution should be applied when discussing the value of ?-ray data based only on 1D explosion models. In light of the limited sensitivity of current ?-ray observatories, the computed theoretical spectra are studied to revisit the issue of whether useful constraints could be obtained for moderately nearby objects. The most useful ?-ray measurements are likely to be of the light curve and time-dependent hardness ratios, but sensitivity higher than currently available, particularly at relatively hard energies (~2-3 MeV), is desirable. © 2008 The Authors. Journal compilation © 2008 RAS.

Time-dependent three-dimensional 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 ?-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 comoving frame. Assuming a photoionization-dominated plasma, the equations of ionization 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. © 2009 RAS.

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 to +42d 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 ~8d before and ~10d 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 Mgii at early epochs and a series of relatively unblended Coii lines at late phases allow us to constrain the regions of the ejecta in which the respective elements are abundant. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

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

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

Multidimensional simulations of radiative transfer in Type Ia supernovae

A three-dimensional Monte Carlo code for modelling radiation transport in Type Ia supernovae is described. In addition to tracking Monte Carlo quanta to follow the emission, scattering and deposition of radiative energy, a scheme involving volume-based Monte Carlo estimators is used to allow properties of the emergent radiation field to be extracted for specific viewing angles in a multidimensional structure. This eliminates the need to compute spectra or light curves by angular binning of emergent quanta. The code is applied to two test problems to illustrate consequences of multidimensional structure on the modelling of light curves. First, elliptical models are used to quantify how large-scale asphericity can introduce angular dependence to light curves. Secondly, a model which incorporates complex structural inhomogeneity, as predicted by modern explosion models, is used to investigate how such structure may affect light-curve properties. © 2006 RAS.

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

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

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 (SNe). 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 prototypical 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 d 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.03Mconsisting 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. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.

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

We present results for a suite of 14 three-dimensional, high-resolution hydrodynamical simulations of delayed-detonation models of Type Ia supernova (SN Ia) explosions. This model suite comprises the first set of three-dimensional SN Ia simulations with detailed isotopic yield information. As such, it may serve as a data base for Chandrasekhar-mass delayed-detonation model nucleosynthetic yields and for deriving synthetic observables such as spectra and light curves. We employ aphysically motivated, stochastic model based on turbulent velocity fluctuations and fuel density to calculate in situ the deflagration-to-detonation transition probabilities. To obtain different strengths of the deflagration phase and thereby different degrees of pre-expansion, we have chosen a sequence of initial models with 1, 3, 5, 10, 20, 40, 100, 150, 200, 300 and 1600 (two different realizations) ignition kernels in a hydrostatic white dwarf with a central density of 2.9 × 10 g cm, as well as one high central density (5.5 × 10 g cm) and one low central density (1.0 × 10 g cm) rendition of the 100 ignition kernel configuration. For each simulation, we determined detailed nucleosynthetic yields by postprocessing10 tracer particles with a 384 nuclide reaction network. All delayed-detonation models result in explosions unbinding thewhite dwarf, producing a range of 56Ni masses from 0.32 to 1.11M. As a general trend, the models predict that the stableneutron-rich iron-group isotopes are not found at the lowest velocities, but rather at intermediate velocities (~3000×10 000 km s) in a shell surrounding a Ni-rich core. The models further predict relatively low-velocity oxygen and carbon, with typical minimum velocities around 4000 and 10 000 km s, respectively. © 2012 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

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 (WDs) - 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 Co production. SN 2011fe is close enough to be followed sufficiently long to study this effect. © © 2012 The American Astronomical Society. All rights reserved.