Impact of stellar winds on the pair-instability supernova rate

TL;DR

This study models the impact of stellar winds on the rate of pair-instability supernovae from very massive stars, showing that realistic wind models significantly lower predicted rates and depend strongly on metallicity.

Contribution

It introduces new stellar evolution models with improved wind prescriptions, demonstrating their effect on supernova rates and the importance of metallicity in stellar fate.

Findings

Predicted PISN rate at z~0 is about 0.1 Gpc$^{-3}$ yr$^{-1}$, much lower than previous models.

Stars with metallicity above 0.002 do not undergo PISN, even at high masses.

Mass loss due to winds at low metallicity prevents the development of the He core needed for PISN.

Abstract

Very massive stars (VMSs, $M_{\star}$ $\geq$ 100 M$_{\odot}$) play a crucial role in several astrophysical processes. At low metallicity, they might collapse directly into black holes, or end their lives as pair-instability supernovae. Recent observational results set an upper limit of $0.7\,{}\mathrm{ yr}^{-1} \,{}\mathrm{ Gpc}^{-3}$ on the rate density of pair-instability supernovae in the nearby Universe. However, most theoretical models predict rates exceeding this limit. Here, we compute new VMS tracks with the MESA code, and use them to analyze the evolution of the (pulsational) pair-instability supernova rate density across cosmic time. We show that stellar wind models accounting for the transition between optically thin and thick winds yield a pair-instability supernova rate $\mathcal{R}_{\mathrm{PISN}}\sim{}0.1$ Gpc$^{-3}$ yr$^{-1}$ at redshift $z\sim{}0$, about two orders of magnitude lower than our previous models. We find that the main contribution to the pair-instability supernova rate comes from stars with metallicity $Z\sim{}0.001-0.002$. Stars with higher metallicities cannot enter the pair-instability supernova regime, even if their zero-age main sequence mass is up to 500 M$_\odot$. The main reason is that VMSs enter the regime for optically thick winds during the main sequence at metallicity as low as $Z\sim{4}\times{}10^{-4}$. This enhances the mass loss rate, quenching the growth of the He core and thus preventing the onset of pair-instability in later evolutionary stages. This result highlights the critical role of mass loss in shaping the final fate of very massive stars and the rate of pair-instability supernovae.