Source author record

Daniel R. van Rossum

Daniel R. van Rossum appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

13 published item(s)

preprint2016arXiv

Fast evolving pair-instability supernova models: evolution, explosion, light curves

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

preprint2016arXiv

Light Curves and Spectra from a Thermonuclear Explosion of a White Dwarf Merger

Double-degenerate (DD) mergers of carbon-oxygen white dwarfs have recently emerged as a leading candidate for normal Type Ia supernovae (SNe Ia). However, many outstanding questions surround DD mergers, including the characteristics of their light curves and spectra. We have recently identified a spiral instability in the post-merger phase of DD mergers and demonstrated that this instability self-consistently leads to detonation in some cases. We call this the spiral merger SN Ia model. Here, we utilize the SuperNu radiative transfer software to calculate three-dimensional synthetic light curves and spectra of the spiral merger simulation with a system mass of 2.1 $M_\odot$ from Kashyap et al. Because of their large system masses, both violent and spiral merger light curves are slowly declining. The spiral merger resembles very slowly declining SNe Ia, including SN 2001ay, and provides a more natural explanation for its observed properties than other SN Ia explosion models. Previous synthetic light curves and spectra of violent DD mergers demonstrate a strong dependence on viewing angle, which is in conflict with observations. Here, we demonstrate that the light curves and spectra of the spiral merger are less sensitive to the viewing angle than violent mergers, in closer agreement with observation. We find that the spatial distribution of 56Ni and IMEs follows a characteristic hourglass shape. We discuss the implications of the asymmetric distribution of 56Ni for the early-time gamma-ray observations of 56Ni from SN 2014J. We suggest that DD mergers that agree with the light curves and spectra of normal SNe Ia will likely require a lower system mass.

preprint2016arXiv

On Measuring the Metallicity of a Type Ia Supernova's Progenitor

In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition scenario with different $^{56}$Ni yields and the W7 simulation to control for differences between explosion models and total yields. While the overall yields of intermediate mass elements (16 $<$ A $\leq$ 40) differ between the three cases, trends in the yields are similar. With increasing metallicity, $^{28}$Si yields remain nearly constant, $^{40}$Ca yields decline, and Ti and $^{54}$Fe yields increase. In the synthetic spectra, we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: a Ti feature around 4200\,Å and a Fe feature around 5200\,Å\@. In all three simulations, their pseudo equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two features may allow differentiation among progenitor metallicities of observed \sneia and potentially help reduce the intrinsic Hubble scatter.

preprint2014arXiv

Pan-chromatic observations of the remarkable nova LMC 2012

We present the results of an intensive multiwavelength campaign on nova LMC 2012. This nova evolved very rapidly in all observed wavelengths. The time to fall two magnitudes in the V band was only 2 days. In X-rays the super soft phase began 13$\pm$5 days after discovery and ended around day 50 after discovery. During the super soft phase, the \Swift/XRT and \Chandra\ spectra were consistent with the underlying white dwarf being very hot, $\sim$ 1 MK, and luminous, $\sim$ 10$^{38}$ erg s$^{-1}$. The UV, optical, and near-IR photometry showed a periodic variation after the initial and rapid fading had ended. Timing analysis revealed a consistent 19.24$\pm$0.03 hr period in all UV, optical, and near-IR bands with amplitudes of $\sim$ 0.3 magnitudes which we associate with the orbital period of the central binary. No periods were detected in the corresponding X-ray data sets. A moderately high inclination system, $i$ = 60$\pm$10$^{\arcdeg}$, was inferred from the early optical emission lines. The {\it HST}/STIS UV spectra were highly unusual with only the \ion{N}{5} (1240Å) line present and superposed on a blue continuum. The lack of emission lines and the observed UV and optical continua from four epochs can be fit with a low mass ejection event, $\sim$ 10$^{-6}$ M$_{\odot}$, from a hot and massive white dwarf near the Chandrasekhar limit. The white dwarf, in turn, significantly illuminated its subgiant companion which provided the bulk of the observed UV/optical continuum emission at the later dates. The inferred extreme white dwarf characteristics and low mass ejection event favor nova LMC 2012 being a recurrent nova of the U Sco subclass.

