Showing 1–49 of 49 results for author: Kerzendorf, W E

Supernova Shocks Cannot Explain the Inflated State of Hypervelocity Runaways from White Dwarf Binaries

Authors: Aakash Bhat, Evan B. Bauer, Rüdiger Pakmor, Ken J. Shen, Ilaria Caiazzo, Abinaya Swaruba Rajamuthukumar, Kareem El-Badry, Wolfgang E. Kerzendorf

Abstract: Recent observations have found a growing number of hypervelocity stars with speeds of $\approx 1500-2500\,$km\,s$^{-1}$ which could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been… ▽ More Recent observations have found a growing number of hypervelocity stars with speeds of $\approx 1500-2500\,$km\,s$^{-1}$ which could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been as white dwarfs filling their Roche lobe. While many simulations exist studying the dynamical phase leading to supernova detonation in these systems, no detailed calculations of the long-term structure of the runaways have yet been performed. We use an existing \textsc{Arepo} hydrodynamical simulation of a supernova in a white dwarf binary as a starting point for the evolution of these stars with the 1 dimensional stellar evolution code MESA. We show that the supernova shock is not enough to inflate the white dwarf over timescales longer than a few thousand years, significantly shorter than the $10^{5-6}$ year lifetimes inferred for observed hypervelocity runaways. Despite experiencing a shock from a supernova less than $\approx 0.02\,R_\odot$ away, our models do not experience significant interior heating, and all contract back to radii around $0.01\,R_\odot$ within about $10^4$\,years. Explaining the observed inflated states requires either an additional source of significant heating or some other physics that is not yet accounted for in the subsequent evolution. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: Submitted to A\&A. 15 pages, 17 figures

Scalable stellar evolution forecasting: Deep learning emulation vs. hierarchical nearest neighbor interpolation

Authors: K. Maltsev, F. R. N. Schneider, F. K. Roepke, A. I. Jordan, G. A. Qadir, W. E. Kerzendorf, K. Riedmiller, P. van der Smagt

Abstract: Many astrophysical applications require efficient yet reliable forecasts of stellar evolution tracks. One example is population synthesis, which generates forward predictions of models for comparison with observations. The majority of state-of-the-art rapid population synthesis methods are based on analytic fitting formulae to stellar evolution tracks that are computationally cheap to sample stati… ▽ More Many astrophysical applications require efficient yet reliable forecasts of stellar evolution tracks. One example is population synthesis, which generates forward predictions of models for comparison with observations. The majority of state-of-the-art rapid population synthesis methods are based on analytic fitting formulae to stellar evolution tracks that are computationally cheap to sample statistically over a continuous parameter range. The computational costs of running detailed stellar evolution codes, such as MESA, over wide and densely sampled parameter grids are prohibitive, while stellar-age based interpolation in-between sparsely sampled grid points leads to intolerably large systematic prediction errors. In this work, we provide two solutions for automated interpolation methods that offer satisfactory trade-off points between cost-efficiency and accuracy. We construct a timescale-adapted evolutionary coordinate and use it in a two-step interpolation scheme that traces the evolution of stars from ZAMS all the way to the end of core helium burning while covering a mass range from ${0.65}$ to $300 \, \mathrm{M_\odot}$. The feedforward neural network regression model (first solution) that we train to predict stellar surface variables can make millions of predictions, sufficiently accurate over the entire parameter space, within tens of seconds on a 4-core CPU. The hierarchical nearest-neighbor interpolation algorithm (second solution) that we hard-code to the same end achieves even higher predictive accuracy, the same algorithm remains applicable to all stellar variables evolved over time, but it is two orders of magnitude slower. Our methodological framework is demonstrated to work on the MIST (Choi et al. 2016) data set. Finally, we discuss the prospective applications of these methods and provide guidelines for generalizing them to higher dimensional parameter spaces. △ Less

Submitted 27 October, 2023; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: Accepted at A&A

Journal ref: A&A 681, A86 (2024)

Ground-based and JWST Observations of SN 2022pul: II. Evidence from Nebular Spectroscopy for a Violent Merger in a Peculiar Type-Ia Supernova

Authors: Lindsey A. Kwok, Matthew R. Siebert, Joel Johansson, Saurabh W. Jha, Stephane Blondin, Luc Dessart, Ryan J. Foley, D. John Hillier, Conor Larison, Ruediger Pakmor, Tea Temim, Jennifer E. Andrews, Katie Auchettl, Carles Badenes, Barnabas Barna, K. Azalee Bostroem, Max J. Brenner Newman, Thomas G. Brink, Maria Jose Bustamante-Rosell, Yssavo Camacho-Neves, Alejandro Clocchiatti, David A. Coulter, Kyle W. Davis, Maxime Deckers, Georgios Dimitriadis , et al. (56 additional authors not shown)

Abstract: We present an analysis of ground-based and JWST observations of SN~2022pul, a peculiar "03fg-like" (or "super-Chandrasekhar") Type Ia supernova (SN Ia), in the nebular phase at 338d post explosion. Our combined spectrum continuously covers 0.4--14 $μ$m and includes the first mid-infrared spectrum of an 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization… ▽ More We present an analysis of ground-based and JWST observations of SN~2022pul, a peculiar "03fg-like" (or "super-Chandrasekhar") Type Ia supernova (SN Ia), in the nebular phase at 338d post explosion. Our combined spectrum continuously covers 0.4--14 $μ$m and includes the first mid-infrared spectrum of an 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. Strong, broad, centrally peaked [Ne II] line at 12.81 $μ$m was previously predicted as a hallmark of "violent merger'' SN Ia models, where dynamical interaction between two sub-$M_{ch}$ white dwarfs (WDs) causes disruption of the lower mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to better reproduce the optical iron emission, and add mass in the innermost region ($< 2000$ km s$^{-1}$) to account for the observed narrow [O I]~$λ\lambda6300$, 6364 emission. A violent WD-WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SN Ia. △ Less

Submitted 23 May, 2024; v1 submitted 23 August, 2023; originally announced August 2023.

Comments: 20 pages, 10 figures, published in ApJ

Journal ref: The Astrophysical Journal, Volume 966, Issue 1, id.135, 18 pp., May 2024

Ground-based and JWST Observations of SN 2022pul: I. Unusual Signatures of Carbon, Oxygen, and Circumstellar Interaction in a Peculiar Type Ia Supernova

Authors: Matthew R. Siebert, Lindsey A. Kwok, Joel Johansson, Saurabh W. Jha, Stéphane Blondin, Luc Dessart, Ryan J. Foley, D. John Hillier, Conor Larison, Rüdiger Pakmor, Tea Temim, Jennifer E. Andrews, Katie Auchettl, Carles Badenes, Barnabas Barna, K. Azalee Bostroem, Max J. Brenner Newman, Thomas G. Brink, María José Bustamante-Rosell, Yssavo Camacho-Neves, Alejandro Clocchiatti, David A. Coulter, Kyle W. Davis, Maxime Deckers, Georgios Dimitriadis , et al. (57 additional authors not shown)

Abstract: Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" S… ▽ More Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon-oxygen rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity ($M_{B}=-18.9$ mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak $B$-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O I] $λλ6300,\ 6364$ (${\rm FWHM} \approx 2{,}000$ km s$^{-1}$), strong, broad emission from [Ca II] $λλ7291,\ 7323$ (${\rm FWHM} \approx 7{,}300$ km s$^{-1}$), and a rapid Fe III to Fe II ionization change. Finally, we present the first-ever optical-to-mid-infrared (MIR) nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with $T \approx 500$ K), combined with prominent [O I] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within carbon/oxygen-rich CSM. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: 23 pages, 11 figures, submitted to ApJ

No plateau observed in late-time near-infrared observations of the underluminous Type Ia supernova 2021qvv

Authors: O. Graur, E. Padilla Gonzalez, J. Burke, M. Deckers, S. W. Jha, L. Galbany, E. Karamenhmetoglu, M. D. Stritzinger, K. Maguire, D. A. Howell, R. Fisher, A. G. Fullard, R. Handberg, D. Hiramatsu, G. Hosseinzadeh, W. E. Kerzendorf, C. McCully, M. Newsome, C. Pellegrino, A. Rest, A. G. Riess, I. R. Seitenzahl, M. M. Shara, K. J. Shen, G. Terreran , et al. (1 additional authors not shown)

