ALBERT

Description: We present detailed spectroscopic analysis of the extraordinarily fast-evolving transient AT2018kzr. The transient’s observed light curve showed a rapid decline rate, comparable to the kilonova AT2017gfo. We calculate a self-consistent sequence of radiative transfer models (using tardis) and determine that the ejecta material is dominated by intermediate-mass elements (O, Mg, and Si), with a photospheric velocity of ∼12 000–14 500 $ m {km}, s^{-1}$. The early spectra have the unusual combination of being blue but dominated by strong Fe ii and Fe iii absorption features. We show that this combination is only possible with a high Fe content (3.5 per cent). This implies a high Fe/(Ni+Co) ratio. Given the short time from the transient’s proposed explosion epoch, the Fe cannot be 56Fe resulting from the decay of radioactive 56Ni synthesized in the explosion. Instead, we propose that this is stable 54Fe, and that the transient is unusually rich in this isotope. We further identify an additional, high-velocity component of ejecta material at ∼20 000–26 000 $ m {km}, s^{-1}$, which is mildly asymmetric and detectable through the Ca ii near-infrared triplet. We discuss our findings with reference to a range of plausible progenitor systems and compare with published theoretical work. We conclude that AT2018kzr is most likely the result of a merger between an ONe white dwarf and a neutron star or black hole. As such, it would be the second plausible candidate with a good spectral sequence for the electromagnetic counterpart of a compact binary merger, after AT2017gfo.

Description: We assemble a sample of 24 hydrogen-poor superluminous supernovae (SLSNe). Parameterizing the light-curve shape through rise and decline time-scales 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 time-scale. 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 mag brighter and have light curves three 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 (6 M ), with a range of 3–30 M . 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.

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

Description: 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 M , hybrid WDs in a binary system may easily approach the Chandrasekhar mass ( M Ch ) by accretion and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in a near- M Ch hybrid WD under the assumption that nuclear burning only occurs in carbon-rich material. Performing hydrodynamics simulations of the explosion and detailed nucleosynthesis post-processing calculations, we find that only 0.014 M of material is ejected while the remainder of the mass stays bound. The ejecta consist predominantly of iron-group elements, O, C, Si and S. We also calculate synthetic observables for our model and find reasonable agreement with the faint Type Iax SN 2008ha. This shows for the first time that deflagrations in near- M Ch WDs can in principle explain the observed diversity of Type Iax supernovae. Leaving behind a near- M Ch bound remnant opens the possibility for recurrent explosions or a subsequent accretion-induced collapse in faint Type Iax SNe, if further accretion episodes occur. From binary population synthesis calculations, we find the rate of hybrid WDs approaching M Ch to be of the order of 1 per cent of the Galactic SN Ia rate.

Description: Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 10 43  erg s –1 . Here, we introduce a simple treatment of clumping, and find that a filling factor of ~0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C iv 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.

Description: Calculations of synthetic spectropolarimetry are one means to test multidimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and 0.9 M 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 2D double-detonation model and a 3D 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 ~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.

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

Description: We present tardis – an open-source code for rapid spectral modelling of supernovae (SNe). Our goal is to develop a tool that is sufficiently fast to allow exploration of the complex parameter spaces of models for SN ejecta. This can be used to analyse the growing number of high-quality SN spectra being obtained by transient surveys. The code uses Monte Carlo methods to obtain a self-consistent description of the plasma state and to compute a synthetic spectrum. It has a modular design to facilitate the implementation of a range of physical approximations that can be compared to assess both accuracy and computational expediency. This will allow users to choose a level of sophistication appropriate for their application. Here, we describe the operation of the code and make comparisons with alternative radiative transfer codes of differing levels of complexity ( syn ++, python and artis ). We then explore the consequence of adopting simple prescriptions for the calculation of atomic excitation, focusing on four species of relevance to Type Ia SN spectra – Si ii , S ii , Mg ii and Ca ii . We also investigate the influence of three methods for treating line interactions on our synthetic spectra and the need for accurate radiative rate estimates in our scheme.

Description: The double-detonation explosion scenario of Type Ia supernovae (SNe Ia) has gained increased support from the SN Ia community as a viable progenitor model, making it a promising candidate alongside the well-known single degenerate and double degenerate scenarios. We present delay times of double-detonation SNe, in which a sub-Chandrasekhar mass carbon–oxygen white dwarf (WD) accretes non-dynamically from a helium-rich companion. One of the main uncertainties in quantifying SN rates from double detonations is the (assumed) retention efficiency of He-rich matter. Therefore, we implement a new prescription for the treatment of accretion/accumulation of He-rich matter on WDs. In addition, we test how the results change depending on which criteria are assumed to lead to a detonation in the helium shell. In comparing the results to our standard case (Ruiter et al.), we find that regardless of the adopted He accretion prescription, the SN rates are reduced by only ~25 per cent if low-mass He shells (0.05 M ) are sufficient to trigger the detonations. If more massive (0.1 M ) shells are needed, the rates decrease by 85 per cent and the delay time distribution is significantly changed in the new accretion model – only SNe with prompt (〈500 Myr) delay times are produced. Since theoretical arguments favour low-mass He shells for normal double-detonation SNe, we conclude that the rates from double detonations are likely to be high, and should not critically depend on the adopted prescription for accretion of He.

Description: The ejected mass distribution of Type Ia supernovae (SNe Ia) 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 SNe Ia, models allowing SNe Ia 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 ej and 56 Ni masses M Ni for a sample of 337 SNe Ia with redshifts z  〈 0.7 used in recent cosmological analyses. We use hierarchical Bayesian inference to reconstruct the joint M ej – M Ni distribution, accounting for measurement errors. The inferred marginal distribution of M ej has a long tail towards sub-Chandrasekhar masses, but cuts off sharply above 1.4 M . Our results imply that 25–50 per cent of normal SNe Ia 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 per cent of all spectroscopically normal SNe Ia. We interpret the SN Ia width–luminosity relation as an underlying relation between M ej and M Ni , and show that the inferred relation is not naturally explained by the predictions of any single known explosion mechanism.