Showing 1–50 of 138 results for author: Sim, S A

Helium features are inconsistent with the spectral evolution of the kilonova AT2017gfo

Authors: Albert Sneppen, Rasmus Damgaard, Darach Watson, Christine E. Collins, Luke Shingles, Stuart A. Sim

Abstract: The spectral features observed in kilonovae (KNe) reveal the elemental composition and the velocity structures of matter ejected from neutron star mergers. In the spectra of the kilonova AT2017gfo, a P Cygni line at about 1$μ$m has been linked to Sr II, providing the first direct evidence of freshly synthesised $r$-process material. An alternative explanation to Sr II was proposed - He I… ▽ More The spectral features observed in kilonovae (KNe) reveal the elemental composition and the velocity structures of matter ejected from neutron star mergers. In the spectra of the kilonova AT2017gfo, a P Cygni line at about 1$μ$m has been linked to Sr II, providing the first direct evidence of freshly synthesised $r$-process material. An alternative explanation to Sr II was proposed - He I $λ1083.3$nm under certain non-local-thermodynamic-equilibrium (NLTE) conditions. A key way to robustly discriminate between these identifications, and indeed other proposed identifications, is to analyse the temporal emergence and evolution of the feature. In this analysis we trace the earliest appearance of the observed feature and detail its spectro-temporal evolution, which we compare with a collisional-radiative model of helium. We show that the 1$μ$m P Cygni line is inconsistent with a He I interpretation both in emergence time and in subsequent spectral evolution. Self-consistent helium masses cannot reproduce the observed feature, due to the diminishing strength of radiative pathways leaving triplet helium. △ Less

Submitted 17 July, 2024; originally announced July 2024.

Comments: Submitted to A&A, 10 pages, 6 figures. Comments are welcome

SN 2023zaw: the low-energy explosion of an ultra-stripped star, with non-radioactive heating

Authors: Thomas Moore, James Gillanders, Matt Nicholl, Mark Huber, Stephen Smartt, Shubham Srivastav, Heloise Stevance, Ting-Wan Chen, Kenneth Chambers, Joseph Anderson, Michael Fulton, Samantha Oates, Charlotte Angus, Giuliano Pignata, Nicolas Erasmus, Hua Gao, Joanna Bulger, Chien-Cheng Lin, Thomas Lowe, Eugene Magnier, Paloma Minguez, Chow-Choong Ngeow, Xinyue Sheng, Stuart A. Sim, Ken Smith , et al. (4 additional authors not shown)

Abstract: Most stripped envelope supernova progenitors are formed through binary interaction, losing hydrogen and/or helium from their outer layers. An emerging class of supernovae with the highest degree of envelope-stripping are thought to be the product of stripping by a NS companion. However, relatively few examples are known and the outcomes of such systems can be diverse and are poorly understood at p… ▽ More Most stripped envelope supernova progenitors are formed through binary interaction, losing hydrogen and/or helium from their outer layers. An emerging class of supernovae with the highest degree of envelope-stripping are thought to be the product of stripping by a NS companion. However, relatively few examples are known and the outcomes of such systems can be diverse and are poorly understood at present. Here, we present spectroscopic observations and high cadence multi-band photometry of SN 2023zaw, a low ejecta mass and rapidly evolving supernova. SN 2023zaw was discovered in a nearby spiral galaxy at D = 39.7 Mpc, with significant Milky Way extinction, $E(B-V) = 0.21$, and significant (but uncertain) host extinction. Bayesian evidence comparison reveals that nickel is not the only power source and an additional energy source is required to explain our observations. Our models suggest an ejecta mass of $M_{\rm ej} \sim 0.07\,\rm M_\odot$ and a synthesised nickel mass of $M_{\rm ej} \sim 0.007\,\rm M_\odot$ is required to explain the explosion. However an additional heating from a magnetar or interaction with circumstellar material is required to power the early light curve. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Including a Luminous Central Remnant in Radiative Transfer Simulations for Type Iax Supernovae

Authors: F. P. Callan, S. A. Sim, C. E. Collins, L. J. Shingles, F. Lach, F. K. Roepke, R. Pakmor, M. Kromer, S. Srivastav

Abstract: Type Iax supernovae (SNe Iax) are proposed to arise from deflagrations of Chandrasekhar mass white dwarfs (WDs). Previous deflagration simulations have achieved good agreement with the light curves and spectra of intermediate-luminosity and bright SNe Iax. However, the model light curves decline too quickly after peak, particularly in red optical and near-infrared (NIR) bands. Deflagration models… ▽ More Type Iax supernovae (SNe Iax) are proposed to arise from deflagrations of Chandrasekhar mass white dwarfs (WDs). Previous deflagration simulations have achieved good agreement with the light curves and spectra of intermediate-luminosity and bright SNe Iax. However, the model light curves decline too quickly after peak, particularly in red optical and near-infrared (NIR) bands. Deflagration models with a variety of ignition configurations do not fully unbind the WD, leaving a remnant polluted with $^{56}\mathrm{Ni}$. Emission from such a remnant may contribute to the luminosity of SNe Iax. Here we investigate the impact of adding a central energy source, assuming instantaneous powering by $^{56}\mathrm{Ni}$ decay in the remnant, in radiative transfer calculations of deflagration models. Including the remnant contribution improves agreement with the light curves of SNe Iax, particularly due to the slower post-maximum decline of the models. Spectroscopic agreement is also improved, with intermediate-luminosity and faint models showing greatest improvement. We adopt the full remnant $^{56}\mathrm{Ni}$ mass predicted for bright models, but good agreement with intermediate-luminosity and faint SNe Iax is only possible for remnant $^{56}\mathrm{Ni}$ masses significantly lower than those predicted. This may indicate that some of the $^{56}\mathrm{Ni}$ decay energy in the remnant does not contribute to the radiative luminosity but instead drives mass ejection, or that escape of energy from the remnant is significantly delayed. Future work should investigate the structure of remnants predicted by deflagration models and the potential roles of winds and delayed energy escape, as well as extend radiative transfer simulations to late times. △ Less

Submitted 19 April, 2024; v1 submitted 22 March, 2024; originally announced March 2024.

Comments: 17 pages, 6 figures. Lightcurves and spectra available at https://hesma.h-its.org

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 530, Issue 2, May 2024, Pages 1457 to 1473

Type Ia supernova explosion models are inherently multidimensional

Authors: R. Pakmor, I. R. Seitenzahl, A. J. Ruiter, S. A. Sim, F. K. Roepke, S. Taubenberger, R. Bieri, S. Blondin

Abstract: Theoretical and observational approaches to settling the important questions surrounding the progenitor systems and the explosion mechanism of normal Type Ia supernovae have thus far failed. With its unique capability to obtain continuous spectra through the near- and mid-infrared, JWST now offers completely new insights into Type Ia supernovae. In particular, observing them in the nebular phase a… ▽ More Theoretical and observational approaches to settling the important questions surrounding the progenitor systems and the explosion mechanism of normal Type Ia supernovae have thus far failed. With its unique capability to obtain continuous spectra through the near- and mid-infrared, JWST now offers completely new insights into Type Ia supernovae. In particular, observing them in the nebular phase allows us to directly see the central ejecta and thereby constrain the explosion mechanism. We aim to understand and quantify differences in the structure and composition of the central ejecta of various Type Ia supernova explosion models. We examined the currently most popular explosion scenarios using self-consistent multidimensional explosion simulations of delayed-detonation and pulsationally assisted, gravitationally confined delayed detonation Chandrasekhar-mass models and double-detonation sub-Chandrasekhar-mass and violent merger models. We find that the distribution of radioactive and stable nickel in the final ejecta, both observable in nebular spectra, are significantly different between different explosion scenarios. Therefore, comparing synthetic nebular spectra with JWST observations should allow us to distinguish between explosion models. We show that the explosion ejecta are inherently multidimensional for all models, and the Chandrasekhar-mass explosions simulated in spherical symmetry in particular lead to a fundamentally unphysical ejecta structure. Moreover, we show that radioactive and stable nickel cover a significant range of densities at a fixed velocity of the homologously expanding ejecta. Any radiation transfer postprocessing has to take these variations into account to obtain faithful synthetic observables; this will likely require multidimensional radiation transport simulations. △ Less

Submitted 26 April, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: 7 pages, 2 figures, accepted by A&A, comments welcome

Journal ref: A&A 686, A227 (2024)

State-of-the-art simulations of line-driven accretion disc winds: realistic radiation-hydrodynamics leads to weaker outflows

Authors: Nick Higginbottom, Nicolas Scepi, Christian Knigge, Knox S. Long, James H. Matthews, Stuart A. Sim

Abstract: Disc winds are a common feature in accreting astrophysical systems on all scales. In active galactic nuclei (AGN) and accreting white dwarfs (AWDs), specifically, radiation pressure mediated by spectral lines is a promising mechanism for driving these outflows. Previous hydrodynamical simulations have largely supported this idea, but relied on highly approximate treatments of ionization and radiat… ▽ More Disc winds are a common feature in accreting astrophysical systems on all scales. In active galactic nuclei (AGN) and accreting white dwarfs (AWDs), specifically, radiation pressure mediated by spectral lines is a promising mechanism for driving these outflows. Previous hydrodynamical simulations have largely supported this idea, but relied on highly approximate treatments of ionization and radiative transfer. Given the sensitivity of line driving to the ionization state and radiation field in the outflow, here we present a new method for carrying out 2.5D radiation-hydrodynamic simulations that takes full account of the frequency-dependent radiative transfer through the wind, the corresponding ionization state and the resulting radiative accelerations. Applying our method to AWDs, we find that it is much harder to drive a powerful line-driven outflow when the interaction between matter and radiation is treated self-consistently. This conclusion is robust to changes in the adopted system parameters. The fundamental difficulty is that discs luminous enough to drive such a wind are also hot enough to over-ionize it. As a result, the mass-loss rates in our simulations are much lower than those found in earlier, more approximate calculations. We also show that the ultraviolet spectra produced by our simulations do not match those observed in AWDs. We conclude that, unless the over-ionization problem can be mitigated (e.g. by sub-grid clumping or a softer-than-expected radiation field), line driving may not be a promising mechanism for powering the outflows from AWDs. These conclusions are likely to have significant implications for disc winds in AGN also. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS. 14 pages, 10 figures + 3 figures in Appendix

