Showing 1–42 of 42 results for author: Even, W

Halted-Pendulum Relaxation: Application to White Dwarf Binary Initial Data

Authors: M. Alexander R. Kaltenborn, Michael J. Falato, Oleg Korobkin, Irina Sagert, Wesley P. Even

Abstract: Studying compact star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, Smoothed Particle Hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacu… ▽ More Studying compact star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, Smoothed Particle Hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, Halted-Pendulum Relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods. △ Less

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

Report number: LA-UR-22-29924

Abundances and Transients from Neutron Star-White Dwarf Mergers

Authors: Mark Alexander Kaltenborn, Chris L. Fryer, Ryan T. Wollaeger, Krzysztof Belczynski, Wesley Even, Chryssa Kouveliotou

Abstract: We systematically investigate the mergers of neutron star-white dwarf binaries from beginning to end, with focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a l… ▽ More We systematically investigate the mergers of neutron star-white dwarf binaries from beginning to end, with focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a large accretion disk around the neutron star. The material evolves according to an advection-dominated accretion disk model with nuclear burning, neutrino-emissions, and disk-surface wind ejection. The extreme dynamics of the entire process has proven difficult for analytic analysis, and thus currently the properties are poorly understood. The outflows from the mergers are iron- and nickel-rich, giving rise to optical and infrared emissions powered from the decay of the radioactive iron-type isotopes, calculated via the SuperNu light-curve code. We find these systems capable of powering bright, yet short-lived, optical transients with the potential to power gamma-ray bursts. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: LA-UR-22-29879

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

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

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

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

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

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

Conservation of Angular Momentum in the Fast Multipole Method

Authors: Oleg Korobkin, Hyun Lim, Irina Sagert, Julien Loiseau, Christopher Mauney, M. Alexander R. Kaltenborn, Bing-Jyun Tsao, Wesley P. Even

