ABSTRACT
The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These ...observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy r-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound–bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion–opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and thus eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion–opacity formalism. The agreement between the line-binned and expansion–opacity results is explained as arising from the similarity in their opacities in the limit of low optical depth, where radiation transport is important in the ejecta. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semirelativistic approaches.
We provide here a significant extension of the NuGrid Set 1 models in mass coverage and towards lower metallicity, adopting the same physics assumptions. The combined data set now includes the ...initial masses MZAMS/M⊙ = 1, 1.65, 2, 3, 4, 5, 6, 7, 12, 15, 20, 25 for Z = 0.02, 0.01, 0.006, 0.001, 0.0001 with α-enhanced composition for the lowest three metallicities. These models are computed with the MESA stellar evolution code and are evolved up to the AGB, the white dwarf stage, or until core collapse. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid MPPNP code. Explosive nucleosynthesis is based on semi-analytic 1D shock models. Metallicity-dependent mass-loss, convective boundary mixing in low- and intermediate-mass models and H and He core burning massive star models are included. Convective O-C shell mergers in some stellar models lead to the strong production of odd-Z elements P, Cl, K, and Sc. In AGB models with hot dredge-up, the convective boundary mixing efficiency is reduced to accommodate for its energetic feedback. In both low-mass and massive star models at the lowest metallicity, H-ingestion events are observed and lead to i-process nucleosynthesis and substantial 15N production. Finally, complete yield data tables, derived data products and online analytic data access are provided.
Gamma-ray bursts (GRBs) are flashes of high-energy radiation arising from energetic cosmic explosions. Bursts of long (greater than two seconds) duration are produced by the core-collapse of massive ...stars
, and those of short (less than two seconds) duration by the merger of compact objects, such as two neutron stars
. A third class of events with hybrid high-energy properties was identified
, but never conclusively linked to a stellar progenitor. The lack of bright supernovae rules out typical core-collapse explosions
, but their distance scales prevent sensitive searches for direct signatures of a progenitor system. Only tentative evidence for a kilonova has been presented
. Here we report observations of the exceptionally bright GRB 211211A, which classify it as a hybrid event and constrain its distance scale to only 346 megaparsecs. Our measurements indicate that its lower-energy (from ultraviolet to near-infrared) counterpart is powered by a luminous (approximately 10
erg per second) kilonova possibly formed in the ejecta of a compact object merger.
A Broad Grid of 2D Kilonova Emission Models Wollaeger, R. T.; Fryer, C. L.; Chase, E. A. ...
Astrophysical journal/The Astrophysical journal,
09/2021, Letnik:
918, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Abstract
Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The ...characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of two-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1
M
⊙
and 0.05 to 0.3
c
, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (
Y
e
) ejecta with a robust
r
-process composition and either a spherical or lobed high-
Y
e
ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code
SuperNu
using a full suite of lanthanide and fourth-row element opacities. We examine the trends of these models with parameter variation, show how they can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.
Context.
As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of ...gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events.
Aims.
In this work, we present new developments in the stellar-remnant formation and related schemes of the current
N
-body evolution program
NBODY7
. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or
BSE
of the
NBODY7
code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised
StarTrack
population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick.
Methods.
All the above newly implemented schemes are also shared with the standalone versions of
SSE
and
BSE
. All these demonstrations are performed with both the updated standalone
BSE
and the updated
NBODY7
/
BSE
.
Results.
When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing.
Conclusions.
Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.
ABSTRACT
Core-collapse explosions of massive stars leave behind neutron stars, with a known diversity that includes the ‘Central Compact Objects’ (CCOs). Typified by the neutron star discovered near ...the centre of the Cas A supernova remnant (SNR), CCOs have been observed to shine only in X-rays. To address their supernova progenitors, we perform a systematic study of SNRs that contain a CCO and display X-ray emission from their shock-heated ejecta. We make use of X-ray data primarily using the Chandra X-ray observatory, complemented with XMM–Newton. This study uses a systematic approach to the analysis of each SNR aimed at addressing the supernova progenitor as well as the explosion properties (energy and ambient density). After fitting for the ejecta abundances estimated from a spatially resolved spectroscopic study, we compare the data to six nucleosynthesis models making predictions on supernova ejecta yields in core-collapse explosions. We find that the explosion models commonly used by the astrophysics community do not match the ejecta yields for any of the SNRs, suggesting additional physics, for example multidimensional explosion models or updated progenitor structures, are required. Overall we find low-mass (≤25 solar masses) progenitors among the massive stars population and low-energy explosions (<1051 ergs). We discuss degeneracies in our model fitting, particularly how altering the explosion energy affects the estimate of the progenitor mass. Our systematic study highlights the need for improving on the theoretical models for nucleosynthesis predictions as well as for sensitive, high-resolution spectroscopy observations to be acquired with next generation X-ray missions.
We present the first three-dimensional smoothed particle hydrodynamics simulations of the induced gravitational collapse scenario of long-duration gamma-ray bursts (GRBs) associated with supernovae ...(SNe). We simulate the SN explosion of a carbon-oxygen core (COcore) forming a binary system with a neutron star (NS) companion. We follow the evolution of the SN ejecta, including their morphological structure, subject to the gravitational field of both the new NS ( NS) formed at the center of the SN and the one of the NS companion. We compute the accretion rate of the SN ejecta onto the NS companion, as well as onto the NS from SN matter fallback. We determine the fate of the binary system for a wide parameter space including different COcore and NS companion masses, orbital periods, and SN explosion geometry and energies. We identify, for selected NS nuclear equations of state, the binary parameters leading the NS companion, by hypercritical accretion, either to the mass-shedding limit or to the secular axisymmetric instability for gravitational collapse to a black hole (BH), or to a more massive, fast-rotating, stable NS. We also assess whether the binary remains gravitationally bound after the SN explosion, hence exploring the space of binary and SN explosion parameters leading to NS-NS and NS-BH binaries. The consequences of our results for the modeling of long GRBs, i.e., X-ray flashes and binary-driven hypernovae, are discussed.
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. The post-bounce temperature and density evolution of zero-age main-sequence 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid 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, and 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 the prospects of observing these radionuclides in supernova remnants.
ABSTRACT We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our ...goal is to provide an internally consistent and comprehensive nuclear production and yield database for applications in areas such as presolar grain studies. Our non-rotating models assume convective boundary mixing (CBM) where it has been adopted before. We include 8 (12) initial masses for Z = 0.01 (0.02). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by simple analytic core-collapse supernova (SN) models with two options for fallback and shock velocities. The explosions show which pre-SN yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impact the light elements and the s and p process. For low- and intermediate-mass models, our stellar yields from H to Bi include the effect of CBM at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. The entirety of our stellar nucleosynthesis profile and time evolution output are available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.
Californium-254 and Kilonova Light Curves Zhu, Y.; Wollaeger, R. T.; Vassh, N. ...
Astrophysical journal. Letters,
08/2018, Letnik:
863, Številka:
2
Journal Article
Recenzirano
Odprti dostop
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, 254Cf, 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 254Cf has the greatest impact to the light curve. Future observations in the mid-infrared, which are bright during this regime, could indicate the production of actinide nucleosynthesis.