Towards inferring the geometry of kilonovae

Authors: Christine E. Collins, Luke J. Shingles, Andreas Bauswein, Stuart A. Sim, Theodoros Soultanis, Vimal Vijayan, Andreas Floers, Oliver Just, Gerrit Leck, Georgios Lioutas, Gabriel Martínez-Pinedo, Albert Sneppen, Darach Watson, Zewei Xiong

Abstract: Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hy… ▽ More Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hydrodynamical merger simulation can be compatible with the observational constraints suggesting a high degree of sphericity in AT2017gfo. Specifically, we determine whether fitting a simple P-Cygni line profile model leads to a value for the photospheric velocity that is consistent with the value obtained from the expanding photosphere method. We would infer that our kilonova simulation is highly spherical at early times, when the spectra resemble a blackbody distribution. The two independently inferred photospheric velocities can be very similar, implying a high degree of sphericity, which can be as spherical as inferred for AT2017gfo, demonstrating that the photosphere can appear spherical even for asymmetrical ejecta. The last-interaction velocities of radiation escaping the simulation show a high degree of sphericity, supporting the inferred symmetry of the photosphere. We find that when the synthetic spectra resemble a blackbody the expanding photosphere method can be used to obtain an accurate luminosity distance (within 4-7 per cent). △ Less

Submitted 23 February, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted by MNRAS

Unprecedented early flux excess in the hybrid 02es-like type Ia supernova 2022ywc indicates interaction with circumstellar material

Authors: Shubham Srivastav, T. Moore, M. Nicholl, M. R. Magee, S. J. Smartt, M. D. Fulton, S. A. Sim, J. M. Pollin, L. Galbany, C. Inserra, A. Kozyreva, Takashi J. Moriya, F. P. Callan, X. Sheng, K. W. Smith, J. S. Sommer, J. P. Anderson, M. Deckers, M. Gromadzki, T. E. Müller-Bravo, G. Pignata, A. Rest, D. R. Young

Abstract: We present optical photometric and spectroscopic observations of the 02es-like type Ia supernova (SN) 2022ywc. The transient occurred in the outskirts of an elliptical host galaxy and showed a striking double-peaked light curve with an early excess feature detected in the ATLAS orange and cyan bands. The early excess is remarkably luminous with an absolute magnitude $\sim -19$, comparable in lumin… ▽ More We present optical photometric and spectroscopic observations of the 02es-like type Ia supernova (SN) 2022ywc. The transient occurred in the outskirts of an elliptical host galaxy and showed a striking double-peaked light curve with an early excess feature detected in the ATLAS orange and cyan bands. The early excess is remarkably luminous with an absolute magnitude $\sim -19$, comparable in luminosity to the subsequent radioactively-driven second peak. The spectra resemble the hybrid 02es-like SN 2016jhr, that is considered to be a helium shell detonation candidate. We investigate different physical mechanisms that could power such a prominent early excess and rule out massive helium shell detonation, surface $^{56}$Ni distribution and ejecta-companion interaction. We conclude that SN ejecta interacting with circumstellar material (CSM) is the most viable scenario. Semi-analytical modelling with MOSFiT indicates that SN ejecta interacting with $\sim 0.05\,$M$_{\odot}$ of CSM at a distance of $\sim 10^{14}$ cm can explain the extraordinary light curve. A double-degenerate scenario may explain the origin of the CSM, either by tidally-stripped material from the secondary white dwarf, or disk-originated matter launched along polar axes following the disruption and accretion of the secondary white dwarf. A non-spherical CSM configuration could suggest that a small fraction of 02es-like events viewed along a favourable line of sight may be expected to display a very conspicuous early excess like SN 2022ywc. △ Less

Submitted 25 September, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: Accepted to ApJL after minor revision

Helium as a signature of the double detonation in Type Ia supernovae

Authors: Christine E. Collins, Stuart A. Sim, Luke. J. Shingles, Sabrina Gronow, Friedrich K. Roepke, Ruediger Pakmor, Ivo R. Seitenzahl, Markus Kromer

Abstract: The double detonation is a widely discussed mechanism to explain Type Ia supernovae from explosions of sub-Chandrasekhar mass white dwarfs. In this scenario, a helium detonation is ignited in a surface helium shell on a carbon/oxygen white dwarf, which leads to a secondary carbon detonation. Explosion simulations predict high abundances of unburnt helium in the ejecta, however, radiative transfer… ▽ More The double detonation is a widely discussed mechanism to explain Type Ia supernovae from explosions of sub-Chandrasekhar mass white dwarfs. In this scenario, a helium detonation is ignited in a surface helium shell on a carbon/oxygen white dwarf, which leads to a secondary carbon detonation. Explosion simulations predict high abundances of unburnt helium in the ejecta, however, radiative transfer simulations have not been able to fully address whether helium spectral features would form. This is because helium can not be sufficiently excited to form spectral features by thermal processes, but can be excited by collisions with non-thermal electrons, which most studies have neglected. We carry out a full non-local thermodynamic equilibrium (non-LTE) radiative transfer simulation for an instance of a double detonation explosion model, and include a non-thermal treatment of fast electrons. We find a clear He I λ 10830 feature which is strongest in the first few days after explosion and becomes weaker with time. Initially this feature is blended with the Mg II λ 10927 feature but over time separates to form a secondary feature to the blue wing of the Mg II λ 10927 feature. We compare our simulation to observations of iPTF13ebh, which showed a similar feature to the blue wing of the Mg II λ 10927 feature, previously identified as C I. Our simulation shows a good match to the evolution of this feature and we identify it as high velocity He I λ 10830. This suggests that He I λ 10830 could be a signature of the double detonation scenario. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 7 pages, accepted by MNRAS

Self-consistent 3D radiative transfer for kilonovae: directional spectra from merger simulations

Authors: Luke J. Shingles, Christine E. Collins, Vimal Vijayan, Andreas Flörs, Oliver Just, Gerrit Leck, Zewei Xiong, Andreas Bauswein, Gabriel Martínez-Pinedo, Stuart A. Sim

Abstract: We present three-dimensional radiative transfer calculations for the ejecta from a neutron star merger that include line-by-line opacities for tens of millions of bound-bound transitions, composition from an r-process nuclear network, and time-dependent thermalization of decay products from individual $α$ and $β^-$ decay reactions. In contrast to expansion opacities and other wavelength-binned tre… ▽ More We present three-dimensional radiative transfer calculations for the ejecta from a neutron star merger that include line-by-line opacities for tens of millions of bound-bound transitions, composition from an r-process nuclear network, and time-dependent thermalization of decay products from individual $α$ and $β^-$ decay reactions. In contrast to expansion opacities and other wavelength-binned treatments, a line-by-line treatment enables us include fluorescence effects and associate spectral features with the emitting and absorbing lines of individual elements. We find variations in the synthetic observables with both the polar and azimuthal viewing angles. The spectra exhibit blended features with strong interactions by Ce III, Sr II, Y II, and Zr II that vary with time and viewing direction. We demonstrate the importance of wavelength-calibration of atomic data using a model with calibrated Sr, Y, and Zr data, and find major differences in the resulting spectra, including a better agreement with AT2017gfo. The synthetic spectra for near-polar inclination show a feature at around 8000 A, similar to AT2017gfo. However, they evolve on a more rapid timescale, likely due to the low ejecta mass (0.005 M$_\odot$) as we take into account only the early ejecta. The comparatively featureless spectra for equatorial observers gives a tentative prediction that future observations of edge-on kilonovae will appear substantially different from AT2017gfo. We also show that 1D models obtained by spherically averaging the 3D ejecta lead to dramatically different direction-integrated luminosities and spectra compared to full 3D calculations. △ Less

Submitted 1 September, 2023; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: 12 pages, 5 figures. Accepted by ApJL

Journal ref: 2023, ApJL, 954, L41

Modelling the spectra of the kilonova AT2017gfo -- II: Beyond the photospheric epochs

Authors: J. H. Gillanders, S. A. Sim, S. J. Smartt, S. Goriely, A. Bauswein

Abstract: Binary neutron star mergers are the first confirmed site of element nucleosynthesis by the rapid neutron-capture process (r-process). The kilonova AT2017gfo is the only electromagnetic counterpart of a neutron star merger spectroscopically observed. We analyse the entire spectral sequence of AT2017gfo (from merger to +10.4 days) and identify seven emission-like features. We confirm that the promin… ▽ More Binary neutron star mergers are the first confirmed site of element nucleosynthesis by the rapid neutron-capture process (r-process). The kilonova AT2017gfo is the only electromagnetic counterpart of a neutron star merger spectroscopically observed. We analyse the entire spectral sequence of AT2017gfo (from merger to +10.4 days) and identify seven emission-like features. We confirm that the prominent 1.08 um feature can be explained by the Sr II near-infrared triplet evolving from a P-Cygni profile through to pure emission. We calculate the expected strength of the [Sr II] doublet and show that its absence requires highly clumped ejecta. Near-infrared features at 1.58 and 2.07 um emerge after three days and become more prominent as the spectra evolve. We model these as optically thick P-Cygni profiles and alternatively as pure emission features (with FWHM = 35600 +/- 6600 km/s), and favour the latter interpretation. The profile of the strong 2.07 um emission feature is best reproduced with two lines, centred at 2.059 and 2.135 um. We search for candidate ions for all prominent features in the spectra. Strong, permitted transitions of La III, Ce III, Gd III, Ra II and Ac I are plausible candidates for the emission features. If any of these features are produced by intrinsically weak, forbidden transitions, we highlight candidate ions spanning the three r-process peaks. The second r-process peak elements Te and I have plausible matches to multiple features. We highlight the need for more detailed and quantitative atomic line transition data. △ Less

Submitted 27 November, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: Updating to the final accepted version in MNRAS. Main text: 20 pages, with 7 figures and 6 tables. Appendices: 9 pages with 3 tables