Abstract: Smoothed particle hydrodynamics (SPH) is positioned as having ideal conservation properties. When properly implemented, conservation of total mass, energy, and both linear and angular momentum is guaranteed exactly, up to machine precision. This is particularly important for some applications in computational astrophysics, such as binary dynamics, mergers, and accretion of compact objects (neutron… ▽ More Smoothed particle hydrodynamics (SPH) is positioned as having ideal conservation properties. When properly implemented, conservation of total mass, energy, and both linear and angular momentum is guaranteed exactly, up to machine precision. This is particularly important for some applications in computational astrophysics, such as binary dynamics, mergers, and accretion of compact objects (neutron stars, black holes, and white dwarfs). However, in astrophysical applications that require the inclusion of gravity, calculating pairwise particle interactions becomes prohibitively expensive. In the Fast Multipole Method (FMM), they are, therefore, replaced with symmetric interactions between distant clusters of particles (contained in the tree nodes) Although such an algorithm is linear momentum-conserving, it introduces spurious torques that violate conservation of angular momentum. We present a modification of FMM that is free of spurious torques and conserves angular momentum explicitly. The new method has practically no computational overhead compared to the standard FMM. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 6 pages, 3 figures, proceedings of 2021 international SPHERIC workshop (virtual), June, 8-11 2021

Report number: LA-UR-21-24198

FleCSPH: The Next Generation FleCSIble Parallel Computational Infrastructure for Smoothed Particle Hydrodynamics

Authors: Julien Loiseau, Hyun Lim, Mark Alexander Kaltenborn, Oleg Korobkin, Christopher M. Mauney, Irina Sagert, Wesley P. Even, Benjamin K. Bergen

Abstract: FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution driv… ▽ More FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution drivers, which can be easily modified and extended to user-specific setups. Data input/output uses the H5part format, compatible with modern visualization software. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Axisymmetric Radiative Transfer Models of Kilonovae

Authors: Oleg Korobkin, Ryan Wollaeger, Christopher Fryer, Aimee L. Hungerford, Stephan Rosswog, Christopher Fontes, Matthew Mumpower, Eve Chase, Wesley Even, Jonah M. Miller, G. Wendell Misch, Jonas Lippuner

Abstract: The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g. in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numeric… ▽ More The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g. in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining," which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent {P Cygni} features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances. △ Less

Submitted 4 March, 2021; v1 submitted 31 March, 2020; originally announced April 2020.

Comments: 23 pages, 18 figures; accepted to ApJ

Report number: LA-UR-20-22607

The Nucleosynthetic Yields of Core-Collapse Supernovae, prospects for the Next Generation of Gamma-Ray Astronomy

Authors: Sydney Andrews, Chris L. Fryer, Samuel W. Jones, Wesley P. Even, Marco Pignatari

Abstract: Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key… ▽ More Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. Post-bounce temperature and density evolution of ZAMS 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid (NuGrid) nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio-isotopes that are of particular interest to the next generation of gamma-ray astronomical observations; 43K, 47Ca, 44Sc, 47Sc, 48V, 48Cr, 51Cr, 52Mn, 59Fe, 56Co, 57Co, 57Ni. These nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43K, 47Ca, 44Sc, 47Sc, and 59Fe, those that are dependent on the progenitor structure are 48V, 51Cr, and 57Co, and those that probe neither are 48Cr, 52Mn, 57Ni, and 56Co. We discuss prospects of observing these radionuclides in supernova remnants. △ Less

Submitted 22 December, 2019; originally announced December 2019.

Comments: 17 pages, 19 figures, submitted to ApJ

Gamma-rays from kilonova: a potential probe of r-process nucleosynthesis

Authors: Oleg Korobkin, Aimee M. Hungerford, Christopher L. Fryer, Matthew R. Mumpower, G. Wendell Misch, Trevor M. Sprouse, Jonas Lippuner, Rebecca Surman, Aaron J. Couture, Peter F. Bloser, Farzane Shirazi, Wesley P. Even, W. Thomas Vestrand, Richard S. Miller

Abstract: The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of $r$-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron s… ▽ More The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of $r$-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. We show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra. △ Less

Submitted 9 April, 2020; v1 submitted 13 May, 2019; originally announced May 2019.

Comments: 14 pages, 11 figures

Report number: LA-UR-19-24338

Journal ref: The Astrophysical Journal, Volume 889, Issue 2, id.168, 2020

Composition Effects on Kilonova Spectra and Light Curves: I

Authors: Wesley Even, Oleg Korobkin, Christopher L. Fryer, Christopher J. Fontes, R. T. Wollaeger, Aimee Hungerford, Jonas Lippuner, Jonah Miller, Matthew R. Mumpower, G. Wendell Misch

Abstract: The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opa… ▽ More The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infra-red wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infra-red in lanthanide opacities have been shown to significantly alter the light-curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. Here we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 micron are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2-1 micron wavelength range and kilonova observations may also be able to constrain these abundances. △ Less

Submitted 30 April, 2019; originally announced April 2019.

Impact of Pulsar and Fallback Sources on Multifrequency Kilonova Models

Authors: Ryan T. Wollaeger, Chris L. Fryer, Christopher J. Fontes, Jonas Lippuner, W. Thomas Vestrand, Matthew R. Mumpower, Oleg Korobkin, Aimee L. Hungerford, Wesley P. Even

Abstract: We explore the impact of pulsar electromagnetic dipole and fallback accretion emission on the luminosity of a suite of kilonova models. The pulsar models are varied over pulsar magnetic field strength, pulsar lifetime, ejecta mass, and elemental abundances; the fallback models are varied over fallback accretion rate and ejecta mass. For the abundances, we use Fe and Nd as representatives of the wi… ▽ More We explore the impact of pulsar electromagnetic dipole and fallback accretion emission on the luminosity of a suite of kilonova models. The pulsar models are varied over pulsar magnetic field strength, pulsar lifetime, ejecta mass, and elemental abundances; the fallback models are varied over fallback accretion rate and ejecta mass. For the abundances, we use Fe and Nd as representatives of the wind and dynamical ejecta, respectively. We simulate radiative transfer in the ejecta in either 1D spherical or 2D cylindrical spatial geometry. For the grid of 1D simulations, the mass fraction of Nd is 0, $10^{-4}$, or $10^{-3}$ and the rest is Fe. Our models that fit the bolometric luminosity of AT 2017gfo (the kilonova associated with the first neutron star merger discovered in gravitational waves, GW170817) do not simultaneously fit the B, V, and I time evolution. However, we find that the trends of the evolution in B and V magnitudes are better matched by the fallback model relative to the pulsar model, implying the time dependence of the remnant source influences the color evolution. Further exploration of the parameter space and model deficiencies is needed before we can describe AT 2017gfo with a remnant source. △ Less

Submitted 31 May, 2019; v1 submitted 11 April, 2019; originally announced April 2019.

Comments: 15 pages, 10 figures, accepted to ApJ. Simulation data available (under "Source term study") at: https://ccsweb.lanl.gov/astro/transient/transients_astro.html

Report number: LA-UR-19-22590 (manuscript), LA-UR-19-22878 (data)

Synthetic Spectra of Pair-Instability Supernovae in 3D

Authors: E. Chatzopoulos, M. S. Gilmer, R. T. Wollaeger, C. Frohlich, W. P. Even

Abstract: Pair-Instability Supernovae (PISNe) may signal the deaths of extremely massive stars in the local Universe or massive primordial stars after the end of the Cosmic Dark Ages. Hydrodynamic simulations of these explosions, performed in 1D, 2D, and 3D geometry, have revealed the strong dependence of mixing in the PISN ejecta on dimensionality. This chemical rearrangement is mainly driven by Rayleigh-T… ▽ More Pair-Instability Supernovae (PISNe) may signal the deaths of extremely massive stars in the local Universe or massive primordial stars after the end of the Cosmic Dark Ages. Hydrodynamic simulations of these explosions, performed in 1D, 2D, and 3D geometry, have revealed the strong dependence of mixing in the PISN ejecta on dimensionality. This chemical rearrangement is mainly driven by Rayleigh-Taylor instabilities that start to grow shortly after the collapse of the carbon-oxygen core. We investigate the effects of such mixing on the spectroscopic evolution of PISNe by post-processing explosion profiles with the radiation diffusion-equilibrium code SNEC and the implicit Monte Carlo-discrete diffusion Monte Carlo (IMC-DDMC) radiation transport code SuperNu. The first 3D radiation transport calculation of a PISN explosion is presented yielding viewing angle-dependent synthetic spectra and lightcurves. We find that while 2D and 3D mixing does not significantly affect the lightcurves of PISNe, their spectroscopic and color evolution is impacted. Strong features of intermediate mass elements dominated by silicon, magnesium and oxygen appear at different phases and reach different intensities depending on the extent of mixing in the silicon/oxygen interface of the PISN ejecta. On the other hand, we do not find a significant dependence of PISN lightcurves and spectra on viewing angle. Our results showcase the capabilities of SuperNu to handle 3D radiation transport and highlight the importance of modeling time-series of spectra in identifying PISNe with future missions. △ Less

Submitted 18 February, 2019; originally announced February 2019.

Comments: 15 pages, 10 figures

$^{60}$Fe in core-collapse supernovae and prospects for X-ray and $γ$-ray detection in supernova remnants

Authors: Samuel Jones, Heiko Moeller, Chris L. Fryer, Christopher J. Fontes, Reto Trappitsch, Wesley P. Even, Aaron Couture, Matthew R. Mumpower, Samar Safi-Harb

Abstract: We investigate $^{60}$Fe in massive stars and core-collapse supernovae focussing on uncertainties that influence its production in 15, 20 and 25 $M_\odot$ stars at solar metallicity. We find that the $^{60}$Fe yield is a monotonic increasing function of the uncertain $^{59}$Fe$(n,γ)^{60}$Fe cross section and that a factor of 10 reduction in the reaction rate results in a factor 8-10 reduction in t… ▽ More We investigate $^{60}$Fe in massive stars and core-collapse supernovae focussing on uncertainties that influence its production in 15, 20 and 25 $M_\odot$ stars at solar metallicity. We find that the $^{60}$Fe yield is a monotonic increasing function of the uncertain $^{59}$Fe$(n,γ)^{60}$Fe cross section and that a factor of 10 reduction in the reaction rate results in a factor 8-10 reduction in the $^{60}$Fe yield; while a factor of 10 increase in the rate increases the yield by a factor 4-7. We find that none of the 189 simulations we have performed are consistent with a core-collapse supernova triggering the formation of the Solar System, and that only models using $^{59}$Fe$(n,γ)^{60}$Fe cross section that is less than or equal to that from NON-SMOKER can reproduce the observed $^{60}$Fe/$^{26}$Al line flux ratio in the diffuse ISM. We examine the prospects of detecting old core-collapse supernova remnants (SNRs) in the Milky Way from their $γ$-ray emission from the decay of $^{60}$Fe, finding that the next generation of gamma-ray missions could be able to discover up to $\sim100$ such old SNRs as well as measure the $^{60}$Fe yields of a handful of known Galactic SNRs. We also predict the X-ray spectrum that is produced by atomic transitions in $^{60}$Co following its ionization by internal conversion and give theoretical X-ray line fluxes as a function of remnant age as well as the Doppler and fine-structure line broadening effects. The X-ray emission presents an interesting prospect for addressing the missing SNR problem with future X-ray missions. △ Less

Submitted 15 February, 2019; originally announced February 2019.

Comments: 25 pages, 25 figures, 1 table. Accepted for publication in MNRAS

Californium-254 and kilonova light curves

Authors: Y. Zhu, R. T. Wollaeger, N. Vassh, R. Surman, T. M. Sprouse, M. R. Mumpower, P. Moller, G. C. McLaughlin, O. Korobkin, T. Kawano, P. J. Jaffke, E. M. Holmbeck, C. L. Fryer, W. P. Even, A. J. Couture, J. Barnes

Abstract: Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on t… ▽ More Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on the order of 15 to 250 days. This is due to the anomalously long half-life of this isotope and the efficiency of fission thermalization compared to other nuclear channels. We estimate the fission fragment yield of this nucleus and outline the astrophysical conditions under which $^{254}$Cf has the greatest impact to the light curve. Future observations in the middle-IR which are bright during this regime could indicate the production of actinide nucleosynthesis. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: 9 pages, 5 figures

Report number: LA-UR-18-25559

The role of dredge-up in double white dwarf mergers

Authors: Jan E. Staff, Brandon Wiggins, Dominic Marcello, Patrick M. Motl, Wesley Even, Christopher L. Fryer, Cody Raskin, Geoffrey C. Clayton, Juhan Frank

Abstract: We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two white dwarfs is $\lesssim1.0~{\rm M_\odot}$. Mergers involving a CO and a He white dwarf have been suggested… ▽ More We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two white dwarfs is $\lesssim1.0~{\rm M_\odot}$. Mergers involving a CO and a He white dwarf have been suggested as a possible formation channel for R Coronae Borealis type stars, and we are interested in testing if such mergers lead to conditions and outcomes in agreement with observations. Even if the conditions during the merger and subsequent nucleosynthesis favor the production of $^{18}{\mathrm O}$, the merger must avoid dredging up large amounts of $^{16}{\mathrm O}$, or else it will be difficult to produce sufficient $^{18}{\mathrm O}$ to explain the oxygen ratio observed to be of order unity. We performed a total of 9 simulations using two different grid-based hydrodynamics codes using fixed and adaptive meshes, and one smooth particle hydrodynamics (SPH) code. We find that in most of the simulations, $>10^{-2}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ is indeed dredged up during the merger. However, in SPH simulations where the accretor is a hybrid He/CO white dwarf with a $\sim 0.1~{\rm M_\odot}$ layer of helium on top, we find that no $^{16}{\mathrm O}$ is being dredged up, while in the $q=0.8$ simulation $<10^{-4}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ has been brought up, making a WD binary consisting of a hybrid CO/He WD and a companion He WD an excellent candidate for the progenitor of RCB stars. △ Less

Submitted 1 June, 2018; originally announced June 2018.

Comments: Accepted for publication in ApJ

Parameterizing the Supernova Engine and its Effects on Remnants and Basic Yields

Authors: Chris L. Fryer, Sydney Andrews, Wesley Even, Alex Heger, Samar Safi-Harb

Abstract: Core-collapse supernova science is now entering an era where engine models are beginning to make both qualitative and quantitative predictions. Although the evidence in support of the convective engine for core-collapse supernova continues to grow, it is difficult to place quantitative constraints on this engine. Some studies have made specific predictions for the remnant distribution from the con… ▽ More Core-collapse supernova science is now entering an era where engine models are beginning to make both qualitative and quantitative predictions. Although the evidence in support of the convective engine for core-collapse supernova continues to grow, it is difficult to place quantitative constraints on this engine. Some studies have made specific predictions for the remnant distribution from the convective engine, but the results differ between different groups. Here we use a broad parameterization for the supernova engine to understand the differences between distinct studies. With this broader set of models, we place error bars on the remnant mass and basic yields from the uncertainties in the explosive engine. We find that, even with only 3 progenitors and a narrow range of explosion energies, we can produce a wide range of remnant masses and nucleosynthetic yields. △ Less

Submitted 9 December, 2017; originally announced December 2017.

Comments: 14 pages, 10 figures, submitted to ApJ

The Emergence of a Lanthanide-Rich Kilonova Following the Merger of Two Neutron Stars

Authors: N. R. Tanvir, A. J. Levan, C. Gonzalez-Fernandez, O. Korobkin, I. Mandel, S. Rosswog, J. Hjorth, P. D'Avanzo, A. S. Fruchter, C. L. Fryer, T. Kangas, B. Milvang-Jensen, S. Rosetti, D. Steeghs, R. T. Wollaeger, Z. Cano, C. M. Copperwheat, S. Covino, V. D'Elia, A. de Ugarte Postigo, P. A. Evans, W. P. Even, S. Fairhurst, R. Figuera Jaimes, C. J. Fontes , et al. (29 additional authors not shown)

Abstract: We report the discovery and monitoring of the near-infrared counterpart (AT2017gfo) of a binary neutron-star merger event detected as a gravitational wave source by Advanced LIGO/Virgo (GW170817) and as a short gamma-ray burst by Fermi/GBM and Integral/SPI-ACS (GRB170817A). The evolution of the transient light is consistent with predictions for the behaviour of a "kilonova/macronova", powered by t… ▽ More We report the discovery and monitoring of the near-infrared counterpart (AT2017gfo) of a binary neutron-star merger event detected as a gravitational wave source by Advanced LIGO/Virgo (GW170817) and as a short gamma-ray burst by Fermi/GBM and Integral/SPI-ACS (GRB170817A). The evolution of the transient light is consistent with predictions for the behaviour of a "kilonova/macronova", powered by the radioactive decay of massive neutron-rich nuclides created via r-process nucleosynthesis in the neutron-star ejecta. In particular, evidence for this scenario is found from broad features seen in Hubble Space Telescope infrared spectroscopy, similar to those predicted for lanthanide dominated ejecta, and the much slower evolution in the near-infrared Ks-band compared to the optical. This indicates that the late-time light is dominated by high-opacity lanthanide-rich ejecta, suggesting nucleosynthesis to the 3rd r-process peak (atomic masses A~195). This discovery confirms that neutron-star mergers produce kilo-/macronovae and that they are at least a major - if not the dominant - site of rapid neutron capture nucleosynthesis in the universe. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Swift and NuSTAR observations of GW170817: detection of a blue kilonova

Authors: P. A. Evans, S. B. Cenko, J. A. Kennea, S. W. K. Emery, N. P. M. Kuin, O. Korobkin, R. T. Wollaeger, C. L. Fryer, K. K. Madsen, F. A. Harrison, Y. Xu, E. Nakar, K. Hotokezaka, A. Lien, S. Campana, S. R. Oates, E. Troja, A. A. Breeveld, F. E. Marshall, S. D. Barthelmy, A. P. Beardmore, D. N. Burrows, G. Cusumano, A. D'Ai, P. D'Avanzo , et al. (34 additional authors not shown)

Abstract: With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star… ▽ More With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow). △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: Science, in press; 56 pages, 12 figures

