skip to main content

Title: Relativistic Coulomb screening in pulsational pair instability supernovae
Context. Pulsational pair-instability supernovae (PPISNe) and pair instability supernovae (PISNe) are the result of a thermonuclear runaway in the presence of a background electron-positron pair plasma. As such, their evolution and resultant black hole masses could possibly be affected by screening corrections due to the electron pair plasma. Aims. The sensitivity of PISNe and PPISNe to relativistic weak screening has been explored. Methods. In this paper a weak screening model that includes effects from relativistic pair production has been developed and applied at temperatures approaching and exceeding the threshold for pair production. This screening model replaces “classical” screening commonly used in astrophysics. Modifications to the weak screening electron Debye length were incorporated in a computationally tractable analytic form. Results. In PPISNe the BH masses were found to increase somewhat at high temperatures, though this increase is small. The BH collapse is also found to occur at earlier times, and the pulsational morphology also changes. In addition to the resultant BH mass, the sensitivity to the screening model of the pulsational period, the pulse structure, the PPISN-to-PISN transition, and the shift in the BH mass gap has been analyzed. The dependence of the composition of the ejected mass was also examined.
; ; ; ; ;
Award ID(s):
2108339 2020275 2110218
Publication Date:
Journal Name:
Astronomy & Astrophysics
Sponsoring Org:
National Science Foundation
More Like this

    We infer the expected detection number of pair instability supernovae (PISNe) during the operation of the Euclid space telescope based on binary population models. Our models reproduce the global maximum at the primary BH mass of ∼9–10 M⊙ and the overall gradient of the primary BH mass distribution in the binary BH merger rate consistent with recent observations. We consider different PISN conditions depending on the 12C(α, γ)16O reaction rate. The fiducial and 3σ models adopt the standard and 3σ smaller reaction rates, respectively. Our fiducial model predicts that Euclid detects several hydrogen-poor PISNe. For the 3σ model, detection of ∼1 hydrogen-poor PISN by Euclid is expected if the stellar mass distribution extends to Mmax  = 600 M⊙, but the expected number becomes significantly smaller if Mmax  = 300 M⊙. We may be able to distinguish the fiducial and 3σ models by the observed PISN rate. This will help us to constrain the origin of binary BHs and the reaction rate, although there remains a degeneracy between Mmax  and the reaction rate. PISN ejecta mass estimates from light curves and spectra obtained by follow-up observations would be important to disentangle the degeneracy.

  2. GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair-instability in particular. We discuss the possibility that stars at low or zero metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime and produce a black hole with a mass in the mass gap by fallback. We present a series of new stellar evolution models at zero and low metallicity computed with the Geneva and MESA stellar evolution codes and compare to existing grids of models. Models with a metallicity in the range 0-0.0004 have three properties which favour higher BH masses as compared to higher metallicity models. These are (i) lower mass-loss rates during the post-MS phase, (ii) a more compact star disfavouring binary interaction and (iii) possible H-He shell interactions which lower the CO core mass. We conclude that it is possible that GW190521 may be the merger of black holes produced directly by massive stars from the first stellar generations. Our models indicate BH masses up to 70-75 Msun. Uncertainties related to convective mixing, mass loss, H-He shell interactions and pair-instability pulsations may increase this limit to ~85 Msun.
  3. Abstract

    Gravitational-wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernovae (PISNe). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections of merging BHs. We useMESAto evolve single, nonrotating, massive helium cores with a metallicity ofZ= 10−5, until they either collapse to form a BH or explode as a PISN, without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the ±3σuncertainty in our high-resolution tabulated12C(α,γ)16O reaction rate probability distribution function. We extensively test temporal and spatial resolutions for resolving the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower-resolution models reach convergence. We establish a new lower edge of the upper mass gap asMlower6014+32Mfrom the ±3σuncertainty inmore »the12C(α,γ)16O rate. We explore the effect of a larger 3αrate on the lower edge of the upper mass gap, findingMlower6918+34M. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog.

    « less
  4. Abstract

    It was recently proposed that exotic particles can trigger a new stellar instability that is analogous to theee+pair instability if they are produced and reach equilibrium in the stellar plasma. In this study, we construct axion instability supernova (AISN) models caused by the new instability to predict their observational signatures. We focus on heavy axion-like particles (ALPs) with masses of ∼400 keV–2 MeV and coupling with photons ofgaγ∼ 10−5GeV−1. It is found that the56Ni mass and the explosion energy are significantly increased by ALPs for a fixed stellar mass. As a result, the peak times of the light curves of AISNe occur earlier than those of standard pair-instability supernovae by 10–20 days when the ALP mass is equal to the electron mass. Also, the event rate of AISNe is 1.7–2.6 times higher than that of pair-instability supernovae, depending on the high mass cutoff of the initial mass function.

  5. Abstract

    We introduce a semiparametric model for the primary mass distribution of binary black holes (BBHs) observed with gravitational waves (GWs) that applies a cubic-spline perturbation to a power law. We apply this model to the 46 BBHs included in the second gravitational-wave transient catalog (GWTC-2). The spline perturbation model recovers a consistent primary mass distribution with previous results, corroborating the existence of a peak at 35M(>97% credibility) found with the Powerlaw+Peakmodel. The peak could be the result of pulsational pair-instability supernovae. The spline perturbation model finds potential signs of additional features in the primary mass distribution at lower masses similar to those previously reported by Tiwari and Fairhurst. However, with fluctuations due to small-number statistics, the simpler Powerlaw+Peakand BrokenPowerlawmodels are both still perfectly consistent with observations. Our semiparametric approach serves as a way to bridge the gap between parametric and nonparametric models to more accurately measure the BBH mass distribution. With larger catalogs we will be able to use this model to resolve possible additional features that could be used to perform cosmological measurements and will build on our understanding of BBH formation, stellar evolution, and nuclear astrophysics.