Abstract
Our understanding of the convective-engine paradigm driving core-collapse supernovae has been used for two decades to predict the remnant mass distribution from stellar collapse. These ...predictions improve as our understanding of this engine increases. In this paper, we review our current understanding of convection (in particular, the growth rate of convection) in stellar collapse and study its effect on the remnant mass distribution. We show how the depth of the mass gap between neutron stars and black holes can help probe this convective growth. We include a study of the effects of stochasticity in both the stellar structure and the convective seeds caused by stellar burning. We study the role of rotation and its effect on the pair-instability mass gap. Under the paradigm limiting stellar rotation to those stars in tight binaries, we determine the effect of rotation on the remnant mass distribution.
ABSTRACT
Understanding the astrophysical phenomena involving compact objects requires an insight about the engine behind the core-collapse supernovae (SNe) and the fate of the stellar collapse of ...massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of the growing populations of the detected black hole–black hole, black hole–neutron star, and neutron star–neutron star mergers. The time-scale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. We adopt new formulas for the relation between the pre-SN star properties and their remnants and check how they impact the population of double compact object (DCO) mergers formed via the isolated binary evolution. The new formulas give one the ability to test a wide spectrum of assumptions on the convection growth time. In particular, the different variants allow for a smooth transition between having a deep mass gap and a remnant mass distribution filled by massive neutron stars and low-mass black holes. We present the distribution of masses, mass ratios, and the local merger rate densities of DCO for the different variants of new formulas and test them together with different approaches to other highly uncertain processes. We find that the mass distribution of DCO mergers is sensitive to the adopted assumption on the SN convection growth time-scale up to $m_1+m_2 \lesssim 35 \, \mathrm{M}_{\odot }$. Between the two extreme tested variants the probability of compact object formation within the mass gap may differ by up to approximately two orders of magnitude.
The mass distribution of neutron stars and stellar-mass black holes provides vital clues into the nature of stellar core collapse and the physical engine responsible for supernova explosions. A ...number of supernova engines have been proposed: neutrino- or oscillation-driven explosions enhanced by early (developing in 10-50 ms) and late-time (developing in 200 ms) convection as well as magnetic field engines (in black hole accretion disks or neutron stars). Using our current understanding of supernova engines, we derive mass distributions of stellar compact remnants. We provide analytic prescriptions for both single-star models (as a function of initial star mass) and for binary-star models-prescriptions for compact object masses for major population synthesis codes. These prescriptions have implications for a range of observations: X-ray binary populations, supernova explosion energies, and gravitational wave sources. We show that advanced gravitational radiation detectors (like LIGO/VIRGO or the Einstein Telescope) will be able to further test the supernova explosion engine models once double black hole inspirals are detected.
The Origin of r-process Elements in the Milky Way Côté, Benoit; Fryer, Chris L.; Belczynski, Krzysztof ...
Astrophysical journal/The Astrophysical journal,
03/2018, Letnik:
855, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site ...between core-collapse supernovae and neutron star-neutron star (NS-NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by quantifying the frequency and yields required to reproduce the amount of europium (Eu) observed in galaxies. Although NS-NS mergers have become popular candidates, their required frequency (or rate) needs to be consistent with that obtained from gravitational wave measurements. Here, we address the first NS-NS merger detected by LIGO/Virgo (GW170817) and its associated gamma-ray burst and analyze their implication for the origin of r-process elements. The range of NS-NS merger rate densities of 320-4740 Gpc−3 yr−1 provided by LIGO/Virgo is remarkably consistent with the range required by GCE to explain the Eu abundances in the Milky Way with NS-NS mergers, assuming the solar r-process abundance pattern for the ejecta. Under the same assumption, this event has produced about 1-5 Earth masses of Eu, and 3-13 Earth masses of gold. When using theoretical calculations to derive Eu yields, constraining the role of NS-NS mergers becomes more challenging because of nuclear astrophysics uncertainties. This is the first study that directly combines nuclear physics uncertainties with GCE calculations. If GW170817 is a representative event, NS-NS mergers can produce Eu in sufficient amounts and are likely to be the main r-process site.
Core-collapse supernova science is now entering an era in which engine models are beginning to make both qualitative and, in some cases, 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 three progenitors and a narrow range of explosion energies, we can produce a wide range of remnant masses and nucleosynthetic yields.
Abstract
The combined detection of a gravitational-wave signal, kilonova, and short gamma-ray burst (sGRB) from GW170817 marked a scientific breakthrough in the field of multimessenger astronomy. But ...even before GW170817, there have been a number of sGRBs with possible associated kilonova detections. In this work, we re-examine these ‘historical’ sGRB afterglows with a combination of state-of-the-art afterglow and kilonova models. This allows us to include optical/near-infrared synchrotron emission produced by the sGRB as well as ultraviolet/optical/near-infrared emission powered by the radioactive decay of r-process elements (i.e. the kilonova). Fitting the light curves, we derive the velocity and the mass distribution as well as the composition of the ejected material. The posteriors on kilonova parameters obtained from the fit were turned into distributions for the peak magnitude of the kilonova emission in different bands and the time at which this peak occurs. From the sGRB with an associated kilonova, we found that the peak magnitude in H bands falls in the range −16.2, −13.1 ($95{{\ \rm per\ cent}}$ of confidence) and occurs within $0.8\!-\!3.6\, \rm d$ after the sGRB prompt emission. In g band instead we obtain a peak magnitude in range −16.8, −12.3 occurring within the first 18 h after the sGRB prompt. From the luminosity distributions of GW170817/AT2017gfo, kilonova candidates GRB130603B, GRB050709, and GRB060614 (with the possible inclusion of GRB150101B, GRB050724A, GRB061201, GRB080905A, GRB150424A, and GRB160821B) and the upper limits from all the other sGRBs not associated with any kilonova detection we obtain for the first time a kilonova luminosity distribution in different bands.
We constrain the properties of massive binaries by comparing radial velocity data on early-type stars in Cygnus OB2 with the expectations of Monte Carlo models. Our comparisons test several popular ...prescriptions for massive binary parameters. We explore a range of true binary fraction, F, a range of power-law slopes, a, describing the distribution of companion masses between the limits q sub(low) and 1, and a range of power-law slopes, beta , describing the distribution of orbital separations between the limits r sub(in) and r sub(out). We also consider distributions of secondary masses described by a Miller-Scalo type IMF and by a two-component IMF that Includes a substantial "twin" population with M sub(2) unk M sub(1). Several seemingly disparate prescriptions for massive binary characteristics can be reconciled by adopting carefully chosen values for F, r sub(in), and r sub(out). We show that binary fractions F < 0.7 are less probable than F greater than or equal to 0.8 for reasonable choices of r sub(in) and r sub(out). Thus, the true binary fraction is high. For F = 1.0 and a distribution of orbital separations near the canonical "Opik's law distribution (i.e., flat; beta = 0), the power-law slope of the mass ratio distribution is alpha = -0.6 to 0.0. For F unk 0.8, alpha is somewhat larger, in the range -0.4 to 1.0. In any case, the secondary star mass function is inconsistent with a Miller-Scalo-like IMF unless the lower end is truncated below similar to 2-4 M unk. In other words, massive stars preferentially have massive companions. The best-fitting models are described by a Salpeter or Miller-Scalo IMF for 60% of secondary star masses with the other similar to 40% of secondaries having M sub(2) unk M sub(1), i.e., "twins." These model parameters simultaneously predict the fraction of Type Ib/c supernovae to be 30%-40% of all core-collapse supernovae, in agreement with recent observational estimates.
The induced gravitational collapse (IGC) paradigm has been successfully applied to the explanation of the concomitance of gamma-ray bursts (GRBs) with supernovae (SNe) Ic. The progenitor is a tight ...binary system composed of a carbon-oxygen (CO) core and a neutron star (NS) companion. The explosion of the SN leads to hypercritical accretion onto the NS companion, which reaches the critical mass, hence inducing its gravitational collapse to a black hole (BH) with consequent emission of the GRB. The first estimates of this process were based on a simplified model of the binary parameters and the Bondi-Hoyle-Lyttleton accretion rate. We present here the first full numerical simulations of the IGC phenomenon. We simulate the core-collapse and SN explosion of CO stars to obtain the density and ejection velocity of the SN ejecta. We follow the hydrodynamic evolution of the accreting material falling into the Bondi-Hoyle surface of the NS all the way up to its incorporation in the NS surface. The simulations go up to BH formation when the NS reaches the critical mass. For appropriate binary parameters, the IGC occurs in short timescales ~10 super(2)-10 super(3) s owing to the combined effective action of the photon trapping and the neutrino cooling near the NS surface. We also show that the IGC scenario leads to a natural explanation for why GRBs are associated only with SNe Ic with totally absent or very little helium.
Probing the origin of r-process elements in the universe represents a multidisciplinary challenge. We review the observational evidence that probes the properties of r-process sites, and address them ...using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contributed to the total mass of r-process elements present in our Galaxy. We found discrepancies with the neutron star (NS-NS) merger scenario. When we assume that they are the only site, the decreasing trend of Eu/Fe at Fe/H > −1 in the disk of the Milky Way cannot be reproduced while accounting for the delay-time distribution (DTD) of coalescence times (∝t−1) derived from short gamma-ray bursts (GRBs) and population synthesis models. Steeper DTD functions (∝t−1.5) or power laws combined with a strong burst of mergers before the onset of supernovae (SNe) Ia can reproduce the Eu/Fe trend, but this scenario is inconsistent with the similar fraction of short GRBs and SNe Ia occurring in early-type galaxies, and it reduces the probability of detecting GW170817 in an early-type galaxy. One solution is to assume an additional production site of Eu that would be active in the early universe, but would fade away with increasing metallicity. If this is correct, this additional site could be responsible for roughly 50% of the Eu production in the early universe before the onset of SNe Ia. Rare classes of supernovae could be this additional r-process source, but hydrodynamic simulations still need to ensure the conditions for a robust r-process pattern.
It is firmly established that the stellar mass distribution is smooth, covering the range 0.1-100 M sub(middot in circle). It is to be expected that the masses of the ensuing compact remnants ...correlate with the masses of their progenitor stars, and thus it is generally thought that the remnant masses should be smoothly distributed from the lightest white dwarfs to the heaviest black holes (BHs). However, this intuitive prediction is not borne out by observed data. In the rapidly growing population of remnants with observationally determined masses, a striking mass gap has emerged at the boundary between neutron stars (NSs) and BHs. The heaviest NSs reach a maximum of two solar masses, while the lightest BHs are at least five solar masses. Over a decade after the discovery, the gap has become a significant challenge to our understanding of compact object formation. We offer new insights into the physical processes that bifurcate the formation of remnants into lower-mass NSs and heavier BHs. Combining the results of stellar modeling with hydrodynamic simulations of supernovae, we both explain the existence of the gap and also put stringent constraints on the inner workings of the supernova explosion mechanism. In particular, we show that core-collapse supernovae are launched within 100-200 ms of the initial stellar collapse, implying that the explosions are driven by instabilities with a rapid (10-20 ms) growth time. Alternatively, if future observations fill in the gap, this will be an indication that these instabilities develop over a longer (>200 ms) timescale.