Light Curves and Spectra from a Unimodal Core-Collapse SuperNova

Authors: Ryan T. Wollaeger, Aimee L. Hungerford, Chris L. Fryer, Allan B. Wollaber, Daniel R. van Rossum, Wesley Even

Abstract: To assess the effectiveness of optical emission as a probe of spatial asymmetry in core-collapse supernovae (CCSNe), we apply the radiative transfer software, \supernu, to a unimodal CCSN model. The \snsph\ radiation-hydrodynamics software was used to simulate an asymmetric explosion of a 16 M$_{\odot}$ ZAMS binary star. The ejecta has 3.36 M$_{\odot}$ with 0.024 M$_{\odot}$ of radioactive… ▽ More To assess the effectiveness of optical emission as a probe of spatial asymmetry in core-collapse supernovae (CCSNe), we apply the radiative transfer software, \supernu, to a unimodal CCSN model. The \snsph\ radiation-hydrodynamics software was used to simulate an asymmetric explosion of a 16 M$_{\odot}$ ZAMS binary star. The ejecta has 3.36 M$_{\odot}$ with 0.024 M$_{\odot}$ of radioactive $^{56}$Ni, with unipolar asymmetry along the z-axis. For 96 discrete angular views, we find the ratio between maximum and minimum peak total luminosities is $\sim$1.36. The brightest light curves emerge from views orthogonal to the z-axis. Multigroup spectra from UV to IR are obtained. We find a shift in wavelength with viewing angle in a near-IR Ca II emission feature, consistent with Ca being mostly in the unimode. We compare emission from the grey gamma-ray transfer in \supernu\ and from the detailed gamma-ray transfer code \maverick. Relative to the optical light curves, the brightness of the gamma-ray emission is more monotonic with respect to viewing angle. UBVRI broad-band light curves are also calculated. Parallel with the unimode, the U and B bands have excess luminosity at $\gtrsim 10$ days post-explosion, due to $^{56}$Ni on the unimode. We compare our CCSN model with SN 2002ap, which is thought to have a similar ejecta morphology. △ Less

Submitted 3 August, 2017; originally announced August 2017.

Comments: Accepted to ApJ - 18 pages, 17 figures