preprint2014arXiv

Radiation Transport for Explosive Outflows: Opacity Regrouping

Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC-DDMC. We explore a multigroup IMC-DDMC scheme that, in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure "opacity regrouping". Opacity regrouping has previously been applied to IMC-DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in \supernu, a code designed to do radiation transport in high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC-DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile-Fleck factor for modified IMC-DDMC. We test \supernu\ using numerical experiments including a quasi-manufactured analytic solution, a simple ten-group problem, and the W7 problem for Type Ia supernovae. We find that the opacity regrouping is necessary to make our IMC-DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the \supernu\ and \phoenix\ radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using \supernu\ is $\sim$10% less than that calculated using \phoenix.

preprint2014arXiv

The Nature of the Si II 6150A, Ca II HK, Ca II IR-triplet, and other Spectral Features in Supernova Type Ia Spectra

Spectra of Type Ia Supernovae (SN Ia) show both continuum-like and absorption-like line features. The pseudo equivalent width (pEW) and Doppler shift of absorption-line-like features, such as Si II 6150A, Ca II HK 3750A, and the Ca II IR-triplet 8150A quantities that are often associated with the optical depth and the velocity of a shell of interacting material, are important tools for interpreting and classifying the SN Ia spectra. In this paper, we examine the nature of spectral features in SN Ia spectra using W7 model spectra and a technique that we call "knock-out" spectra. We show that the P-Cyg profiles of Si II 6150A and many other features are largely emission dominated, rather than absorption dominated, and that the concepts of "absorption line" and "continuum" are therefore not adequate to describe the nature of spectral features in SNe Ia. Apparent absorption features, like Si II 6150A, are frequently just coincidental troughs between two (or more) uncorrelated emission features. In this situation, the pEW measured between these emission peaks is little related to the true strength of the presumed absorption feature. Furthermore, using the same knock-out technique, we demonstrate how spectral features overlap each other at different times after the explosion. This overlap distorts individual line profiles and affects measured absorption-line velocities. With synthetic SN Ia ejecta stratifications that are tuned to match specific observations, the method presented in this paper can in principle be used to quantify the effects of specific line opacities on measured line velocity, line strength, and line blending, and improve the interpretation and informative value of observed SN Ia spectra.

preprint2014arXiv

Three-dimensional Simulations of Pure Deflagration Models for Thermonuclear Supernovae

We present a systematic study of the pure deflagration model of Type Ia supernovae using three-dimensional, high-resolution, full-star hydrodynamical simulations, nucleosynthetic yields calculated using Lagrangian tracer particles, and light curves calculated using radiation transport. We evaluate the simulations by comparing their predicted light curves with many observed SNe Ia using the SALT2 data-driven model and find that the simulations may correspond to under-luminous SNe Iax. We explore the effects of the initial conditions on our results by varying the number of randomly selected ignition points from 63 to 3500, and the radius of the centered sphere they are confined in from 128 to 384 km. We find that the rate of nuclear burning depends on the number of ignition points at early times, the density of ignition points at intermediate times, and the radius of the confining sphere at late times. The results depend primarily on the number of ignition points, but we do not expect this to be the case in general. The simulations with few ignition points release more nuclear energy $E_{\mathrm{nuc}}$, have larger kinetic energies $E_{\mathrm{K}}$, and produce more $^{56}$Ni than those with many ignition points, and differ in the distribution of $^{56}$Ni, Si, and C/O in the ejecta. For these reasons, the simulations with few ignition points exhibit higher peak B-band absolute magnitudes $M_\mathrm{B}$ and light curves that rise and decline more quickly; their $M_\mathrm{B}$ and light curves resemble those of under-luminous SNe Iax, while those for simulations with many ignition points are not.

preprint2013arXiv

Comparing the light curves of simulated Type Ia Supernovae with observations using data-driven models

