
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39227590
51863
10.1038/s41467-024-51863-z
Article
Binary progenitor systems for Type Ic supernovae
http://orcid.org/0000-0002-3148-1359
Solar Martín martin.solar@amu.edu.pl

1
Michałowski Michał J. michal.michalowski@amu.edu.pl

1
http://orcid.org/0000-0003-1440-9061
Nadolny Jakub 1
http://orcid.org/0000-0002-1296-6887
Galbany Lluís 23
Hjorth Jens 4
http://orcid.org/0000-0002-7464-498X
Zapartas Emmanouil 5
http://orcid.org/0000-0003-1546-6615
Sollerman Jesper 6
http://orcid.org/0000-0001-9162-2371
Hunt Leslie 7
Klose Sylvio 8
http://orcid.org/0000-0001-5785-1154
Koprowski Maciej 9
Leśniewska Aleksandra 14
http://orcid.org/0009-0006-7826-2416
Małkowski Michał 1
Nicuesa Guelbenzu Ana M. 8
Ryzhov Oleh 1
http://orcid.org/0000-0003-2354-3238
Savaglio Sandra 101112
Schady Patricia 13
http://orcid.org/0000-0001-6797-1889
Schulze Steve 14
de Ugarte Postigo Antonio 15
Vergani Susanna D. 1617
http://orcid.org/0000-0002-4465-8264
Watson Darach 1819
Wróblewski Radosław 1
1 https://ror.org/04g6bbq64 grid.5633.3 0000 0001 2097 3545 Astronomical Observatory Institute, Faculty of Physics, Adam Mickiewicz University, Poznań, Poland
2 https://ror.org/01vbgty78 grid.450286.d 0000 0004 1793 4897 Institute of Space Sciences (ICE, CSIC), Campus UAB, Barcelona, Spain
3 https://ror.org/00k6njn28 grid.435450.3 0000 0004 1784 9780 Institut d’Estudis Espacials de Catalunya (IEEC), Barcelona, Spain
4 https://ror.org/035b05819 grid.5254.6 0000 0001 0674 042X DARK, Niels Bohr Institute, University of Copenhagen, Copenhagen N, Denmark
5 grid.4834.b 0000 0004 0635 685X Institute of Astrophysics, FORTH, Heraklion, Greece
6 grid.10548.38 0000 0004 1936 9377 The Oskar Klein Centre, Department of Astronomy, Stockholm University, Albanova University Center, Stockholm, Sweden
7 https://ror.org/00fbze943 grid.426239.8 0000 0000 9176 4495 INAF—Osservatorio Astrofisico di Arcetri, Firenze, Italy
8 https://ror.org/05jyz4092 grid.440503.6 0000 0004 0646 0278 Thüringer Landessternwarte Tautenburg, Tautenburg, Germany
9 https://ror.org/0102mm775 grid.5374.5 0000 0001 0943 6490 Institute of Astronomy, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Toruń, Poland
10 https://ror.org/02rc97e94 grid.7778.f 0000 0004 1937 0319 Department of physics, University of Calabria, Arcavacata di Rende (CS), Italy
11 https://ror.org/00gqsp710 INAF—Osservatorio di Astrofisica e Scienza dello Spazio, Bologna, Italy
12 https://ror.org/049jf1a25 grid.463190.9 0000 0004 0648 0236 INFN—Laboratori Nazionali di Frascati, Frascati, Italy
13 https://ror.org/002h8g185 grid.7340.0 0000 0001 2162 1699 Department of Physics, University of Bath, Bath, UK
14 grid.10548.38 0000 0004 1936 9377 The Oskar Klein Centre, Department of Physics, Stockholm University, Albanova University Center, Stockholm, Sweden
15 grid.460782.f 0000 0004 4910 6551 Artemis, Observatoire de la Cotê d’Azur, Université Cotê d’Azur, Nice, France
16 grid.440907.e 0000 0004 1784 3645 GEPI, Observatoire de Paris, PSL University, CNRS, Meudon, France
17 https://ror.org/022bnxw24 grid.435813.8 0000 0001 0540 8249 Institut d’Astrophysique de Paris, UMR 7095, CNRS-SU, Paris, France
18 grid.5254.6 0000 0001 0674 042X The Cosmic Dawn Center (DAWN), Copenhagen, Denmark
19 https://ror.org/035b05819 grid.5254.6 0000 0001 0674 042X Niels Bohr Institute, University of Copenhagen, Copenhagen N, Denmark
3 9 2024
3 9 2024
2024
15 766726 8 2023
16 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Core-collapse supernovae are explosions of massive stars at the end of their evolution. They are responsible for metal production and for halting star formation, having a significant impact on galaxy evolution. The details of these processes depend on the nature of supernova progenitors, but it is unclear if Type Ic supernovae (without hydrogen or helium lines in their spectra) originate from core-collapses of very massive stars (>30 M⊙) or from less massive stars in binary systems. Here we show that Type II (with hydrogen lines) and Ic supernovae are located in environments with similar molecular gas densities, therefore their progenitors have comparable lifetimes and initial masses. This supports a binary interaction for most Type Ic supernova progenitors, which explains the lack of hydrogen and helium lines. This finding can be implemented in sub-grid prescriptions in numerical cosmological simulations to improve the feedback and chemical mixing.

Type Ic supernovae (SNe) originate either from the core collapse of very massive stars or from less massive stars in binary systems. Here, the authors show that progenitors of Type II and Ic SNe have comparable lifetimes and initial masses, which supports binary interaction for most Type Ic SNe progenitors.

Subject terms

Stellar evolution
Interstellar medium
501100004281 Narodowe Centrum Nauki (National Science Centre) 2018/30/E/ST9/00208 issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Massive stars (>8 M⊙) have a significant impact on the interstellar medium (ISM) by regulating it through stellar winds, ionising radiation, and supernova (SN) explosions. SNe contribute to the origin of heavy elements, a still poorly understood aspect in cosmology1. The details of SN feedback and metal production depend primarily on which stars explode as which type of SNe because the mechanism of explosion and element production yield is different for each of them. SN feedback models that take into consideration progenitor stars and mechanisms of explosions are essential to improve simulations.

The connection between SN types and their progenitors is of great importance for galaxy evolution and cosmology. This can be directly studied by identifying progenitor stars on pre-explosion images, with follow-up observations that confirm their disappearance. However, to date, only 23 (18 Type II, four Type IIb, and one Type Ib) core-collapse SN (CCSN) progenitor stars have been confirmed to disappear in post-explosion images2. No Type Ic SN (without hydrogen or helium in the spectrum) progenitor has been confirmed in this way, with SN 2017ein as the only candidate, but with a wide range of derived progenitor masses3–6. Hence, most of our knowledge on the nature of Type Ic SN progenitors is based on photometric and spectroscopic observations after the explosion.

Type Ic SNe likely originates either from core-collapses of very massive stars or from less massive stars in binary systems. The mechanism responsible for the lack of hydrogen or helium lines is still a subject of debate. In the case of the very massive star model, the hydrogen and helium envelope is lost due to stellar winds from the star itself7,8 (e.g., radiatively driven winds, episodic mass-loss, or rapid rotation). In the binary model, the companion radically affects the evolution of the progenitor due to the mass exchange9–11. Non-detections of Type Ic SN progenitors in pre-explosion images and relative rates of different types of SNe suggest that most of Type Ic SN progenitors are binaries with initial masses  <  20 M⊙12,13.

Another way to address the relation between progenitor stars and resulting SNe is to investigate the molecular gas properties at the explosion location. Molecular gas at the locations of SNe of different types was recently investigated at a spatial resolution comparable to giant molecular clouds (GMCs)14, using the Atacama Large Millimetre Array (ALMA) carbon monoxide 2-1 line transition [CO(2-1)] observations from the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)15,16 survey. Their sample consisted of a total of 59 SNe: 12 thermonuclear (Type Ia SNe), 32 Type II SNe, eight stripped-envelope SNe (SESNe, hereafter, Type Ib, Ic or Ib/c), and seven unclassified. They found that Type Ia and II SNe are associated with little or no molecular gas emission, while SESNe and unclassified SNe mostly show strong molecular gas emission. They concluded that there is a clear dependence of the type of SN and the molecular gas environment, however, their conclusions are drawn based on a low sample size for SESNe and, thus are not statistically significant.

In this work, our goal is to constrain lifetimes and initial masses of Type Ic SN progenitors. To this end, we compare the molecular gas densities at the positions of Type II and Ic SNe. By targeting a large sample of SNe we aim to uncover their nature. This statistical approach offers strong constraints on the overall progenitor characteristics of different SN populations but does not provide a strong constraint on individual SN progenitor properties. We report a statistically significant study to do so with spatial resolution comparable to the GMC sizes. This is an important factor because molecular hydrogen column surface densities and lifetimes of GMCs can only be measured accurately if the resolution at least matches the cloud sizes17,18.