Photometry and spectroscopy of the Type Icn supernova 2021ckj: The diverse properties of the ejecta and circumstellar matter of Type Icn SNe

Authors: T. Nagao, H. Kuncarayakti, K. Maeda, T. Moore, A. Pastorello, S. Mattila, K. Uno, S. J. Smartt, S. A. Sim, L. Ferrari, L. Tomasella, J. P. Anderson, T. -W. Chen, L. Galbany, H. Gao, M. Gromadzki, C. P. Gutiérrez, C. Inserra, E. Kankare, E. A. Magnier, T. E. Müller-Bravo, A. Reguitti, D. R. Young

Abstract: We present photometric and spectroscopic observations of the Type Icn supernova (SN) 2021ckj. Spectral modeling of SN 2021ckj reveals that its composition is dominated by oxygen, carbon and iron group elements, and the photospheric velocity at peak is ~10000 km/s. From the light curve (LC) modeling applied to SNe 2021ckj, 2019hgp, and 2021csp, we find that the ejecta and CSM properties of Type Icn… ▽ More We present photometric and spectroscopic observations of the Type Icn supernova (SN) 2021ckj. Spectral modeling of SN 2021ckj reveals that its composition is dominated by oxygen, carbon and iron group elements, and the photospheric velocity at peak is ~10000 km/s. From the light curve (LC) modeling applied to SNe 2021ckj, 2019hgp, and 2021csp, we find that the ejecta and CSM properties of Type Icn SNe are diverse. SNe 2021ckj and 2021csp likely have two ejecta components (an aspherical high-energy component and a spherical standard-energy component) with a roughly spherical CSM, while SN 2019hgp can be explained by a spherical ejecta-CSM interaction alone. The ejecta of SNe 2021ckj and 2021csp have larger energy per ejecta mass than the ejecta of SN 2019hgp. The density distribution of the CSM is similar in these three SNe, and is comparable to those of Type Ibn SNe. This may imply that the mass-loss mechanism is common between Type Icn (and also Type Ibn) SNe. The CSM masses of SN 2021ckj and SN 2021csp are higher than that of SN 2019hgp, although all these values are within the diversity seen in Type Ibn SNe. The early spectrum of SN 2021ckj shows narrow emission lines from C II and C III, without a clear absorption component, in contrast with that observed in SN 2021csp. The similarity of the emission components of these lines implies that the emitting regions of SNe 2021ckj and 2021csp have similar ionization states, and thus suggests that they have similar properties of the ejecta and CSM, which is inferred also from the LC modeling. Taking into account the difference in the strength of the absorption features, this heterogeneity may be attributed to viewing angle effects in otherwise common aspherical ejecta. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 13 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 673, A27 (2023)

The optical light curve of GRB 221009A: the afterglow and the emerging supernova

Authors: M. D. Fulton, S. J. Smartt, L. Rhodes, M. E. Huber, A. V. Villar, T. Moore, S. Srivastav, A. S. B. Schultz, K. C. Chambers, L. Izzo, J. Hjorth, T. -W. Chen, M. Nicholl, R. J. Foley, A. Rest, K. W. Smith, D. R. Young, S. A. Sim, J. Bright, Y. Zenati, T. de Boer, J. Bulger, J. Fairlamb, H. Gao, C. -C. Lin , et al. (24 additional authors not shown)

Abstract: We present extensive optical photometry of the afterglow of GRB~221009A. Our data cover $0.9 - 59.9$\,days from the time of \textit{Swift} and \textit{Fermi} GRB detections. Photometry in $rizy$-band filters was collected primarily with Pan-STARRS and supplemented by multiple 1- to 4-meter imaging facilities. We analyzed the Swift X-ray data of the afterglow and found a single decline rate power-l… ▽ More We present extensive optical photometry of the afterglow of GRB~221009A. Our data cover $0.9 - 59.9$\,days from the time of \textit{Swift} and \textit{Fermi} GRB detections. Photometry in $rizy$-band filters was collected primarily with Pan-STARRS and supplemented by multiple 1- to 4-meter imaging facilities. We analyzed the Swift X-ray data of the afterglow and found a single decline rate power-law $f(t) \propto t^{-1.556\pm0.002}$ best describes the light curve. In addition to the high foreground Milky Way dust extinction along this line of sight, the data favour additional extinction to consistently model the optical to X-ray flux with optically thin synchrotron emission. We fit the X-ray-derived power-law to the optical light curve and find good agreement with the measured data up to $5-6$\,days. Thereafter we find a flux excess in the $riy$ bands which peaks in the observer frame at $\sim20$\,days. This excess shares similar light curve profiles to the type Ic broad-lined supernovae SN~2016jca and SN~2017iuk once corrected for the GRB redshift of $z=0.151$ and arbitrarily scaled. This may be representative of a supernova emerging from the declining afterglow. We measure rest-frame absolute peak AB magnitudes of $M_g=-19.8\pm0.6$ and $M_r=-19.4\pm0.3$ and $M_z=-20.1\pm0.3$. If this is an SN component, then Bayesian modelling of the excess flux would imply explosion parameters of $M_{\rm ej}=7.1^{+2.4}_{-1.7}$ M$_{\odot}$, $M_{\rm Ni}=1.0^{+0.6}_{-0.4}$ M$_{\odot}$, and $v_{\rm ej}=33,900^{+5,900}_{-5,700} kms^{-1}$, for the ejecta mass, nickel mass and ejecta velocity respectively, inferring an explosion energy of $E_{\rm kin}\simeq 2.6-9.0\times10^{52}$ ergs. △ Less

Submitted 23 March, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: Submitted to ApJL on 23rd January 2023, for consideration for publication in the special issue on GRB 221009A. Accepted on 2nd March 2023. The results of this paper are under press embargo until 28th March 2023. 15 pages, 7 figures, 1 table

The luminous type Ia supernova 2022ilv and its early excess emission

Authors: Shubham Srivastav, S. J. Smartt, M. E. Huber, G. Dimitriadis, K. C. Chambers, Michael D. Fulton, Thomas Moore, F. P. Callan, James H. Gillanders, K. Maguire, M. Nicholl, Luke J. Shingles, S. A. Sim, K. W. Smith, J. P. Anderson, Thomas de Boer, Ting-Wan Chen, Hua Gao, D. R. Young

Abstract: We present observations and analysis of the host-less and luminous type Ia supernova 2022ilv, illustrating it is part of the 2003fg-like family, often referred to as super-Chandrasekhar (Ia-SC) explosions. The ATLAS light curve shows evidence of a short-lived, pulse-like early excess, similar to that detected in another luminous type Ia supernova (SN 2020hvf). The light curve is broad and the earl… ▽ More We present observations and analysis of the host-less and luminous type Ia supernova 2022ilv, illustrating it is part of the 2003fg-like family, often referred to as super-Chandrasekhar (Ia-SC) explosions. The ATLAS light curve shows evidence of a short-lived, pulse-like early excess, similar to that detected in another luminous type Ia supernova (SN 2020hvf). The light curve is broad and the early spectra are remarkably similar to SN 2009dc. Adopting a redshift of $z=0.026 \pm 0.005$ for SN 2022ilv based on spectral matching, our model light curve requires a large $^{56}$Ni mass in the range $0.7-1.5$ M$_{\odot}$, and a large ejecta mass in the range $1.6-2.3$ M$_{\odot}$. The early excess can be explained by fast-moving SN ejecta interacting with a thin, dense shell of circumstellar material close to the progenitor ($\sim 10^{13}$ cm), a few hours after the explosion. This may be realised in a double-degenerate scenario, wherein a white dwarf merger is preceded by ejection of a small amount ($\sim 10^{-3}-10^{-2}$ M$_{\odot}$) of hydrogen and helium-poor tidally stripped material. A deep pre-explosion Pan-STARRS1 stack indicates no host galaxy to a limiting magnitude of $r \sim 24.5$. This implies a surprisingly faint limit for any host of $M_r \gtrsim -11$, providing further evidence that these types of explosion occur predominantly in low-metallicity environments. △ Less

Submitted 22 January, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: Accepted to ApJL after minor revision

StaNdaRT: A repository of standardized test models and outputs for supernova radiative transfer

Authors: Stéphane Blondin, Sergei Blinnikov, Fionntan P. Callan, Christine E. Collins, Luc Dessart, Wesley Even, Andreas Flörs, Andrew G. Fullard, D. John Hillier, Anders Jerkstrand, Daniel Kasen, Boaz Katz, Wolfgang Kerzendorf, Alexandra Kozyreva, Jack O'Brien, Ezequiel A. Pássaro, Nathaniel Roth, Ken J. Shen, Luke Shingles, Stuart A. Sim, Jaladh Singhal, Isaac G. Smith, Elena Sorokina, Victor P. Utrobin, Christian Vogl , et al. (4 additional authors not shown)

Abstract: We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M… ▽ More We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M$_\odot$) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time, are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe in detail the test models, radiative-transfer codes and output formats and provide access to the repository. We present example results of several key diagnostic features. △ Less

Submitted 15 April, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: Accepted for publication in A&A. 27 pages, 12 figures (v4: updated to match published version). The ejecta models and output files from the simulations are available at https://github.com/sn-rad-trans/data1

Journal ref: A&A 668, A163 (2022)

3D radiative transfer kilonova modelling for binary neutron star merger simulations

Authors: Christine E. Collins, Andreas Bauswein, Stuart A. Sim, Vimal Vijayan, Gabriel Martínez-Pinedo, Oliver Just, Luke J. Shingles, Markus Kromer