Abstract: Near-infrared (NIR) observations of normal Type Ia supernovae (SNe Ia) obtained between 150 to 500 d past maximum light reveal the existence of an extended plateau. Here, we present observations of the underluminous, 1991bg-like SN 2021qvv. Early, ground-based optical and NIR observations show that SN 2021qvv is similar to SN 2006mr, making it one of the dimmest, fastest-evolving 1991bg-like SNe t… ▽ More Near-infrared (NIR) observations of normal Type Ia supernovae (SNe Ia) obtained between 150 to 500 d past maximum light reveal the existence of an extended plateau. Here, we present observations of the underluminous, 1991bg-like SN 2021qvv. Early, ground-based optical and NIR observations show that SN 2021qvv is similar to SN 2006mr, making it one of the dimmest, fastest-evolving 1991bg-like SNe to date. Late-time (170-250 d) Hubble Space Telescope observations of SN 2021qvv reveal no sign of a plateau. An extrapolation of these observations backwards to earlier-phase NIR observations of SN 2006mr suggests the complete absence of a NIR plateau, at least out to 250 d. This absence may be due to a higher ionization state of the ejecta, as predicted by certain sub-Chandrasekhar-mass detonation models, or to the lower temperatures of the ejecta of 1991bg-like SNe, relative to normal SNe Ia, which might preclude their becoming fluorescent and shifting ultraviolet light into the NIR. This suggestion can be tested by acquiring NIR imaging of a sample of 1991bg-like SNe that covers the entire range from slowly-evolving to fast-evolving events ($0.2 \lesssim s_\mathrm{BV} \lesssim 0.6$). A detection of the NIR plateau in slower-evolving, hotter 1991bg-like SNe would provide further evidence that these SNe exist along a continuum with normal SNe Ia. Theoretical progenitor and explosion scenarios would then have to match the observed properties of both SN Ia subtypes. △ Less

Submitted 26 September, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS (14 pages, 9 figures)

1991T-Like Type Ia Supernovae as an Extension of the Normal Population

Authors: John T. O'Brien, Wolfgang E. Kerzendorf, Andrew Fullard, Reudiger Pakmor, Johannes Buchner, Christian Vogl, Nutan Chen, Patrick van der Smagt, Marc Williamson, Jaladh Singhal

Abstract: Type Ia supernovae remain poorly understood despite decades of investigation. Massive computationally intensive hydrodynamic simulations have been developed and run to model an ever-growing number of proposed progenitor channels. Further complicating the matter, a large number of sub-types of Type Ia supernovae have been identified in recent decades. Due to the massive computational load required,… ▽ More Type Ia supernovae remain poorly understood despite decades of investigation. Massive computationally intensive hydrodynamic simulations have been developed and run to model an ever-growing number of proposed progenitor channels. Further complicating the matter, a large number of sub-types of Type Ia supernovae have been identified in recent decades. Due to the massive computational load required, inference of the internal structure of Type Ia supernovae ejecta directly from observations using simulations has previously been computationally intractable. However, deep-learning emulators for radiation transport simulations have alleviated such barriers. We perform abundance tomography on 40 Type Ia supernovae from optical spectra using the radiative transfer code TARDIS accelerated by the probabilistic DALEK deep-learning emulator. We apply a parametric model of potential ejecta structures to comparatively investigate abundance distributions and internal ionization fractions of intermediate-mass elements between normal and 1991T-like Type Ia supernovae. Our inference shows that 1991T-like Type Ia supernovae are under-abundant in the typical intermediate mass elements that heavily contribute to the spectral line formation seen in normal Type Ia supernovae at early times. Additionally, we find that the intermediate-mass elements present in 1991T-like Type Ia supernovae are highly ionized compared to those in the normal Type Ia population. Finally, we conclude that the transition between normal and 1991T-like Type Ia supernovae appears to be continuous observationally and that the observed differences come out of a combination of both abundance and ionization fractions in these supernovae populations. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 15 Pages, 6 Figures, Article, Submitted to ApJ

Over 500 Days in the Life of the Photosphere of the Type Iax Supernova SN 2014dt

Authors: Yssavo Camacho-Neves, Saurabh W. Jha, Barnabas Barna, Mi Dai, Alexei V. Filippenko, Ryan J. Foley, Griffin Hosseinzadeh, D. Andrew Howell, Joel Johansson, Patrick Kelly, Wolfgang E. Kerzendorf, Lindsey A. Kwok, Conor Larison, Mark R. Magee, Curtis McCully, John T. O'Brien, Yen-Chen Pan, Viraj Pandya, Jaladh Singhal, Benjamin E. Stahl, Tamás Szalai, Meredith Wieber, Marc Williamson

Abstract: Type Iax supernovae (SN Iax) are the largest known class of peculiar white dwarf supernovae, distinct from normal Type Ia supernovae (SN Ia). The unique properties of SN Iax, especially their strong photospheric lines out to extremely late times, allow us to model their optical spectra and derive physical parameters for the long-lasting photosphere. We present an extensive spectral timeseries, inc… ▽ More Type Iax supernovae (SN Iax) are the largest known class of peculiar white dwarf supernovae, distinct from normal Type Ia supernovae (SN Ia). The unique properties of SN Iax, especially their strong photospheric lines out to extremely late times, allow us to model their optical spectra and derive physical parameters for the long-lasting photosphere. We present an extensive spectral timeseries, including 21 new spectra, of SN Iax 2014dt from +11 to +562 days after maximum light. We are able to reproduce the entire timeseries with a self-consistent, nearly unaltered deflagration explosion model from Fink et al. (2014) using TARDIS, an open-source radiative transfer code (Kerzendorf & Sim 2014; Kerzendorf et al. 2023). We find that the photospheric velocity of SN 2014dt slows its evolution between +64 and +148 days, which closely overlaps the phase when we see SN 2014dt diverge from the normal spectral evolution of SN Ia (+90 to +150 days). The photospheric velocity at these epochs, ~400$-$1000 km s$^{-1}$, may demarcate a boundary within the ejecta below which the physics of SN Iax and normal SN Ia differ. Our results suggest that SN 2014dt is consistent with a weak deflagration explosion model that leaves behind a bound remnant and drives an optically thick, quasi-steady-state wind creating the photospheric lines at late times. The data also suggest that this wind may weaken at epochs past +450 days, perhaps indicating a radioactive power source that has decayed away. △ Less

Submitted 1 May, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: Accepted to ApJ, 22 pages, 8 figures, 3 tables

Journal ref: 2023, ApJ, 951, 67

A JWST Near- and Mid-Infrared Nebular Spectrum of the Type Ia Supernova 2021aefx

Authors: Lindsey A. Kwok, Saurabh W. Jha, Tea Temim, Ori D. Fox, Conor Larison, Yssavo Camacho-Neves, Max J. Brenner Newman, Justin D. R. Pierel, Ryan J. Foley, Jennifer E. Andrews, Carles Badenes, Barnabas Barna, K. Azalee Bostroem, Maxime Deckers, Andreas Flors, Peter Garnavich, Melissa L. Graham, Or Graur, Griffin Hosseinzadeh, D. Andrew Howell, John P. Hughes, Joel Johansson, Sarah Kendrew, Wolfgang E. Kerzendorf, Keiichi Maeda , et al. (33 additional authors not shown)

Abstract: We present JWST near- and mid-infrared spectroscopic observations of the nearby normal Type Ia supernova SN 2021aefx in the nebular phase at $+255$ days past maximum light. Our Near Infrared Spectrograph (NIRSpec) and Mid Infrared Instrument (MIRI) observations, combined with ground-based optical data from the South African Large Telescope (SALT), constitute the first complete optical $+$ NIR $+$… ▽ More We present JWST near- and mid-infrared spectroscopic observations of the nearby normal Type Ia supernova SN 2021aefx in the nebular phase at $+255$ days past maximum light. Our Near Infrared Spectrograph (NIRSpec) and Mid Infrared Instrument (MIRI) observations, combined with ground-based optical data from the South African Large Telescope (SALT), constitute the first complete optical $+$ NIR $+$ MIR nebular SN Ia spectrum covering 0.3$-$14 $μ$m. This spectrum unveils the previously unobserved 2.5$-$5 $μ$m region, revealing strong nebular iron and stable nickel emission, indicative of high-density burning that can constrain the progenitor mass. The data show a significant improvement in sensitivity and resolution compared to previous Spitzer MIR data. We identify numerous NIR and MIR nebular emission lines from iron-group elements and as well as lines from the intermediate-mass element argon. The argon lines extend to higher velocities than the iron-group elements, suggesting stratified ejecta that are a hallmark of delayed-detonation or double-detonation SN Ia models. We present fits to simple geometric line profiles to features beyond 1.2 $μ$m and find that most lines are consistent with Gaussian or spherical emission distributions, while the [Ar III] 8.99 $μ$m line has a distinctively flat-topped profile indicating a thick spherical shell of emission. Using our line profile fits, we investigate the emissivity structure of SN 2021aefx and measure kinematic properties. Continued observations of SN 2021aefx and other SNe Ia with JWST will be transformative to the study of SN Ia composition, ionization structure, density, and temperature, and will provide important constraints on SN Ia progenitor and explosion models. △ Less

Submitted 10 February, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: published in ApJ Letters, 17 pages, 12 figures

Journal ref: ApJL, Volume 944 L3, 2023

The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

Authors: The Astropy Collaboration, Adrian M. Price-Whelan, Pey Lian Lim, Nicholas Earl, Nathaniel Starkman, Larry Bradley, David L. Shupe, Aarya A. Patil, Lia Corrales, C. E. Brasseur, Maximilian Nöthe, Axel Donath, Erik Tollerud, Brett M. Morris, Adam Ginsburg, Eero Vaher, Benjamin A. Weaver, James Tocknell, William Jamieson, Marten H. van Kerkwijk, Thomas P. Robitaille, Bruce Merry, Matteo Bachetti, H. Moritz Günther, Thomas L. Aldcroft , et al. (111 additional authors not shown)

