Pair-instability Supernovae Research Articles

SDSS J0018-0939: A Clear Signature of Sub-Chandrasekhar Mass Type 1a Supernova

Very metal-poor (VMP) stars ([Fe/H] ≤ −2) that have subsolar values of [X/Fe] for α elements such as Mg, Si, and Ca, are referred to as α-poor VMP stars. They are quite rare among VMP stars and are thought to have formed from gas enriched predominantly by a single Type Ia supernova (SN 1a) in contrast to most VMP stars, which are α-enhanced and usually associated with core-collapse supernovae. The observed abundance pattern in such stars can provide a direct way to probe the nucleosynthesis in individual SN 1a. Although the abundance patterns in some α-poor VMP stars have been shown to be consistent with SN 1a ejecta, a clear nucleosynthetic signature for SN 1a resulting from the explosion of a near Chandrasekhar mass (near-M Ch) or a sub-Chandrasekhar mass (sub-M Ch) white dwarf, has not been unambiguously detected. We perform a detailed analysis of various formation channels of VMP stars and find that the α-poor VMP star SDSS J0018-0939, which was earlier reported as a star with potential pair-instability supernova origin, provides a near-smoking-gun signature of a sub-M Ch SN 1a resulting from He detonation. We find that compared to other α-poor VMP stars that were previously identified with SN 1a, SDSS J0018-0939 is the only star that has a clear and unambiguous signature of SN 1a. Interestingly, our results are consistent with constraints on SN1 a from recent galactic chemical evolution studies that indicate that sub-M Ch SN 1a accounts for ∼50%–75% of all SN 1a and is possibly the dominant channel in the early Galaxy.

Origin of LAMOST J1010+2358 Revisited

Signature from Pop III massive stars of 140– that end their lives as pair-instability supernovae (PISNe) are expected to be seen in very metal-poor (VMP) stars of . Although thousands of VMP stars have been discovered, the identification of a VMP star with a PISN signature has been elusive. Recently, the VMP star LAMOST J1010+2358 was claimed to be the first star with a clear PISN signature. A subsequent study showed that ejecta from low-mass core-collapse supernovae (CCSNe) can also fit the abundance pattern equally well and additional elements such as C and Al are required to differentiate the two sources. Follow-up observations of LAMOST J1010+2358 by two independent groups were able to detect both C and Al. Additionally, key odd elements such as Na and Sc were also detected whose abundances were found to be higher than the upper limits found in the original detection. We perform a detailed analysis of the newly observed abundance patterns by exploring various possible formation channels for VMP stars. We find that purely low-mass CCSN ejecta as well as the combination of CCSN and Type 1a SN ejecta can provide an excellent fit to the newly observed abundance pattern. Our results confirm earlier analysis that the newly observed abundance pattern is peculiar but has no signatures of PISN.

The effect of stellar rotation on black hole mass and spin

Abstract The gravitational wave signature of a binary black hole (BBH) merger is dependent on its component mass and spin. If such black holes originate from rapidly rotating progenitors, the large angular momentum reserve in the star could drive a collapsar-like supernova explosion, hence substantially impacting these characteristics of the black holes in the binary. To examine the effect of stellar rotation on the resulting black hole mass and spin, we conduct a 1D general relativistic study of the end phase of the stellar collapse. We find that the resulting black hole mass at times differs significantly from the previously assumed values. We quantify the dependence of the black hole spin magnitude on the hydrodynamics of the accretion flow, providing analytical relations for calculating the mass and spin based on the progenitor’s pre-collapse properties. Depending on the nature of the accretion flow, our findings have implications for the black hole upper mass gap resulting from pair-instability supernovae, the maximum mass of a maximally rotating stellar black hole, and the maximum effective spin of a BBH formed in a tidally locked helium star - black hole binary.

The rates and host galaxies of pair-instability supernovae through cosmic time: predictions from BPASS and IllustrisTNG