Abstract: The detection of GW170817 and the accompanying electromagnetic counterpart, AT2017gfo, have provided an important set of observational constraints for theoretical models of neutron star mergers, nucleosynthesis, and radiative transfer for kilonovae. We apply the 3D Monte Carlo radiative transfer code ARTIS to produce synthetic light curves of the dynamical ejecta from a neutron star merger, which… ▽ More The detection of GW170817 and the accompanying electromagnetic counterpart, AT2017gfo, have provided an important set of observational constraints for theoretical models of neutron star mergers, nucleosynthesis, and radiative transfer for kilonovae. We apply the 3D Monte Carlo radiative transfer code ARTIS to produce synthetic light curves of the dynamical ejecta from a neutron star merger, which has been modelled with 3D smooth-particle hydrodynamics (SPH) and included neutrino interactions. Nucleosynthesis calculations provide the energy released from radioactive decays of r-process nuclei, and radiation transport is performed using grey opacities given as functions of the electron fraction. We present line-of-sight dependent bolometric light curves, and find the emission along polar lines of sight to be up to a factor of ~2 brighter than along equatorial lines of sight. Instead of a distinct emission peak, our bolometric light curve exhibits a monotonic decline, characterised by a shoulder at the time when the bulk ejecta becomes optically thin. We show approximate band light curves based on radiation temperatures and compare these to the observations of AT2017gfo. We find that the rapidly declining temperatures lead to a blue to red colour evolution similar to that shown by AT2017gfo. We also investigate the impact of an additional, spherically symmetric secular ejecta component, and we find that the early light curve remains nearly unaffected, while after about 1 day the emission is strongly enhanced and dominated by the secular ejecta, leading to the shift of the shoulder from 1-2 to 6-10 days. △ Less

Submitted 23 February, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: Accepted by MNRAS

Double detonations: variations in Type Ia supernovae due to different core and He shell masses -- II: synthetic observables

Authors: Christine E. Collins, Sabrina Gronow, Stuart A. Sim, Friedrich K. Roepke

Abstract: Double detonations of sub-Chandrasekhar mass white dwarfs are a promising explosion scenario for Type Ia supernovae, whereby a detonation in a surface helium shell triggers a secondary detonation in a carbon-oxygen core. Recent work has shown that low mass helium shell models reproduce observations of normal SNe Ia. We present 3D radiative transfer simulations for a suite of 3D simulations of the… ▽ More Double detonations of sub-Chandrasekhar mass white dwarfs are a promising explosion scenario for Type Ia supernovae, whereby a detonation in a surface helium shell triggers a secondary detonation in a carbon-oxygen core. Recent work has shown that low mass helium shell models reproduce observations of normal SNe Ia. We present 3D radiative transfer simulations for a suite of 3D simulations of the double detonation explosion scenario for a range of shell and core masses. We find light curves broadly able to reproduce the faint end of the width-luminosity relation shown by SNe Ia, however, we find that all of our models show extremely red colours, not observed in normal SNe Ia. This includes our lowest mass helium shell model. We find clear Ti II absorption features in the model spectra, which would lead to classification as peculiar SNe Ia, as well as line blanketing in some lines of sight by singly ionised Cr and Fe-peak elements. Our radiative transfer simulations show that these explosion models remain promising to explain peculiar SNe Ia. Future full non-LTE simulations may improve the agreement of these explosion models with observations of normal SNe Ia. △ Less

Submitted 5 July, 2023; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: Spectra available at https://doi.org/10.5281/zenodo.7997388

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

Authors: G. A. Matzeu, M. Lieu, M. T. Costa, J. N. Reeves, V. Braito, M. Dadina, E. Nardini, P. G. Boorman, M. L. Parker, S. A. Sim, D. Barret, E. Kammoun, R. Middei, M. Giustini, M. Brusa, J. Pérez Cabrera, S. Marchesi

Abstract: 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 s… ▽ More 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. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 20 pages, 17 figures, MNRAS accepted for publication

Modelling the ionisation state of Type Ia supernovae in the nebular-phase

Authors: Luke J. Shingles, Andreas Flörs, Stuart A. Sim, Christine E. Collins, Friedrich K. Roepke, Ivo R. Seitenzahl, Ken J. Shen

Abstract: The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a me… ▽ More The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a method to determine the collisional excitation conditions from [Fe II] line ratios independently from the ionisation state and find that it cannot be applied to highly-ionised models due to the influence of recombination cascades on Fe$^{+}$ level populations. When the ionisation state is artificially lowered, the line ratios (and excitation conditions) are too similar to distinguish between explosion scenarios. We investigate changes to the treatment of non-thermal energy deposition as a way to reconcile over-ionised theoretical models with observations and find that a simple work function approximation provides closer agreement with the data for sub-Mch models than a detailed Spencer-Fano treatment with widely-used cross section data. To quantify the magnitude of additional heating processes that would be required to sufficiently reduce ionisation from fast leptons, we artificially boost the rate of energy loss to free electrons. We find that the equivalent of as much as an eight times increase to the plasma loss rate would be needed to reconcile the sub-Mch model with observed spectra. Future studies could distinguish between reductions in the non-thermal ionisation rates and increased recombination rates, such as by clumping. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 13 pages. Accepted by MNRAS

On the fate of the secondary white dwarf in double-degenerate double-detonation Type Ia supernovae

Authors: R. Pakmor, F. P. Callan, C. E. Collins, S. E. de Mink, A. Holas, W. E. Kerzendorf, M. Kromer, P. G. Neunteufel, John T. O'Brien, F. K. Roepke, A. J. Ruiter, I. R. Seitenzahl, Luke J. Shingles, S. A. Sim, S. Taubenberger

Abstract: The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and exp… ▽ More The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and explodes it. We investigate the fate of the secondary white dwarf, focusing on changes of the ejecta and observables of the explosion if the secondary explodes as well rather than survives. We simulate a binary system of a $1.05\,\mathrm{M_\odot}$ and a $0.7\,\mathrm{M_\odot}$ carbon-oxygen white dwarf with $0.03\,\mathrm{M_\odot}$ helium shells each. We follow the system self-consistently from inspiral to ignition, through the explosion, to synthetic observables. We confirm that the primary white dwarf explodes self-consistently. The helium detonation around the secondary white dwarf, however, fails to ignite a carbon detonation. We restart the simulation igniting the carbon detonation in the secondary white dwarf by hand and compare the ejecta and observables of both explosions. We find that the outer ejecta at $v~>~15\,000$\,km\,s$^{-1}$ are indistinguishable. Light curves and spectra are very similar until $\sim~40\,\mathrm{d}$ after explosion and the ejecta are much more spherical than violent merger models. The inner ejecta differ significantly slowing down the decline rate of the bolometric light curve after maximum of the model with a secondary explosion by $\sim20$ per cent. We expect future synthetic 3D nebular spectra to confirm or rule out either model. △ Less

Submitted 25 October, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 12 pages, 7 figures, accepted to MNRAS, comments welcome

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

Authors: J. H. Gillanders, S. J. Smartt, S. A. Sim, A. Bauswein, S. Goriely

Abstract: 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… ▽ More 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. △ Less

Submitted 3 May, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: 20 pages, 14 figures and 3 tables. Accepted for publication in MNRAS

Optical line spectra of tidal disruption events from reprocessing in optically thick outflows

Authors: Edward J. Parkinson, Christian Knigge, James H. Matthews, Knox S. Long, Nick Higginbottom, Stuart A. Sim, Samuel W. Mangham

Abstract: A significant number of tidal disruption events (TDEs) radiate primarily at optical and ultraviolet (UV) wavelengths, with only weak soft X-ray components. One model for this optical excess proposes that thermal X-ray emission from a compact accretion disc is reprocessed to longer wavelengths by an optically thick envelope. Here, we explore this reprocessing scenario in the context of an optically… ▽ More A significant number of tidal disruption events (TDEs) radiate primarily at optical and ultraviolet (UV) wavelengths, with only weak soft X-ray components. One model for this optical excess proposes that thermal X-ray emission from a compact accretion disc is reprocessed to longer wavelengths by an optically thick envelope. Here, we explore this reprocessing scenario in the context of an optically thick accretion disc wind. Using state-of-the-art Monte Carlo radiative transfer and ionization software, we produce synthetic UV and optical spectra for wind and disc-hosting TDEs. Our models are inspired by observations, spanning a realistic range of accretion rates and wind kinematics. We find that such outflows can efficiently reprocess the disc emission and produce the broad Balmer and helium recombination features commonly seen in TDEs and exhibit asymmetric red wings. Moreover, the characteristic colour temperature of the reprocessed spectral energy distribution (SED) is much lower than that of the accretion disc. We show explicitly how changes in black hole mass, accretion rate and wind properties affect the observed broadband SED and line spectrum. In general, slower, denser winds tend to reprocess more radiation and produce stronger Balmer emission. Most of the outflows we consider are too highly ionized to produce UV absorption features, but this is sensitive to the input SED. For example, truncating the inner disc at just 4 $R_{ISCO}$ lowers the wind ionization state sufficiently to produce UV absorption features for sight lines looking into the wind △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: Accepted to MNRAS. 19 pages, 11 figures

Models of pulsationally assisted gravitationally confined detonations with different ignition conditions

Authors: F. Lach, F. P. Callan, S. A. Sim, F. K. Roepke

Abstract: 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… ▽ More 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. △ Less

Submitted 19 March, 2022; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: 15 pages, 12 figures, published in A&A, corrected reference

Journal ref: A&A 659, A27 (2022)

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

Authors: Shubham Srivastav, S. J. Smartt, M. E. Huber, K. C. Chambers, C. R. Angus, T. -W. Chen, F. P. Callan, J. H. Gillanders, O. R. McBrien, S. A. Sim, M. Fulton, J. Hjorth, K. W. Smith, D. R. Young, K. Auchettl, J. P. Anderson, G. Pignata, T. J. L. de Boer, C. -C. Lin, E. A. Magnier

Abstract: 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… ▽ More 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. △ Less

Submitted 18 January, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: Accepted to MNRAS after minor revision

An analysis of the spectroscopic signatures of layering in the ejecta of type Iax supernovae

Authors: M. R. Magee, J. H. Gillanders, K. Maguire, S. A. Sim, F. P. Callan