Results

In order to investigate the environments of a significant number of SNe at high resolution, we initiated the ALMA CO SN (ACOS) survey, obtaining CO(2-1) observations of the locations of 16 Type Ic SNe. Together with the PHANGS survey this results in a sample of 63 SNe: 12 Type Ia, 30 Type II, and 21 Type Ic SNe. These CO(2-1) observations have a spatial resolution of ~ 100 pc, similar to sizes of GMCs. The spatial resolution and the large sample allow us to study the immediate environments in which the SNe exploded. Our main conclusion is from the comparison of Type II and Ic SNe, whereas Type Ia SNe are shown only to contrast the different progenitor natures. See 'Methods', subsections ALMA CO SN survey, PHANGS–ALMA data, and Supernova sample for detailed information about the SNe and their host galaxies. Supplementary Data 1 lists the information for the SNe used in this work.

As an example, Fig. 1 shows the molecular hydrogen gas surface density (Σmol) map of NGC 4254 (M 99, a typical PHANGS–ALMA galaxy) with its four CCSNe; three Type II (SN 1967H, SN 1972Q, and SN 1986I) and one Type Ic (SN 2014L) SNe, plus the location of SN 2009bb (Type Ic SN) hosted in NGC 3278 (observed by ACOS). Σmol is computed from the CO(2-1) line intensity, see methods, subsection CO-to-Σmol for the description used.Fig. 1 Distribution of molecular gas in galaxy host and SN locations.

SN 1967H, SN 1972Q and SN 1986l hosted in NGC 4254 (a), and SN 2009bb hosted in NGC 3278 (b). Colour-coded Σmol intensity is represented in logarithmic scale. Pixels without signal are masked and shown as white. Red dots represent SN locations. North is up and East is left.

The Σmol value for a given SN was calculated in two ways. First, we measured it at the exact pixel of the SN location and this is denoted as “SN location”. However, the SN location might not be the exact site where the progenitor star formed. The true location of the formation of the SN progenitor star could be shifted due to astrometric displacement and/or peculiar motion of the progenitor system with respect to the parent GMC. To take into account these effects, together with the maximum size of GMCs, we also measured the maximum value within a radius of 200 pc centred at the SN position (see methods, subsection 200 pc regions). To assess an average level of Σmol for our host galaxies, we also measured Σmol for 103 random pixels in each PHANGS galaxy and the respective maximum value within a radius of 200 pc centred at these random pixels. The ACOS sample was not included in the random pixel calculations because these observations cover the SN positions but not the whole galaxy, so it is not possible to obtain a representative measurement of Σmol for these galaxies.

Figure 2 shows the empirical cumulative distribution function (eCDF) of Σmol in the environments of SNe with their respective confidence intervals calculated using the Dvoretzky–Kiefer–Wolfowitz inequality19 at a significance level of 68%. If a given Σmol measurement has a significance lower than 2σ (non-detection), where σ is the noise at the pixel location, then the 2σ upper limit is shown on the eCDF. This happens for 66% of the random pixels (15077 out of 23000), 67% of Type Ia (eight out of 12), 37% Type II (11 out of 30) and 43% Type Ic SNe (nine out of 21). These upper limits are not used for any calculations, only for the visualisation of the eCDFs. Supplementary Data 2 and 3 lists the Σmol environments of SNe and random locations in galaxy hosts, respectively. Random pixel locations extend towards lower values when compared to the SN positions, indicating that the SN locations have a strong connection with their GMC parents or larger molecular structures. This is an expected result for CCSNe, since they are associated with massive progenitor stars and therefore they are expected to explode close to their birthplaces in spiral arms where molecular gas densities are higher20.Fig. 2 Σmol eCDFs for SN locations.

Random locations, Type Ia, Type II, and Type Ic SNe are represented by black, blue, red, and grey lines, respectively. Upper limits (2σ) are used in case of non-detections. The shaded areas represent confidence intervals at 1σ.

To quantify if our samples are drawn from different parent GMC populations, we performed a Kolmogorov—Smirnov (KS) test for each combination. KS statistics and p-values are shown in Table 1 (and for the 200 pc regions in Supplementary Table 1). Comparing the locations of CCSNe and random pixels, the low p-values indicate that their distributions are different, suggesting high densities of molecular gas environments for the sites where SNe was observed. On the other hand, high p-values (> 0.05) for SN combinations indicate that it is not possible to reject the null hypothesis that the Σmol at the SN positions are from the same distribution.Table 1 KS statistics (and p-value in parenthesis) for different SN groups (SN locations)

	Random	Ia	II	Ic	
Random	0.0 (1.0)	3.1e-01 (1.6e-01)	4.5e-01 (4.7e-06)	4.1e-01 (1.0e-03)	
Ia		0.0 (1.0)	4.0e-01 (1.0e-01)	2.4e-01 (7.1e-01)	
II			0.0 (1.0)	2.5e-01 (3.7e-01)	
Ic				0.0 (1.0)	

In order to obtain confidence interval ranges of Σmol for each SN type and random locations, we performed Monte Carlo simulations. Each measurement of an individual SN was perturbed according to its uncertainty, which was assumed to have a Normal distribution. We created 104 of such perturbed sets, each time calculating its median value. The uncertainty was adopted to be at 16% and 84% of the simulated distribution. We show the results of these simulations in Fig. 3. The medians and 1σ confidence intervals for SN locations and 200 pc regions are summarised in Supplementary Table 2.Fig. 3 Median values and 1σ confidence intervals (using 104 Monte Carlo simulations) of Σmol for SN locations.

Random locations, Type Ia, Type II, and Type Ic SNe are represented by black, blue, red, and grey error bars, respectively.

Median molecular gas densities increase from the values measured for the random pixels (4.47−0.04+0.05M⊙pc−2), through Type Ia SNe (6.93−2.36+3.70M⊙pc−2), to Type II and Ic SNe (20.15−2.46+3.38 and 20.62−4.88+4.28M⊙pc−2, respectively). The results are shown in Fig. 3. The random locations and Type Ia SNe have much lower median molecular gas densities than CCSNe. At the positions of Type II and Ic SNe we obtained similar molecular gas densities within the 1σ confidence levels.

Under the assumption of single very massive star progenitors for Type Ic SNe (for masses of 30–100 M⊙ with lifetimes of 7–3 Myr, respectively21), it is expected that the respective parent GMC would not have been dispersed before the SN explosion due to a short progenitor lifetime, and the progenitors would not have enough time to shift away from their birthplaces significantly. Therefore, if progenitors of Type Ic SNe were very massive stars, then the distribution of Σmol at their positions should be shifted toward higher values compared to that of Type II SNe. This is because the lower masses of progenitors of Type II SNe imply longer lifetimes, and therefore more time for the parent clouds to disperse and for the progenitors to move away. This scenario is not supported by our results. In the alternative scenario, the progenitors of Type II SNe evolve as single stars (or in wide binaries in which their hydrogen layers are not affected) and those of Type Ic SNe are similarly massive stars that evolve in binary systems with a companion being responsible for removing the external layers of hydrogen and helium21. Then the progenitor masses, and therefore lifetimes, of both types are similar. Thereby, their distributions of Σmol should be comparable, which is consistent with our data. See methods, subsection Timing the SN progenitor lifetime with molecular gas information for the justification of using the molecular gas densities to constrain the stellar population age.

We note that at lower gas densities the distribution of Type Ic SNe is slightly lower than that of Type II SNe (see Fig. 2). This may indicate that the lifetimes of some of the Type Ic progenitors were increased by the binary interaction21. However, the statistical significance of this difference is too low to draw any definitive conclusions.

In order to assess the maximum difference between the lifetimes of progenitors of Type II and Ic SNe, we assumed that the molecular gas density at the SN progenitor positions, Σmol,SN, decreases exponentially with time as1 Σmol,SN=Σ0e−tSN/τGMC,

where Σ0 is a normalisation constant, tSN is the lifetime of the SN progenitor and τGMC is the lifetime of the GMC. In order to constrain the lifetimes of Type Ic SN progenitors, we made use of our measured ratio of Σmol,SN for Type Ic and II SNe (Σmol,Ic/Σmol,II=1.02−0.27+0.27) and a measured value of the GMC lifetime of τGMC=16−5+6Myr22.