Report number: LANL Report Number: LA-UR-16-20746

Impact of ejecta morphology and composition on the electromagnetic signatures of neutron star mergers

Authors: Ryan T. Wollaeger, Oleg Korobkin, Christopher J. Fontes, Stephan K. Rosswog, Wesley P. Even, Christopher L. Fryer, Jesper Sollerman, Aimee L. Hungerford, Daniel R. van Rossum, Allan B. Wollaber

Abstract: The electromagnetic transients accompanying compact binary mergers ($γ$-ray bursts, afterglows and 'macronovae') are crucial to pinpoint the sky location of gravitational wave sources. Macronovae are caused by the radioactivity from freshly synthesised heavy elements, e.g. from dynamic ejecta and various types of winds. We study macronova signatures by using multi-dimensional radiative transfer ca… ▽ More The electromagnetic transients accompanying compact binary mergers ($γ$-ray bursts, afterglows and 'macronovae') are crucial to pinpoint the sky location of gravitational wave sources. Macronovae are caused by the radioactivity from freshly synthesised heavy elements, e.g. from dynamic ejecta and various types of winds. We study macronova signatures by using multi-dimensional radiative transfer calculations. We employ the radiative transfer code SuperNu and state-of-the-art LTE opacities for a few representative elements from the wind and dynamical ejecta (Cr, Pd, Se, Te, Br, Zr, Sm, Ce, Nd, U) to calculate synthetic light curves and spectra for a range of ejecta morphologies. The radioactive power of the resulting macronova is calculated with the detailed input of decay products. We assess the detection prospects for our most complex models, based on the portion of viewing angles that are sufficiently bright, at different cosmological redshifts ($z$). The brighter emission from the wind is unobscured by the lanthanides (or actinides) in some of the models, permitting non-zero detection probabilities for redshifts up to $z=0.07$. We also find the nuclear mass model and the resulting radioactive heating rate are crucial for the detectability. While for the most pessimistic heating rate (from the FRDM model) no reasonable increase in the ejecta mass or velocity, or wind mass or velocity, can possibly make the light curves agree with the observed nIR excess after GRB130603B, a more optimistic heating rate (from the Duflo-Zuker model) leads to good agreement. We conclude that future reliable macronova observations would constrain nuclear heating rates, and consequently help constrain nuclear mass models. △ Less

Submitted 19 May, 2017; originally announced May 2017.

Comments: 34 pages, submitted to MNRAS

Report number: LA-UR-17-24109

A Comparison of Grid-Based and SPH Binary Mass-Transfer and Merger Simulations

Authors: Patrick M Motl, Juhan Frank, Jan Staff, Geoffrey C Clayton, Christopher L Fryer, Wesley Even, Steven Diehl, Joel E Tohline

Abstract: Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different nu… ▽ More Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms --- a finite-volume "grid" code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations an atmosphere forms around the accretor which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code. △ Less

Submitted 12 February, 2017; originally announced February 2017.

Comments: 55 pages, 41 figures, accepted for publication in ApJS

The Supernovae Analysis Application (SNAP)

Authors: Amanda J. Bayless, Chris L. Fryer, Brandon Wiggins, Wesley Even, Ryan Wollaeger, Janie de la Rosa, Peter W. A. Roming, Lucy Frey, Patrick A. Young, Rob Thorpe, Luke Powell, Rachel Landers, Heather D. Persson, Rebecca Hay

Abstract: The SuperNovae Analysis aPplication (SNAP) is a new tool for the analysis of SN observations and validation of SN models. SNAP consists of an open source relational database with (a) observational light curve, (b) theoretical light curve, and (c) correlation table sets, statistical comparison software, and a web interface available to the community. The theoretical models are intended to span a gr… ▽ More The SuperNovae Analysis aPplication (SNAP) is a new tool for the analysis of SN observations and validation of SN models. SNAP consists of an open source relational database with (a) observational light curve, (b) theoretical light curve, and (c) correlation table sets, statistical comparison software, and a web interface available to the community. The theoretical models are intended to span a gridded range of parameter space. The goal is to have users to upload new SN models or new SN observations and run the comparison software to determine correlations via the web site. There are looming problems on the horizon that SNAP begins to solve. Namely, large surveys will discover thousands of SNe annually. Frequently, the parameter space of a new SN event is unbounded. SNAP will be a resource to constrain parameters and determine if an event needs follow-up without spending resources to create new light curve models from scratch. Secondly, there is not a rapidly available, systematic way to determine degeneracies between parameters or even what physics is needed to model a realistic SNe. The correlations made within the SNAP system begin to solve these problems. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: Submitted to AAS publishing, 22 pages, 8 figures

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

Authors: Alexandra Kozyreva, Matthew Gilmer, Raphael Hirschi, Carla Frohlich, Sergey Blinnikov, Ryan T. Wollaeger, Ulrich M. Noebauer, Daniel R. van Rossum, Alexander Heger, Wesley P. Even, Roni Waldman, Alexey Tolstov, Emmanouil Chatzopoulos, Elena Sorokina

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

Submitted 4 October, 2016; originally announced October 2016.

Comments: 14 pages, 1 table, 19 figures. Accepted for publication in MNRAS

Extreme Supernova Models for the Superluminous Transient ASASSN-15lh

Authors: E. Chatzopoulos, J. C. Wheeler, J. Vinko, A. P. Nagy, B. K. Wiggins, W. P. Even

Abstract: The recent discovery of the unprecedentedly superluminous transient ASASSN-15lh (or SN 2015L) with its UV-bright secondary peak challenges all the power-input models that have been proposed for superluminous supernovae. Here we examine some of the few viable interpretations of ASASSN-15lh in the context of a stellar explosion, involving combinations of one or more power inputs. We model the lightc… ▽ More The recent discovery of the unprecedentedly superluminous transient ASASSN-15lh (or SN 2015L) with its UV-bright secondary peak challenges all the power-input models that have been proposed for superluminous supernovae. Here we examine some of the few viable interpretations of ASASSN-15lh in the context of a stellar explosion, involving combinations of one or more power inputs. We model the lightcurve of ASASSN-15lh with a hybrid model that includes contributions from magnetar spin-down energy and hydrogen-poor circumstellar interaction. We also investigate models of pure circumstellar interaction with a massive hydrogen-deficient shell and discuss the lack of interaction features in the observed spectra. We find that, as a supernova ASASSN-15lh can be best modeled by the energetic core-collapse of a ~40 Msun star interacting with a hydrogen-poor shell of ~20 Msun. The circumstellar shell and progenitor mass are consistent with a rapidly rotating pulsational pair-instability supernova progenitor as required for strong interaction following the final supernova explosion. Additional energy injection by a magnetar with initial period of 1-2 ms and magnetic field of 0.1-1 x 10^14 G may supply the excess luminosity required to overcome the deficit in single-component models, but this requires more fine-tuning and extreme parameters for the magnetar, as well as the assumption of efficient conversion of magnetar energy into radiation. We thus favor a single-input model where the reverse shock formed in a strong SN ejecta-CSM interaction following a very powerful core-collapse SN explosion can supply the luminosity needed to reproduce the late-time UV-bright plateau. △ Less

Submitted 13 June, 2016; v1 submitted 22 March, 2016; originally announced March 2016.

Comments: 8 pages, 3 figures

Assessing the Observability of Hypernovae and Pair-Instability Supernovae in the Early Universe

Authors: Brandon K. Wiggins, Joseph M. Smidt, Daniel J. Whalen, Wesley P. Even, Victor Migenes, Chris L. Fryer