Abstract: Investigations of some type Iax supernovae have led to the suggestion that their ejecta must be layered to some degree. Such an ejecta structure has been argued as inconsistent with the well-mixed composition predicted by pure deflagrations. Based on explosion models, we create toy models in which the ejecta are artificially stratified and progressively mixed until a uniform composition is obtaine… ▽ More Investigations of some type Iax supernovae have led to the suggestion that their ejecta must be layered to some degree. Such an ejecta structure has been argued as inconsistent with the well-mixed composition predicted by pure deflagrations. Based on explosion models, we create toy models in which the ejecta are artificially stratified and progressively mixed until a uniform composition is obtained. We find that models that are heavily mixed, containing burned and unburned material at all velocities, produce reasonably good agreement with SN 2012Z, for which a layered structure has been suggested. We also discuss how existing ejecta compositions determined for type Iax supernovae do not necessarily contradict pure deflagration models and may be consistent with a steeper density profile. We investigate previous claims that differences in line profile shapes may be due to strong blending, by presenting a series of models with different plasma states. These models indicate that blending could indeed explain differences in the observed profiles. Alternatively, stratification could also explain such differences, however all of our models indicate that this does not necessarily require stratification in abundance. Sufficient stratification in ionisation state can be achieved even for a well-mixed model. Based on our analysis, we demonstrate that there is insufficient evidence to suggest the ejecta of type Iax supernovae must be layered and therefore argue the pure deflagration scenario is not ruled out, even for the brightest type Iax supernovae. Our analysis does not indicate the ejecta cannot be layered to some degree, but observations within days of explosion are necessary to determine the extent to which the outer ejecta could be layered. △ Less

Submitted 23 October, 2021; originally announced October 2021.

Comments: 19 pages, 12 figures, 8 tables. Accepted for publication in MNRAS

Type Iax supernovae from deflagrations in Chandrasekhar mass white dwarfs

Authors: F. Lach, F. P. Callan, D. Bubeck, F. K. Roepke, S. A. Sim, M. Schrauth, S. T. Ohlmann, M. Kromer

Abstract: 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… ▽ More 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. △ Less

Submitted 3 December, 2021; v1 submitted 7 September, 2021; originally announced September 2021.

Comments: 27 pages, 17 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 658, A179 (2022)

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

Authors: S. Huber, S. H. Suyu, D. Ghoshdastidar, S. Taubenberger, V. Bonvin, J. H. H. Chan, M. Kromer, U. M. Noebauer, S. A. Sim, L. Leal-Taixé

Abstract: 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… ▽ More 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. △ Less

Submitted 21 December, 2021; v1 submitted 5 August, 2021; originally announced August 2021.

Comments: 25 pages, 28 figures; accepted for publication in A&A

Journal ref: A&A 658, A157 (2022)

Metallicity-dependent nucleosynthetic yields of Type Ia supernovae originating from double detonations of sub-M$_{\text{Ch}}$ white dwarfs

Authors: Sabrina Gronow, Benoit Cote, Florian Lach, Ivo R. Seitenzahl, Christine E. Collins, Stuart A. Sim, Friedrich K. Roepke

Abstract: Double detonations in sub-Chandrasekhar mass carbon-oxygen white dwarfs with helium shell are a potential explosion mechanism for a Type Ia supernova (SNe Ia). It comprises a shell detonation and subsequent core detonation. The focus of our study is on the effect of the progenitor metallicity on the nucleosynthetic yields. For this, we compute and analyse a set of eleven different models with vary… ▽ More Double detonations in sub-Chandrasekhar mass carbon-oxygen white dwarfs with helium shell are a potential explosion mechanism for a Type Ia supernova (SNe Ia). It comprises a shell detonation and subsequent core detonation. The focus of our study is on the effect of the progenitor metallicity on the nucleosynthetic yields. For this, we compute and analyse a set of eleven different models with varying core and shell masses at four different metallicities each. This results in a total of 44 models at metallicities between 0.01$Z_\odot$ and 3$Z_\odot$. Our models show a strong impact of the metallicity in the high density regime. The presence of $^{22}$Ne causes a neutron-excess which shifts the production from $^{56}$Ni to stable isotopes such as $^{54}$Fe and $^{58}$Ni in the $α$-rich freeze-out regime. The isotopes of the metallicity implementation further serve as seed nuclei for additional reactions in the shell detonation. Most significantly, the production of $^{55}$Mn increases with metallicity confirming the results of previous work. A comparison of elemental ratios relative to iron shows a relatively good match to solar values for some models. Super-solar values are reached for Mn at 3$Z_\odot$ and solar values in some models at $Z_\odot$. This indicates that the required contribution of SNe Ia originating from Chandrasekhar mass WDs can be lower than estimated in orevious work to reach solar values of [Mn/Fe] at [Fe/H]$=0$. Our galactic chemical evolution models suggest that SNe Ia from sub-Chandrasekhar mass white dwarfs, along with core-collapse supernovae, could account for more than 80% of the solar Mn abundance. Using metallicity-dependent SN Ia yields helps to reproduce the upward trend of [Mn/Fe] as a function of metallicity for the solar neighborhood. These chemical evolution predictions, however, depend on the massive star yields adopted in the calculations. △ Less

Submitted 1 September, 2021; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: accepted for publication by A&A

Journal ref: A&A 656, A94 (2021)

Double detonations of sub-Chandrasekhar mass CO white dwarfs: Can different core and He shell masses explain variations of Type Ia supernovae?

Authors: Sabrina Gronow, Christine E. Collins, Stuart A. Sim, Friedrich K. Roepke

Abstract: Sub-Chandrasekhar mass carbon-oxygen white dwarfs (CO WDs) with a surface helium (He) shell have been proposed as progenitors of Type Ia supernovae (SNe Ia). If true, the resulting thermonuclear explosions should be able to account for at least some of the range of SNe Ia observables. To study this, we conduct a parameter study based on 3D simulations of double detonations in CO WDs with a He shel… ▽ More Sub-Chandrasekhar mass carbon-oxygen white dwarfs (CO WDs) with a surface helium (He) shell have been proposed as progenitors of Type Ia supernovae (SNe Ia). If true, the resulting thermonuclear explosions should be able to account for at least some of the range of SNe Ia observables. To study this, we conduct a parameter study based on 3D simulations of double detonations in CO WDs with a He shell, assuming different core and shell masses. An admixture of C to the shell and solar metallicity are included in the models. The hydrodynamic simulations are carried out using the AREPO code. This allows us to follow the He shell detonation with high numerical resolution, and improves the reliability of predicted nucleosynthetic shell detonation yields. The addition of C to the shell leads to a lower production of 56Ni while including solar metallicity increases the production of IMEs. The production of higher mass elements is further shifted to stable isotopes at solar metallicity. Moreover, we find different core detonation ignition mechanisms depending on the core and shell mass configuration. This has an influence on the ejecta structure. We present the bolometric light curves predicted from our explosion simulations using the radiative transfer code ARTIS, and make comparisons with SNe Ia data. The bolometric light curves of our models show a range of brightnesses, able to account for sub-luminous to normal SNe Ia. We show the model bolometric width-luminosity relation compared to data for a range of viewing angles. We find that, on average, our brighter models lie within the observed data. The ejecta asymmetries produce a wide distribution of observables, which might account for outliers in the data. However, the models overestimate the extent of this compared to data. We also find the bolometric decline rate over 40 days appears systematically faster than data. (abridged) △ Less

Submitted 12 February, 2021; originally announced February 2021.

Comments: submitted to A&A, second revision after referee's comments

Journal ref: A&A 649, A155 (2021)

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

Authors: M. R. Magee, K. Maguire, R. Kotak, S. A. Sim

Abstract: 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… ▽ More 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. △ Less

Submitted 24 January, 2021; originally announced January 2021.

Comments: 21 pages, 25 figures, 1 table. Accepted for publication in MNRAS. Model light curves available at https://github.com/MarkMageeAstro/TURTLS-Light-curves

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

Authors: S. Huber, S. H. Suyu, U. M. Noebauer, J. H. H. Chan, M. Kromer, S. A. Sim, D. Sluse, S. Taubenberger

Abstract: 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 w… ▽ More 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$. △ Less

Submitted 20 January, 2021; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 17 pages, 15 figures; accepted for publication in A&A

Journal ref: A&A 646, A110 (2021)

AT2018kzr: the merger of an oxygen-neon white dwarf and a neutron star or black hole

Authors: James H. Gillanders, Stuart A. Sim, Stephen J. Smartt

Abstract: We present detailed spectroscopic analysis of the extraordinarily fast-evolving transient AT2018kzr. The transient's observed lightcurve 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 photosp… ▽ More We present detailed spectroscopic analysis of the extraordinarily fast-evolving transient AT2018kzr. The transient's observed lightcurve 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 $\sim$12000-14500km/s. The early spectra have the unusual combination of being blue but dominated by strong FeII and FeIII absorption features. We show that this combination is only possible with a high Fe content (3.5%). This implies a high Fe/(Ni+Co) ratio. Given the short time from the transient's proposed explosion epoch, the Fe cannot be $^{56}$Fe resulting from the decay of radioactive $^{56}$Ni synthesised in the explosion. Instead, we propose that this is stable $^{54}$Fe, and that the transient is unusually rich in this isotope. We further identify an additional, high-velocity component of ejecta material at $\sim$20000-26000km/s, which is mildly asymmetric and detectable through the CaII NIR 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. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: 17 pages, 11 figures, 5 tables

Accretion Disc Winds in Tidal Disruption Events: Ultraviolet Spectral Lines as Orientation Indicators

Authors: Edward J. Parkinson, Christian Knigge, Knox S. Long, James H. Matthews, Nick Higginbottom, Stuart A. Sim, Henrietta A. Hewitt

