The Progenitor of the S147 Supernova Remnant

By analyzing Gaia DR3 data to model the local stellar population around supernova remnant S147, the study identifies a high-mass progenitor range of $21.5M_{\odot}$to$41.1M_{\odot}$ for the supernova, based on the ages of the most luminous nearby stars including the unbound binary companion HD 37424.

The Gist

Imagine the night sky as a giant, cosmic crime scene. A star exploded long ago, leaving behind a glowing, expanding cloud of debris called a Supernova Remnant (SNR). In the center of this cloud, S147, sits a rapidly spinning "dead" star called a pulsar (PSR J0538+2817).

The big mystery astronomers have been trying to solve is: What kind of star exploded? Was it a lightweight star that barely made the cut, or a massive, heavyweight giant?

This paper is like a team of cosmic detectives using a new set of high-tech magnifying glasses (data from the Gaia satellite) to look at the neighborhood where the explosion happened. They aren't just looking at the crime scene; they are looking at the survivors and the neighbors to figure out who the victim was.

Here is the story of their investigation, broken down simply:

1. The "Runaway" Suspect

The first clue was a star named HD 37424. It's a bright, blue star that is currently flying away from the center of the explosion at high speed.

• The Analogy: Imagine a car crash where one car is totaled (the supernova) and the other car (HD 37424) is thrown off the road, speeding away.
• The Theory: Scientists believe these two stars were originally a couple, orbiting each other. When the heavier partner exploded, the explosion was so violent it blew the couple apart, sending the survivor flying. This makes HD 37424 the "unbound companion" of the star that died.

2. The Neighborhood Watch

To find out how big the dead star was, the astronomers didn't just look at the survivor. They looked at the entire neighborhood (a cylinder of space about 360 light-years long and 200 light-years wide) surrounding the explosion.

They used data from the Gaia satellite (a space telescope that maps the positions and distances of billions of stars) to pick out 439 stars that live in this specific neighborhood.

• The Analogy: Think of it like walking into a high school cafeteria to figure out how old the graduating class was. You don't just look at the valedictorian; you look at the whole room. If the room is full of 18-year-olds, you know the class is graduating now. If it's full of 10-year-olds, you know the explosion happened much later.

3. The "Time Machine" (Color and Brightness)

Stars change color and brightness as they age, just like people change from energetic toddlers to gray-haired elders.

• The Method: The team plotted these 439 stars on a graph (a Color-Magnitude Diagram). It's like a "family tree" where the position of a star tells you its age and mass.
• The Twist: They had to account for the fact that some stars are actually binary systems (two stars orbiting each other that look like one bright dot). They ran two different computer simulations: one assuming all stars are single, and one assuming they are pairs.

4. The Verdict: A Heavyweight Champion

After crunching the numbers, the results pointed to a very specific conclusion:

• The Survivor's Weight: The runaway star (HD 37424) weighs about 13.5 times the mass of our Sun.
• The Victim's Weight: For the victim to have exploded before its partner (HD 37424), the victim had to be heavier than the partner.
• The Statistical Clue: The "neighborhood" analysis showed that the most likely age for the explosion corresponds to a star that was 21.5 to 41.1 times the mass of our Sun.

In simple terms: The explosion wasn't caused by a "medium-sized" star. It was caused by a cosmic heavyweight. The star that exploded was likely more than twice as heavy as the star that survived.

5. Why This Matters

This is a big deal for two reasons:

1. It's a Rare Find: This is one of the very few places in our galaxy where we can clearly see a star that was "kicked out" of a binary system by a supernova. It's like finding a smoking gun.
2. The Method Works: The astronomers used a technique usually reserved for distant galaxies (where we can't see individual stars easily) and applied it right here in our own backyard. They proved that by studying the "crowd" of stars around a supernova, we can statistically deduce the size of the star that died, even if we can't see the dead star itself.

The Bottom Line

The paper concludes that the star that created the S147 supernova remnant was a massive, short-lived giant (21–41 times the mass of the Sun). It lived fast, died young, and blew its partner (HD 37424) out of the neighborhood, leaving behind a spinning neutron star and a cloud of debris that we can still see today.

It's a story of a cosmic breakup that ended in a spectacular explosion, solved by looking at the stars that were left behind.

Technical Summary

Here is a detailed technical summary of the paper "The Progenitor of the S147 Supernova Remnant" by Cruz-Cruz & Kochanek (2025).

1. Problem Statement

Understanding the progenitor masses of core-collapse supernovae (ccSN) is critical for constraining massive star evolution and supernova mechanisms. While methods exist to estimate progenitor masses in external galaxies by modeling the local stellar population's Color-Magnitude Diagram (CMD), applying this to Galactic supernova remnants (SNRs) has been historically difficult due to:

• Distance Uncertainties: Galactic SNRs often lack precise distance measurements.
• Projection Effects: Difficulty distinguishing stars physically associated with the SNR from foreground/background contamination.
• Binary Complexity: Massive stars frequently exist in binaries; interactions can alter evolutionary paths and explosion outcomes, complicating single-star models.

The S147 SNR is a unique case study because it contains the pulsar PSR J0538+2817 and a likely unbound binary companion, HD 37424. It is the only strong Galactic candidate for a binary system unbound by a core-collapse supernova. However, the mass of the progenitor star that created S147 remains debated, with previous estimates ranging widely (e.g., $20-35 M_\odot$).

2. Methodology

The authors employ a statistical stellar population analysis using Gaia Data Release 3 (DR3) to reconstruct the star formation history (SFH) of the region surrounding S147.

A. Stellar Sample Selection

• Volume Definition: A truncated cylinder centered on S147 (RA: $84.75^\circ$, DEC:$27.83^\circ$) with a projected radius of 100 pc.
• Distance Constraints: Based on recent distance measurements (Kochanek et al. 2024) of $1.37$kpc and the parallax of HD 37424, the authors selected stars with parallaxes$0.614 < \varpi < 0.787$mas. This corresponds to a line-of-sight depth of$\approx 360$ pc.
• Magnitude Cuts: Stars with $G < 15$ (absolute magnitude $M_G \le 4.3$) and extinction-corrected colors $-1.0 < (B_P - R_P) < 3.5$ were selected.
• Final Sample: 439 stars with $-8 < M_G < 0$.

B. Individual Star Modeling

• The 12 most luminous stars ($M_G < -3$) were individually modeled using Spectral Energy Distributions (SEDs).
• Tools: SED fitting was performed using DUSTY within a Markov Chain Monte Carlo (MCMC) framework, utilizing MARCS and Castelli & Kurucz model atmospheres.
• Data: Multi-wavelength fluxes from UV (Thompson et al., Wesselius et al.) to mid-IR (2MASS, ALLWISE) were used.
• Binary Systems: Specific attention was paid to known binaries (e.g., HD 37366, ET Tau, HD 37424) to determine primary masses and ages.

C. Population Density Modeling
The authors constructed two distinct stellar density models to fit the observed CMD:

1. Single-Star Model: Stars drawn from a Salpeter Initial Mass Function (IMF, $x=2.35$) between $1-100 M_\odot$.
2. Non-Interacting Binary Model: Primary masses drawn from the IMF; secondary masses drawn uniformly from the same isochrone ($1 M_\odot \le M_2 \le M_1$). Fluxes were summed to simulate binary CMD positions.

• Isochrones: Solar metallicity PARSEC isochrones were used, spanning ages from $10^{6.3}$to$10^{10.1}$ years in 0.3 dex bins.
• Likelihood Analysis: An MCMC approach (emcee) was used to optimize the star formation rate (SFR) in each age bin. The likelihood function was based on a Poisson probability of star counts in CMD bins, renormalized to match the observed total star count ($N_\star = 439$).
• Progenitor Probability: The probability of the SNR progenitor originating from a specific age bin was calculated based on the fraction of massive stars ($M > M_{min}$) that would have died within the last $\delta t = 10^5$ years.

3. Key Contributions

• Application of CMD Analysis to Galactic SNR: Successfully applied the external galaxy SFH method to a Galactic SNR, leveraging precise Gaia DR3 parallaxes and a newly derived distance to S147.
• Refined Progenitor Mass Constraints: Provided a statistically robust mass range for the S147 progenitor, moving beyond simple upper/lower limits derived from kinematics.
• Binary Population Synthesis: Demonstrated the impact of including non-interacting binary populations in CMD modeling, showing that binaries increase the number of young, luminous stars in the model.
• Validation of the Unbound Companion Scenario: Confirmed that the age of the unbound companion (HD 37424) is consistent with the derived age of the progenitor, supporting the hypothesis that S147 resulted from a binary system disrupted by the supernova.

4. Key Results

A. Individual Star Properties

• HD 37424 (Unbound Companion): Estimated mass $\approx 13.5 M_\odot$, age $\approx 10^{6.8}$ yr.
• HD 37367: Second most luminous single star, mass $\approx 14.3 M_\odot$.
• HD 37366: Most luminous star (O9.5 IV spectroscopic binary), primary mass $\approx 20.9 M_\odot$.
• ET Tau: Eclipsing binary (B8), primary mass $\approx 16.7 M_\odot$.

B. Progenitor Mass Distribution
The analysis of the 439-star population yielded the following probability distributions for the progenitor mass:

• Preferred Mass Range: Both single-star and binary models strongly favor a high-mass progenitor in the range $21.5 M_\odot - 41.1 M_\odot$.
• Probability Peaks:
  • In the Single-Star Model, the age bin $10^{6.6} - 10^{6.9}$yr (corresponding to$21.5-41.1 M_\odot$) has a differential probability roughly **4 times higher** than bins for lower masses ($\le 21.5 M_\odot$) or higher masses ($\ge 41.1 M_\odot$).
  • In the Binary Model, the same mass range remains the most probable, though the peak is slightly less pronounced (approx. 1.4 times higher than the highest mass bin, but still $\sim 3$ times higher than lower mass bins).
• Integral Probability: The age range $10^{6.6} - 10^{6.9}$yr encompasses$\sim 83%$of the probability for the single-star model and$\sim 64%$ for the binary model.
• Consistency: The derived progenitor mass ($> 21.5 M_\odot$) is consistent with the requirement that the progenitor must have been more massive than its companion (HD 37424, $\sim 13.5 M_\odot$) to explode first and unbind the system.

C. Star Formation History

• The local stellar population shows evidence of a recent burst of star formation (within the last $10^7$ years), which is necessary to produce the observed population of massive O and B stars.

5. Significance and Implications

• Progenitor Mass Constraints: The study provides one of the most precise mass estimates for a Galactic ccSN progenitor, narrowing the range to $21.5-41.1 M_\odot$. This supports theories that high-mass stars are the primary drivers of core-collapse events in such environments.
• Binary Evolution: The results validate the "unbound binary" scenario for S147. The progenitor was likely a massive star in a binary system with HD 37424; the supernova explosion disrupted the orbit, leaving HD 37424 as a runaway star and the pulsar as the compact remnant.
• Methodological Limitations: The authors caution that while powerful, this method relies heavily on the assumption that the local stellar population is dominated by a single starburst. Projection effects and contamination can obscure the true SFH. They suggest that for Galactic SNRs, a comparative statistical approach (comparing SNR regions to control fields without SNRs) may be more robust for ensemble studies.
• Future Directions: The success of this method on S147, aided by precise distance measurements via high-velocity absorption lines, paves the way for a statistical survey of Galactic SNRs to better understand the demographics of supernova progenitors and binary interactions.