Abstract: The era of the universe's first (Population III) stars is essentially unconstrained by observation. Ultra-luminous and massive stars from this time altered the chemistry of the cosmos, provided the radiative scaffolding to support the formation of the first protogalaxies, and facilitated the creation and growth of now-supermassive black holes. Unfortunately, because these stars lie literally at th… ▽ More The era of the universe's first (Population III) stars is essentially unconstrained by observation. Ultra-luminous and massive stars from this time altered the chemistry of the cosmos, provided the radiative scaffolding to support the formation of the first protogalaxies, and facilitated the creation and growth of now-supermassive black holes. Unfortunately, because these stars lie literally at the edge of the observable universe, they will remain beyond the reach of even the next generation of telescopes such as the James Webb Space Telescope and the Thirty-Meter Telescope. In this paper, we provide a primer to supernovae modeling and the first stars to make our discussion accessible to those new to or outside our field. We review recent work of the Los Alamos Supernova Light Curve Project and Brigham Young University to explore the possibility of probing this era through observations of the spectacular deaths of the first stars. We find that many such brilliant supernova explosions will be observable as far back as $\sim 99$% of the universe's current age, tracing primordial star formation rates and the locations of their protogalaxies on the sky. The observation of Population III supernovae will be among the most spectacular discoveries in observational astronomy in the coming decade. △ Less

Submitted 27 March, 2015; originally announced March 2015.

Report number: LA-UR-14-24307

Finding the First Cosmic Explosions. IV. 90 - 140 M$_{\odot}$ Pair-Instability Supernovae

Authors: Joseph Smidt, Daniel J. Whalen, E. Chatzopoulos, Brandon K. Wiggins, Ke-Jung Chen, Alexandra Kozyreva, Wesley Even

Abstract: Population III stars that die as pair-instability supernovae are usually thought to fall in the mass range of 140 - 260 M$_{\odot}$. But several lines of work have now shown that rotation can build up the He cores needed to encounter the pair instability at stellar masses as low as 90 $_{\odot}$. Depending on the slope of the initial mass function of Population III stars, there could be 4 - 5 time… ▽ More Population III stars that die as pair-instability supernovae are usually thought to fall in the mass range of 140 - 260 M$_{\odot}$. But several lines of work have now shown that rotation can build up the He cores needed to encounter the pair instability at stellar masses as low as 90 $_{\odot}$. Depending on the slope of the initial mass function of Population III stars, there could be 4 - 5 times as many stars from 90 - 140 $_{\odot}$ in the primordial universe than in the usually accepted range. We present numerical simulations of the pair-instability explosions of such stars performed with the MESA, FLASH and RAGE codes. We find that they will be visible to supernova factories such as Pan-STARRS and LSST in the optical out to z $\sim$ 1 - 2 and to JWST and the 30 m-class telescopes in the NIR out to $z \sim$ 7 - 10. Such explosions will thus probe the stellar populations of the first galaxies and cosmic star formation rates in the era of cosmological reionization. These supernovae are also easily distinguished from more massive pair-instability explosions, underscoring the fact that there is far greater variety to the light curves of these events than previously understood. △ Less

Submitted 19 November, 2014; originally announced November 2014.

Comments: 12 pages. 8 figures

Report number: LA-UR-14-28622

Observational Constraints of Stellar Collapse: Diagnostic Probes of Nature's Extreme Matter Experiment

Authors: C. L. Fryer, W. Even, B. W. Grefenstette, T. -W. Wong

Abstract: Supernovae are Nature's high-energy, high density laboratory experiments, reaching densities in excess of nuclear densities and temperatures above 10MeV. Astronomers have built up a suite of diagnostics to study these supernovae. If we can utilize these diagnostics, and tie them together with a theoretical understanding of supernova physics, we can use these cosmic explosions to study the nature o… ▽ More Supernovae are Nature's high-energy, high density laboratory experiments, reaching densities in excess of nuclear densities and temperatures above 10MeV. Astronomers have built up a suite of diagnostics to study these supernovae. If we can utilize these diagnostics, and tie them together with a theoretical understanding of supernova physics, we can use these cosmic explosions to study the nature of matter at these extreme densities and temperatures. Capitalizing on these diagnostics will require understanding a wide range of additional physics. Here we review the diagnostics and the physics needed to use them to learn about the supernova engine, and ultimate nuclear physics. △ Less

Submitted 14 March, 2014; originally announced March 2014.

Comments: 17 pages, 13 figures, accepted by AIP Advances

Report number: LA-UR-14-20200

Population III Hypernovae

Authors: Joseph Smidt, Daniel J. Whalen, Wesley Even, Brandon K. Wiggins, Jarrett L. Johnson, Chris L. Fryer

Abstract: Population III supernovae have been of growing interest of late for their potential to directly probe the properties of the first stars, particularly the most energetic events that are visible near the edge of the observable universe. But until now, hypernovae, the unusually energetic Type Ib/c supernovae that are sometimes associated with gamma-ray bursts, have been overlooked as cosmic beacons a… ▽ More Population III supernovae have been of growing interest of late for their potential to directly probe the properties of the first stars, particularly the most energetic events that are visible near the edge of the observable universe. But until now, hypernovae, the unusually energetic Type Ib/c supernovae that are sometimes associated with gamma-ray bursts, have been overlooked as cosmic beacons at the highest redshifts. In this, the latest of a series of studies on Population III supernovae, we present numerical simulations of 25 - 50 M$_{\odot}$ hypernovae and their light curves done with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible at z = 10 - 15 to the James Webb Space Telescope (JWST) and z = 4 - 5 to the Wide-Field Infrared Survey Telescope (WFIRST), tracing star formation rates in the first galaxies and at the end of cosmological reionization. If, however, the hypernova crashes into a dense shell ejected by its progenitor, it is expected that a superluminous event will occur that may be seen at z ~ 20, in the first generation of stars. △ Less

Submitted 14 October, 2015; v1 submitted 22 January, 2014; originally announced January 2014.

Comments: 15 pages, 9 figures, accepted by ApJ

Report number: LA-UR-13-23214

The Effects on Supernova Shock Breakout and Swift Light Curves Due to the Mass of the Hydrogen-Rich Envelope

Authors: Amanda J. Bayless, Wesley Even, Lucille H. Frey, Chris L. Fryer, Peter W. A. Roming, Patrick A. Young

Abstract: Mass loss remains one of the primary uncertainties in stellar evolution. In the most massive stars, mass loss dictates the circumstellar medium and can significantly alter the fate of the star. Mass loss is caused by a variety of wind mechanisms and also through binary interactions. Supernovae are excellent probes of this mass loss, both the circumstellar material and the reduced mass of the hydro… ▽ More Mass loss remains one of the primary uncertainties in stellar evolution. In the most massive stars, mass loss dictates the circumstellar medium and can significantly alter the fate of the star. Mass loss is caused by a variety of wind mechanisms and also through binary interactions. Supernovae are excellent probes of this mass loss, both the circumstellar material and the reduced mass of the hydrogen-rich envelope. In this paper, we focus on the effects of reducing the hydrogen-envelope mass on the supernova light curve, studying both the shock breakout and peak light curve emission for a wide variety of mass loss scenarios. Even though the trends of this mass loss will be masked somewhat by variations caused by different progenitors, explosion energies, and circumstellar media, these trends have significant effects on the supernova light-curves that should be seen in supernova surveys. We conclude with a comparison of our results to a few key observations. △ Less