Abstract: Some tidal disruption events (TDEs) exhibit blueshifted broad absorption lines (BALs) in their rest-frame ultraviolet (UV) spectra, while others display broad emission lines (BELs). Similar phenomenology is observed in quasars and accreting white dwarfs, where it can be interpreted as an orientation effect associated with line formation in an accretion disc wind.We propose and explore a similar un… ▽ More Some tidal disruption events (TDEs) exhibit blueshifted broad absorption lines (BALs) in their rest-frame ultraviolet (UV) spectra, while others display broad emission lines (BELs). Similar phenomenology is observed in quasars and accreting white dwarfs, where it can be interpreted as an orientation effect associated with line formation in an accretion disc wind.We propose and explore a similar unification scheme for TDEs. We present synthetic UV spectra for disc and wind-hosting TDEs, produced by a state-of-the-art Monte Carlo ionization and radiative transfer code. Our models cover a wide range of disc wind geometries and kinematics. Such winds naturally reproduce both BALs and BELs. In general, sight lines looking into the wind cone preferentially produce BALs, while other orientations preferentially produce BELs. We also study the effect of wind clumping and CNO-processed abundances on the observed spectra. Clumpy winds tend to produce stronger UV emission and absorption lines, because clumping increases both the emission measure and the abundances of the relevant ionic species, the latter by reducing the ionization state of the outflow. The main effect of adopting CNO-processed abundances is a weakening of C~{\sc iv}~1550~Å~ and an enhancement of N \textsc{v}~1240~Å~ in the spectra. We conclude that line formation in an accretion disc wind is a promising mechanism for explaining the diverse UV spectra of TDEs. If this is correct, the relative number of BAL and BEL TDEs can be used to estimate the covering factor of the outflow. The models in this work are publicly available online and upon request. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: 16 pages, 9 figures

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

Authors: M. Bulla, Z. -W. Liu, F. K. Roepke, S. A. Sim, M. Fink, M. Kromer, R. Pakmor, I. R. Seitenzahl

Abstract: 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 depend… ▽ More 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. △ Less

Submitted 27 February, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 7 pages, 3 figures; accepted for publication in A&A; metadata update

Report number: NORDITA 2020-019

Journal ref: A&A 635, A179 (2020)

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

Authors: S. H. Suyu, S. Huber, R. Cañameras, M. Kromer, S. Schuldt, S. Taubenberger, A. Yıldırım, V. Bonvin, J. H. H. Chan, F. Courbin, U. Nöbauer, S. A. Sim, D. Sluse

Abstract: 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… ▽ More 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. △ Less

Submitted 11 November, 2020; v1 submitted 19 February, 2020; originally announced February 2020.

Comments: 15 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 644, A162 (2020)

SNe Ia from double detonations: Impact of core-shell mixing on the carbon ignition mechanism

Authors: Sabrina Gronow, Christine Collins, Sebastian T. Ohlmann, Ruediger Pakmor, Markus Kromer, Ivo R. Seitenzahl, Stuart A. Sim, Friedrich K. Roepke

Abstract: Sub-Chandrasekhar mass white dwarfs accreting a helium shell on a carbon-oxygen core are potential progenitors of normal Type Ia supernovae. This work focuses on the details of the onset of the carbon detonation in the double detonation sub-Chandrasekhar model. In order to simulate the influence of core-shell mixing on the carbon ignition mechanism, the helium shell and its detonation are followed… ▽ More Sub-Chandrasekhar mass white dwarfs accreting a helium shell on a carbon-oxygen core are potential progenitors of normal Type Ia supernovae. This work focuses on the details of the onset of the carbon detonation in the double detonation sub-Chandrasekhar model. In order to simulate the influence of core-shell mixing on the carbon ignition mechanism, the helium shell and its detonation are followed with an increased resolution compared to the rest of the star treating the propagation of the detonation wave more accurately. This significantly improves the predictions of the nucleosynthetic yields from the helium burning. The simulations were carried out with the AREPO code. A carbon-oxygen core with a helium shell was set up in one dimension and mapped to three dimensions. We ensured the stability of the white dwarf with a relaxation step before the hydrodynamic detonation simulation started. Synthetic observables were calculated with the radiative transfer code ARTIS. An ignition mechanism of the carbon detonation was observed, which received little attention before. In this "scissors mechanism", the impact the helium detonation wave has on unburnt material when converging opposite to its ignition spot is strong enough to ignite a carbon detonation. This is possible in a carbon enriched transition region between the core and shell. The detonation mechanism is found to be sensitive to details of the core-shell transition and our models illustrate the need to consider core-shell mixing taking place during the accretion process. Even though the detonation ignition mechanism differs form the converging shock mechanism, the differences in the synthetic observables are not significant. Though they do not fit observations better than previous simulations, they illustrate the need for multi-dimensional simulations. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 635, A169 (2020)

The lowest of the low: discovery of SN 2019gsc and the nature of faint Iax supernovae

Authors: Shubham Srivastav, Stephen J. Smartt, Giorgos Leloudas, Mark E. Huber, Ken Chambers, Daniele B. Malesani, Jens Hjorth, James H. Gillanders, A. Schultz, Stuart A. Sim, Katie Auchettl, Johan P. U. Fynbo, Christa Gall, Owen R. McBrien, Armin Rest, Ken W. Smith, Radoslaw Wojtak, David R. Young

Abstract: We present the discovery and optical follow-up of the faintest supernova-like transient known. The event (SN 2019gsc) was discovered in a star-forming host at 53\,Mpc by ATLAS. A detailed multi-colour light curve was gathered with Pan-STARRS1 and follow-up spectroscopy was obtained with the NOT and Gemini-North. The spectra near maximum light show narrow features at low velocities of 3000 to 4000… ▽ More We present the discovery and optical follow-up of the faintest supernova-like transient known. The event (SN 2019gsc) was discovered in a star-forming host at 53\,Mpc by ATLAS. A detailed multi-colour light curve was gathered with Pan-STARRS1 and follow-up spectroscopy was obtained with the NOT and Gemini-North. The spectra near maximum light show narrow features at low velocities of 3000 to 4000 km s$^{-1}$, similar to the extremely low luminosity SNe 2010ae and 2008ha, and the light curve displays a similar fast decline (\dmr $0.91 \pm 0.10$ mag). SNe 2010ae and 2008ha have been classified as type Iax supernovae, and together the three either make up a distinct physical class of their own or are at the extreme low luminosity end of this diverse supernova population. The bolometric light curve is consistent with a low kinetic energy of explosion ($E_{\rm k} \sim 10^{49}$ erg s$^{-1}$), a modest ejected mass ($M_{\rm ej} \sim 0.2$ \msol) and radioactive powering by $^{56}$Ni ($M_{\rm Ni} \sim 2 \times 10^{-3}$ \msol). The spectra are quite well reproduced with radiative transfer models (TARDIS) and a composition dominated by carbon, oxygen, magnesium, silicon and sulphur. Remarkably, all three of these extreme Iax events are in similar low-metallicity star-forming environments. The combination of the observational constraints for all three may be best explained by deflagrations of near $M_{\rm Ch}$ hybrid carbon-oxygen-neon white dwarfs which have short evolutionary pathways to formation. △ Less

Submitted 18 February, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: 12 pages, 5 figures, accepted to ApJL, minor changes to submitted version

Thermal and radiation driving can produce observable disk winds in hard-state X-ray binaries

Authors: Nick Higginbottom, Christian Knigge, Stuart A. Sim, Knox S. Long, James H. Matthews, Henrietta A. Hewitt, Edward J. Parkinson, Sam W. Mangham

Abstract: X-ray signatures of outflowing gas have been detected in several accreting black-hole binaries, always in the soft state. A key question raised by these observations is whether these winds might also exist in the hard state. Here, we carry out the first full-frequency radiation hydrodynamic simulations of luminous ($\rm{L = 0.5 \, L_{\mathrm{Edd}}}$) black-hole X-ray binary systems in both the har… ▽ More X-ray signatures of outflowing gas have been detected in several accreting black-hole binaries, always in the soft state. A key question raised by these observations is whether these winds might also exist in the hard state. Here, we carry out the first full-frequency radiation hydrodynamic simulations of luminous ($\rm{L = 0.5 \, L_{\mathrm{Edd}}}$) black-hole X-ray binary systems in both the hard and the soft state, with realistic spectral energy distributions (SEDs). Our simulations are designed to describe X-ray transients near the peak of their outburst, just before and after the hard-to-soft state transition. At these luminosities, it is essential to include radiation driving, and we include not only electron scattering, but also photoelectric and line interactions. We find powerful outflows with $\rm{\dot{M}_{wind} \simeq 2 \,\dot{M}_{acc}}$ are driven by thermal and radiation pressure in both hard and soft states. The hard-state wind is significantly faster and carries approximately 20 times as much kinetic energy as the soft-state wind. However, in the hard state the wind is more ionized, and so weaker X-ray absorption lines are seen over a narrower range of viewing angles. Nevertheless, for inclinations $\gtrsim 80^{\circ}$, blue-shifted wind-formed Fe XXV and Fe XXVI features should be observable even in the hard state. Given that the data required to detect these lines currently exist for only a single system in a {\em luminous} hard state -- the peculiar GRS~1915+105 -- we urge the acquisition of new observations to test this prediction. The new generation of X-ray spectrometers should be able to resolve the velocity structure. △ Less

Submitted 27 January, 2020; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: Accepted for publication by MNRAS

Stratified disc wind models for the AGN broad-line region: ultraviolet, optical and X-ray properties

Authors: James H. Matthews, Christian Knigge, Nick Higginbottom, Knox S. Long, Stuart A. Sim, Samuel W. Mangham, Edward J. Parkinson, Henrietta A. Hewitt