ABSTRACT Pair-instability supernovae (PISNe) have long been predicted to be the final fates of near-zero-metallicity very massive stars ($Z \lt Z_\odot /3$, M$_\mathrm{ZAMS} \gtrsim 140\, \text{M}_\odot$). However, no definite PISN has been observed to date, leaving theoretical modelling validation open. To investigate the observability of these explosive transients, we combine detailed stellar evolution models for PISNe formation, computed from the binary population and spectral synthesis code suite, bpass, with the star formation history of all individual computational elements in the Illustris-TNG simulation. This allows us to compute comic PISN rates and predict their host galaxy properties. Of particular importance is that IllustrisTNG galaxies do not have uniform metallicities throughout, with metal-enriched galaxies often harbouring metal-poor pockets of gas where PISN progenitors may form. Accounting for the chemical inhomogeneities within these galaxies, we find that the peak redshift of PISNe formation is $z=3.5$ instead of the value of $z=6$ when ignoring chemical inhomogeneities within galaxies. Furthermore, the rate increases by an order of magnitude from 1.9 to 29 PISN Gpc$^{-3}$ yr$^{-1}$ at $z=0$, if the chemical inhomogeneities are considered. Using state-of-the-art theoretical PISN light curves, we find an observed rate of 13.8 (1.2) visible PISNe per year for the Euclid-Deep survey, or 83 (7.3) over the 6-yr lifetime of the mission when considering chemically inhomogeneous (homogenous) systems. Interestingly, only 12 per cent of helium PISN progenitors are sufficiently massive to power a superluminous supernova event, which can potentially explain why PISN identification in time-domain surveys remains elusive and progress requires dedicated strategies.

The evolution of accreting population III stars at 10−6–103 M⊙ yr−1

Context. The first stars formed over five orders of magnitude in mass by accretion in primordial dark matter halos. Aims. We study the evolution of massive, very massive and supermassive primordial (Pop III) stars over nine orders of magnitude in accretion rate. Methods. We use the stellar evolution code GENEC to evolve accreting Pop III stars from 10−6–103 M⊙ yr−1 and study how these rates determine final masses. The stars are evolved until either the end central Si burning or they encounter the general relativistic instability (GRI). We also examine how metallicity affects the evolution of the star at one accretion rate. Results. At rates below ∼2.5 × 10−5 M⊙ yr−1 the final mass of the star falls below that required for pair-instability supernovae. The minimum rate required to produce black holes with masses above 250 M⊙ is ∼5 × 10−5 M⊙ yr−1, well within the range of infall rates found in numerical simulations of halos that cool via H2, ≲10−3 M⊙ yr−1. At rates of 5 × 10−5 M⊙ yr−1 to 4 × 10−2 M⊙ yr−1, like those expected for halos cooling by both H2 and Lyα, the star collapses after Si burning. At higher accretion rates the GRI triggers the collapse of the star during central H burning. Stars that grow at above these rates are cool red hypergiants with effective temperatures log(Teff) = 3.8 and luminosities that can reach 1010.5 L⊙. At accretion rates of 100–1000 M⊙ yr−1 the gas encounters the general relativistic instability prior to the onset of central hydrogen burning and collapses to a black hole with a mass of ∼106 M⊙ without ever having become a star. Conclusions. Our models corroborate previous studies of Pop III stellar evolution with and without hydrodynamics over separate, smaller ranges in accretion rate. They also reveal for the first time the critical transition rate in accretion above which catastrophic baryon collapse, like that which can occur during galaxy collisions in the high-redshift Universe, produces supermassive black holes via dark collapse.

The cosmic rate of pair-instability supernovae

Abstract Pair-instability supernovae (PISNe) have crucial implications for many astrophysical topics, including the search for very massive stars, the black hole mass spectrum, and galaxy chemical enrichment. To this end, we need to understand where PISNe are across cosmic time, and what are their favourable galactic environments. We present a new determination of the PISN rate as a function of redshift, obtained by combining up-to-date stellar evolution tracks from the PARSEC and FRANEC codes, with an up-to-date semi-empirical determination of the star formation rate and metallicity evolution of star-forming galaxies throughout cosmic history. We find the PISN rate to exhibit a huge dependence on the model assumptions, including the criterion to identify stars unstable to pair production, and the upper limit of the stellar initial mass function. Remarkably, the interplay between the maximum metallicity at which stars explode as PISNe, and the dispersion of the galaxy metallicity distribution, dominates the uncertainties, causing a ∼ seven-orders-of-magnitude PISN rate range. Furthermore, we show a comparison with the core-collapse supernova rate, and study the properties of the favourable PISN host galaxies. According to our results, the main contribution to the PISN rate comes from metallicities between ∼10−3 and 10−2, against the common assumption that views very-low-metallicity, Population III stars as exclusive or dominant PISN progenitors. The strong dependencies we find offer the opportunity to constrain stellar and galaxy evolution models based on possible future (or the lack of) PISN observations.