Abstract: The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as… ▽ More The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 43 pages, 5 figures. To appear in ApJ. The author list has two parts: the authors that made significant contributions to the writing and/or coordination of the paper, followed by maintainers of and contributors to the Astropy Project. The position in the author list does not correspond to contributions to the Astropy Project as a whole

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

New mass estimates for massive binary systems: a probabilistic approach using polarimetric radiative transfer

Authors: Andrew G. Fullard, John T. O'Brien, Wolfgang E. Kerzendorf, Manisha Shrestha, Jennifer L. Hoffman, Richard Ignace, Patrick van der Smagt

Abstract: Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational wave sources such as binary black holes or neutron stars. For Wolf-Rayet stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimate… ▽ More Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational wave sources such as binary black holes or neutron stars. For Wolf-Rayet stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimates from binary systems which have spectroscopic $M \sin i$ measurements. Orbitally-phased polarization signals can encode the inclination angle of binary systems, where the Wolf-Rayet winds act as scattering regions. We investigated four Wolf-Rayet + O star binary systems, WR 42, WR 79, WR 127, and WR 153, with publicly available phased polarization data to estimate their masses. To avoid the biases present in analytic models of polarization while retaining computational expediency, we used a Monte Carlo radiative transfer model accurately emulated by a neural network. We used the emulated model to investigate the posterior distribution of parameters of our four systems. Our mass estimates calculated from the estimated inclination angles put strong constraints on existing mass estimates for three of the systems, and disagrees with the existing mass estimates for WR 153. We recommend a concerted effort to obtain polarization observations that can be used to estimate the masses of Wolf-Rayet binary systems and increase our understanding of their evolutionary paths. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: 24 pages, 13 figures, accepted to ApJ

Probabilistic Reconstruction of Type Ia Supernova SN 2002bo

Authors: John T. O'Brien, Wolfgang E. Kerzendorf, Andrew Fullard, Marc Williamson, Ruediger Pakmor, Johannes Buchner, Stephan Hachinger, Christian Vogl, James H. Gillanders, Andreas Floers, Patrick van der Smagt

Abstract: Manual fits to spectral times series of Type Ia supernovae have provided a method of reconstructing the explosion from a parametric model but due to lack of information about model uncertainties or parameter degeneracies direct comparison between theory and observation is difficult. In order to mitigate this important problem we present a new way to probabilistically reconstruct the outer ejecta o… ▽ More Manual fits to spectral times series of Type Ia supernovae have provided a method of reconstructing the explosion from a parametric model but due to lack of information about model uncertainties or parameter degeneracies direct comparison between theory and observation is difficult. In order to mitigate this important problem we present a new way to probabilistically reconstruct the outer ejecta of the normal Type Ia supernova SN 2002bo. A single epoch spectrum, taken 10 days before maximum light, is fit by a 13-parameter model describing the elemental composition of the ejecta and the explosion physics (density, temperature, velocity, and explosion epoch). Model evaluation is performed through the application of a novel rapid spectral synthesis technique in which the radiative transfer code, TARDIS, is accelerated by a machine-learning framework. Analysis of the posterior distribution reveals a complex and degenerate parameter space and allows direct comparison to various hydrodynamic models. Our analysis favors detonation over deflagration scenarios and we find that our technique offers a novel way to compare simulation to observation. △ Less

Submitted 12 October, 2021; v1 submitted 17 May, 2021; originally announced May 2021.

Comments: 12 Pages, 4 Figures, submitted to AAS Journals, comments and constructive criticism welcome

HOLISMOKES -- V. Microlensing of type II supernovae and time-delay inference through spectroscopic phase retrieval

Authors: J. Bayer, S. Huber, C. Vogl, S. H. Suyu, S. Taubenberger, D. Sluse, J. H. H. Chan, W. E. Kerzendorf

Abstract: We investigate strongly gravitationally lensed type II supernovae (LSNe II) for time-delay cosmography incorporating microlensing effects, which expands on previous microlensing studies of type Ia supernovae (SNe Ia). We use the radiative-transfer code ${\rm \small TARDIS}$ to recreate five spectra of the prototypical SN 1999em at different times within the plateau phase of the light curve. The mi… ▽ More We investigate strongly gravitationally lensed type II supernovae (LSNe II) for time-delay cosmography incorporating microlensing effects, which expands on previous microlensing studies of type Ia supernovae (SNe Ia). We use the radiative-transfer code ${\rm \small TARDIS}$ to recreate five spectra of the prototypical SN 1999em at different times within the plateau phase of the light curve. The microlensing-induced deformations of the spectra and light curves are calculated by placing the SN into magnification maps generated with the code ${\rm \small GERLUMPH}$. We study the impact of microlensing on the color curves and find that there is no strong influence on them during the investigated time interval of the plateau phase. The color curves are only weakly affected by microlensing due to the almost achromatic behavior of the intensity profiles. However, the lack of non-linear structure in the color curves makes time-delay measurements difficult given the possible presence of differential dust extinction. Therefore, we further investigate SN phase inference through spectral absorption lines under the influence of microlensing and Gaussian noise. As the spectral features shift to longer wavelengths with progressing time after explosion, the measured wavelength of a specific absorption line provides information on the epoch of the SN. The comparison between retrieved epochs of two observed lensing images then gives the time delay of the images. We find that the phase retrieval method using spectral features yields accurate delays with uncertainties $\small {\lesssim}$2 days, making it a promising approach. △ Less

Submitted 8 September, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 19 pages, 17 figures

Journal ref: A&A 653, A29 (2021)

Production of very light elements and strontium in the early ejecta of neutron star mergers

Authors: A. Perego, D. Vescovi, A. Fiore, L. Chiesa, C. Vogl, S. Benetti, S. Bernuzzi, M. Branchesi, S. Cristallo, E. Cappellaro, A. Flors, W. E. Kerzendorf, D. Radice

Abstract: We study the production of very light elements ($Z < 20$) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition dependent nuclear equations o… ▽ More We study the production of very light elements ($Z < 20$) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition dependent nuclear equations of state and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, $\sim 0.5-2 \times 10^{-6} M_{\odot}$ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while $\sim 1.5-11 \times 10^{-6} M_{\odot}$ of Helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r-process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including non local thermodynamics equilibrium effects. The latter could be crucial to observe He lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long lived remnant) the one required to explain early spectral features in the kilonova of GW170817. △ Less

Submitted 8 November, 2021; v1 submitted 18 September, 2020; originally announced September 2020.

Comments: 21 pages, 9 figures, Accepted for publication in The Astrophysical Journal

Journal ref: ApJ 925, 22 (2022)

Dalek -- a deep-learning emulator for TARDIS

Authors: Wolfgang E. Kerzendorf, Christian Vogl, Johannes Buchner, Gabriella Contardo, Marc Williamson, Patrick van der Smagt

Abstract: Supernova spectral time series contain a wealth of information about the progenitor and explosion process of these energetic events. The modeling of these data requires the exploration of very high dimensional posterior probabilities with expensive radiative transfer codes. Even modest parametrizations of supernovae contain more than ten parameters and a detailed exploration demands at least sever… ▽ More Supernova spectral time series contain a wealth of information about the progenitor and explosion process of these energetic events. The modeling of these data requires the exploration of very high dimensional posterior probabilities with expensive radiative transfer codes. Even modest parametrizations of supernovae contain more than ten parameters and a detailed exploration demands at least several million function evaluations. Physically realistic models require at least tens of CPU minutes per evaluation putting a detailed reconstruction of the explosion out of reach of traditional methodology. The advent of widely available libraries for the training of neural networks combined with their ability to approximate almost arbitrary functions with high precision allows for a new approach to this problem. Instead of evaluating the radiative transfer model itself, one can build a neural network proxy trained on the simulations but evaluating orders of magnitude faster. Such a framework is called an emulator or surrogate model. In this work, we present an emulator for the TARDIS supernova radiative transfer code applied to Type Ia supernova spectra. We show that we can train an emulator for this problem given a modest training set of a hundred thousand spectra (easily calculable on modern supercomputers). The results show an accuracy on the percent level (that are dominated by the Monte Carlo nature of TARDIS and not the emulator) with a speedup of several orders of magnitude. This method has a much broader set of applications and is not limited to the presented problem. △ Less

Submitted 3 July, 2020; originally announced July 2020.

Comments: 6 pages;5 figures submitted to AAS Journals. Constructive Criticism invited

Distributed peer review enhanced with natural language processing and machine learning

Authors: Wolfgang E. Kerzendorf, Ferdinando Patat, Dominic Bordelon, Glenn van de Ven, Tyler A. Pritchard

