-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
$^{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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
2D Implosion Simulations with a Kinetic Particle Code
Authors:
Irina Sagert,
Wesley P. Even,
Terrance T. Strother
Abstract:
We perform two-dimensional (2D) implosion simulations using a Monte Carlo kinetic particle code. The paper is motivated by the importance of non-equilibrium effects in inertial confinement fusion (ICF) capsule implosions. These cannot be fully captured by hydrodynamic simulations while kinetic methods, as the one presented in this study, are able to describe continuum and rarefied regimes within o…
▽ More
We perform two-dimensional (2D) implosion simulations using a Monte Carlo kinetic particle code. The paper is motivated by the importance of non-equilibrium effects in inertial confinement fusion (ICF) capsule implosions. These cannot be fully captured by hydrodynamic simulations while kinetic methods, as the one presented in this study, are able to describe continuum and rarefied regimes within one approach. In the past, our code has been verified via traditional shock wave and fluid instability simulations. In the present work, we focus on setups that are closer to applications in ICF. We perform simple 2D disk implosion simulations using one particle species. The obtained results are compared to simulations using the hydrodynamics code RAGE. In a first study, the implosions are powered by energy deposition in the outer layers of the disk. We test the impact of the particle mean-free-path and find that while the width of the implosion shock broadens, its location as a function of time remains very similar. In a second study, we focus on the formation of fluid instabilities from induced perturbations. We find good agreement with hydrodynamic studies regarding the location of the shock, the implosion dynamics and the formation of imposed fluid instabilities. Differences are due to the higher resolution of RAGE and statistical noise in the kinetic simulations. The current studies represent proofs-of-principle and serve as a stepping stone for more sophisticated ICF simulations in the future.
△ Less
Submitted 30 May, 2017; v1 submitted 14 November, 2016;
originally announced November 2016.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.