STELLA Lightcurves of Energetic Pair-instability Supernovae in the Context of SN2018ibb

SN2018ibb is a recently observed hydrogen-poor superluminous supernova that appears to be powered by the decay of 30 M ⊙ of radioactive nickel. This supernova has been suggested to show hybrid signatures of a pair-instability supernova and an interacting supernova. In a previous paper, we found that rotating, metal-enriched pair-instability supernova progenitors appeared to check both of these boxes. In this paper, we model the lightcurves of the pair-instability supernovae using STELLA. We find that the STELLA models can explain the overall shape of the bolometric lightcurve of SN2018ibb, though not specific morphological features such as the luminosity peak or the bump at roughly 300 days after the peak. We also estimate the contribution from interaction and find that with relatively low wind velocities, the circumstellar medium originating from the stellar winds is consistent with the evidence for interaction in the spectra. The observed values of the photosphere velocity in the 100 days after peak luminosity are similar to the STELLA models, but the deceleration is lower. This leads to the biggest inconsistency, which is the blackbody temperature of SN2018ibb being much hotter than any of the STELLA models. We note that this high temperature (and the flat velocity) may be difficult to reconcile with the long rise time of SN2018ibb, but nevertheless conclude that if it is accurate, this discrepancy represents a challenge for SN2018ibb being a robust PISN candidate. This result is noteworthy given the lack of other scenarios for this supernova.

LAMOST J1010+2358 is not a Pair-Instability Supernova Relic

The discovery of a star formed out of pair-instability supernova ejecta would have massive implications for the Population III star initial mass function and the existence of stars over 100 Msun, but none have yet been found. Recently, the star LAMOST J1010+2358 was claimed to be a star that formed out of gas enriched by a pair-instability supernova. We present a non-LTE abundance analysis of a new high-resolution Keck/HIRES spectrum of J1010+2358. We determined the carbon and aluminum abundances needed to definitively distinguish between enrichment by a pair-instability and core-collapse supernova. Our new analysis demonstrates that J1010+2358 does not have the unique abundance pattern of a a pair-instability supernova, but was instead enriched by the ejecta of a low mass core-collapse supernova. Thus, there are still no known stars displaying unambiguous signatures of pair-instability supernovae.

Primary and secondary source of energy in the superluminous supernova 2018ibb

We examined a possible pair-instability origin for the superluminous supernova 2018ibb. As the base model, we used a non-rotating stellar model with an initial mass of 250 at about 1/15 solar metallicity. We considered three versions of the model as input for radiative transfer simulations done with the stella and artis codes: with 25 of 56Ni, with 34 of 56Ni, and a chemically mixed case with 34 of 56Ni. We present light curves and spectra in comparison to the observed data of SN\,2018ibb, and conclude that the pair-instability supernova model with 34 of 56Ni explains the broadband light curves reasonably well between $-100$ and 250 days around the peak. Our synthetic spectra have many similarities with the observed spectra. The luminosity excess in the light curves and the blue-flux excess in the spectra can be explained by an additional energy source, which may be an interaction of the supernova ejecta with circumstellar matter. We discuss possible mechanisms that could have caused the circumstellar matter to be ejected in the decades before the pair-instability explosion.

On the Pair-instability Supernova Origin of J1010+2358* * Based on observations made with the Very Large Telescope (VLT) of the ESO at the La Silla Paranal Observatory under the program ID 112.26WJ.001

The first (Population III) stars formed only out of H and He and were likely more massive than present-day stars. Massive Population III stars in the range 140–260 M ⊙ are predicted to end their lives as pair-instability supernovae (PISNe), enriching the environment with a unique abundance pattern, with high ratios of odd to even elements. Recently, the most promising candidate for a pure descendant of a zero-metallicity massive PISN (260 M ⊙) was discovered by the LAMOST survey, the star J1010+2358. However, key elements to verify the high PISN contribution, C and Al, were missing from the analysis. To rectify this, we obtained and analyzed a high-resolution Very Large Telescope/UVES spectrum, correcting for 3D and/or non-local thermodynamic equilibrium effects. Our measurements of both C and Al give much higher values (∼1 dex) than expected from a 260 M ⊙ PISN. Furthermore, we find significant discrepancies with the previous analysis and therefore a much less pronounced odd–even pattern. Our results show that J1010+2358 cannot be a pure descendant of a 260 M ⊙ PISN. Instead, we find that the best-fit model consists of a 13 M ⊙ Population II core-collapse supernova combined with a Population III supernova. Alternative, less favored solutions ( χ2/χbest2≈2.3 ) include a 50% contribution from a 260 M ⊙ PISN or a 40% contribution from a Population III Type Ia supernova. Ultimately, J1010+2358 is certainly a unique star giving insights into the earliest chemical enrichment; however, this star is not a pure PISN descendant.