Submitted 24 March, 2015; v1 submitted 17 January, 2014; originally announced January 2014.

Comments: Accepted by ApJ

Pair-Instability Supernovae in the Local Universe

Authors: Daniel J. Whalen, Wesley Even, Joseph Smidt, Alexander Heger, Raphael Hirschi, Norhasliza Yusof, Massimo Stiavelli, Chris L. Fryer, Ke-Jung Chen, Candace C. Joggerst

Abstract: The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500… ▽ More The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500 M$_{\odot}$ Z $\sim$ 0.1 - 0.4 Z$_{\odot}$ stars. These stars lose up to 80% of their mass to strong line-driven winds and explode as bare He cores. We find that their light curves and spectra are quite different from those of Population III pair-instability explosions, which therefore cannot be used as templates for low-redshift events. Although non-zero metallicity pair-instability supernovae are generally dimmer than their Population III counterparts, in some cases they will be bright enough to be detected at the earliest epochs at which they can occur, the formation of the first galaxies at $z \sim$ 10 - 15. Others can masquerade as dim, short duration supernovae that are only visible in the local universe and that under the right conditions could be hidden in a wide variety of supernova classes. We also report for the first time that some pair-instability explosions can create black holes with masses of $\sim$ 100 M$_{\odot}$. △ Less

Submitted 30 November, 2014; v1 submitted 18 December, 2013; originally announced December 2013.

Comments: 16 pages, 15 figures, submitted to ApJ, revised in response to comments by the referee

Report number: LA-UR-13-29505

Finding the First Cosmic Explosions. III. Pulsational Pair-Instability Supernovae

Authors: Daniel J. Whalen, Joseph Smidt, Wesley Even, S. E. Woosley, Alexander Heger, Massimo Stiavelli, Chris L. Fryer

Abstract: Population III supernovae have been the focus of growing attention because of their potential to directly probe the properties of the first stars, particularly the most energetic events that can be seen at the edge of the observable universe. But until now pair-pulsation supernovae, in which explosive thermonuclear burning in massive stars fails to unbind them but can eject their outer layers into… ▽ More Population III supernovae have been the focus of growing attention because of their potential to directly probe the properties of the first stars, particularly the most energetic events that can be seen at the edge of the observable universe. But until now pair-pulsation supernovae, in which explosive thermonuclear burning in massive stars fails to unbind them but can eject their outer layers into space, have been overlooked as cosmic beacons at the earliest redshifts. These shells can later collide and, like Type IIn supernovae, produce superluminous events in the UV at high redshifts that could be detected in the near infrared today. We present numerical simulations of a 110 M$_{\odot}$ pair-pulsation explosion done with the Los Alamos radiation hydrodynamics code RAGE. We find that collisions between consecutive pair pulsations are visible in the near infrared out to z $\sim$ 15 - 20 and can probe the earliest stellar populations at cosmic dawn. △ Less

Submitted 24 December, 2013; v1 submitted 5 November, 2013; originally announced November 2013.

Comments: 9 pages, 4 figures, accepted by ApJ

Report number: LA-UR-13-28360

The Biggest Explosions in the Universe. II

Authors: Daniel J. Whalen, Jarrett L. Johnson, Joseph Smidt, Alexander Heger, Wesley Even, Chris L. Fryer

Abstract: One of the leading contenders for the origin of supermassive black holes at $z \gtrsim$ 7 is catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 15. In this scenario, a few protogalaxies form in the presence of strong Lyman-Werner UV backgrounds that quench H$_2$ formation in their constituent halos, preventing them from forming stars or blowing heavy elements into the intergalacti… ▽ More One of the leading contenders for the origin of supermassive black holes at $z \gtrsim$ 7 is catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 15. In this scenario, a few protogalaxies form in the presence of strong Lyman-Werner UV backgrounds that quench H$_2$ formation in their constituent halos, preventing them from forming stars or blowing heavy elements into the intergalactic medium prior to formation. At masses of 10$^ 8$ \Ms\ and virial temperatures of 10$^4$ K, gas in these halos rapidly cools by H lines, in some cases forming 10$^4$ - 10$^6$ \Ms\ Pop III stars and, a short time later, the seeds of supermassive black holes. Instead of collapsing directly to black holes some of these stars died in the most energetic thermonuclear explosions in the universe. We have modeled the explosions of such stars in the dense cores of line-cooled protogalaxies in the presence of cosmological flows. In stark contrast to the explosions in diffuse regions in previous simulations, these SNe briefly engulf the protogalaxy but then collapse back into its dark matter potential. Fallback drives turbulence that efficiently distributes metals throughout the interior of the halo and fuels the rapid growth of nascent black holes at its center. The accompanying starburst and x-ray emission from these line-cooled galaxies easily distinguish them from more slowly evolving neighbors and might reveal the birthplaces of supermassive black holes on the sky. △ Less

Submitted 6 September, 2013; v1 submitted 14 August, 2013; originally announced August 2013.

Comments: 9 pages, 8 figures, accepted by ApJ

Report number: LA-UR-13-26263

SN 2010jl: Optical to hard X-ray observations reveal an explosion embedded in a ten solar mass cocoon

Authors: E. O. Ofek, A. Zoglauer, S. E. Boggs, N. M. Barriere, S. P. Reynolds, C. L. Fryer, F. A. Harrison, S. B. Cenko, S. R. Kulkarni, A. Gal-Yam, I. Arcavi, E. Bellm, J. S. Bloom, F. Christensen, W. W. Craig, W. Even, A. V. Filippenko, B. Grefenstette, C. J. Hailey, R. Laher, K. Madsen, E. Nakar, P. E. Nugent, D. Stern, M. Sullivan , et al. (2 additional authors not shown)

Abstract: (Abridged) Some supernovae (SNe) may be powered by the interaction of the SN ejecta with a large amount of circumstellar matter (CSM). Here we outline a method to measure the mass of the optically thick CSM around such SNe. We present observations of SN2010jl, including the first detection of a SN using NuSTAR. The total radiated luminosity of SN2010jl is extreme, at least 9e50 erg. By modeling th… ▽ More (Abridged) Some supernovae (SNe) may be powered by the interaction of the SN ejecta with a large amount of circumstellar matter (CSM). Here we outline a method to measure the mass of the optically thick CSM around such SNe. We present observations of SN2010jl, including the first detection of a SN using NuSTAR. The total radiated luminosity of SN2010jl is extreme, at least 9e50 erg. By modeling the visible-light data, we robustly show that the mass of the circumstellar material within ~1e16 cm of the progenitor was in excess of 10 solar masses, likely ejected tens of years prior to the SN explosion. Our modeling suggests that the shock velocity during shock breakout was ~6000 km/s, decelerating to ~2600 km/s about two years after maximum light. Our late-time NuSTAR+XMM spectra of the SN presumably provide the first direct measurement of SN shock velocity two years after the SN maximum light -- measured to be in the range of 2000 to 4500 km/s if the ions and electrons are in equilibrium, and >~2000 km/s if they are not in equilibrium. This measurement is in agreement with the shock velocity predicted by our modeling of the optical data. We also show that the mean radial density distribution of the CSM roughly follows an r^-2 law. A possible explanation for the massive CSM with a wind-like profile is that they are the result of multiple pulsational pair instability events prior to the SN explosion, separated from each other by years. △ Less

Submitted 8 July, 2013; originally announced July 2013.

Comments: Submitted to ApJ, 17 pages, 13 figures