Abstract: While ancient scientists often had patrons to fund their work, peer review of proposals for the allocation of resources is a foundation of modern science. A very common method is that proposals are evaluated by a small panel of experts (due to logistics and funding limitations) nominated by the grant-giving institutions. The expert panel process introduces several issues - most notably: 1) biases… ▽ More While ancient scientists often had patrons to fund their work, peer review of proposals for the allocation of resources is a foundation of modern science. A very common method is that proposals are evaluated by a small panel of experts (due to logistics and funding limitations) nominated by the grant-giving institutions. The expert panel process introduces several issues - most notably: 1) biases introduced in the selection of the panel. 2) experts have to read a very large number of proposals. Distributed Peer Review promises to alleviate several of the described problems by distributing the task of reviewing among the proposers. Each proposer is given a limited number of proposals to review and rank. We present the result of an experiment running a machine-learning enhanced distributed peer review process for allocation of telescope time at the European Southern Observatory. In this work, we show that the distributed peer review is statistically the same as a `traditional' panel, that our machine learning algorithm can predict expertise of reviewers with a high success rate, and we find that seniority and reviewer expertise have an influence on review quality. The general experience has been overwhelmingly praised from the participating community (using an anonymous feedback mechanism). △ Less

Submitted 8 April, 2020; originally announced April 2020.

Comments: 20 pages including supplementary information. Constructive criticism more than welcome. Accepted and published in Nature Astronomy

Journal ref: Nature Astronomy 2020

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)

Search for Surviving Companions of Progenitors of Young LMC Type Ia Supernova Remnants

Authors: Chuan-Jui Li, Wolfgang E. Kerzendorf, You-Hua Chu, Ting-Wan Chen, Tuan Do, Robert A. Gruendl, Abigail Holmes, Ryoko Ishioka, Bruno Leibundgut, Kuo-Chuan Pan, Paul M. Ricker, Daniel Weisz

Abstract: We have used two methods to search for surviving companions of Type Ia supernova progenitors in three Balmer-dominated supernova remnants (SNRs) in the Large Magellanic Cloud: 0519-69.0, 0505-67.9 (DEM L71), and 0548-70.4. In the first method, we use the Hubble Space Telescope photometric measurements of stars to construct color-magnitude diagrams (CMDs), and compare positions of stars in the CMDs… ▽ More We have used two methods to search for surviving companions of Type Ia supernova progenitors in three Balmer-dominated supernova remnants (SNRs) in the Large Magellanic Cloud: 0519-69.0, 0505-67.9 (DEM L71), and 0548-70.4. In the first method, we use the Hubble Space Telescope photometric measurements of stars to construct color-magnitude diagrams (CMDs), and compare positions of stars in the CMDs with those expected from theoretical post-impact evolution of surviving main sequence or helium star companions. No obvious candidates of surviving companion are identified in this photometric search. Future models for surviving red giant companions or with different explosion mechanisms are needed for thorough comparisons with these observations in order to make more definitive conclusions. In the second method, we use Multi-Unit Spectroscopic Explorer (MUSE) observations of 0519-69.0 and DEM L71 to carry out spectroscopic analyses of stars in order to use large peculiar radial velocities as diagnostics of surviving companions. We find a star in 0519-69.0 and a star in DEM L71 moving at radial velocities of 182 $\pm$ 0 km s$^{-1}$ and 213 $\pm$ 0 km s$^{-1}$, more than 2.5$σ$ from the mean radial velocity of the underlying stellar population, 264 km s$^{-1}$ and 270 km s$^{-1}$, respectively. These stars need higher-quality spectra to investigate their abundances and rotation velocities to determine whether they are indeed surviving companions of the SN progenitors. △ Less

Submitted 12 November, 2021; v1 submitted 2 October, 2019; originally announced October 2019.

Comments: 25 pages, 9 tables, 22 figures, accepted for publication in The Astrophysical Journal (ApJ)

Journal ref: The Astrophysical Journal, Volume 886, Number 2 (2019)

Relativistic redshift of the star S0-2 orbiting the Galactic center supermassive black hole

Authors: Tuan Do, Aurelien Hees, Andrea Ghez, Gregory D. Martinez, Devin S. Chu, Siyao Jia, Shoko Sakai, Jessica R. Lu, Abhimat K. Gautam, Kelly Kosmo O'Neil, Eric E. Becklin, Mark R. Morris, Keith Matthews, Shogo Nishiyama, Randy Campbell, Samantha Chappell, Zhuo Chen, Anna Ciurlo, Arezu Dehghanfar, Eulalia Gallego-Cano, Wolfgang E. Kerzendorf, James E. Lyke, Smadar Naoz, Hiromi Saida, Rainer Schödel , et al. (4 additional authors not shown)

Abstract: General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach… ▽ More General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $Υ$. Our result, $Υ=0.88 \pm 0.17$, is consistent with General Relativity ($Υ=1$) and excludes a Newtonian model ($Υ=0$ ) with a statistical significance of 5 $σ$. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 78 pages, 21 figures, accepted to Science, astrometry and radial velocity measurements, and posterior chains available here, http://www.astro.ucla.edu/~ghezgroup/gc/gr/data.zip

Envisioning the next decade of Galactic Center science: a laboratory for the study of the physics and astrophysics of supermassive black holes

Authors: Tuan Do, Andrea Ghez, Jessica R. Lu, Mark Morris, Matthew Hosek Jr., Aurelien Hees, Smadar Naoz, Anna Ciurlo, Philip J. Armitage, Rachael L Beaton, Eric Becklin, Andrea Bellini, Rory O. Bentley, Joss Bland-Hawthorn, Sukanya Chakrabarti, Zhuo Chen, Devin S. Chu, Arezu Dehghanfar, Charles F. Gammie, Abhimat K. Gautam, Reinhard Genzel, Jenny Greene, Daryl Haggard, Joseph Hora, Wolfgang E. Kerzendorf , et al. (16 additional authors not shown)

Abstract: As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct desc… ▽ More As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 9 pages, 4 figures, submitted for Astro2020 White Paper

Multi-Messenger Astronomy with Extremely Large Telescopes

Authors: Ryan Chornock, Philip S. Cowperthwaite, Raffaella Margutti, Dan Milisavljevic, Kate D. Alexander, Igor Andreoni, Iair Arcavi, Adriano Baldeschi, Jennifer Barnes, Eric Bellm, Paz Beniamini, Edo Berger, Christopher P. L. Berry, Federica Bianco, Peter K. Blanchard, Joshua S. Bloom, Sarah Burke-Spolaor, Eric Burns, Dario Carbone, S. Bradley Cenko, Deanne Coppejans, Alessandra Corsi, Michael Coughlin, Maria R. Drout, Tarraneh Eftekhari , et al. (60 additional authors not shown)

Abstract: The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment hi… ▽ More The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: White paper submitted to the Astro2020 Decadal Survey

Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s

Authors: R. J. Foley, K. D. Alexander, I. Andreoni, I. Arcavi, K. Auchettl, J. Barnes, G. Baym, E. C. Bellm, A. M. Beloborodov, N. Blagorodnova, J. P. Blakeslee, P. R. Brady, M. Branchesi, J. S. Brown, N. Butler, M. Cantiello, R. Chornock, D. O. Cook, J. Cooke, D. L. Coppejans, A. Corsi, S. M. Couch, M. W. Coughlin, D. A. Coulter, P. S. Cowperthwaite , et al. (88 additional authors not shown)

Abstract: As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse s… ▽ More As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submitted to Astro2020

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)

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)

Limits on stable iron in Type$\,$Ia supernovae from NIR spectroscopy

Authors: A. Flörs, J. Spyromilio, K. Maguire, S. Taubenberger, W. E. Kerzendorf, S. Dhawan

Abstract: We obtained optical and near-infrared spectra of Type$\,$Ia supernovae (SNe$\,$Ia) at epochs ranging from 224 to 496 days after the explosion. The spectra show emission lines from forbidden transitions of singly ionised iron and cobalt atoms. We used non-local thermodynamic equilibrium (NLTE) modelling of the first and second ionisation stages of iron, nickel, and cobalt to fit the spectra using a… ▽ More We obtained optical and near-infrared spectra of Type$\,$Ia supernovae (SNe$\,$Ia) at epochs ranging from 224 to 496 days after the explosion. The spectra show emission lines from forbidden transitions of singly ionised iron and cobalt atoms. We used non-local thermodynamic equilibrium (NLTE) modelling of the first and second ionisation stages of iron, nickel, and cobalt to fit the spectra using a sampling algorithm allowing us to probe a broad parameter space. We derive velocity shifts, line widths, and abundance ratios for iron and cobalt. The measured line widths and velocity shifts of the singly ionised ions suggest a shared emitting region. Our data are fully compatible with radioactive $^{56}$Ni decay as the origin for cobalt and iron. We compare the measured abundance ratios of iron and cobalt to theoretical predictions of various SN$\,$Ia explosion models. These models include, in addition to $^{56}$Ni, different amounts of $^{57}$Ni and stable $^{54,56}$Fe. We can exclude models that produced only $^{54,56}$Fe or only $^{57}$Ni in addition to $^{56}$Ni. If we consider a model that has $^{56}$Ni, $^{57}$Ni, and $^{54,56}$Fe then our data imply that these ratios are $^{54,56}$Fe / $^{56}$Ni $=0.272\pm0.086$ and $^{57}$Ni / $^{56}$Ni $=0.032\pm0.011$. △ Less

Submitted 25 October, 2018; originally announced October 2018.

Comments: 10 pages, 7 figures, Accepted for publication in A&A

Journal ref: A&A 620, A200 (2018)

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

Three Hypervelocity White Dwarfs in Gaia DR2: Evidence for Dynamically Driven Double-Degenerate Double-Detonation Type Ia Supernovae

Authors: Ken J. Shen, Douglas Boubert, Boris T. Gänsicke, Saurabh W. Jha, Jennifer E. Andrews, Laura Chomiuk, Ryan J. Foley, Morgan Fraser, Mariusz Gromadzki, James Guillochon, Marissa M. Kotze, Kate Maguire, Matthew R. Siebert, Nathan Smith, Jay Strader, Carles Badenes, Wolfgang E. Kerzendorf, Detlev Koester, Markus Kromer, Broxton Miles, Rüdiger Pakmor, Josiah Schwab, Odette Toloza, Silvia Toonen, Dean M. Townsley , et al. (1 additional authors not shown)

Abstract: Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 k… ▽ More Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia. △ Less

Submitted 18 July, 2018; v1 submitted 30 April, 2018; originally announced April 2018.

Comments: Accepted for publication in ApJ. Minor corrections for clarity. D6 spectra are available as ancillary data files

Tycho-B: an unlikely companion for SN 1572

Authors: W. E. Kerzendorf, Knox S. Long, P. Frank Winkler, Tuan Do

Abstract: If some or all Type Ia SNe arise from accretion onto a massive WD from a companion, then the companion will remain in some form after the SN explosion. Tycho-B is an unusual, relatively hot star along the line of sight to Tycho's SNR -- conclusively shown to be a Type Ia -- and has been suggested as such a companion. If the interior of Tycho's SNR contains unshocked Fe, and if Tycho-B is either wi… ▽ More If some or all Type Ia SNe arise from accretion onto a massive WD from a companion, then the companion will remain in some form after the SN explosion. Tycho-B is an unusual, relatively hot star along the line of sight to Tycho's SNR -- conclusively shown to be a Type Ia -- and has been suggested as such a companion. If the interior of Tycho's SNR contains unshocked Fe, and if Tycho-B is either within the SNR shell or in the background, then one might hope to see evidence of this in the UV spectrum. Such is the case for SN\,1006, where spectra of the background Schweizer-Middleditch star, as well as two AGNs, show broad absorption lines of Fe II. To test this idea, we have used STIS on HST to obtain a UV spectrum of Tycho-B. The observed spectrum, however, shows no evidence of Fe II absorption. Furthermore, a luminosity distance estimate using UV and optical spectra of Tycho-B suggests that the star is consistent with a foreground interloper. We conclude either that Tycho B is nearer than Tycho's SNR, or that all of the Fe in the interior of Tycho's SNR is more highly ionized. △ Less

Submitted 24 July, 2018; v1 submitted 20 March, 2018; originally announced March 2018.

Comments: 9 pages, 7 figures, submitted to MNRAS and comments (particularly constructive criticisms) are always welcome

Modelling the early time behaviour of type Ia supernovae: effects of the $^{56}$Ni distribution

Authors: M. R. Magee, S. A. Sim, R. Kotak, W. E. Kerzendorf

Abstract: Recent studies have demonstrated the diversity in type Ia supernovae (SNe Ia) at early times and highlighted a need for a better understanding of the explosion physics as manifested by observations soon after explosion. To this end, we present a Monte Carlo code designed to model the light curves of radioactively driven, hydrogen-free transients from explosion to approximately maximum light. In th… ▽ More Recent studies have demonstrated the diversity in type Ia supernovae (SNe Ia) at early times and highlighted a need for a better understanding of the explosion physics as manifested by observations soon after explosion. To this end, we present a Monte Carlo code designed to model the light curves of radioactively driven, hydrogen-free transients from explosion to approximately maximum light. In this initial study, we have used a parametrised description of the ejecta in SNe Ia, and performed a parameter study of the effects of the $^{56}$Ni distribution on the observed colours and light curves for a fixed $^{56}$Ni mass of 0.6 $M_\odot$. For a given density profile, we find that models with $^{56}$Ni extending throughout the entirety of the ejecta are typically brighter and bluer shortly after explosion. Additionally, the shape of the density profile itself also plays an important role in determining the shape, rise time, and colours of observed light curves. We find that the multi-band light curves of at least one SNe Ia (SN 2009ig) are inconsistent with less extended $^{56}$Ni distributions, but show good agreement with models that incorporate $^{56}$Ni throughout the entire ejecta. We further demonstrate that comparisons with full $UVOIR$ colour light curves are powerful tools in discriminating various $^{56}$Ni distributions, and hence explosion models. △ Less

Submitted 21 August, 2018; v1 submitted 12 March, 2018; originally announced March 2018.

Comments: 14 pages, 8 figures, 2 tables. Minor changes in notation to match published version in Astronomy & Astrophysics

Journal ref: A&A 614, A115 (2018)

The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package

Authors: The Astropy Collaboration, A. M. Price-Whelan, B. M. Sipőcz, H. M. Günther, P. L. Lim, S. M. Crawford, S. Conseil, D. L. Shupe, M. W. Craig, N. Dencheva, A. Ginsburg, J. T. VanderPlas, L. D. Bradley, D. Pérez-Suárez, M. de Val-Borro, T. L. Aldcroft, K. L. Cruz, T. P. Robitaille, E. J. Tollerud, C. Ardelean, T. Babej, M. Bachetti, A. V. Bakanov, S. P. Bamford, G. Barentsen , et al. (112 additional authors not shown)

Abstract: The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy p… ▽ More The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project. △ Less

Submitted 16 January, 2018; v1 submitted 8 January, 2018; originally announced January 2018.

Comments: Minor changes to author list and title. Comments and feedback welcome through the paper source repository: https://github.com/astropy/astropy-v2.0-paper For more information about Astropy, see http://www.astropy.org/

No surviving stellar companion for Cassiopeia A

Authors: Wolfgang E. Kerzendorf, Tuan Do, Selma E. de Mink, Ylva Götberg, Dan Millisaljevic, Emmanouil Zapartas, Mathieu Renzo, Stephen Justham, Philipp Podsiadlowski, Robert A. Fesen

Abstract: Massive stars in binaries can give rise to extreme phenomena such as X-ray binaries and gravitational wave sources after one or both stars end their lives as core-collapse supernovae. Stars in close orbit around a stellar or compact companion are expected to explode as "stripped-envelope supernovae", showing no (Type Ib/c) or little (Type IIb) signs of hydrogen in the spectra, because hydrogen-ric… ▽ More Massive stars in binaries can give rise to extreme phenomena such as X-ray binaries and gravitational wave sources after one or both stars end their lives as core-collapse supernovae. Stars in close orbit around a stellar or compact companion are expected to explode as "stripped-envelope supernovae", showing no (Type Ib/c) or little (Type IIb) signs of hydrogen in the spectra, because hydrogen-rich progenitors are too large to fit. The physical processes responsible for the stripping process and the fate of the companion are still very poorly understood. Aiming to find new clues, we investigate Cas~A, which is a very young ($\sim$340 \,yr) and near ($\sim$3.4\,kpc) remnant of a core collapse supernova. Cas~A has been subject to several searches for possible companions, all unsuccessfully. We present new measurements of the proper motions and photometry of stars in the vicinity based on deep HST ACS/WFC and WFC3-IR data. We identify stellar sources that are close enough in projection, but using their proper motions we show that none are compatible with being at the location of center at the time of explosion, in agreement with earlier findings. Our photometric measurements allow us to place much deeper (order of magnitude) upper limits on the brightness of possible undetected companions. We systematically compare them with model predictions for a wide variety of scenarios. We can confidently rule out the presence of any stellar companion of any reasonable mass and age (main sequence, pre main sequence or stripped). Although this finding is not in direct conflict with model predictions, it does rule out what many considered to be the most likely evolutionary scenario for type IIb supernova. (abstract abbreviated) △ Less

Submitted 5 December, 2018; v1 submitted 31 October, 2017; originally announced November 2017.

Comments: 13 pages, accepted by A&A

Journal ref: A&A 623, A34 (2019)

A kilonova as the electromagnetic counterpart to a gravitational-wave source

Authors: S. J. Smartt, T. -W. Chen, A. Jerkstrand, M. Coughlin, E. Kankare, S. A. Sim, M. Fraser, C. Inserra, K. Maguire, K. C. Chambers, M. E. Huber, T. Kruhler, G. Leloudas, M. Magee, L. J. Shingles, K. W. Smith, D. R. Young, J. Tonry, R. Kotak, A. Gal-Yam, J. D. Lyman, D. S. Homan, C. Agliozzo, J. P. Anderson, C. R. Angus C. Ashall , et al. (96 additional authors not shown)

Abstract: Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a r… ▽ More Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements. △ Less

Submitted 17 October, 2017; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: Nature, in press, DOI 10.1038/nature24303. Data files will be made available at http://www.pessto.org

A Search for a Surviving White Dwarf Companion in SN 1006

Authors: W. E. Kerzendorf, G. Strampelli, K. J. Shen, J. Schwab, R. Pakmor, T. Do, J. Buchner, A. Rest

Abstract: Multiple channels have been proposed to produce Type Ia supernovae, with many scenarios suggesting that the exploding white dwarf accretes from a binary companion pre-explosion. In almost all cases, theory suggests that this companion will survive. However, no such companion has been unambiguously identified in ancient supernova remnants -- possibly falsifying the accretion scenario. Existing surv… ▽ More Multiple channels have been proposed to produce Type Ia supernovae, with many scenarios suggesting that the exploding white dwarf accretes from a binary companion pre-explosion. In almost all cases, theory suggests that this companion will survive. However, no such companion has been unambiguously identified in ancient supernova remnants -- possibly falsifying the accretion scenario. Existing surveys, however, have only looked for stars as faint as $\approx 0.1 L_\odot$ and thus would have missed a surviving white dwarf companion. In this work, we present very deep DECAM imaging $(u, g, r, z)$ of the Type Ia supernova remnant SN 1006 specifically to search for a potential surviving white dwarf companion. We find no object within the inner third of the SN 1006 remnant that is consistent with a relatively young cooling white dwarf. We find that if there is a companion white dwarf, it must have formed long ago and cooled undisturbed for $> 10^8$ yr to be consistent with the redder objects in our sample. We conclude that our findings are consistent with the complete destruction of the secondary (such as in a merger) or an unexpectedly cool and thus very dim surviving companion white dwarf. △ Less

Submitted 19 September, 2017; originally announced September 2017.

Comments: 9 pages, 5 figure, submitted to MNRAS - comments welcome

Spectral sequences of Type Ia supernovae. I. Connecting normal and sub-luminous SN Ia and the presence of unburned carbon

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

Abstract: Type Ia supernovae are generally agreed to arise from thermonuclear explosions of carbon-oxygen white dwarfs. The actual path to explosion, however, remains elusive, with numerous plausible parent systems and explosion mechanisms suggested. Observationally, type Ia supernovae have multiple subclasses, distinguished by their lightcurves and spectra. This raises the question whether these reflect th… ▽ More Type Ia supernovae are generally agreed to arise from thermonuclear explosions of carbon-oxygen white dwarfs. The actual path to explosion, however, remains elusive, with numerous plausible parent systems and explosion mechanisms suggested. Observationally, type Ia supernovae have multiple subclasses, distinguished by their lightcurves and spectra. This raises the question whether these reflect that multiple mechanisms occur in nature, or instead that explosions have a large but continuous range of physical properties. We revisit the idea that normal and 91bg-like supernovae can be understood as part of a spectral sequence, in which changes in temperature dominate. Specifically, we find that a single ejecta structure is sufficient to provide reasonable fits of both the normal type Ia supernova SN~2011fe and the 91bg-like SN~2005bl, provided that the luminosity and thus temperature of the ejecta are adjusted appropriately. This suggests that the outer layers of the ejecta are similar, thus providing some support of a common explosion mechanism. Our spectral sequence also helps to shed light on the conditions under which carbon can be detected in pre-maximum SN~Ia spectra -- we find that emission from iron can "fill in" the carbon trough in cool SN~Ia. This may indicate that the outer layers of the ejecta of events in which carbon is detected are relatively metal poor compared to events where carbon is not detected. △ Less

Submitted 26 July, 2017; originally announced July 2017.

Abundance Tomography of Type Iax SN 2011ay with TARDIS

Authors: Barnabas Barna, Tamas Szalai, Markus Kromer, Wolfgang E. Kerzendorf, Jozsef Vinko, Jeffrey M. Silverman, G. H. Marion, J. Craig Wheeler

Abstract: We present a detailed spectral analysis of Type Iax SN 2011ay. Our spectra cover epochs between -3 and +19 days with respect to the maximum light in B-band. This time range allows us to employ a so-called abundance tomography technique. The synthetic spectral fitting was made with the 1D Monte Carlo radiative transfer code TARDIS. In this paper, we describe our method to fit multiple epochs with a… ▽ More We present a detailed spectral analysis of Type Iax SN 2011ay. Our spectra cover epochs between -3 and +19 days with respect to the maximum light in B-band. This time range allows us to employ a so-called abundance tomography technique. The synthetic spectral fitting was made with the 1D Monte Carlo radiative transfer code TARDIS. In this paper, we describe our method to fit multiple epochs with a self-consistent, stratified atmospheric model. We compare our results to previously published SYN++ models and the predictions of different explosion scenarios. Using a fixed density profile (exponential fit of W7), we find that a uniform abundance model cannot reproduce the spectral features before maximum light because of the emergence of excessively strong Fe lines. In our best-fit TARDIS model, we find an abundance profile that separated into two different regimes: a well-mixed region under 10,000 km s$^{-1}$ and a stratified region with decreasing IGE abundances above 10,000 km s$^{-1}$. Based on a detailed comparative analysis, our conclusion is that the available pure deflagration models cannot fully explain either the observed properties of SN 2011ay or the results of our TARDIS modeling. Further examinations are necessary to find an adequate explanation for the origin of this object. △ Less

Submitted 25 July, 2017; originally announced July 2017.

Comments: 11 pages, 13 figures, 2 tables; accepted for publication in MNRAS

Extremely late photometry of SN~2011fe

Authors: W. E. Kerzendorf, C. McCully, S. Taubenberger, A. Jerkstrand, I. Seitenzahl, A. J. Ruiter, J. Spyromilio, K. S. Long, C. Fransson

Abstract: Type Ia supernovae are widely accepted to be the outcomes of thermonuclear explosions in white dwarf stars. However, many details of these explosions remain uncertain (e.g. the mass, ignition mechanism, and flame speed). Theory predicts that at very late times (beyond 1000 d) it might be possible to distinguish between explosion models. Few very nearby supernovae can be observed that long after th… ▽ More Type Ia supernovae are widely accepted to be the outcomes of thermonuclear explosions in white dwarf stars. However, many details of these explosions remain uncertain (e.g. the mass, ignition mechanism, and flame speed). Theory predicts that at very late times (beyond 1000 d) it might be possible to distinguish between explosion models. Few very nearby supernovae can be observed that long after the explosion. The Type Ia supernova SN 2011fe located in M101 and along a line of sight with negligible extinction, provides us with the once-in-a-lifetime chance to obtain measurements that may distinguish between theoretical models. In this work, we present the analysis of photometric data of SN 2011fe taken between 900 and 1600 days after explosion with Gemini and HST. At these extremely late epochs theory suggests that the light curve shape might be used to measure isotopic abundances which is a useful model discriminant. However, we show in this work that there are several currently not well constrained physical processes introducing large systematic uncertainties to the isotopic abundance measurement. We conclude that without further detailed knowledge of the physical processes at this late stage one cannot reliably exclude any models on the basis of this dataset. △ Less

Submitted 20 September, 2017; v1 submitted 5 June, 2017; originally announced June 2017.

Comments: accepted by MNRAS, comments still welcome

Bayesian modelling of uncertainties of Monte Carlo radiative-transfer simulations

Authors: Frederik Beaujean, Hans C. Eggers, Wolfgang E. Kerzendorf

Abstract: One of the big challenges in astrophysics is the comparison of complex simulations to observations. As many codes do not directly generate observables (e.g. hydrodynamic simulations), the last step in the modelling process is often a radiative-transfer treatment. For this step, the community relies increasingly on Monte Carlo radiative transfer due to the ease of implementation and scalability wit… ▽ More One of the big challenges in astrophysics is the comparison of complex simulations to observations. As many codes do not directly generate observables (e.g. hydrodynamic simulations), the last step in the modelling process is often a radiative-transfer treatment. For this step, the community relies increasingly on Monte Carlo radiative transfer due to the ease of implementation and scalability with computing power. We consider simulations in which the number of photon packets is Poisson distributed, while the weight assigned to a single photon packet follows any distribution of choice. We show how to estimate the statistical uncertainty of the sum of weights in each bin from the output of a single radiative-transfer simulation. Our Bayesian approach produces a posterior distribution that is valid for any number of packets in a bin, even zero packets, and is easy to implement in practice. Our analytic results for large number of packets show that we generalise existing methods that are valid only in limiting cases. The statistical problem considered here appears in identical form in a wide range of Monte Carlo simulations including particle physics and importance sampling. It is particularly powerful in extracting information when the available data are sparse or quantities are small. △ Less

Submitted 7 May, 2018; v1 submitted 29 May, 2017; originally announced May 2017.

Comments: 8 pages, 2 figures. Match affiliations and abstract with journal version

Knowledge discovery through text-based similarity searches for astronomy literature

Authors: W. E. Kerzendorf

Abstract: The increase in the number of researchers coupled with the ease of publishing and distribution of scientific papers (due to technological advancements) has resulted in a dramatic increase in astronomy literature. This has likely led to the predicament that the body of the literature is too large for traditional human consumption and that related and crucial knowledge is not discovered by researche… ▽ More The increase in the number of researchers coupled with the ease of publishing and distribution of scientific papers (due to technological advancements) has resulted in a dramatic increase in astronomy literature. This has likely led to the predicament that the body of the literature is too large for traditional human consumption and that related and crucial knowledge is not discovered by researchers. In addition to the increased production of astronomical literature, recent decades have also brought several advancements in computational linguistics. Especially, the machine-aided processing of literature dissemination might make it possible to convert this stream of papers into a coherent knowledge set. In this paper, we present the application of computational linguistics techniques to astronomy literature. In particular, we developed a tool that will find similar articles purely based on text content from an input paper. We find that our technique performs robustly in comparison with other tools recommending articles given a reference paper (known as recommender system). Our novel tool shows the great power in combining computational linguistics with astronomy literature and suggests that additional research in this endeavor will likely produce even better tools that will help researchers cope with the vast amounts of knowledge being produced. △ Less

Submitted 20 August, 2018; v1 submitted 16 May, 2017; originally announced May 2017.

Comments: 6 pages, 5 figures, subm., comments welcome

No evidence for younger stellar generations within the intermediate age massive clusters NGC 1783, NGC 1696 and NGC 411

Authors: I. Cabrera-Ziri, F. Niederhofer, N. Bastian, M. Rejkuba, E. Balbinot, W. E. Kerzendorf, S. S. Larsen, A. D. Mackey, E. Dalessandro, A. Mucciarelli, C. Charbonnel, M. Hilker, M. Gieles, V. Hénault-Brunet

Abstract: Recently, Li et al. (2016) claimed to have found evidence for multiple generations of stars in the intermediate age clusters NGC 1783, NGC 1696 and NGC 411 in the Large and Small Magellanic Clouds (LMC/SMC). Here we show that these young stellar populations are present in the field regions around these clusters and are not likely associated with the clusters themselves. Using the same datasets, we… ▽ More Recently, Li et al. (2016) claimed to have found evidence for multiple generations of stars in the intermediate age clusters NGC 1783, NGC 1696 and NGC 411 in the Large and Small Magellanic Clouds (LMC/SMC). Here we show that these young stellar populations are present in the field regions around these clusters and are not likely associated with the clusters themselves. Using the same datasets, we find that the background subtraction method adopted by the authors does not adequately remove contaminating stars in the small number Poisson limit. Hence, we conclude that their results do not provide evidence of young generations of stars within these clusters. △ Less

Submitted 20 April, 2016; originally announced April 2016.

Comments: 7 pages, 8 figures, accepted for publication in MNRAS

Spectroscopy of the Type Ia supernova 2011fe past 1000 days

Authors: S. Taubenberger, N. Elias-Rosa, W. E. Kerzendorf, S. Hachinger, J. Spyromilio, C. Fransson, M. Kromer, A. J. Ruiter, I. R. Seitenzahl, S. Benetti, E. Cappellaro, A. Pastorello, M. Turatto, A. Marchetti

Abstract: In this letter we present an optical spectrum of SN 2011fe taken 1034 d after the explosion, several hundred days later than any other spectrum of a Type Ia supernova (disregarding light-echo spectra and local-group remnants). The spectrum is still dominated by broad emission features, with no trace of a light echo or interaction of the supernova ejecta with surrounding interstellar material. Comp… ▽ More In this letter we present an optical spectrum of SN 2011fe taken 1034 d after the explosion, several hundred days later than any other spectrum of a Type Ia supernova (disregarding light-echo spectra and local-group remnants). The spectrum is still dominated by broad emission features, with no trace of a light echo or interaction of the supernova ejecta with surrounding interstellar material. Comparing this extremely late spectrum to an earlier one taken 331 d after the explosion, we find that the most prominent feature at 331 d - [Fe III] emission around 4700 A - has entirely faded away, suggesting a significant change in the ionisation state. Instead, [Fe II] lines are probably responsible for most of the emission at 1034 d. An emission feature at 6300-6400 A has newly developed at 1034 d, which we tentatively identify with Fe I λ6359, [Fe I] λλ6231, 6394 or [O I] λλ6300, 6364. Interestingly, the features in the 1034-d spectrum seem to be collectively redshifted, a phenomenon that we currently have no convincing explanation for. We discuss the implications of our findings for explosion models, but conclude that sophisticated spectral modelling is required for any firm statement. △ Less

Submitted 15 December, 2014; v1 submitted 27 November, 2014; originally announced November 2014.

Comments: 6 pages, 3 figures, MNRAS Letters accepted

Very late photometry of SN 2011fe

Authors: W. E. Kerzendorf, S. Taubenberger, I. R. Seitenzahl, A. J. Ruiter

Abstract: The Type Ia supernova SN 2011fe is one of the closest supernovae of the past decades. Due to its proximity and low dust extinction, this object provides a very rare opportunity to study the extremely late time evolution ($>900$ days) of thermonuclear supernovae. In this Letter, we present our photometric data of SN 2011fe taken at an unprecedented late epoch of $\approx 930$ days with GMOS-N mount… ▽ More The Type Ia supernova SN 2011fe is one of the closest supernovae of the past decades. Due to its proximity and low dust extinction, this object provides a very rare opportunity to study the extremely late time evolution ($>900$ days) of thermonuclear supernovae. In this Letter, we present our photometric data of SN 2011fe taken at an unprecedented late epoch of $\approx 930$ days with GMOS-N mounted on the Gemini North telescope ($g=23.43 \pm 0.28$, $r=24.14 \pm 0.14$, $i=23.91 \pm 0.18$, and $z=23.90 \pm 0.17$) to study the energy production and retention in the ejecta of SN 2011fe. Together with previous measurements by other groups, our result suggests that the optical supernova light curve can still be explained by the full thermalization of the decay positrons of Co56. This is in spite of theoretical predicted effects (e.g. infrared catastrophe, positron escape, and dust) that advocate a substantial energy redistribution and/or loss via various processes that result in a more rapid dimming at these very late epochs. △ Less

Submitted 12 November, 2014; v1 submitted 23 June, 2014; originally announced June 2014.

Comments: 8 pages, 3 figures, accepted by ApJ Letters

A spectral synthesis code for rapid modelling of supernovae

Authors: Wolfgang E. Kerzendorf, Stuart A. Sim

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

Submitted 21 January, 2014; originally announced January 2014.

Comments: 18 pages, 16 Figures, accepted for publication in MNRAS

A reconnaissance of the possible donor stars to the Kepler supernova

Authors: Wolfgang E. Kerzendorf, Michael Childress, Julia Scharwaechter, Tuan Do, Brian P. Schmidt

Abstract: The identity of Type Ia supernova progenitors remains a mystery, with various lines of evidence pointing towards either accretion from a non-degenerate companion, or the rapid merger of two degenerate stars leading to the thermonuclear destruction of a white dwarf. In this paper we spectroscopically scrutinize 24 of the brightest stars residing in the central 38" x 38" of the SN 1604 (Kepler) supe… ▽ More The identity of Type Ia supernova progenitors remains a mystery, with various lines of evidence pointing towards either accretion from a non-degenerate companion, or the rapid merger of two degenerate stars leading to the thermonuclear destruction of a white dwarf. In this paper we spectroscopically scrutinize 24 of the brightest stars residing in the central 38" x 38" of the SN 1604 (Kepler) supernova remnant to search for a possible surviving companion star. We can rule out, with high certainty, a red giant companion star - a progenitor indicated by some models of the supernova remnant. Furthermore, we find no star that exhibits properties uniquely consistent with those expected of a donor star down to L>10Lsun. While the distribution of star properties towards the remnant are consistent with unrelated stars, we identify the most promising candidates for further astrometric and spectroscopic follow-up. Such a program would either discover the donor star, or place strong limits on progenitor systems to luminosities with L<<Lsun. △ Less

Submitted 23 September, 2013; originally announced September 2013.

Comments: accepted by ApJ

Astropy: A Community Python Package for Astronomy

Authors: The Astropy Collaboration, Thomas P. Robitaille, Erik J. Tollerud, Perry Greenfield, Michael Droettboom, Erik Bray, Tom Aldcroft, Matt Davis, Adam Ginsburg, Adrian M. Price-Whelan, Wolfgang E. Kerzendorf, Alexander Conley, Neil Crighton, Kyle Barbary, Demitri Muna, Henry Ferguson, Frédéric Grollier, Madhura M. Parikh, Prasanth H. Nair, Hans M. Günther, Christoph Deil, Julien Woillez, Simon Conseil, Roban Kramer, James E. H. Turner , et al. (20 additional authors not shown)

Abstract: We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, ph… ▽ More We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions. △ Less

Submitted 23 July, 2013; originally announced July 2013.

Comments: Accepted for publication in Astronomy and Astrophysics. For more information about Astropy, visit http://www.astropy.org

The rotation of surviving companion stars after type Ia supernova explosions in the WD+MS scenario

Authors: Zheng-Wei Liu, Ruediger Pakmor, Friedrich K. Roepke, Philipp Edelmann, Wolfgang Hillebrandt, Wolfgang E. Kerzendorf, Bo Wang, Zhan-Wen Han

Abstract: In the SD scenario of SNe Ia the companion survives the SN explosion and thus should be visible near the center of the SN remnant and may show some unusual features. A promising approach to test progenitor models of SNe Ia is to search for the companion in SNRs. Here we present the results of 3D hydrodynamics simulations of the interaction between the SN Ia blast wave and a MS companion taking int… ▽ More In the SD scenario of SNe Ia the companion survives the SN explosion and thus should be visible near the center of the SN remnant and may show some unusual features. A promising approach to test progenitor models of SNe Ia is to search for the companion in SNRs. Here we present the results of 3D hydrodynamics simulations of the interaction between the SN Ia blast wave and a MS companion taking into consideration its orbital motion and spin. The primary goal of this work is to investigate the rotation of surviving companions after SN Ia explosions in the WD+MS scenario. We use Eggleton's code including the optically thick accretion wind model to obtain realistic models of companions. The impact of the SN blast wave on these companions is followed in 3D hydrodynamic simulations employing the SPH code GADGET3. We find that the rotation of the companion does not significantly affect the amount of stripped mass and the kick velocity caused by the SN impact. However, in our simulations, the rotational velocity of the companion is significantly reduced to about 14% to 32% of its pre-explosion value due to the expansion of the companion and the fact that 55%-89% of the initial angular momentum is carried away by the stripped matter. Compared with the observed rotational velocity of the presumed companion star of Tycho's SN, Tycho G, of 6 km/s the final rotational velocity we obtain is still higher by at least a factor of two. Whether this difference is significant, and may cast doubts on the suggestion that Tycho G is the companion of SN 1572, has to be investigated in future studies. Based on binary population synthesis results we present, for the first time, the expected distribution of rotational velocities of companions after the explosion which may provide useful information for the identification of the surviving companion in observational searches in other historical SNRs. △ Less

Submitted 11 March, 2013; originally announced March 2013.

Comments: 13 pages, 15 figures, accepted for publication by Astronomy and Astrophysics

A High-Resolution Spectroscopic Search for the Remaining Donor for Tycho's Supernova

Authors: Wolfgang E. Kerzendorf, David Yong, Brian P. Schmidt, Joshua D. Simon, C. Simon Jeffery, Jay Anderson, Philipp Podsiadlowski, Avishay Gal-Yam, Jeffrey M. Silverman, Alexei V. Filippenko, Ken'ichi Nomoto, Simon J. Murphy, Michael S. Bessell, Kim A. Venn, Ryan J. Foley

Abstract: In this paper, we report on our analysis using Hubble Space Telescope astrometry and Keck-I HIRES spectroscopy of the central six stars of Tycho's supernova remnant (SN 1572). With these data, we measured the proper motions, radial velocities, rotational velocities, and chemical abundances of these objects. Regarding the chemical abundances, we do not confirm the unusu- ally high [Ni/Fe] ratio pre… ▽ More In this paper, we report on our analysis using Hubble Space Telescope astrometry and Keck-I HIRES spectroscopy of the central six stars of Tycho's supernova remnant (SN 1572). With these data, we measured the proper motions, radial velocities, rotational velocities, and chemical abundances of these objects. Regarding the chemical abundances, we do not confirm the unusu- ally high [Ni/Fe] ratio previously reported for Tycho-G. Rather, we find that for all metrics in all stars, none exhibit the characteristics expected from traditional SN Ia single-degenerate-scenario calculations. The only possible exception is Tycho-B, a rare, metal-poor A-type star; however, we are unable to find a suitable scenario for it. Thus, we suggest that SN 1572 cannot be explained by the standard single-degenerate model. △ Less

Submitted 9 October, 2012; originally announced October 2012.

Comments: 34 pages, 11 Figures, revised and resubmitted to ApJ

Hunting for the progenitor of SN 1006: High resolution spectroscopic search with the FLAMES instrument

Authors: Wolfgang E. Kerzendorf, Brian P. Schmidt, John B. Laird, Philipp Podsiadlowski, Michael S. Bessell

Abstract: Type Ia supernovae play a significant role in the evolution of the Universe and have a wide range of applications. It is widely believed that these events are the thermonuclear explosions of carbon-oxygen white dwarfs close to the Chandrasekhar mass (1.38 M\odot). However, CO white dwarfs are born with masses much below the Chandrasekhar limit and thus require mass accretion to become Type Ia supe… ▽ More Type Ia supernovae play a significant role in the evolution of the Universe and have a wide range of applications. It is widely believed that these events are the thermonuclear explosions of carbon-oxygen white dwarfs close to the Chandrasekhar mass (1.38 M\odot). However, CO white dwarfs are born with masses much below the Chandrasekhar limit and thus require mass accretion to become Type Ia supernovae. There are two main scenarios for accretion. First, the merger of two white dwarfs and, second, a stable mass accretion from a companion star. According to predictions, this companion star (also referred to as donor star) survives the explosion and thus should be visible in the center of Type Ia remnants. In this paper we scrutinize the central stars (79 in total) of the SN 1006 remnant to search for the surviving donor star as predicted by this scenario. We find no star consistent with the traditional accretion scenario in SN1006. △ Less

Submitted 18 July, 2012; originally announced July 2012.

Comments: 11 pages, accepted by ApJ

Abundance stratification in Type Ia supernovae - III. The normal SN 2003du

Authors: Masaomi Tanaka, Paolo A. Mazzali, Vallery Stanishev, Immanuel Maurer, Wolfgang E. Kerzendorf, Ken'ichi Nomoto

Abstract: The element abundance distributions in the ejecta of Type Ia supernova (SN) is studied by modelling a time series of optical spectra of SN 2003du until ~1 year after the explosion. Since SN 2003du is a very normal Type Ia SN both photometrically and spectroscopically, the abundance distribution derived for it can be considered as representative of normal Type Ia SNe. We find that the innermost lay… ▽ More The element abundance distributions in the ejecta of Type Ia supernova (SN) is studied by modelling a time series of optical spectra of SN 2003du until ~1 year after the explosion. Since SN 2003du is a very normal Type Ia SN both photometrically and spectroscopically, the abundance distribution derived for it can be considered as representative of normal Type Ia SNe. We find that the innermost layers are dominated by stable Fe-group elements, with a total mass of ~ 0.2 Msun, which are synthesized through electron capture. Above the core of stable elements there are thick 56Ni-rich layers. The total mass of 56Ni is 0.65 Msun. The Si- and S-rich layers are located above the 56Ni-rich layers. The dominant element in the outermost layers (M_r > 1.1 Msun, v > 13000 km/s) is O, with a small amount of Si. Little unburned C remains, with an upper limit of 0.016 Msun. The element distributions in the ejecta are moderately mixed, but not fully mixed as seen in three-dimensional deflagration models. △ Less

Submitted 18 August, 2010; originally announced August 2010.

Comments: 15 pages, 11 figures, Accepted for publication in MNRAS

Subaru high-resolution spectroscopy of Star G in the Tycho supernova remnant

Authors: W. E. Kerzendorf, B. P. Schmidt, M. Asplund, K. Nomoto, Ph. Podsiadlowski, A. Frebel, R. A. Fesen, D. Yong

Abstract: It is widely believed that Type Ia supernovae (SN Ia) originate in binary systems where a white dwarf accretes material from a companion star until its mass approaches the Chandrasekhar mass and carbon is ignited in the white dwarf's core. This scenario predicts that the donor star should survive the supernova explosion, providing an opportunity to understand the progenitors of Type Ia supernova… ▽ More It is widely believed that Type Ia supernovae (SN Ia) originate in binary systems where a white dwarf accretes material from a companion star until its mass approaches the Chandrasekhar mass and carbon is ignited in the white dwarf's core. This scenario predicts that the donor star should survive the supernova explosion, providing an opportunity to understand the progenitors of Type Ia supernovae.In this paper we argue that rotation is a generic signature expected of most non-giant donor stars that is easily measurable. \citep{2004Natur.431.1069R} examined stars in the center of the remnant of SN 1572 (Tycho's SN) and showed evidence that a subgiant star (Star G by their naming convention) near the remnant's centre was the system's donor star. We present high-resolution (R \simeq 40000) spectra taken with the High Dispersion Spectrograph on Subaru of this candidate donor star and measure the star's radial velocity as $79\pm 2$ \kms with respect to the LSR and put an upper limit on the star's rotation of 7.5 \kms. In addition, by comparing images that were taken in 1970 and 2004, we measure the proper motion of Star G to be $μ_l = -1.6 \pm 2.1$ \masyr and $μ_b = -2.7 \pm 1.6$ \masyr. We demonstrate that all of the measured properties of Star G presented in this paper are consistent with those of a star in the direction of Tycho's SN that is not associated with the supernova event. However, we discuss an unlikely, but still viable scenario for Star G to be the donor star, and suggest further observations that might be able to confirm or refute it. △ Less

Submitted 4 June, 2009; originally announced June 2009.

Comments: Accepted for publication in ApJ

Journal ref: Astrophys.J.701:1665-1672,2009