Chemical Diagnostics to Unveil Environments Enriched by First Stars

Unveiling the chemical fingerprints of the first (Population III, hereafter Pop III) stars is crucial for indirectly studying their properties and probing their massive nature. In particular, very massive Pop III stars explode as energetic pair-instability supernovae (PISNe), allowing their chemical products to escape in the diffuse medium around galaxies, opening the possibility to observe their fingerprints in distant gas clouds. Recently, three z > 6.3 absorbers with abundances consistent with an enrichment from PISNe have been observed with JWST. In this Letter, we present novel chemical diagnostics to uncover environments mainly imprinted by PISNe. Furthermore, we revise the JWST low-resolution measurements by analyzing the publicly available high-resolution X-Shooter spectra for two of these systems. Our results reconcile the chemical abundances of these absorbers with those from literature, which are found to be consistent with an enrichment dominated (>50% metals) by normal Pop II SNe. We show the power of our novel diagnostics in isolating environments uniquely enriched by PISNe from those mainly polluted by other Pop III and Pop II SNe. When the subsequent enrichment from Pop II SNe is included, however, we find that the abundances of PISN-dominated environments partially overlap with those predominantly enriched by other Pop III and Pop II SNe. We dub these areas confusion regions. Yet, the odd–even abundance ratios [Mg,Si/Al] are extremely effective in pinpointing PISN-dominated environments and allowed us to uncover, for the first time, an absorber consistent with a combined enrichment by a PISN and another Pop III SN for all the six measured elements.

Impacts of the 12C(α, γ)16O reaction rate on 56Ni nucleosynthesis in pair-instability supernovae

ABSTRACT Nuclear reactions are key to our understanding of stellar evolution, particularly the $^{12}{\rm C}(\alpha ,\gamma)^{16}{\rm O}\,$ rate, which is known to significantly influence the lower and upper ends of the black hole (BH) mass distribution due to pair-instability supernovae (PISNe). However, these reaction rates have not been sufficiently determined. We use the mesa stellar evolution code to explore the impact of uncertainty in the $^{12}{\rm C}(\alpha ,\gamma)^{16}{\rm O}\,$ rate on PISN explosions, focusing on nucleosynthesis and explosion energy by considering the high resolution of the initial mass. Our findings show that the mass of synthesized radioactive nickel (56Ni) and the explosion energy increase with $^{12}{\rm C}(\alpha ,\gamma)^{16}{\rm O}\,$ rate for the same initial mass, except in the high-mass edge region. With a high (about twice the starlib standard value) rate, the maximum amount of nickel produced falls below 70 M⊙, while with a low rate (about half of the standard value) it increases up to 83.9 M⊙. These results highlight that carbon ‘preheating’ plays a crucial role in PISNe by determining core concentration when a star initiates expansion. Our results also suggest that the onset of the expansion, which means the end of compression, competes with collapse caused by helium photodisintegration, and the maximum mass that can lead to an explosion depends on the $^{12}{\rm C}(\alpha ,\gamma)^{16}{\rm O}\,$ reaction rate.

Dark matter annihilation and pair-instability supernovae

Dark matter annihilation and pair-instability supernovae

Contribution of very massive stars to the sulfur abundance in star-forming galaxies: Role of pair-instability supernovae

Context. Recent work presented increasing evidence of high non-constant S/O abundance ratios observed in star-forming metal-poor galaxies that deviated from the constant canonical S/O across a wide range of O/H abundances. Similar peculiarly high Fe/O ratios have also recently been detected. Aims. We investigate whether these high S/O ratios at low metallicities could be explained when the process of pair-instability supernovae (PISN) in chemical modelling is included, through which a similar behaviour of the Fe/O ratios was reproduced successfully. Methods. We used chemical evolution models that considered the stages of PISN in the previously published yields and adopted a suitable initial mass function (IMF) to characterize this evolutionary stage appropriately. Results. The peculiarly high values and the behaviour of the observed S/O versus O/H relation can be reproduced when the ejecta of very massive stars that go through the process of PISN are taken into account. Additionally, a bimodal top-heavy IMF and an initial strong burst of star formation are required to reach the reported high S/O values. Conclusions. We show that the role of very massive stars going through the process of PISN should be taken into account to explain the chemical enrichment of sulphur and oxygen in metal-poor star-forming regions.