This characteristic cloud evolution lifetime is in agreement with theoretical works and simulations18,22–27 and while the exact value of the assumed average GMC lifetime influences these calculations, it does not change the interpretation when lifetimes of two SN types are compared. Assuming that progenitors of Type II and Ic SNe are born in GMCs with similar average initial conditions, i.e. that average Σ0 is the same for both, from eq. (1) it is possible to calculate the difference between the SN progenitor lifetimes as2 tII−tIc=τGMClnΣmol,Ic/Σmol,II=0.37−4.26+4.27Myr.

A zero-age main sequence (ZAMS) mass Minit,SN for Type II SN progenitors of Minit,II=10.66−0.20+0.20M⊙ (obtained by averaging nine SNe with pre-explosion detection28 and confirmed by the disappearance in post-explosion images2) yields a lifetime of tII = 25.22−0.80+0.80Myr. Using this tII in eq. (2) results in a lifetime for Type Ic SN progenitors of tIc=24.9−4.3+4.3Myr and a ZAMS mass of Minit,Ic = 10.90−1.20+1.20M⊙. On the other hand, if we assume a typical mass for red supergiant progenitors for type II SNe29 of Minit,II=15−1+1M⊙, then we obtain Minit,Ic = 15.3−3.2+3.2M⊙. This also means that if Type II SN progenitors include rare examples of very massive stars, so can Type Ic SN progenitors30. To account for signal unrelated to SNe, as the first step we also subtracted the random values of Σmol from those of SNe, which resulted in tII−tIc=0.47−5.45+5.47Myr, tIc=24.8−5.5+5.5Myr, and Minit,Ic = 10.9−1.5+1.5M⊙, indistinguishable from the original results.

Another effect to take into account is that SN progenitors may be runaway stars, which are moving away from their parent clusters with significant velocities. Maximum velocities for OB runaway stars are ~ 30 km s−1 or ~ 30 pc/Myr31, so they would need ~ 3 Myr to cross a GMC. Replacing τGMC with the effective crossing timescale τcro, Eq. (1) will then be smaller, including both the movement of the progenitors and the cloud dispersal. Using τcro = 3 Myr, Eq. (2) yields an even smaller lifetime difference between Type II and Ic SNe of ~ 0.06 Myr, what makes our conclusions even stronger.

If most of Type Ic SN progenitors were very massive stars with masses around 30 M⊙, then for their lifetime of 7 Myr, Eq. (1) results in molecular gas densities a factor of 4 higher than those at the positions of Type II SNe, which is not in agreement with our results.

Discussion

Our findings indicate that the binary interaction model (mass transfer due to a companion) is the main mechanism extracting outer layers for most Type Ic SN progenitors. However, we remark that we do not reject the possibility that strong stellar winds of a high-mass star can blow away the layers of hydrogen and helium leading to an explosion as a Type Ic SN, individual Type Ic SNe can be due to very massive star progenitors32. Indeed our calculations show that the 1σ range of the accepted fraction of low-mass stars in the Type Ic SNe population is 66–100% (see methods, subsection Statistical significance of the sample). Another possible binary progenitor scenario includes mergers or accretion of significant mass from a companion so that an initially low-mass star becomes massive and luminous enough to launch strong winds shedding the outer envelope before explosion.

Our results are consistent with low measurements of Type Ic SN progenitor masses from light curve modelling33, the comparison of the SN rates12, the modelling of emission lines at the positions of SNe Type Ib/c29, and the direct observational evidence of a binary system for Type Ic SN 2022jli34. Moreover, in our Galaxy, ~ 70% of massive O-type stars are formed in close binary systems and are expected to experience mass transfer during their lifetime35. All these works support our conclusion about binary systems as progenitors of Type Ic SNe.

On the other hand, our low progenitor masses for Type Ic SNe are inconsistent with some other measurements. The 56Ni yields of Type Ic SNe are five times higher than those of Type II SNe36. However, this may not necessarily imply higher progenitor masses, but different explosion mechanisms or energy sources37. If Type Ic SN explosions are asymmetric or not predominantly powered by the radioactive decay of 56Ni, then the nickel mass will be overestimated for them.

Moreover, using PHANGS data, ref. 14 reported that SESNe (Type IIb, Ib, Ib/c, and Ic SNe) are associated with stronger CO(2-1) emission than those for Type Ia and II SNe. Their sample included 59 SNe in 28 galaxies, with the Type Ic sample size being the main difference compared to this work (they included only five of them, see methods, subsection Statistical significance of the sample for a statistical comparison).

SESNe (including Type Ic SNe) are associated with higher Hα intensities than Type II SNe, suggesting younger star-forming regions and higher progenitor masses38–44. However, these results are complicated by the fact that the Hα emission disappears on a timescale of several Myr only for isolated star-forming regions, whereas for larger complexes (several hundreds of pc, as probed by these observations), this timescale may be as long as 20 Myr45. Low ages for Type Ic SN locations were also obtained from population fitting of stars in the SN vicinity46,47, but the ages may be underestimated due to blending of stars47. Moreover, stronger associations of Type Ic SNe with Hα and UV-bright stars may also be explained within a scenario in which their lower-mass progenitors prefer regions of high stellar density or more top-heavy initial mass function (IMF), which increase the binarity fraction, but also the Hα and UV emission. High-resolution multi-wavelength observations are needed to alleviate this tension.

Broad-lined Type Ic (Type Ic-BL) SNe and gamma-ray bursts, both known to be connected with very massive progenitor stars48, show twice the amount of synthesised 56Ni than regular Type Ic SNe49. Our sample contains only four Type Ic-BL SNe and when they are excluded from the Type Ic SN sample, the results remain the same (at the 1σ level). Similarly, when Type IIn and IIn/LBV (2 in total) are removed from the Type II SN sample, the results do not change.

The lack of statistically significant difference between Type Ia and II SN positions could be due to a low number of statistics for the former. Subtracting the median molecular density at the random positions from that of Type Ia SNe yields a non-detection: Σmol=1.93−2.36+3.70M⊙pc−2. Using a corresponding 2σ upper limit in Eq. (2) gives a lower limit on the lifetime difference between the Type Ia and II SN progenitors of tIa − tII > 12 Myr, consistent with the expected higher lifetimes of Type Ia SN progenitors.

Implementation of our findings into cosmological simulations will have important impact on our understanding of the feedback processes. The progenitor type can be implemented in sub-grid prescriptions in numerical cosmological simulations50–52 with regards to SN feedback and chemical yields into the ISM. Moreover, whether the progenitors of Type Ic SNe are binary systems or very massive stars changes their contribution to the formation of heavy elements, one of the key aspects of stellar and galaxy evolution1. This includes carbon, one of the most fundamental elements in the Universe and building block of life. At solar metallicity, within the Local Universe (including the Milky Way), stellar winds and SN explosions from binary–stripped stars are found to produce twice more 12C than similar single stars53. The fourth most abundant element in the Solar system is carbon (after hydrogen, helium, and oxygen) and a significant contribution could be produced from Type Ic SNe. This method can also be applied to larger samples divided into different properties (e.g., explosion characteristics, host galaxy types, environmental metallicities, etc) and more rare events to learn about their nature.

Methods

ALMA CO SN survey

ACOS (ALMA ID 2021.1.00099.S, P.I. M.J.M.) consists in observations of the J = 2 → 1 transition of the 12CO line in the environment of 16 Type Ic SNe with an angular spatial resolution range of 0.4–1.1″ so that the physical resolution is around 50–100 pc. The selection criteria were the observability with ALMA, i.e. declination < 20°, and distances < 55 Mpc (redshift z < 0.013) to allow the detection of individual GMCs in a reasonable observing time. We excluded seven hosts which are edge-on, for which projection effects would make it difficult to measure the gas surface density at the SN position. The distances, masses and luminosities of these galaxies are comparable to PHANGS–ALMA galaxies.

PHANGS–ALMA data

PHANGS survey provides CO(2-1) line observations using ALMA for 74 galaxies in the Local Universe (<20 Mpc), which mostly are face-on (i < 75o). The typical resolution was ~ 2″, corresponding to ~ 100 pc comparable to the sizes of GMCs. We used calibrated data from ref. 15 for the PHANGS–ALMA galaxy sample16. Further information about these procedures can be found in refs. 15,16.

Supernova sample