We propose a robust, quantitative method to compare the synthetic light curves of a Type Ia Supernova (SNIa) explosion model with a large set of observed SNeIa, and derive a figure of merit for the explosion model's agreement with observations. The synthetic light curves are fit with the data-driven model SALT2 which returns values for stretch, color, and magnitude at peak brightness, as well as a goodness-of-fit parameter. Each fit is performed multiple times with different choices of filter bands and epoch range in order to quantify the systematic uncertainty on the fitted parameters. We use a parametric population model for the distribution of observed SNIa parameters from large surveys, and extend it to represent red, dim, and bright outliers found in a low-redshift SNIa data set. We discuss the potential uncertainties of this population model and find it to be reliable given the current uncertainties on cosmological parameters. Using our population model, we assign each set of fitted parameters a likelihood of being observed in nature, and a figure of merit based on this likelihood. We define a second figure of merit based on the quality of the light curve fit, and combine the two measures into an overall figure of merit for each explosion model. We compute figures of merit for a variety of 1D, 2D and 3D explosion models and show that our evaluation method allows meaningful inferences across a wide range of light curve quality and fitted parameters.

preprint2013arXiv

Radiation Transport for Explosive Outflows: A Multigroup Hybrid Monte Carlo Method

We explore Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) for radiation transport in high-velocity outflows with structured opacity. The IMC method is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking Monte Carlo particles through optically thick materials. DDMC accelerates IMC in diffusive domains. Abdikamalov extended IMC and DDMC to multigroup, velocity-dependent transport with the intent of modeling neutrino dynamics in core-collapse supernovae. Densmore has also formulated a multifrequency extension to the originally grey DDMC method. We rigorously formulate IMC and DDMC over a high-velocity Lagrangian grid for possible application to photon transport in the post-explosion phase of Type Ia supernovae. This formulation includes an analysis that yields an additional factor in the standard IMC-to-DDMC spatial interface condition. To our knowledge the new boundary condition is distinct from others presented in prior DDMC literature. The method is suitable for a variety of opacity distributions and may be applied to semi-relativistic radiation transport in simple fluids and geometries. Additionally, we test the code, called SuperNu, using an analytic solution having static material, as well as with a manufactured solution for moving material with structured opacities. Finally, we demonstrate with a simple source and 10 group logarithmic wavelength grid that IMC-DDMC performs better than pure IMC in terms of accuracy and speed when there are large disparities between the magnitudes of opacities in adjacent groups. We also present and test our implementation of the new boundary condition.

preprint2012arXiv

A Public Set of Synthetic Spectra from Expanding Atmospheres for X-Ray Novae. I. Solar Abundances

X-ray grating observations have revealed great detail in the spectra of Novae in the Super Soft Source (SSS) phase. Notable features in the SSS spectra are blue-shifted absorption lines, P-Cygni line profiles, and the absence of strong ionization edges, all of which are indicators of an expanding atmosphere. We present, and make publicly available, a set of 672 wind-type (WT) synthetic spectra, obtained from the expanding NLTE SSS models introduced in Van Rossum 2010 with the PHOENIX stellar atmosphere code. The set presented in this paper is limited to solar abundances with the aim to focus on the basic model parameters and their effect on the spectra, providing the basis upon which abundance effects can be studied using a much bigger non-solar set in the next paper in this series. We fit the WT spectra to the five grating spectra taken in the SSS phase of nova V4743 Sgr 2003 as an example application of the WT models. Within the limits of solar abundances we demonstrate that the following parameters are constrained by the data (in order of decreasing accuracy): column density N_H, bolometric luminosity L_bol, effective temperature T_eff, white dwarf radius R, wind asymptotic velocity v_inf, and the mass-loss rate M_dot. The models are also sensitive to the assumed white dwarf mass M_wd but the effect on the spectra can largely be compensated by the other model parameters. The WT spectra with solar abundances fit the data better than abundance optimized hydro-static models.

preprint2012arXiv

Failed-Detonation Supernovae: Sub-Luminous Low-Velocity Ia Supernovae and Their Kicked Remnant White Dwarfs with Iron-Rich Cores