Abstract: The origin, geometry and kinematics of the broad line region (BLR) gas in quasars and active galactic nuclei (AGN) are uncertain. We demonstrate that clumpy biconical disc winds illuminated by an AGN continuum can produce BLR-like spectra. We first use a simple toy model to illustrate that disc winds make quite good BLR candidates, because they are self-shielded flows and can cover a large portion… ▽ More The origin, geometry and kinematics of the broad line region (BLR) gas in quasars and active galactic nuclei (AGN) are uncertain. We demonstrate that clumpy biconical disc winds illuminated by an AGN continuum can produce BLR-like spectra. We first use a simple toy model to illustrate that disc winds make quite good BLR candidates, because they are self-shielded flows and can cover a large portion of the ionizing flux-density ($φ_H$-$n_H$) plane. We then conduct Monte Carlo radiative transfer and photoionization calculations, which fully account for self-shielding and multiple scattering in a non-spherical geometry. The emergent model spectra show broad emission lines with equivalent widths and line ratios comparable to those observed in AGN, provided that the wind has a volume filling factor of $f_V\lesssim0.1$. Similar emission line spectra are produced for a variety of wind geometries (polar or equatorial) and for launch radii that differ by an order of magnitude. The line emission arises almost exclusively from plasma travelling below the escape velocity, implying that `failed winds' are important BLR candidates. The behaviour of a line-emitting wind (and possibly any `smooth flow' BLR model) is similar to that of the locally optimally-emitting cloud (LOC) model originally proposed by Baldwin et al (1995), except that the gradients in ionization state and temperature are large-scale and continuous, rather than within or between distinct clouds. Our models also produce UV absorption lines and X-ray absorption features, and the stratified ionization structure can partially explain the different classes of broad absorption line quasars. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: Accepted to MNRAS. 20 pages, 13 figures plus appendix. Supplementary plots, some simulation data and demo notebooks can be found at https://github.com/jhmatthews/windy-blr-2020

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

Authors: M. R. Magee, K. Maguire, R. Kotak, S. A. Sim, J. H. Gillanders, S. J. Prentice, K. Skillen

Abstract: 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 mod… ▽ More 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. △ Less

Submitted 16 December, 2019; originally announced December 2019.

Comments: 31 pages, 32 figures, 2 tables. Accepted for publication in A&A. Model light curves available at https://github.com/MarkMageeAstro/TURTLS-Light-curves

Journal ref: A&A 634, A37 (2020)

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

Authors: L. J. Shingles, S. A. Sim, M. Kromer, K. Maguire, M. Bulla, C. Collins, C. P. Ballance, A. S. Michel, C. A. Ramsbottom, F. K. Roepke, I. R. Seitenzahl, N. B. Tyndall

Abstract: 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-the… ▽ More 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. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 15 pages, 12 figures. Accepted for publication by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 492, Issue 2, February 2020, Pages 2029-2043

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

Authors: C. Vogl, W. E. Kerzendorf, S. A. Sim, U. M. Noebauer, S. Lietzau, W. Hillebrandt

Abstract: 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… ▽ More 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. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Comments: 18 pages, 13 figures, 3 tables, submitted to A&A

Journal ref: A&A 633, A88 (2020)

SN2018kzr: a rapidly declining transient from the destruction of a white dwarf

Authors: Owen R. McBrien, Stephen J. Smartt, Ting-Wan Chen, Cosimo Inserra, James H. Gillanders, Stuart A. Sim, Anders Jerkstrand, Armin Rest, Stefano Valenti, Rupak Roy, Mariusz Gromadzki, Stefan Taubenberger, Andreas Flörs, Mark E. Huber, Ken C. Chambers, Avishay Gal-Yam, David R. Young, Matt Nicholl, Erkki Kankare, Ken W. Smith, Kate Maguire, Ilya Mandel, Simon Prentice, Ósmar Rodríguez, Jonathon Pineda Garcia , et al. (9 additional authors not shown)

Abstract: We present SN2018kzr, the fastest declining supernova-like transient, second only to the kilonova, AT2017gfo. SN2018kzr is characterized by a peak magnitude of $M_r = -17.98$, peak bolometric luminosity of ${\sim} 1.4 \times 10^{43}$erg s$^{\mathrm{-1}}$ and a rapid decline rate of $0.48 \pm 0.03$ mag day$^{\textrm{-1}}$ in the $r$ band. The bolometric luminosity evolves too quickly to be explaine… ▽ More We present SN2018kzr, the fastest declining supernova-like transient, second only to the kilonova, AT2017gfo. SN2018kzr is characterized by a peak magnitude of $M_r = -17.98$, peak bolometric luminosity of ${\sim} 1.4 \times 10^{43}$erg s$^{\mathrm{-1}}$ and a rapid decline rate of $0.48 \pm 0.03$ mag day$^{\textrm{-1}}$ in the $r$ band. The bolometric luminosity evolves too quickly to be explained by pure $^{\mathrm{56}}$Ni heating, necessitating the inclusion of an alternative powering source. Incorporating the spin-down of a magnetized neutron star adequately describes the lightcurve and we estimate a small ejecta mass of $M_\mathrm{ej} = 0.10 \pm 0.05$ $\textrm{M}_{\odot}$. Our spectral modelling suggests the ejecta is composed of intermediate mass elements including O, Si and Mg and trace amounts of Fe-peak elements, which disfavours a binary neutron star merger. We discuss three explosion scenarios for SN2018kzr, given the low ejecta mass, intermediate mass element composition and the high likelihood of additional powering - core collapse of an ultra-stripped progenitor, the accretion induced collapse of a white dwarf and the merger of a white dwarf and neutron star. The requirement for an alternative input energy source favours either the accretion induced collapse with magnetar powering or a white dwarf - neutron star merger with energy from disk wind shocks. △ Less

Submitted 19 November, 2019; v1 submitted 10 September, 2019; originally announced September 2019.

Comments: 11 pages, 5 figures, 3 tables Accepted by ApJL on 2019 October 15

Monte Carlo Radiative Transfer

Authors: Ulrich M. Noebauer, Stuart A. Sim

Abstract: The theory and numerical modelling of radiation processes and radiative transfer play a key role in astrophysics: they provide the link between the physical properties of an object and the radiation it emits. In the modern era of increasingly high-quality observational data and sophisticated physical theories, development and exploitation of a variety of approaches to the modelling of radiative tr… ▽ More The theory and numerical modelling of radiation processes and radiative transfer play a key role in astrophysics: they provide the link between the physical properties of an object and the radiation it emits. In the modern era of increasingly high-quality observational data and sophisticated physical theories, development and exploitation of a variety of approaches to the modelling of radiative transfer is needed. In this article, we focus on one remarkably versatile approach: Monte Carlo Radiative Transfer (MCRT). We describe the principles behind this approach, and highlight the relative ease with which they can (and have) been implemented for application to a range of astrophysical problems. All MCRT methods have in common a need to consider the adverse consequences of Monte Carlo noise in simulation results. We overview a range of methods used to suppress this noise and comment on their relative merits for a variety of applications. We conclude with a brief review of specific applications for which MCRT methods are currently popular and comment on the prospects for future developments. △ Less

Submitted 23 July, 2019; originally announced July 2019.

Comments: 113 pages, 26 figures, invited review for Living Reviews in Computational Astrophysics

Journal ref: Living Rev Comput Astrophys (2019) 5: 1

Do Reverberation Mapping Analyses Provide an Accurate Picture of the Broad Line Region?

Authors: S. W. Mangham, C. Knigge, P. Williams, Keith Horne, A. Pancoast, J. H. Matthews, K. S. Long, S. A. Sim, N. Higginbottom

Abstract: Reverberation mapping (RM) is a powerful approach for determining the nature of the broad-line region (BLR) in active galactic nuclei. However, inferring physical BLR properties from an observed spectroscopic time series is a difficult inverse problem. Here, we present a blind test of two widely used RM methods: MEMEcho (developed by Horne) and CARAMEL (developed by Pancoast and collaborators). Th… ▽ More Reverberation mapping (RM) is a powerful approach for determining the nature of the broad-line region (BLR) in active galactic nuclei. However, inferring physical BLR properties from an observed spectroscopic time series is a difficult inverse problem. Here, we present a blind test of two widely used RM methods: MEMEcho (developed by Horne) and CARAMEL (developed by Pancoast and collaborators). The test data are simulated spectroscopic time series that track the H$α$ emission line response to an empirical continuum light curve. The underlying BLR model is a rotating, biconical accretion disc wind, and the synthetic spectra are generated via self-consistent ionization and radiative transfer simulations. We generate two mock data sets, representing Seyfert galaxies and QSOs. The Seyfert model produces a largely *negative* response, which neither method can recover. However, both fail $``gracefully''$, neither generating spurious results. For the QSO model both CARAMEL and expert interpretation of MEMEcho's output both capture the broadly annular, rotation-dominated nature of the line-forming region, though MEMEcho analysis overestimates its size by 50%, but CARAMEL is unable to distinguish between additional inflow and outflow components. Despite fitting individual spectra well, the CARAMEL velocity-delay maps and RMS line profiles are strongly inconsistent with the input data. Finally, since the H$α$ line-forming region is rotation dominated, neither method recovers the disc wind nature of the underlying BLR model. Thus considerable care is required when interpreting the results of RM analyses in terms of physical models. △ Less

Submitted 26 June, 2019; originally announced June 2019.

Comments: 22 pages, 24 figures

Report number: MN-19-0268-MJ.R1

Spectral sequences of Type Ia supernovae. II. Carbon as a diagnostic tool for explosion mechanisms

Authors: E. Heringer, M. H. van Kerkwijk, S. A. Sim, W. E. Kerzendorf, Melissa L. Graham

Abstract: How an otherwise inert carbon-oxygen white dwarf can be made to explode as a Type Ia supernova remains unknown. A promising test of theoretical models is to constrain the distribution of material that is left unburned, in particular of carbon. So far, most investigations used line identification codes to detect carbon in the ejecta, a method that cannot be readily compared against model prediction… ▽ More How an otherwise inert carbon-oxygen white dwarf can be made to explode as a Type Ia supernova remains unknown. A promising test of theoretical models is to constrain the distribution of material that is left unburned, in particular of carbon. So far, most investigations used line identification codes to detect carbon in the ejecta, a method that cannot be readily compared against model predictions because it requires assumed opacities and temperatures. Here, we instead use tomographic techniques to investigate the amount of carbon in the inner layers of SN~2011fe, starting from the previously published tomographic analysis of Mazzali et al. From the presence of the carbon feature in the optical at early epochs and its disappearance later on, we derive an average carbon mass fraction between 0.001 and 0.05 for velocities in the range $13500 \lesssim v \lesssim 16000\ \rm{km\ s^{-1}}$, and an upper limit of 0.005 inside that region. Based on our models and the assumed density profile, only small amounts of carbon should be in the neutral state, too little to be responsible for features seen in near-infrared spectra that were previously identified as due to neutral carbon; We discuss possible reasons for this discrepancy. We compare our results against a suite of explosion models, although uncertainties in both the models and our simulations make it difficult to draw definitive conclusions. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Journal ref: The Astrophysical Journal, Volume 871, Number 2 (2019, Feb. 5)