1100 days in the life of the supernova 2018ibb

Stars with zero-age main sequence masses between 140 and 260 M⊙ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN 2018ibb is a hydrogen-poor SLSN at z = 0.166 that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the near-infrared (NIR) with 2–10 m class telescopes. SN 2018ibb radiated > 3 × 1051 erg during its evolution, and its bolometric light curve reached > 2 × 1044 erg s−1 at its peak. The long-lasting rise of > 93 rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source (56Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions as to their photometric and spectroscopic properties. SN 2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, and potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25–44 M⊙ of freshly nucleosynthesised 56Ni, pointing to the explosion of a metal-poor star with a helium core mass of 120–130 M⊙ at the time of death. This interpretation is also supported by the tentative detection of [Co II] λ 1.025 μm, which has never been observed in any other PISN candidate or SLSN before. We observe a significant excess in the blue part of the optical spectrum during the nebular phase, which is in tension with predictions of existing PISN models. However, we have compelling observational evidence for an eruptive mass-loss episode of the progenitor of SN 2018ibb shortly before the explosion, and our dataset reveals that the interaction of the SN ejecta with this oxygen-rich circumstellar material contributed to the observed emission. That may explain this specific discrepancy with PISN models. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN 2018ibb by far the best candidate for being a PISN, to date.

The Gaia RVS benchmark stars

Context. The Gaia satellite has already provided the astronomical community with three data releases, and the Radial Velocity Spectrometer (RVS) on board Gaia has provided the radial velocity for 33 million stars. Aims. When deriving the radial velocity from the RVS spectra, several stars are measured to have large values. To verify the credibility of these measurements, we selected some bright stars with the modulus of radial velocity in excess of 500 km s−1 to be observed with SOPHIE at OHP and UVES at VLT. This paper is devoted to investigating the chemical composition of the stars observed with UVES. Methods. We derived atmospheric parameters using Gaia photometry and parallaxes, and we performed a chemical analysis using the MyGIsFOS code. Results. We find that the sample consists of metal-poor stars, although none have extremely low metallicities. The abundance patterns match what has been found in other samples of metal-poor stars selected irrespective of their radial velocities. We highlight the presence of three stars with low Cu and Zn abundances that are likely descendants of pair-instability supernovae. Two stars are apparently younger than 1 Ga, and their masses exceed twice the turn-off mass of metal-poor populations. This makes it unlikely that they are blue stragglers because it would imply they formed from triple or multiple systems. We suggest instead that they are young metal-poor stars accreted from a dwarf galaxy. Finally, we find that the star RVS721 is associated with the Gjoll stream, which itself is associated with the Globular Cluster NGC 3201.

EMPRESS. XIII. Chemical Enrichment of Young Galaxies Near and Far at z ∼ 0 and 4–10: Fe/O, Ar/O, S/O, and N/O Measurements with a Comparison of Chemical Evolution Models

We present gas-phase elemental abundance ratios of thirteen local extremely metal-poor galaxies (EMPGs), including our new Keck/LRIS spectroscopy determinations together with 33 James Webb Space Telescope z ∼ 4–10 star-forming galaxies in the literature, and compare chemical evolution models. We develop chemical evolution models with the yields of core-collapse supernovae (CCSNe), Type Ia SNe, hypernovae (HNe), and pair-instability supernovae (PISNe), and compare the EMPGs and high-z galaxies in conjunction with dust depletion contributions. We find that high Fe/O values of EMPGs can (cannot) be explained by PISN metal enrichments (CCSN/HN enrichments even with the mixing-and-fallback mechanism enhancing iron abundance), while the observed Ar/O and S/O values are much smaller than the predictions of the PISN models. The abundance ratios of EMPGs can be explained by the combination of Type Ia SNe and CCSNe/HNe whose inner layers of argon and sulfur mostly fallback, which are comparable to the Sculptor stellar chemical abundance distribution, suggesting that early chemical enrichment has taken place in the EMPGs. Comparing our chemical evolution models with the star-forming galaxies at z ∼ 4–10, we find that the Ar/O and S/O ratios of the high-z galaxies are comparable to those of the CCSN/HN models, while the majority of high-z galaxies do not have constraints good enough to rule out contributions from PISNe. The high N/O ratio recently reported in GN-z11 cannot be explained even by rotating PISNe, but could be reproduced by the winds of rotating Wolf–Rayet stars that end up as a direct collapse.

True Pair-instability Supernova Descendant: Implications for the First Stars’ Mass Distribution

The initial mass function (IMF) of the first Population III (Pop III) stars remains a persistent mystery. Their predicted massive nature implies the existence of stars exploding as pair-instability supernovae (PISNe), but no observational evidence had been found. Now, the LAMOST survey claims to have discovered a pure PISN descendant, J1010+2358, at [Fe/H] = − 2.4. Here we confirm that a massive 250–260 M ⊙ PISN is needed to reproduce the abundance pattern of J1010+2358. However, the PISN contribution can be as low as 10%, since key elements are missing to discriminate between scenarios. We investigate the implications of this discovery for the Pop III IMF, by statistical comparison with the predictions of our cosmological galaxy formation model, NEFERTITI. First, we show that the nondetection of mono-enriched PISN descendants at [Fe/H] < − 2.5 allows us to exclude (i) a flat IMF at a 90% confidence level; and (ii) a Larson-type IMF with characteristic mass m ch/M ⊙ > 191.16x − 132.44, where x is the slope, at a 75% confidence level. Second, we show that if J1010+2358 has only inherited <70% of its metals from a massive PISN, no further constraints can be put on the Pop III IMF. If, instead, J1010+2358 will be confirmed to be a nearly pure (>90%) PISN descendant, it will offer strong and complementary constraints on the Pop III IMF, excluding the steepest and bottom-heaviest IMFs: m ch/M ⊙ < 143.21x − 225.94. Our work shows that even a single detection of a pure PISN descendant can be crucial to our understanding of the mass distribution of the first stars.

Rapid Chemical Enrichment by Intermittent Star Formation in GN-z11

We interpret the peculiar supersolar nitrogen abundance recently reported by the James Webb Space Telescope observations for GN-z11 (z = 10.6) using our state-of-the-art chemical evolution models. The observed CNO ratios can be successfully reproduced—independently of the adopted initial mass function, nucleosynthesis yields, and presence of supermassive (>1000M ⊙) stars—if the galaxy has undergone an intermittent star formation history with a quiescent phase lasting ∼100 Myr, separating two strong starbursts. Immediately after the second burst, Wolf–Rayet stars (up to 120M ⊙) become the dominant enrichment source, also temporarily (<1 Myr) enhancing particular elements (N, F, Na, and Al) and isotopes (13C and 18O). Alternative explanations involving (i) single burst models, also including very massive stars and/or pair-instability supernovae, or (ii) pre-enrichment scenarios fail to match the data. Feedback-regulated, intermittent star formation might be common in early systems. Elemental abundances can be used to test this hypothesis and to get new insights on nuclear and stellar astrophysics.

Two of a Kind: Comparing Big and Small Black Holes in Binaries with Gravitational Waves

When modeling the population of merging binary black holes, analyses have generally focused on characterizing the distribution of primary (i.e., more massive) black holes in the binary, while using simplistic prescriptions for the distribution of secondary masses. However, the secondary mass distribution and its relationship to the primary mass distribution provide a fundamental observational constraint on the formation history of coalescing binary black holes. If both black holes experience similar stellar evolutionary processes prior to collapse, as might be expected in dynamical formation channels, the primary and secondary mass distributions would show similar features. If they follow distinct evolutionary pathways (for example, due to binary interactions that break symmetry between the initially more massive and less massive stars), their mass distributions may differ. We present the first analysis of the binary black hole population that explicitly fits for the secondary mass distribution. We find that the data is consistent with a ∼30 M ⊙ peak existing only in the distribution of secondary rather than primary masses. This would have major implications for our understanding of the formation of these binaries. Alternatively, the data is consistent with the peak existing in both component mass distributions, a possibility not included in most previous studies. In either case, the peak is observed at 31.4−2.6+2.3M⊙ , which is shifted lower than the value obtained in previous analyses of the marginal primary mass distribution, placing this feature in further tension with expectations from a pulsational pair-instability supernova pileup.