Type Ia supernovae (SNe Ia) originate from the thermonuclear explosions of carbon-oxygen (C-O) white dwarfs (WDs). The single-degenerate scenario is a well-explored model of SNe Ia where unstable thermonuclear burning initiates in an accreting, Chandrasekhar-mass WD and forms an advancing flame. By several proposed physical processes the rising, burning material triggers a detonation, which subsequently consumes and unbinds the WD. However, if a detonation is not triggered and the deflagration is too weak to unbind the star, a completely different scenario unfolds. We explore the failure of the Gravitationally-Confined Detonation (GCD) mechanism of SNe Ia, and demonstrate through 2D and 3D simulations the properties of failed-detonation SNe. We show that failed-detonation SNe expel a few 0.1 solar masses of burned and partially-burned material and that a fraction of the material falls back onto the WD, polluting the remnant WD with intermediate-mass and iron-group elements, that likely segregate to the core forming an WD whose core is iron rich. The remaining material is asymmetrically ejected at velocities comparable to the escape velocity from the WD, and in response, the WD is kicked to velocities of a few hundred km/s. These kicks may unbind the binary and eject a runaway/hyper-velocity WD. Although the energy and ejected mass of the failed-detonation SN are a fraction of typical thermonuclear SNe, they are likely to appear as sub-luminous low-velocity SNe Ia. Such failed detonations might therefore explain or are related to the observed branch of peculiar SNe Ia, such as the family of low-velocity sub-luminous SNe (SN 2002cx/SN 2008ha-like SNe).

preprint2012arXiv

Massive NLTE models for X-ray novae with PHOENIX

X-ray grating spectra provide the confirmation of continued mass loss from novae in the super-soft source (SSS) phase of the outburst. In this work expanding nova atmosphere models are developed and used to study the effect of mass loss on the SSS spectra. The very high temperatures combined with high expansion velocities and large radial extension make nova in the SSS phase very interesting but also difficult objects to model. The radiation transport code PHOENIX was applied to SSS novae before, but careful analysis of the old results has revealed a number of problems which lead to new methods and improvements to the code: 1) an improved NLTE module (a new opacity formalism, rate matrix solver, global iteration scheme, and temperature correction method); 2) a new hybrid hydrostatic-dynamic nova atmosphere setup; 3) the models are treated in pure NLTE (no LTE approximation for any opacity). With the new framework a modest amount of models (limited by computation time) are calculated. These show: 1) systematic behaviour for various atmospheric conditions, 2) the effect of expansion on the model spectrum is significant, and 3) the spectra are sensitive to the details of the atmospheric structure. The models are compared to the ten well-exposed grating spectra presently available: 5x V4743 Sgr, 3x RS Oph, and 2x V2491 Cyg. Although the models are on a coarse grid they do match the observations surprisingly well. Also, hydrostatic models are computed. The reproduction of the data is clearly inferior to the expanding models and, more importantly, their interpretation with hydrostatic models leads to conclusions opposite to those from expanding models. The models enable the derivation of accurate constraints on the physical conditions deep in the nova atmosphere that are revealed only in the SSS phase.

preprint2012arXiv

Radiation Energy-Balance Method for Calculating the Time Evolution of Type Ia Supernovae During the Post-Explosion Phase

A new method is presented for calculating the time evolution of spherically symmetric Type Ia Supernova in the post-explosion phase, enabling light curves and spectra to be simulated in a physically self-consistent way. The commonly exploited radiative equilibrium, that is in essence a /gas energy balance/ condition, is unsuitable for this purpose for important physical and numerical reasons. Firstly, the RE depends on the heating and cooling rates of the gas by the radiation field, two quantities that almost completely cancel and are very difficult to calculate accurately. Secondly, the internal energy of the gas is only a tiny fraction of the total energy in the system (the vast majority of the energy resides in the radiation field), so that the vast majority of the energy is neglected in solving for the energy balance. The method presented in this paper, based on the /radiation energy balance/, addresses the bulk of the energy, does not depend on the heating/cooling rates, guarantees an accurate run of the bolometric luminosity over time while bringing the gas temperatures into consistence with the radiation field. We have implemented the method in the stellar atmosphere code PHOENIX and applied it to the classical W7 model. The results illustrate the importance of each of the four physical contributions to the energy balance as a function of time. The simulated spectra and light curves for W7 show good resemblance to the observations, which demonstrates what can be done using PHOENIX with the REB method.