Detecting the signatures of helium in type Iax supernovae

Authors: M. R. Magee, S. A. Sim, R. Kotak, K. Maguire, A. Boyle

Abstract: Recent studies have argued that the progenitor system of type Iax supernovae must consist of a carbon-oxygen white dwarf accreting from a helium star companion. Based on existing explosion models invoking the pure deflagration of carbon-oxygen white dwarfs, we investigate the likelihood of producing spectral features due to helium in type Iax supernovae. From this scenario, we select those explosi… ▽ More Recent studies have argued that the progenitor system of type Iax supernovae must consist of a carbon-oxygen white dwarf accreting from a helium star companion. Based on existing explosion models invoking the pure deflagration of carbon-oxygen white dwarfs, we investigate the likelihood of producing spectral features due to helium in type Iax supernovae. From this scenario, we select those explosion models producing ejecta and $^{56}$Ni masses that are broadly consistent with those estimated for type Iax supernovae (0.014 - 0.478~$M_{\odot}$ and $\sim0.003$ - 0.183~$M_{\odot}$, respectively). To this end, we present a series of models of varying luminosities ($-18.4 \lesssim M_{\rm{V}} \lesssim -14.5$~mag) with helium abundances accounting for up to $\sim$36\% of the ejecta mass, and covering a range of epochs beginning a few days before B$-$band maximum to approximately two weeks after maximum. We find that the best opportunity for detecting \ion{He}{i} features is at near-infrared wavelengths, and in the post-maximum spectra of the fainter members of this class. We show that the optical spectrum of SN~2007J is potentially consistent with a large helium content (a few 10$^{-2}~M_{\odot}$), but argue that current models of accretion and material stripping from a companion struggle to produce compatible scenarios. We also investigate the presence of helium in all objects with near-infrared spectra. We show that SNe~2005hk, 2012Z, and 2015H contain either no helium or their helium abundances are constrained to much lower values ($\lesssim$10$^{-3}~M_{\odot}$). Our results demonstrate the differences in helium content among type Iax supernovae, perhaps pointing to different progenitor channels. Either SN~2007J is an outlier in terms of its progenitor system, or it is not a true member of the type Iax supernova class. △ Less

Submitted 20 December, 2018; originally announced December 2018.

Comments: 15 pages, 12 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 622, A102 (2019)

Progenitor candidate for the type II-P supernova SN 2018aoq in NGC 4151

Authors: D. O'Neill, R. Kotak, M. Fraser, S. A. Sim, S. Benetti, S. J. Smartt, S. Mattila, C. Ashall, E. Callis, N. Elias-Rosa, M. Gromadzki, S. J. Prentice

Abstract: We present our findings based on pre- and post-explosion data of the type II-Plateau SN 2018aoq that exploded in NGC 4151. As distance estimates to NGC 4151 vary by an order of magnitude, we utilised the well-known correlation between ejecta velocity and plateau brightness, i.e. the standard candle method, to obtain a distance of 18.2$\pm$1.2 Mpc, which is in very good agreement with measurements… ▽ More We present our findings based on pre- and post-explosion data of the type II-Plateau SN 2018aoq that exploded in NGC 4151. As distance estimates to NGC 4151 vary by an order of magnitude, we utilised the well-known correlation between ejecta velocity and plateau brightness, i.e. the standard candle method, to obtain a distance of 18.2$\pm$1.2 Mpc, which is in very good agreement with measurements based on geometric methods. The above distance implies a mid-plateau absolute magnitude of $M_{V}^{50}=-15.76\pm$0.14 suggesting that it is of intermediate brightness when compared to IIP SNe such as SN 2005cs at the faint end, and more typical events such as SN 1999em. This is further supported by relatively low expansion velocities (Fe II $λ$5169 $\sim$3000 km s$^{-1}$ at +42 d). Using archival HST/WFC3 imaging data, we find a point source coincident with the supernova position in the F350LP, F555W, F814W, and F160W filters. This source shows no significant variability over the $\sim$2 month time span of the data. From fits to the spectral energy distribution of the candidate progenitor, we find $\log\left(L/L_\odot\right)\sim 4.7$ and $T_{\mathrm{eff}}\sim 3.5$ kK, implying an M-type red supergiant progenitor. From comparisons to single and binary star models, we find that both favour the explosion of a star with a zero-age main sequence mass of $\sim$$10 M_\odot$. △ Less

Submitted 5 February, 2019; v1 submitted 12 December, 2018; originally announced December 2018.

Comments: 9 pages, 8 figures, accepted for publication in Astronomy & Astrophysics (11/12/18)

Journal ref: A&A 622, L1 (2019)

Spectral modeling of type II supernovae. I. Dilution factors

Authors: C. Vogl, S. A. Sim, U. M. Noebauer, W. E. Kerzendorf, W. Hillebrandt

Abstract: We present substantial extensions to the Monte Carlo radiative transfer code TARDIS to perform spectral synthesis for type II supernovae. By incorporating a non-LTE ionization and excitation treatment for hydrogen, a full account of free-free and bound-free processes, a self-consistent determination of the thermal state and by improving the handling of relativistic effects, the improved code versi… ▽ More We present substantial extensions to the Monte Carlo radiative transfer code TARDIS to perform spectral synthesis for type II supernovae. By incorporating a non-LTE ionization and excitation treatment for hydrogen, a full account of free-free and bound-free processes, a self-consistent determination of the thermal state and by improving the handling of relativistic effects, the improved code version includes the necessary physics to perform spectral synthesis for type II supernovae to high precision as required for the reliable inference of supernova properties. We demonstrate the capabilities of the extended version of TARDIS by calculating synthetic spectra for the prototypical type II supernova SN1999em and by deriving a new and independent set of dilution factors for the expanding photosphere method. We have investigated in detail the dependence of the dilution factors on photospheric properties and, for the first time, on changes in metallicity. We also compare our results with two previously published sets of dilution factors by Eastman et al. (1996) and by Dessart & Hillier (2005), and discuss the potential sources of the discrepancies between studies. △ Less

Submitted 6 November, 2018; originally announced November 2018.

Comments: 16 pages, 12 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A 621, A29 (2019)

Type Iax SNe as a few-parameter family

Authors: Barnabas Barna, Tamas Szalai, Wolfgang E. Kerzendorf, Markus Kromer, Stuart A. Sim, Mark R. Magee, Bruno Leibundgut

Abstract: We present direct spectroscopic modeling of five Type Iax supernovae (SNe) with the one dimensional Monte Carlo radiative transfer code TARDIS. The abundance tomography technique is used to map the chemical structure and physical properties of the SN atmosphere. Through via fitting of multiple spectral epochs with self-consistent ejecta models, we can then constrain the location of some elements w… ▽ More We present direct spectroscopic modeling of five Type Iax supernovae (SNe) with the one dimensional Monte Carlo radiative transfer code TARDIS. The abundance tomography technique is used to map the chemical structure and physical properties of the SN atmosphere. Through via fitting of multiple spectral epochs with self-consistent ejecta models, we can then constrain the location of some elements within the ejecta. The synthetic spectra of the best-fit models are able to reproduce the flux continuum and the main absorption features in the whole sample. We find that the mass fractions of IGEs and IMEs show a decreasing trend toward the outer regions of the atmospheres using density profiles similar to those of deflagration models in the literature. Oxygen is the only element, which could be dominant at higher velocities. The stratified abundance structure contradicts the well-mixed chemical profiles predicted by pure deflagration models. Based on the derived densities and abundances, a template model atmosphere is created for the SN Iax class and compared to the observed spectra. Free parameters are the scaling of the density profile, the velocity shift of the abundance template, and the peak luminosity. The results of this test support the idea that all SNe Iax can be described by a similar internal structure, which argues for a common origin of this class of explosions. △ Less

Submitted 1 August, 2018; originally announced August 2018.

Comments: 21 pages, 7 tables, 16 figures, accepted by MNRAS

Thermonuclear explosions of rapidly differentially rotating white dwarfs: Candidates for superluminous Type Ia supernovae?

Authors: M. Fink, M. Kromer, W. Hillebrandt, F. K. Roepke, R. Pakmor, I. R. Seitenzahl, S. A. Sim

Abstract: The observed sub-class of "superluminous" Type Ia supernovae lacks a convincing theoretical explanation. If the emission of such objects were powered exclusively by radioactive decay of 56Ni formed in the explosion, a progenitor mass close to or even above the Chandrasekhar limit for a non-rotating white dwarf star would be required. Masses significantly exceeding this limit can be supported by di… ▽ More The observed sub-class of "superluminous" Type Ia supernovae lacks a convincing theoretical explanation. If the emission of such objects were powered exclusively by radioactive decay of 56Ni formed in the explosion, a progenitor mass close to or even above the Chandrasekhar limit for a non-rotating white dwarf star would be required. Masses significantly exceeding this limit can be supported by differential rotation. We, therefore, explore explosions and predict observables for various scenarios resulting from differentially rotating carbon-oxygen white dwarfs close to their respective limit of stability. Specifically, we have investigated a prompt detonation model, detonations following an initial deflagration phase ("delayed detonation" models), and a pure deflagration model. In postprocessing steps, we performed nucleosynthesis and three-dimensional radiative transfer calculations, that allow us, for the first time, to consistently derive synthetic observables from our models. We find that all explosion scenarios involving detonations produce very bright events. The observables predicted for them, however, are inconsistent with any known subclass of Type Ia supernovae. Pure deflagrations resemble 2002cx-like supernovae and may contribute to this class. We discuss implications of our findings for the explosion mechanism and for the existence of differentially rotating white dwarfs as supernova progenitors. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Comments: 12 pages, 9 figures, 2 tables, accepted for publication in A&A. Model data are available from the Heidelberg Supernova Model Archive (HESMA) at https://hesma.h-its.org/

Journal ref: A&A 618, A124 (2018)