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.  more » « less
Award ID(s):
2108339 2020275 2110218
NSF-PAR ID:
10318487
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
659
ISSN:
0004-6361
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Current observations of binary black hole (BBH) merger events show support for a feature in the primary BH-mass distribution at $\sim \, 35 \ \mathrm{M}_{\odot }$, previously interpreted as a signature of pulsational pair-instability supernovae (PPISNe). Such supernovae are expected to map a wide range of pre-supernova carbon–oxygen (CO) core masses to a narrow range of BH masses, producing a peak in the BH mass distribution. However, recent numerical simulations place the mass location of this peak above $50 \ \mathrm{M}_{\odot }$. Motivated by uncertainties in the progenitor’s evolution and explosion mechanism, we explore how modifying the distribution of BH masses resulting from PPISN affects the populations of gravitational-wave (GW) and electromagnetic (EM) transients. To this end, we simulate populations of isolated BBH systems and combine them with cosmic star formation rates. Our results are the first cosmological BBH-merger predictions made using the binary_c rapid population synthesis framework. We find that our fiducial model does not match the observed GW peak. We can only explain the $35 \ \mathrm{M}_{\odot }$ peak with PPISNe by shifting the expected CO core-mass range for PPISN downwards by $\sim {}15 \ \mathrm{M}_{\odot }$. Apart from being in tension with state-of-the art stellar models, we also find that this is likely in tension with the observed rate of hydrogen-less super-luminous supernovae. Conversely, shifting the mass range upward, based on recent stellar models, leads to a predicted third peak in the BH mass function at $\sim {}64 \ \mathrm{M}_{\odot }$. Thus we conclude that the $\sim {}35 \ \mathrm{M}_{\odot }$ feature is unlikely to be related to PPISN.

     
    more » « less
  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. 
    more » « less
  3. 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.

     
    more » « less
  4. 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.

     
    more » « less
  5. Context. The possible existence of warm ( T eff  ∼ 19 000 K) pulsating DA white dwarf (WD) stars, hotter than ZZ Ceti stars, was predicted in theoretical studies more than 30 yr ago. These studies reported the occurrence of g -mode pulsational instabilities due to the κ mechanism acting in the partial ionization zone of He below the H envelope in models of DA WDs with very thin H envelopes ( M H / M ⋆  ≲ 10 −10 ). However, to date, no pulsating warm DA WD has been discovered, despite the varied theoretical and observational evidence suggesting that a fraction of WDs should be formed with a range of very low H content. Aims. We re-examine the pulsational predictions for such WDs on the basis of new full evolutionary sequences. We analyze all the warm DAs observed by the TESS satellite up to Sector 9 in order to search for the possible pulsational signal. Methods. We computed WD evolutionary sequences of masses 0.58 and 0.80 M ⊙ with H content in the range −14.5 ≲ log( M H / M ⋆ )≲ − 10, appropriate for the study of pulsational instability of warm DA WDs. Initial models were extracted from progenitors that were evolved through very late thermal pulses on the early cooling branch. We use LPCODE stellar code into which we have incorporated a new full-implicit treatment of time-dependent element diffusion to precisely model the H–He transition zone in evolving WD models with very low H content. The nonadiabatic pulsations of our warm DA WD models were computed in the effective temperature range of 30 000 − 10 000 K, focusing on ℓ = 1 g modes with periods in the range 50 − 1500 s. Results. We find that traces of H surviving the very late thermal pulse float to the surface, eventually forming thin, growing pure H envelopes and rather extended H–He transition zones. We find that such extended transition zones inhibit the excitation of g modes due to partial ionization of He below the H envelope. Only in the cases where the H–He transition is assumed much more abrupt than predicted by diffusion do models exhibit pulsational instability. In this case, instabilities are found only in WD models with H envelopes in the range of −14.5 ≲ log( M H / M ⋆ )≲ − 10 and at effective temperatures higher than those typical for ZZ Ceti stars, in agreement with previous studies. None of the 36 warm DAs observed so far by TESS satellite are found to pulsate. Conclusions. Our study suggests that the nondetection of pulsating warm DAs, if WDs with very thin H envelopes do exist, could be attributed to the presence of a smooth and extended H–He transition zone. This could be considered as indirect proof that element diffusion indeed operates in the interior of WDs. 
    more » « less