The Supernova that Destroyed a Protogalaxy: Prompt Chemical Enrichment and Supermassive Black Hole Growth

Authors: Daniel J. Whalen, Jarrett J. Johnson, Joseph Smidt, Avery Meiksin, Alexander Heger, Wesley Even, Chris L. Fryer

Abstract: The first primitive galaxies formed from accretion and mergers by z ~ 15, and were primarily responsible for cosmological reionization and the chemical enrichment of the early cosmos. But a few of these galaxies may have formed in the presence of strong Lyman-Werner UV fluxes that sterilized them of H_2, preventing them from forming stars or expelling heavy elements into the IGM prior to assembly.… ▽ More The first primitive galaxies formed from accretion and mergers by z ~ 15, and were primarily responsible for cosmological reionization and the chemical enrichment of the early cosmos. But a few of these galaxies may have formed in the presence of strong Lyman-Werner UV fluxes that sterilized them of H_2, preventing them from forming stars or expelling heavy elements into the IGM prior to assembly. At masses of 10^8 Ms and virial temperatures of 10^4 K, these halos began to rapidly cool by atomic lines, perhaps forming 10^4 - 10^6 Ms Pop III stars and, later, the seeds of supermassive black holes. We have modeled the explosion of a supermassive Pop III star in the dense core of a line-cooled protogalaxy with the ZEUS-MP code. We find that the supernova (SN) expands to a radius of ~ 1 kpc, briefly engulfing the entire galaxy, but then collapses back into the potential well of the dark matter. Fallback fully mixes the interior of the protogalaxy with metals, igniting a violent starburst and fueling the rapid growth of a massive black hole at its center. The starburst would populate the protogalaxy with stars in greater numbers and at higher metallicities than in more slowly-evolving, nearby halos. The SN remnant becomes a strong synchrotron source that can be observed with eVLA and eMERLIN and has a unique signature that easily distinguishes it from less energetic SN remnants. Such explosions, and their attendant starbursts, may well have marked the birthplaces of supermassive black holes on the sky. △ Less

Submitted 9 July, 2013; v1 submitted 29 May, 2013; originally announced May 2013.

Comments: 11 pages, 10 figures, accepted for publication in ApJ

Report number: LA-UR-13-23474

The Biggest Explosions in the Universe

Authors: Jarrett L. Johnson, Daniel J. Whalen, Wesley Even, Chris L. Fryer, Alex Heger, Joseph Smidt, Ke-Jung Chen

Abstract: Supermassive primordial stars are expected to form in a small fraction of massive protogalaxies in the early universe, and are generally conceived of as the progenitors of the seeds of supermassive black holes (BHs). Supermassive stars with masses of ~55,000 M_Sun, however, have been found to explode and completely disrupt in a supernova (SN) with an energy of up to ~10^55 erg instead of collapsin… ▽ More Supermassive primordial stars are expected to form in a small fraction of massive protogalaxies in the early universe, and are generally conceived of as the progenitors of the seeds of supermassive black holes (BHs). Supermassive stars with masses of ~55,000 M_Sun, however, have been found to explode and completely disrupt in a supernova (SN) with an energy of up to ~10^55 erg instead of collapsing to a BH. Such events, ~10,000 times more energetic than typical SNe today, would be among the biggest explosions in the history of the universe. Here we present a simulation of such a SN in two stages. Using the RAGE radiation hydrodynamics code we first evolve the explosion from an early stage through the breakout of the shock from the surface of the star until the blast wave has propagated out to several parsecs from the explosion site, which lies deep within an atomic cooling dark matter (DM) halo at z ~ 15. Then, using the GADGET cosmological hydrodynamics code we evolve the explosion out to several kiloparsecs from the explosion site, far into the low-density intergalactic medium. The host DM halo, with a total mass of 4 x 10^7 M_Sun, much more massive than typical primordial star-forming halos, is completely evacuated of high density gas after < 10 Myr, although dense metal-enriched gas recollapses into the halo, where it will likely form second-generation stars with metallicities of ~ 0.05 Z_Sun after > 70 Myr. The chemical signature of supermassive star explosions may be found in such long-lived second-generation stars today. △ Less

Submitted 1 August, 2013; v1 submitted 16 April, 2013; originally announced April 2013.

Comments: 9 pages, 7 figures; accepted for publication in ApJ; a short movie of the explosion can be found at http://www.as.utexas.edu/~jljohnson/SMS_SN.mp4

Report number: LA-UR-13-22490

Illuminating the Primeval Universe with Type IIn Supernovae

Authors: Daniel J. Whalen, Wesley Even, C. C. Lovekin, Chris L. Fryer, Massimo Stiavelli, P. W. A. Roming, Jeff Cooke, T. A. Pritchard, Daniel E. Holz, Cynthia Knight

Abstract: The detection of Pop III supernovae could directly probe the primordial IMF for the first time, unveiling the properties of the first galaxies, early chemical enrichment and reionization, and the seeds of supermassive black holes. Growing evidence that some Pop III stars were less massive than 100 solar masses may complicate prospects for their detection, because even though they would have been m… ▽ More The detection of Pop III supernovae could directly probe the primordial IMF for the first time, unveiling the properties of the first galaxies, early chemical enrichment and reionization, and the seeds of supermassive black holes. Growing evidence that some Pop III stars were less massive than 100 solar masses may complicate prospects for their detection, because even though they would have been more plentiful they would have died as core-collapse supernovae, with far less luminosity than pair-instability explosions. This picture greatly improves if the SN shock collides with a dense circumstellar shell ejected during a prior violent LBV type eruption. Such collisions can turn even dim SNe into extremely bright ones whose luminosities can rival those of pair-instability SNe. We present simulations of Pop III Type IIn SN light curves and spectra performed with the Los Alamos RAGE and SPECTRUM codes. Taking into account Lyman-alpha absorption in the early universe and cosmological redshifting, we find that 40 solar mass Pop III Type IIn SNe will be visible out to z ~ 20 with JWST and out to z ~ 7 with WFIRST. Thus, even low mass Pop III SNe can be used to probe the primeval universe. △ Less

Submitted 13 March, 2013; v1 submitted 2 February, 2013; originally announced February 2013.

Comments: 17 pages, 11 figures, revised in response to the referee, submitted to ApJ

Report number: LA-UR-12-26001

Finding the First Cosmic Explosions I: Pair-Instability Supernovae

Authors: Daniel J. Whalen, Wesley Even, Lucille H. Frey, Joseph Smidt, Jarrett L. Johnson, C. C. Lovekin, Chris L. Fryer, Massimo Stiavelli, Daniel E. Holz, Alexander Heger, S. E. Woosley, Aimee L. Hungerford

Abstract: The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova ex… ▽ More The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova explosions, which may soon be detected by a new generation of NIR observatories such as JWST and WFIRST. We present light curves and spectra for Pop III pair-instability supernovae calculated with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan-STARRS will find them at z ~ 7 - 8. The discovery of these ancient explosions will probe the first stellar populations and reveal the existence of primitive galaxies that might not otherwise have been detected. △ Less

Submitted 5 September, 2013; v1 submitted 21 November, 2012; originally announced November 2012.

Comments: 17 pages, 12 figures, accepted by ApJ

Report number: LA-UR-12-25732

Supermassive Population III Supernovae and the Birth of the First Quasars

Authors: Daniel J. Whalen, Wesley Even, Joseph Smidt, Alexander Heger, K. -J. Chen, Chris L. Fryer, Massimo Stiavelli, Hao Xu, Candace C. Joggerst