Our SN sample was compiled from the Open Catalogue for Supernova (https://github.com/astrocatalogs/supernovae/) in April 202254. The SN compilation, as designated in the catalogue, consists of Type Ia (Ia, Ia-02cx, and IaPec), Type II (II, IIP, IIn and IIn/LBV), and Type Ic (Ic and Ic-BL). We did not consider Type I, Ib/c, IIb or SNe without classification because they cannot be categorised as any of our well-defined thermonuclear and CCSNe. Only two Type Ib SNe exploded in PHANGS galaxies, so they were also excluded from the analysis. We obtained a total of 63 SNe hosted in 39 galaxies (either single or multiple SNe located in one galaxy): 16 Type Ic SNe within ACOS and 47 SNe in 23 PHANGS galaxies. Out of these 47 SNe, 12 were thermonuclear Type Ia, 30 Type II and five Type Ic. We note that the sample may be biased by observational limitations, so it can miss dust-embedded SNe and those in the outskirts of the PHANGS galaxies not covered by the ALMA data. We also do not constrain SNe that do emerge from very massive stars that have not been able to escape from their parent GMCs but are so extinguished that they do not make it into an optical sample. This can only be addressed by an infra-red-selected SN sample. We also do not constrain the properties of stars which collapse directly to black holes (BHs)55, because they do not feature in our SN samples. However, their exclusion does not affect significantly our results (see methods, subsection Statistical significance of the sample).

CO-to-Σmol conversion

The CO(2-1) velocity-integrated intensities from moment-0 maps were converted to Σmol using the following equation (eq. 10 from ref. 16):3 Σmol=αCO1−0R21−1ICO(2-1)cosi,

where αCO1−0 is the CO(1-0) conversion factor, R21 is the CO(2-1)-to-CO(1-0) line ratio, i is the inclination angle of the galaxy, and ICO(2-1) is the line-integrated CO(2-1) intensity. We adopt a Galactic CO-to-H2 conversion factor of αCO1−0=5M⊙pc−2(Kkms−1)−1, the same as in ref. 56, and a line ratio of R21 = 0.5, from ref. 57. Inclination angles were taken from ref. 58 for the PHANGS sample, and from the Hyperleda (http://leda.univ-lyon1.fr/) galaxy database59 for ACOS galaxy hosts.

200 pc regions

An SN explosion could be located away from the centre of a cloud. In order to have a better understanding of the parent GMCs, the maximum value Σmol in a circumference within a radius of 200 pc centred on the SN position was also calculated and denoted as “200 pc region”. In the Milky Way, this radius is comparable to the maximum GMC size60. Supplementary Fig. 1 shows the Σmol eCDFs in such 200 pc regions, with a clear shift to higher densities compared with SN locations, as expected. The two-sample KS test from Supplementary Table 1 shows high probabilities that each of the location pairs is drawn from the same distribution. The median and 1σ values obtained via Monte Carlo simulations are shown in Supplementary Fig. 2. The fact that Type II SNe reach molecular gas densities higher than Type Ic SNe strengthens our conclusions.

Cosmological model

We use the nine-year Wilkinson Microwave Anisotropy Probe cosmological model61 with parameters H0 = 69.32 km s−1 Mpc−1, ΩΛ = 0.7134, and Ωm = 0.2865. Redshift values were used only to compute the 200 pc region for each SN and have no influence on the physical interpretation of the results.

Timing the SN progenitor lifetime with molecular gas observations

The use of molecular gas density as a tracer of stellar age is justified by the strong correlation between the age distribution and the cluster-GMC distance62–64. Moreover, there is a close association between the birth environment (i.e. GMC separation) and age of the cluster, measured by the equivalent width (EW) of Hα [EW(Hα)]65. Finally, the analysis of stellar associations and the CO(2-1) emission revealed that the percentage of overlap between the regions of stellar associations and GMCs is ~ 60%66.

In order to test if there is a correlation between molecular gas densities and stellar ages in the PHANGS sample, a pixel-to-pixel comparison was computed for Σmol and EW(Hα), a proxy for age. The Hα maps were obtained from the Multi Unit Spectroscopic Explorer (MUSE)67. The continuum maps were collected from the Wide Field Imager (La Silla’s 2.2m MPG/ESO telescope)68 and also available in the PHANGS–MUSE dataset.

Supplementary Fig. 3 shows the relation of Σmol and EW(Hα) for every pixel of our galaxy sample, i.e. 11 PHANGS galaxies (NGC 628, NGC 1087, NGC 1365, NGC 1385, NGC 1433, NGC 1566, NGC 1672, NGC 3627, NGC 4254, NGC 4303, and NGC 4321), with both ALMA and MUSE data, and hosting at least one SN from our sample. There is a clear correlation with pixels with lower EW(Hα) (older) having lower molecular gas density, The scatter is significant, but we take the scatter of this magnitude into account in our significance test below (and this justifies the need of a sample of the order of a few tens of SNe).

Statistical significance of the sample

To assess the statistical significance of our results with respect to the sample size, we generated 104 sets of synthetic parent GMC densities for 30 Type II SNe (as in our data) and a variable number of very massive stars in order to test if we can distinguish them. We have done it in three ways, first starting from the measured gas densities of Type II SNe (method 1), second starting from the measured gas density distribution in PHANGS galaxies (method 2), and last from lifetimes of binary systems from a numerical model (method 3).

For the former case, in order to have a realistic distribution of GMC densities we need to remove the outliers of Σmol,II data because their high values do not correspond to densities of single GMCs (as we intend to probe), but the accumulation of GMCs along the line-of-sight towards to galaxy centres, where indeed, all identified outliers are located. We obtained the first, second, and third quartiles of Σmol,II (Q1, Q2, and Q3, respectively) and considered outliers as values lower than Q1 − 1.5 ∗ IQR or higher than Q3 + 1.5 ∗ IQR, where IQR = Q3 − Q1 is the interquartile range. After removing outliers (292, 486, 515, 1442, and 5599 M⊙pc−2, higher than Q3 + 1.5 ∗ IQR = 157 M⊙ pc−2), we found an analytical function which best reproduces the distribution of Type II SN Σmol locations by fitting around ~80 different distributions69. The best function was a generalised normal continuous random distribution f(x,β)=β2Γ(1/β)e−∣x∣β, where x is a real number, β > 0 is the shape parameter, and Γ is a gamma function. The fitted parameters were β = 0.51, centred at 11.7 with a scale of 3.52. From this distribution, we constructed two different synthetic distributions corresponding to Type II SNe and very massive stars to assess our ability to distinguish them. For Type II SNe we randomly drew from the function fitted above. For the very massive stars, we made use of Eq. (2) to derive their median Σmol,massive=Σmol,IIe(tII−tmassive)/τGMC=4Σmol,II and drew from a similar function scaled by this factor. In this calculation we assumed the initial mass for Type II SN progenitor of Minit,II = 11 M⊙, corresponding to a lifetime of tII = 25 Myr and an initial mass of Minit,massive = 30 M⊙, corresponding to a lifetime of tmassive = 3 Myr. Finally, we assumed τGMC = 16 Myr.

In the second method, for each SN we drew a random progenitor age from a normal distribution of 25 ± 5 Myr and 3 ± 1 Myr for Type II SNe and very massive stars, respectively, and the lifetime of the GMC of τGMC = 16 ± 5 Myr. We also drew an initial GMC gas density from a distribution with a mean value 0.5 dex higher than the observed distribution70 (because the initial densities were higher for all of the observed clouds) and the same width, so that log(Σ0/M⊙pc−2)=2.0±0.5. The value of this parameter has no influence on the results, as this is only a normalisation and was chosen so that the median of the simulated distribution for Type II SNe is consistent with the observed value. Then we evolved the clouds as an exponential decay to calculate the surface densities at the time of the SN explosions (Eq. (1)).

For the third test, in order to take into account the effect of binarity in a simplified way, we drew samples of Type II SN progenitors and massive stars from the ranges of 6 M⊙ – Mthresh and Mthresh – 100 M⊙, respectively, for Mthresh = 15, 20, 25, and 30 M⊙, weighting with the Kroupa IMF. We randomly assigned an age according to the age probability distribution of SN progenitors for a given initial mass according to the models of ref. 21. Prior mass exchange of the progenitor with its companion leads in general to a longer lifetime. The two mass ranges represent tentative progenitors of Type II SNe and stripped-envelope SNe, although various binary scenarios violate this threshold. In a way, this test takes into account the change of lifetimes due to binarity, without accounting for a possible change in the SN type due to it.

For each method and for each simulated pair of sets (Type II SNe and very massive stars), we performed the KS test in order to check if we could reject the incorrect-by-design null hypothesis that they are drawn from the same distribution. Supplementary Fig. 4 shows the percentages of p-values below 0.05 (to reject the null hypothesis) and 0.37 (measured value from Table 1) for 104 Monte Carlo simulations from a KS two-sample test between the distributions of the 30 random values of Type II SNe and the massive stars constructed above, as a function of sample size for such massive stars. With the sample size of 21, as in our sample of Type Ic SNe, in these simulations, in  ~96% of the cases we obtained the p-value lower than 0.05 (and in 99.9% of cases lower than the measured value of 0.37). This means that we have statistical significance to correctly reject the null hypothesis and if Type Ic SNe were very massive stars, then we would obtain a lower p-value than we measured for virtually all the cases, so our data have enough statistical significance to rule out the very massive star hypothesis.

We also tested how the data can constrain a mixed Type Ic SN population, by analysing the fraction of the simulations with higher p-values than measured when we replaced some of the massive stars by lower mass progenitors in the same range as we assumed for Type II SNe. The 1σ range (68% of the simulated samples having a p-value higher than measured) of the accepted fraction of low-mass stars in the Type Ic SN population is 66–100%. Hence, only a third of the Type Ic SNe could be very massive stars, so that we could still measure the high p-value.

Moreover, in method 2, instead of drawing ages from normal distributions, we also drew masses according to the Kroupa IMFs71 and calculated their lifetimes according to the relationship of ref. 21. In this case, for all the values of Mthresh listed above, the number of the simulated samples having p-values lower than the measured value decreased from 99.9% to 97–98%. Finally, none of these calculations was significantly affected by the exclusion of the mass ranges for which no SNe are expected, due to a direct collapse into BHs, i.e. within the ranges 22–25 and 27–60 M⊙55. If this is taken into account the significance increases by 1–2% due to making the difference between the Type II SNe and massive stars more pronounced.

Lifetime–initial mass relation

We converted the ZAMS masses to lifetimes using the lifetime–initial mass relation for single stars from ref. 21 (see their Fig. 1).

Supplementary information

Supplementary information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51863-z.

Acknowledgements

M.S., M.J.M., J.N., and A.L. acknowledge the support of the National Science Centre, Poland through the SONATA BIS grant 2018/30/E/ST9/00208. M.J.M. acknowledges the support of the Polish National Agency for Academic Exchange Bekker grant BPN/BEK/2022/1/00110. L.G. acknowledges financial support from the Spanish Ministerio de Ciencia e Innovación (MCIN), the Agencia Estatal de Investigación (AEI) 10.13039/501100011033, and the European Social Fund (ESF) ”Investing in your future” under the 2019 Ramón y Cajal program RYC2019-027683-I and the PID2020-115253GA-I00 HOSTFLOWS project, from Centro Superior de Investigaciones Científicas (CSIC) under the PIE project 20215AT016, and the programme Unidad de Excelencia María de Maeztu CEX2020-001058-M. This work was supported by research grants (VIL16599, VIL54489) from VILLUM FONDEN. EZ acknowledges support from the Hellenic Foundation for Research and Innovation under the “3rd Call for H.F.R.I. Research Projects to Support Post-Doctoral Researchers” (Project No: 7933). O.R. acknowledge the support of the National Science Centre, Poland through the grant 2022/01/4/ST9/00037. This research was funded in whole or in part by the National Science Centre, Poland (grant numbers: 2020/39/D/ST9/03078 and 2021/41/N/ST9/02662). We acknowledge David Alex Kann who passed away before the submission of this manuscript and contributed to the writing of the observing proposal and interpretation of the data. The Cosmic Dawn Center (DAWN) is funded by the Danish National Research Foundation under grant DNRF140. D.W. is co-funded by the European Union (ERC, HEAVYMETAL, 101071865). Views and opinions expressed are, however, those of the authors only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. Based on observations taken as part of the PHANGS–MUSE large programme67. This research has made use of the services of the ESO Science Archive Facility. This research has made use of the services of the ESO Science Archive Facility. Science data products from the ESO archive may be distributed by third parties, and disseminated via other services, according to the terms of the Creative Commons Attribution 4.0 International license (https://creativecommons.org/licences/by/4.0/). Credit to the ESO origin of the data must be acknowledged, and the file headers preserved. This work is based [in part] on observations made with the Spitzer Space Telescope, which was operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. We acknowledge the usage of the HyperLeda database (http://leda.univ-lyon1.fr).

Author contributions

M.S. performed most of the analysis and led writing of the manuscript. M.J.M. conceived the idea, led the ALMA proposal 2021.1.00099.S on which this work is based, and performed two of the significance tests. M.S., M.J.M., and J.N. coordinated the project. L.G., J.H., E.Z., and J.S. provided significant contributions to the interpretation of the data. L.G., J.H., L.H., S.K., M.K., A.L., M.M., A.M.N.G., Sandra S., P.S., Steve S., J.S., A.dUP., S.D.V., and D.W. contributed to the writing of the observing proposal. R.W. supported the data analysis and improved the text. M.M. and O.R. compiled the SN list. All the authors contributed to writing the manuscript.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

ACOS imaging is available from https://almascience.eso.org/aq/under the proposal ID 2021.1.00099.S (P.I. M.J.M.). PHANGS–ALMA and PHANGS–MUSE imaging are available from https://www.canfar.net/storage/vault/list/phangs/RELEASES/SNe information (name, astrometric positions, type, and redshift) were collected from https://github.com/astrocatalogs/supernovae/to the date of April 2022. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2012.1.00650.S, ADS/JAO.ALMA#2013.1.01161.S, ADS/JAO.ALMA#2015.1.00121.S, ADS/JAO.ALMA#2015.1.00925.S, ADS/JAO.ALMA#2015.1.00956.S, ADS/JAO.ALMA#2016.1.00386.S, ADS/JAO.ALMA#2017.1.00392.S, ADS/JAO.ALMA#2017.1.00886.L, ADS/JAO.ALMA#2018.1.01651.S, ADS/JAO.ALMA#2021.1.00099.S Based on data products created from observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme(s) 1100.B-0651, 095.C-0473, and 094.C-0623 (PHANGS–MUSE; PI Schinnerer), as well as 094.B-0321 (MAGNUM; PI Marconi), 099.B-0242, 0100.B-0116, 098.B-0551 (MAD; PI Carollo) and 097.B-0640 (TIMER; PI Gadotti). Source data are provided with this paper. The Supplementary Data 1, 2, and 3 generated in this study have been deposited in the Zenodo database under accession code 10.5281/zenodo.13152457. Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Arcones A Thielemann F-K Origin of the elements Astron. Astrophys. Rev. 2023 31 1 10.1007/s00159-022-00146-x
Arcones, A. & Thielemann, F.-K. Origin of the elements. Astron. Astrophys. Rev. 31, 1 (2023).10.1007/s00159-022-00146-x
2. Van Dyk SD The disappearances of six supernova progenitors Mon. Not. R. Astron. Soc. 2023 519 471 482
Van Dyk, S. D. et al. The disappearances of six supernova progenitors. Mon. Not. R. Astron. Soc. 519, 471–482 (2023).
3. Kilpatrick CD A potential progenitor for the Type Ic supernova 2017ein Mon. Not. R. Astron. Soc. 2018 480 2072 2084 10.1093/mnras/sty2022
Kilpatrick, C. D. et al. A potential progenitor for the Type Ic supernova 2017ein. Mon. Not. R. Astron. Soc. 480, 2072–2084 (2018).10.1093/mnras/sty2022
4. Van Dyk SD SN 2017ein and the possible first identification of a Type Ic supernova progenitor Astrophys. J. 2018 860 90 10.3847/1538-4357/aac32c
Van Dyk, S. D. et al. SN 2017ein and the possible first identification of a Type Ic supernova progenitor. Astrophys. J. 860, 90 (2018).10.3847/1538-4357/aac32c
5. Xiang D Observations of SN 2017ein reveal shock breakout emission and a massive progenitor star for a Type Ic supernova Astrophys. J. 2019 871 176 10.3847/1538-4357/aaf8b0
Xiang, D. et al. Observations of SN 2017ein reveal shock breakout emission and a massive progenitor star for a Type Ic supernova. Astrophys. J. 871, 176 (2019).10.3847/1538-4357/aaf8b0
6. Teffs JJ Prentice SJ Mazzali PA Ashall C Observations and spectral modelling of the narrow-lined Type Ic SN 2017ein Mon. Not. R. Astron. Soc. 2021 502 3829 3842 10.1093/mnras/stab258
Teffs, J. J., Prentice, S. J., Mazzali, P. A. & Ashall, C. Observations and spectral modelling of the narrow-lined Type Ic SN 2017ein. Mon. Not. R. Astron. Soc. 502, 3829–3842 (2021).10.1093/mnras/stab258
7. Begelman MC Sarazin CL SN 1985f: death of a Wolf-Rayet star Astrophys. J. Lett. 1986 302 L59 10.1086/184637
Begelman, M. C. & Sarazin, C. L. SN 1985f: death of a Wolf-Rayet star. Astrophys. J. Lett. 302, L59 (1986).10.1086/184637
8. Woosley SE Langer N Weaver TA The presupernova evolution and explosion of helium stars that experience mass loss Astrophys. J. 1995 448 315 10.1086/175963
Woosley, S. E., Langer, N. & Weaver, T. A. The presupernova evolution and explosion of helium stars that experience mass loss. Astrophys. J. 448, 315 (1995).10.1086/175963
9. Podsiadlowski P Joss PC Hsu JJL Presupernova evolution in massive interacting binaries Astrophys. J. 1992 391 246 10.1086/171341
Podsiadlowski, P., Joss, P. C. & Hsu, J. J. L. Presupernova evolution in massive interacting binaries. Astrophys. J. 391, 246 (1992).10.1086/171341
10. Yoon SC Woosley SE Langer N Type Ib/c supernovae in binary systems. I. Evolution and properties of the progenitor stars Astrophys. J. 2010 725 940 954 10.1088/0004-637X/725/1/940
Yoon, S. C., Woosley, S. E. & Langer, N. Type Ib/c supernovae in binary systems. I. Evolution and properties of the progenitor stars. Astrophys. J. 725, 940–954 (2010).10.1088/0004-637X/725/1/940
11. Zapartas E Predicting the presence of companions for stripped-envelope supernovae: the case of the broad-lined Type Ic SN 2002ap Astrophys. J. 2017 842 125 10.3847/1538-4357/aa7467
Zapartas, E. et al. Predicting the presence of companions for stripped-envelope supernovae: the case of the broad-lined Type Ic SN 2002ap. Astrophys. J. 842, 125 (2017).10.3847/1538-4357/aa7467
12. Smartt SJ Eldridge JJ Crockett RM Maund JR The death of massive stars—I. Observational constraints on the progenitors of Type II-P supernovae Mon. Not. R. Astron. Soc. 2009 395 1409 1437 10.1111/j.1365-2966.2009.14506.x
Smartt, S. J., Eldridge, J. J., Crockett, R. M. & Maund, J. R. The death of massive stars—I. Observational constraints on the progenitors of Type II-P supernovae. Mon. Not. R. Astron. Soc. 395, 1409–1437 (2009).10.1111/j.1365-2966.2009.14506.x
13. Eldridge JJ Fraser M Smartt SJ Maund JR Crockett RM The death of massive stars—II. Observational constraints on the progenitors of Type Ibc supernovae Mon. Not. R. Astron. Soc. 2013 436 774 795 10.1093/mnras/stt1612
Eldridge, J. J., Fraser, M., Smartt, S. J., Maund, J. R. & Crockett, R. M. The death of massive stars—II. Observational constraints on the progenitors of Type Ibc supernovae. Mon. Not. R. Astron. Soc. 436, 774–795 (2013).10.1093/mnras/stt1612
14. Mayker Chen N Comparing the locations of supernovae to CO (2-1) emission in their host galaxies Astrophys. J. 2023 944 110 10.3847/1538-4357/acab00
Mayker Chen, N. et al. Comparing the locations of supernovae to CO (2-1) emission in their host galaxies. Astrophys. J. 944, 110 (2023).10.3847/1538-4357/acab00
15. Leroy AK PHANGS-ALMA data processing and pipeline Astrophys. J. Suppl. Ser. 2021 255 19 10.3847/1538-4365/abec80
Leroy, A. K. et al. PHANGS-ALMA data processing and pipeline. Astrophys. J. Suppl. Ser. 255, 19 (2021).10.3847/1538-4365/abec80
16. Leroy AK PHANGS-ALMA: arcsecond CO(2-1) imaging of nearby star-forming galaxies Astrophys. J. Suppl. Ser. 2021 257 43 10.3847/1538-4365/ac17f3
Leroy, A. K. et al. PHANGS-ALMA: arcsecond CO(2-1) imaging of nearby star-forming galaxies. Astrophys. J. Suppl. Ser. 257, 43 (2021).10.3847/1538-4365/ac17f3
17. Schinnerer E The gas-star formation cycle in nearby star-forming galaxies. I. Assessment of multi-scale variations Astrophys. J. 2019 887 49 10.3847/1538-4357/ab50c2
Schinnerer, E. et al. The gas-star formation cycle in nearby star-forming galaxies. I. Assessment of multi-scale variations. Astrophys. J. 887, 49 (2019).10.3847/1538-4357/ab50c2
18. Chevance M The lifecycle of molecular clouds in nearby star-forming disc galaxies Mon. Not. R. Astron. Soc. 2020 493 2872 2909 10.1093/mnras/stz3525
Chevance, M. et al. The lifecycle of molecular clouds in nearby star-forming disc galaxies. Mon. Not. R. Astron. Soc. 493, 2872–2909 (2020).10.1093/mnras/stz3525
19. Dvoretzky A Kiefer J Wolfowitz J Asymptotic minimax character of the sample distribution function and of the classical multinomial estimator Ann. Math. Stat. 1956 27 642 – 669 10.1214/aoms/1177728174
Dvoretzky, A., Kiefer, J. & Wolfowitz, J. Asymptotic minimax character of the sample distribution function and of the classical multinomial estimator. Ann. Math. Stat. 27, 642 – 669 (1956).10.1214/aoms/1177728174
20. Querejeta M Stellar structures, molecular gas, and star formation across the PHANGS sample of nearby galaxies Astron. Astrophys. 2021 656 A133 10.1051/0004-6361/202140695
Querejeta, M. et al. Stellar structures, molecular gas, and star formation across the PHANGS sample of nearby galaxies. Astron. Astrophys. 656, A133 (2021).10.1051/0004-6361/202140695
21. Zapartas E Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction Astron. Astrophys. 2017 601 A29 10.1051/0004-6361/201629685
Zapartas, E. et al. Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction. Astron. Astrophys. 601, A29 (2017).10.1051/0004-6361/201629685
22. Kim J Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies Mon. Not. R. Astron. Soc. 2022 516 3006 3028 10.1093/mnras/stac2339
Kim, J. et al. Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies. Mon. Not. R. Astron. Soc. 516, 3006–3028 (2022).10.1093/mnras/stac2339
23. Dobbs CL Pringle JE The exciting lives of giant molecular clouds Mon. Not. R. Astron. Soc. 2013 432 653 667 10.1093/mnras/stt508
Dobbs, C. L. & Pringle, J. E. The exciting lives of giant molecular clouds. Mon. Not. R. Astron. Soc. 432, 653–667 (2013).10.1093/mnras/stt508
24. Benincasa SM Live fast, die young: GMC lifetimes in the FIRE cosmological simulations of Milky Way mass galaxies Mon. Not. R. Astron. Soc. 2020 497 3993 3999 10.1093/mnras/staa2116
Benincasa, S. M. et al. Live fast, die young: GMC lifetimes in the FIRE cosmological simulations of Milky Way mass galaxies. Mon. Not. R. Astron. Soc. 497, 3993–3999 (2020).10.1093/mnras/staa2116
25. Lancaster L Ostriker EC Kim J-G Kim C-G Star formation regulation and self-pollution by stellar wind feedback Astrophys. J. Lett. 2021 922 L3 10.3847/2041-8213/ac3333
Lancaster, L., Ostriker, E. C., Kim, J.-G. & Kim, C.-G. Star formation regulation and self-pollution by stellar wind feedback. Astrophys. J. Lett. 922, L3 (2021).10.3847/2041-8213/ac3333
26. Jeffreson SMR A scaling relation for the molecular cloud lifetime in Milky Way-like galaxies Mon. Not. R. Astron. Soc. 2021 505 1678 1698 10.1093/mnras/stab1293 34099958
Jeffreson, S. M. R. et al. A scaling relation for the molecular cloud lifetime in Milky Way-like galaxies. Mon. Not. R. Astron. Soc. 505, 1678–1698 (2021).34099958 10.1093/mnras/stab1293
27. Semenov VA Kravtsov AV Gnedin NY Spatial decorrelation of young stars and dense gas as a probe of the star formation-feedback cycle in galaxies Astrophys. J. 2021 918 13 10.3847/1538-4357/ac0a77
Semenov, V. A., Kravtsov, A. V. & Gnedin, N. Y. Spatial decorrelation of young stars and dense gas as a probe of the star formation-feedback cycle in galaxies. Astrophys. J. 918, 13 (2021).10.3847/1538-4357/ac0a77
28. Davies B Beasor ER The initial masses of the red supergiant progenitors to Type II supernovae Mon. Not. R. Astron. Soc. 2018 474 2116 2128 10.1093/mnras/stx2734
Davies, B. & Beasor, E. R. The initial masses of the red supergiant progenitors to Type II supernovae. Mon. Not. R. Astron. Soc. 474, 2116–2128 (2018).10.1093/mnras/stx2734
29. Xiao L Galbany L Eldridge JJ Stanway ER Core-collapse supernovae ages and metallicities from emission-line diagnostics of nearby stellar populations Mon. Not. R. Astron. Soc. 2019 482 384 401 10.1093/mnras/sty2557
Xiao, L., Galbany, L., Eldridge, J. J. & Stanway, E. R. Core-collapse supernovae ages and metallicities from emission-line diagnostics of nearby stellar populations. Mon. Not. R. Astron. Soc. 482, 384–401 (2019).10.1093/mnras/sty2557
30. Karamehmetoglu E A population of Type Ibc supernovae with massive progenitors. Broad lightcurves not uncommon in (i)PTF Astron. Astrophys. 2023 678 A87 10.1051/0004-6361/202245231
Karamehmetoglu, E. et al. A population of Type Ibc supernovae with massive progenitors. Broad lightcurves not uncommon in (i)PTF. Astron. Astrophys. 678, A87 (2023).10.1051/0004-6361/202245231
31. Renzo M Massive runaway and walkaway stars. A study of the kinematical imprints of the physical processes governing the evolution and explosion of their binary progenitors Astron. Astrophys. 2019 624 A66 10.1051/0004-6361/201833297
Renzo, M. et al. Massive runaway and walkaway stars. A study of the kinematical imprints of the physical processes governing the evolution and explosion of their binary progenitors. Astron. Astrophys. 624, A66 (2019).10.1051/0004-6361/201833297
32. Yoon S-C Evolutionary models for Type Ib/c supernova progenitors Publ. Astron. Soc. Aust. 2015 32 e015 10.1017/pasa.2015.16
Yoon, S.-C. Evolutionary models for Type Ib/c supernova progenitors. Publ. Astron. Soc. Aust. 32, e015 (2015).10.1017/pasa.2015.16
33. Barbarino C Type Ic supernovae from the (intermediate) palomar transient factory Astron. Astrophys. 2021 651 A81 10.1051/0004-6361/202038890
Barbarino, C. et al. Type Ic supernovae from the (intermediate) palomar transient factory. Astron. Astrophys. 651, A81 (2021).10.1051/0004-6361/202038890
34. Chen P A 12.4-day periodicity in a close binary system after a supernova Nature 2024 625 253 258 10.1038/s41586-023-06787-x 38200292
Chen, P. et al. A 12.4-day periodicity in a close binary system after a supernova. Nature 625, 253–258 (2024).38200292 10.1038/s41586-023-06787-x
35. Sana H Binary interaction dominates the evolution of massive stars Science 2012 337 444 10.1126/science.1223344 22837522
Sana, H. et al. Binary interaction dominates the evolution of massive stars. Science 337, 444 (2012).22837522 10.1126/science.1223344
36. Anderson JP A meta-analysis of core-collapse supernova 56Ni masses Astron. Astrophys. 2019 628 A7 10.1051/0004-6361/201935027
Anderson, J. P. A meta-analysis of core-collapse supernova 56Ni masses. Astron. Astrophys. 628, A7 (2019).10.1051/0004-6361/201935027
37. Rodriguez O Nakar E Maoz D Stripped-envelope supernova light curves argue for central engine activity Nature 2024 628 733 735 10.1038/s41586-024-07262-x 38632408
Rodriguez, O., Nakar, E. & Maoz, D. Stripped-envelope supernova light curves argue for central engine activity. Nature 628, 733–735 (2024).38632408 10.1038/s41586-024-07262-x
38. Kelly PL Kirshner RP Pahre M Long γ-Ray bursts and Type Ic core-collapse supernovae have similar locations in hosts Astrophys. J. 2008 687 1201 1207 10.1086/591925
Kelly, P. L., Kirshner, R. P. & Pahre, M. Long γ-Ray bursts and Type Ic core-collapse supernovae have similar locations in hosts. Astrophys. J. 687, 1201–1207 (2008).10.1086/591925
39. Leloudas G Do Wolf-Rayet stars have similar locations in hosts as type Ib/c supernovae and long gamma-ray bursts? Astron. Astrophys. 2010 518 A29 10.1051/0004-6361/200913753
Leloudas, G. et al. Do Wolf-Rayet stars have similar locations in hosts as type Ib/c supernovae and long gamma-ray bursts? Astron. Astrophys. 518, A29 (2010).10.1051/0004-6361/200913753
40. Anderson JP Habergham SM James PA Hamuy M Progenitor mass constraints for core-collapse supernovae from correlations with host galaxy star formation Mon. Not. R. Astron. Soc. 2012 424 1372 1391 10.1111/j.1365-2966.2012.21324.x
Anderson, J. P., Habergham, S. M., James, P. A. & Hamuy, M. Progenitor mass constraints for core-collapse supernovae from correlations with host galaxy star formation. Mon. Not. R. Astron. Soc. 424, 1372–1391 (2012).10.1111/j.1365-2966.2012.21324.x
41. Anderson JP James PA Habergham SM Galbany L Kuncarayakti H Statistical studies of supernova environments Publ. Astron. Soc. Aust. 2015 32 e019 10.1017/pasa.2015.19
Anderson, J. P., James, P. A., Habergham, S. M., Galbany, L. & Kuncarayakti, H. Statistical studies of supernova environments. Publ. Astron. Soc. Aust. 32, e019 (2015).10.1017/pasa.2015.19
42. Kelly PL Kirshner RP Core-collapse supernovae and host galaxy stellar populations Astrophys. J. 2012 759 107 10.1088/0004-637X/759/2/107
Kelly, P. L. & Kirshner, R. P. Core-collapse supernovae and host galaxy stellar populations. Astrophys. J. 759, 107 (2012).10.1088/0004-637X/759/2/107
43. Galbany L Molecular gas in supernova local environments unveiled by EDGE Mon. Not. R. Astron. Soc. 2017 468 628 644 10.1093/mnras/stx367
Galbany, L. et al. Molecular gas in supernova local environments unveiled by EDGE. Mon. Not. R. Astron. Soc. 468, 628–644 (2017).10.1093/mnras/stx367
44. Galbany L PISCO: The PMAS/PPak integral-field supernova hosts compilation Astrophys. J. 2018 855 107 10.3847/1538-4357/aaaf20
Galbany, L. et al. PISCO: The PMAS/PPak integral-field supernova hosts compilation. Astrophys. J. 855, 107 (2018).10.3847/1538-4357/aaaf20
45. Crowther PA On the association between core-collapse supernovae and H ii regions Mon. Not. R. Astron. Soc. 2013 428 1927 1943 10.1093/mnras/sts145
Crowther, P. A. On the association between core-collapse supernovae and H ii regions. Mon. Not. R. Astron. Soc. 428, 1927–1943 (2013).10.1093/mnras/sts145
46. Maund JR The very young resolved stellar populations around stripped-envelope supernovae Mon. Not. R. Astron. Soc. 2018 476 2629 2663 10.1093/mnras/sty093
Maund, J. R. The very young resolved stellar populations around stripped-envelope supernovae. Mon. Not. R. Astron. Soc. 476, 2629–2663 (2018).10.1093/mnras/sty093
47. Sun N-C Maund JR Crowther PA A UV census of the environments of stripped-envelope supernovae Mon. Not. R. Astron. Soc. 2023 521 2860 2873 10.1093/mnras/stad690
Sun, N.-C., Maund, J. R. & Crowther, P. A. A UV census of the environments of stripped-envelope supernovae. Mon. Not. R. Astron. Soc. 521, 2860–2873 (2023).10.1093/mnras/stad690
48. Japelj J Host galaxies of SNe Ic-BL with and without long gamma-ray bursts Astron. Astrophys. 2018 617 A105 10.1051/0004-6361/201833209
Japelj, J. et al. Host galaxies of SNe Ic-BL with and without long gamma-ray bursts. Astron. Astrophys. 617, A105 (2018).10.1051/0004-6361/201833209
49. Prentice SJ The bolometric light curves and physical parameters of stripped-envelope supernovae Mon. Not. R. Astron. Soc. 2016 458 2973 3002 10.1093/mnras/stw299
Prentice, S. J. et al. The bolometric light curves and physical parameters of stripped-envelope supernovae. Mon. Not. R. Astron. Soc. 458, 2973–3002 (2016).10.1093/mnras/stw299
50. Vogelsberger M Marinacci F Torrey P Puchwein E Cosmological simulations of galaxy formation Nat. Rev. Phys. 2020 2 42 66 10.1038/s42254-019-0127-2
Vogelsberger, M., Marinacci, F., Torrey, P. & Puchwein, E. Cosmological simulations of galaxy formation. Nat. Rev. Phys. 2, 42–66 (2020).10.1038/s42254-019-0127-2
51. Chaikin, E. et al. A thermal-kinetic subgrid model for supernova feedback in simulations of galaxy formation. Mon. Not. R. Astron. Soc. 523, 3709 (2023).
52. Chaikin E The importance of the way in which supernova energy is distributed around young stellar populations in simulations of galaxies Mon. Not. R. Astron. Soc. 2022 514 249 264 10.1093/mnras/stac1132
Chaikin, E. et al. The importance of the way in which supernova energy is distributed around young stellar populations in simulations of galaxies. Mon. Not. R. Astron. Soc. 514, 249–264 (2022).10.1093/mnras/stac1132
53. Farmer R Laplace E de Mink SE Justham S The cosmic carbon footprint of massive stars stripped in binary systems Astrophys. J. 2021 923 214 10.3847/1538-4357/ac2f44
Farmer, R., Laplace, E., de Mink, S. E. & Justham, S. The cosmic carbon footprint of massive stars stripped in binary systems. Astrophys. J. 923, 214 (2021).10.3847/1538-4357/ac2f44
54. Guillochon J Parrent J Kelley LZ Margutti R An open catalog for supernova data Astrophys. J. 2017 835 64 10.3847/1538-4357/835/1/64
Guillochon, J., Parrent, J., Kelley, L. Z. & Margutti, R. An open catalog for supernova data. Astrophys. J. 835, 64 (2017).10.3847/1538-4357/835/1/64
55. Sukhbold T Ertl T Woosley SE Brown JM Janka HT Core-collapse supernovae from 9 to 120 solar masses based on neutrino-powered explosions Astrophys. J. 2016 821 38 10.3847/0004-637X/821/1/38
Sukhbold, T., Ertl, T., Woosley, S. E., Brown, J. M. & Janka, H. T. Core-collapse supernovae from 9 to 120 solar masses based on neutrino-powered explosions. Astrophys. J. 821, 38 (2016).10.3847/0004-637X/821/1/38
56. Michałowski MJ Gotkiewicz N Hjorth J Kamphuis P Connection of supernovae 2002ap, 2003gd, 2013ej, and 2019krl in M 74 with atomic gas accretion and spiral structure Astron. Astrophys. 2020 638 A47 10.1051/0004-6361/202037692
Michałowski, M. J., Gotkiewicz, N., Hjorth, J. & Kamphuis, P. Connection of supernovae 2002ap, 2003gd, 2013ej, and 2019krl in M 74 with atomic gas accretion and spiral structure. Astron. Astrophys. 638, A47 (2020).10.1051/0004-6361/202037692
57. Carilli CL Walter F Cool gas in high-redshift galaxies Annu. Rev. Astron. Astrophys. 2013 51 105 161 10.1146/annurev-astro-082812-140953
Carilli, C. L. & Walter, F. Cool gas in high-redshift galaxies. Annu. Rev. Astron. Astrophys. 51, 105–161 (2013).10.1146/annurev-astro-082812-140953
58. Lang P PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution Astrophys. J. 2020 897 122 10.3847/1538-4357/ab9953
Lang, P. et al. PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution. Astrophys. J. 897, 122 (2020).10.3847/1538-4357/ab9953
59. Makarov D Prugniel P Terekhova N Courtois H Vauglin I HyperLEDA. III. The catalogue of extragalactic distances Astron. Astrophys. 2014 570 A13 10.1051/0004-6361/201423496
Makarov, D., Prugniel, P., Terekhova, N., Courtois, H. & Vauglin, I. HyperLEDA. III. The catalogue of extragalactic distances. Astron. Astrophys. 570, A13 (2014).10.1051/0004-6361/201423496
60. Miville-Deschênes M-A Murray N Lee EJ Physical properties of molecular clouds for the entire Milky Way disk Astrophys. J. 2017 834 57 10.3847/1538-4357/834/1/57
Miville-Deschênes, M.-A., Murray, N. & Lee, E. J. Physical properties of molecular clouds for the entire Milky Way disk. Astrophys. J. 834, 57 (2017).10.3847/1538-4357/834/1/57
61. Hinshaw G Nine-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: cosmological parameter results Astrophys. J. Suppl. Ser. 2013 208 19 10.1088/0067-0049/208/2/19
Hinshaw, G. et al. Nine-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: cosmological parameter results. Astrophys. J. Suppl. Ser. 208, 19 (2013).10.1088/0067-0049/208/2/19
62. Grasha K Connecting young star clusters to CO molecular gas in NGC 7793 with ALMA-LEGUS Mon. Not. R. Astron. Soc. 2018 481 1016 1027 10.1093/mnras/sty2154
Grasha, K. et al. Connecting young star clusters to CO molecular gas in NGC 7793 with ALMA-LEGUS. Mon. Not. R. Astron. Soc. 481, 1016–1027 (2018).10.1093/mnras/sty2154
63. Turner JA PHANGS: constraining star formation time-scales using the spatial correlations of star clusters and giant molecular clouds Mon. Not. R. Astron. Soc. 2022 516 4612 4626 10.1093/mnras/stac2559
Turner, J. A. et al. PHANGS: constraining star formation time-scales using the spatial correlations of star clusters and giant molecular clouds. Mon. Not. R. Astron. Soc. 516, 4612–4626 (2022).10.1093/mnras/stac2559
64. Peltonen J Clusters, clouds, and correlations: relating young clusters to giant molecular clouds in M33 and M31 Mon. Not. R. Astron. Soc. 2023 522 6137 6149 10.1093/mnras/stad1430
Peltonen, J. et al. Clusters, clouds, and correlations: relating young clusters to giant molecular clouds in M33 and M31. Mon. Not. R. Astron. Soc. 522, 6137–6149 (2023).10.1093/mnras/stad1430
65. Scheuermann F Stellar associations powering H II regions—I. Defining an evolutionary sequence Mon. Not. R. Astron. Soc. 2023 522 2369 2383 10.1093/mnras/stad878
Scheuermann, F. et al. Stellar associations powering H II regions—I. Defining an evolutionary sequence. Mon. Not. R. Astron. Soc. 522, 2369–2383 (2023).10.1093/mnras/stad878
66. Larson KL Multiscale stellar associations across the star formation hierarchy in PHANGS-HST nearby galaxies: methodology and properties Mon. Not. R. Astron. Soc. 2023 523 6061 6081 10.1093/mnras/stad1600
Larson, K. L. et al. Multiscale stellar associations across the star formation hierarchy in PHANGS-HST nearby galaxies: methodology and properties. Mon. Not. R. Astron. Soc. 523, 6061–6081 (2023).10.1093/mnras/stad1600
67. Emsellem E The PHANGS-MUSE survey. Probing the chemo-dynamical evolution of disc galaxies Astron. Astrophys. 2022 659 A191 10.1051/0004-6361/202141727
Emsellem, E. et al. The PHANGS-MUSE survey. Probing the chemo-dynamical evolution of disc galaxies. Astron. Astrophys. 659, A191 (2022).10.1051/0004-6361/202141727
68. Baade D The Wide Field Imager at the 2.2-m MPG/ESO telescope: first views with a 67-million-facette eye Messenger 1999 95 15 16
Baade, D. et al. The Wide Field Imager at the 2.2-m MPG/ESO telescope: first views with a 67-million-facette eye. Messenger 95, 15–16 (1999).
69. Virtanen P SciPy 1.0: fundamental algorithms for scientific computing in Python Nat. Methods 2020 17 261 272 10.1038/s41592-019-0686-2 32015543
Virtanen, P. et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat. Methods 17, 261–272 (2020).32015543 10.1038/s41592-019-0686-2
70. Sun J Cloud-scale molecular gas properties in 15 nearby galaxies Astrophys. J. 2018 860 172 10.3847/1538-4357/aac326
Sun, J. et al. Cloud-scale molecular gas properties in 15 nearby galaxies. Astrophys. J. 860, 172 (2018).10.3847/1538-4357/aac326
71. Kroupa P On the variation of the initial mass function Mon. Not. R. Astron. Soc. 2001 322 231 246 10.1046/j.1365-8711.2001.04022.x
Kroupa, P. On the variation of the initial mass function. Mon. Not. R. Astron. Soc. 322, 231–246 (2001).10.1046/j.1365-8711.2001.04022.x