Abstract: The existence of supermassive black holes as early as z ~ 7 is one of the great unsolved problems in cosmological structure formation. One leading theory argues that they are born during catastrophic baryon collapse in z ~ 15 protogalaxies in strong Lyman-Werner UV backgrounds. Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to directly collapse to 10^4… ▽ More The existence of supermassive black holes as early as z ~ 7 is one of the great unsolved problems in cosmological structure formation. One leading theory argues that they are born during catastrophic baryon collapse in z ~ 15 protogalaxies in strong Lyman-Werner UV backgrounds. Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to directly collapse to 10^4 - 10^5 solar-mass black holes. We have now discovered that some of these fragments can instead become supermassive stars that eventually explode as pair-instability supernovae with energies of ~ 10^55 erg, the most energetic explosions in the universe. We have calculated light curves and spectra for supermassive Pop III PI SNe with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible in NIR all-sky surveys by Euclid out to z ~ 10 - 15 and by WFIRST and WISH out to z ~ 15 - 20, perhaps revealing the birthplaces of the first quasars. △ Less

Submitted 12 September, 2013; v1 submitted 8 November, 2012; originally announced November 2012.

Comments: 13 pages, 11 figures, accepted by ApJ

Report number: LA-UR-12-25778

The Long-Lived UV "Plateau" of SN 2012aw

Authors: Amanda J. Bayless, Tyler Pritchard, Peter W. A. Roming, Paul Kuin, Peter J. Brown, Maria Teresa Botticella, Massimo Dall'Ora, Lucille H. Frey, Wesley Even, Chris L. Fryer, Justyn R. Maund, Morgan Fraser

Abstract: Observations with the Swift UVOT have unambiguously uncovered for the first time a long-lived, UV "plateau" in a Type II-P supernova (SN). Although this flattening in slope is hinted at in a few other SNe, due to its proximity and minimal line-of-sight extinction, SN 2012aw has afforded the first opportunity to clearly observe this UV plateau. The observations of SN 2012aw revealed all Swift UV… ▽ More Observations with the Swift UVOT have unambiguously uncovered for the first time a long-lived, UV "plateau" in a Type II-P supernova (SN). Although this flattening in slope is hinted at in a few other SNe, due to its proximity and minimal line-of-sight extinction, SN 2012aw has afforded the first opportunity to clearly observe this UV plateau. The observations of SN 2012aw revealed all Swift UV and u-band lightcurves initially declined rapidly, but 27 days after explosion the light curves flattened. Some possible sources of the UV plateau are: the same thermal process that cause the optical plateau, heating from radioactive decay, or a combination of both processes. △ Less

Submitted 15 January, 2013; v1 submitted 19 October, 2012; originally announced October 2012.

Comments: In press at ApJL. 5 pages. 4 color figures

Seeing the First Supernovae at the Edge of the Universe with JWST

Authors: Daniel J. Whalen, Chris L. Fryer, Daniel E. Holz, Alexander Heger, S. E. Woosley, Massimo Stiavelli, Wesley Even, Lucille L. Frey

Abstract: The first stars ended the cosmic Dark Ages and created the first heavy elements necessary for the formation of planets and life. The properties of these stars remain uncertain, and it may be decades before individual Pop III stars are directly observed. Their masses, however, can be inferred from their supernova explosions, which may soon be found in both deep-field surveys by JWST and in all-sky… ▽ More The first stars ended the cosmic Dark Ages and created the first heavy elements necessary for the formation of planets and life. The properties of these stars remain uncertain, and it may be decades before individual Pop III stars are directly observed. Their masses, however, can be inferred from their supernova explosions, which may soon be found in both deep-field surveys by JWST and in all-sky surveys by WFIRST. We have performed radiation hydrodynamical simulations of the near infrared signals of Pop III pair-instability supernovae in realistic circumstellar environments with Lyman absorption by the neutral intergalactic medium. We find that JWST and WFIRST will detect these explosions out to z ~ 30 and 20, respectively, unveiling the first generation of stars in the universe. △ Less

Submitted 22 November, 2012; v1 submitted 16 September, 2012; originally announced September 2012.

Comments: 5 pages, 4 figures, accepted by ApJL

Report number: LA-UR 10-06844

Do R Coronae Borealis Stars Form from Double White Dwarf Mergers?

Authors: Jan. E. Staff, Athira Menon, Falk Herwig, Wesley Even, Chris L. Fryer, Patrick M. Motl, Tom Geballe, Marco Pignatari, Geoffrey C. Clayton, Joel E. Tohline

Abstract: A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the resu… ▽ More A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars. △ Less

Submitted 2 August, 2012; originally announced August 2012.

Comments: 42 pages, 19 figures. Accepted for publication in the Astrophysical Journal

The Los Alamos Supernova Light Curve Project: Computational Methods

Authors: Lucille H. Frey, Wesley Even, Daniel J. Whalen, Chris L. Fryer, Aimee L. Hungerford, Christopher J. Fontes, James Colgan

Abstract: We have entered the era of explosive transient astronomy, in which upcoming real-time surveys like the Large Synoptic Survey Telescope (LSST), the Palomar Transient Factory (PTF) and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) will detect supernovae in unprecedented numbers. Future telescopes such as the James Webb Space Telescope may discover supernovae from the earliest sta… ▽ More We have entered the era of explosive transient astronomy, in which upcoming real-time surveys like the Large Synoptic Survey Telescope (LSST), the Palomar Transient Factory (PTF) and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) will detect supernovae in unprecedented numbers. Future telescopes such as the James Webb Space Telescope may discover supernovae from the earliest stars in the universe and reveal their masses. The observational signatures of these astrophysical transients are the key to unveiling their central engines, the environments in which they occur, and to what precision they will pinpoint cosmic acceleration and the nature of dark energy. We present a new method for modeling supernova light curves and spectra with the radiation hydrodynamics code RAGE coupled with detailed monochromatic opacities in the SPECTRUM code. We include a suite of tests that demonstrate how the improved physics is indispensable to modeling shock breakout and light curves. △ Less

Submitted 25 January, 2013; v1 submitted 26 March, 2012; originally announced March 2012.

Comments: 18 pages, 19 figures, published in ApJ Supplements

Report number: LA-UR 12-01363

Journal ref: ApJS, 204, 16, 2013

Constructing Synchronously Rotating Double White Dwarf Binaries

Authors: Wesley Even, Joel E. Tohline

Abstract: We have developed a self-consistent-field technique similar to the one described by Hachisu, Eriguchi, & Nomoto (1986b) that can be used to construct detailed force-balanced models of synchronously rotating, double white dwarf (DWD) binaries that have a wide range of total masses, mass ratios, and separations. In addition to providing a computational tool that can be used to provide quiet initia… ▽ More We have developed a self-consistent-field technique similar to the one described by Hachisu, Eriguchi, & Nomoto (1986b) that can be used to construct detailed force-balanced models of synchronously rotating, double white dwarf (DWD) binaries that have a wide range of total masses, mass ratios, and separations. In addition to providing a computational tool that can be used to provide quiet initial starts for dynamical studies of the onset of mass transfer in DWD systems, we show that this SCF technique can be used to construct model sequences that mimic the last portion of the detached inspiral phase of DWD binary evolutions, and semi-detached model sequences that mimic a phase of conservative mass transfer. △ Less

Submitted 14 August, 2009; originally announced August 2009.

Comments: 51 pages, 10 figures, submitted for publication in Astrophysical Journal Supplement Series