Galactic chemical evolution with the short-lived isotopes Mn-53, Fe-60, Hf-182, and Pu-244

Basic Information

• Paper ID: 2510.13070
• Title: Galactic chemical evolution with the short-lived isotopes Mn-53, Fe-60, Hf-182, and Pu-244
• Authors: B. Wehmeyer, A. Yagüe López, B. Côté, M. K. Pető, C. Kobayashi, M. Lugaro
• Classification: astro-ph.GA (Galactic Astrophysics)
• Conference Publication: EPJ Web of Conferences 297, 01016 (2024) - ISNA23 Conference
• DOI: https://doi.org/10.1051/epjconf/202429701016

Abstract

This study implements a three-dimensional galactic chemical evolution (GCE) model to track the propagation of four short-lived radioactive isotopes (SLRs) throughout the Milky Way: ⁵³Mn from Type Ia supernovae, ⁶⁰Fe from core-collapse supernovae, ¹⁸²Hf from intermediate-mass stars, and ²⁴⁴Pu from neutron star mergers. The research compares the galactic chemical evolution of these SLRs with recent detections in deep-sea sediments, finding that despite originating from different astrophysical sites, these SLRs frequently arrive at Earth simultaneously.

Research Background and Motivation

Core Research Questions

The fundamental problem addressed by this study is understanding the propagation mechanisms of short-lived radioactive isotopes throughout the Milky Way and how they converge at Earth from different astrophysical sources.

Research Significance

1. Temporal Information Value: Short-lived isotopes (with half-lives on the Myr timescale) provide additional temporal information about astrophysical processes
2. Nucleosynthesis Constraints: Utilizing the exponential nature of radioactive decay, unknown parameters can be inferred from three parameters: production yield, observed quantity, and time interval
3. Multi-Source Convergence Phenomenon: Explaining why isotopes from different sources appear simultaneously at Earth

Observational Methodology Innovation

Traditional SLR observation methods include:

• Measuring isotopic excesses in meteorites to study the early solar system
• Detecting current distributions and production sources through γ-ray observations

This study introduces a third method: analyzing SLR deposition records in deep-sea sediments.

Core Contributions

1. Established the first three-dimensional GCE model: Simultaneously tracking four short-lived isotopes from different sources propagating through the Milky Way
2. Discovered the "shock wave surfing" mechanism: Explaining how SLRs from different sources propagate cooperatively through supernova shock waves
3. Predicted the existence of ¹⁸²Hf: Based on model predictions, the fourth SLR should be detectable in deep-sea sediments
4. Validated deep-sea detection data: Model results are consistent with observational data for ⁵³Mn, ⁶⁰Fe, and ²⁴⁴Pu

Methodology Details

Task Definition

Input: Initial conditions of the Milky Way, star formation rate, nucleosynthesis yield data Output: Spatiotemporal density distribution evolution of four SLRs Constraints: Radioactive decay law, hydrodynamic equations, stellar evolution timescales

Model Architecture

Three-Dimensional GCE Model Parameters

• Spatial Scale: Cubic domain with 2 kpc edge length, periodic boundary conditions
• Spatial Resolution: (50 pc)² subgrid
• Time Step: 1 Myr
• Time Window: 13.309-13.332 Gy (approximately 23 Myr)

Core Computational Workflow

Each time step includes the following calculations:

1. Gas Inflow: Following a prescription of early rise and late exponential decline
2. Star Formation: Using Schmidt law (power-law index α=1.5) and Salpeter initial mass function (slope -2.35, mass range 0.1-50 M⊙)
3. Stellar Lifetime: Adopting the Geneva group equation:

   log(t) = (3.79 + 0.24Z) − (3.10 + 0.35Z)log(M) + (0.74 + 0.11Z)log²(M)
   

4. Binary System Probability:
   • SNIa probability: P_SNIa = 6×10⁻³
   • NSM probability: P_NSM = 0.04

Treatment Schemes for Different Stellar Objects

Intermediate-Mass Stars (IMSs)

• Evolutionary Endpoint: Planetary nebulae
• ¹⁸²Hf Yield: Based on yield tables from Lugaro et al.
• Mass Return: Outer material returned to local subgrid

Core-Collapse Supernovae (CCSNe)

• Explosion Energy: ~1 Bethe
• Swept-Up Mass: 5×10⁴ M⊙
• ⁶⁰Fe Yield: Using yield data from Limongi & Chieffi
• Shock Wave Model: Sedov-Taylor theory

Type Ia Supernovae

• ⁵³Mn Yield: 10⁻⁴ M⊙
• Swept-Up Mass: 5×10⁴ M⊙
• Yield Model: CDD2 model

Neutron Star Mergers (NSMs)

• Merger Timescale: t_coal = 10⁸ years
• ²⁴⁴Pu Yield: 10⁻⁸ M⊙
• Total Ejected Mass: 10⁻² M⊙
• Mass Fraction: X_244 = 10⁻⁶

Technical Innovations

1. Multi-Source Isotope Coupling: First simultaneous treatment of four SLRs from different sources in a single model
2. Shock Wave Propagation Mechanism: Detailed modeling of gas sweeping and isotope transport from explosive events
3. Spatiotemporal Correlation Analysis: Revealing cooperative propagation patterns of isotopes from different sources

Experimental Setup

Data Sources

• Deep-Sea Detection Data: From observations by Wallner et al. 11,12,14 and Korschinek et al.
• Nucleosynthesis Yields: Latest data from multiple theoretical calculation groups
• Stellar Evolution Models: Stellar evolution tracks from the Geneva group

Model Validation Strategy

• Shape-Fitting Priority: Focusing on the shape of SLR density variations rather than absolute values
• Temporal Axis Adjustment: Time-shifting deep-sea data to optimize fitting
• Statistical Analysis: Computing 68%, 95%, and 100% statistical fluctuation intervals

Experimental Results

Main Findings

SLR Density Evolution Characteristics

Three primary density variation patterns are observed from Figure 1:

1. Radioactive Decay: At t≈13.330 Gy, ⁶⁰Fe and ¹⁸²Hf show diagonal downward decay trends
2. Explosive Clearing: Sharp density drops (spanning several orders of magnitude) caused by nearby explosive events clearing gas and SLR content
3. Cooperative Enhancement: At t=13.314 Gy, ⁵³Mn, ¹⁸²Hf, and ²⁴⁴Pu simultaneously exhibit jump-like growth

"Shock Wave Surfing" Mechanism

The study reveals the key mechanism for SLR propagation:

• Initial Enrichment: NSM explosion produces ²⁴⁴Pu-enriched shell
• Secondary Acceleration: Subsequent CCSN explosion within the enriched shell produces ⁶⁰Fe
• Cooperative Transport: CCSN kinetic energy propels the shell containing both SLRs outward
• Simultaneous Arrival: Enabling SLRs from different sources to reach observation locations simultaneously

Statistical Results

• Median Density: Yellow line shows median density of each SLR across all subgrids
• Statistical Fluctuations: Black (68%), dark gray (95%), and light gray (100%) regions show density statistical distribution
• Best Fit: Green line represents the subgrid evolution with best match to deep-sea detection data

Model Predictions

Based on model results, the study predicts that ¹⁸²Hf should be detectable in deep-sea sediments, providing theoretical guidance for future observations.

Related Work

GCE Research History

• Early Work: Audouze & Tinsley (1976) established the foundational framework for GCE
• Three-Dimensional Models: Previous work by Wehmeyer et al. 15,16 laid the model foundation for this study
• Nucleosynthesis Research: Recent advances by Kobayashi et al. 5 and Prantzos et al. 6

SLR Observation Techniques

• Meteorite Studies: Traditional method for early solar system SLR research
• γ-Ray Astronomy: Contemporary SLR distribution observations by Diehl et al. 8-10
• Deep-Sea Archaeology: Deep-sea sediment SLR detection pioneered by Wallner et al.

Conclusions and Discussion

Main Conclusions

1. Multi-Source Convergence Mechanism: Confirms that SLRs from different astrophysical sources can arrive at Earth simultaneously through shock wave propagation mechanisms
2. Propagation Timescale: SLR propagation timescales in the interstellar medium are comparable to their half-lives, enabling cooperative transport
3. Observational Verification: The three-dimensional GCE model successfully reproduces the SLR deposition patterns detected in deep-sea sediments

Limitations

1. Spatial Scale Constraints: The 2 kpc simulation volume may not fully represent the complex structure of the Milky Way
2. Temporal Resolution: The 1 Myr time step may miss certain rapid evolutionary processes
3. Simplified Assumptions:
   • Probabilistic treatment of binary evolution neglects detailed stellar evolution
   • Fixed merger timescale (10⁸ years) may not apply to all systems
   • Non-realistic nature of periodic boundary conditions

Future Directions

1. Larger Spatial Scales: Extending simulations to galactic-scale dimensions
2. Finer Temporal Resolution: Improving temporal resolution to capture rapid processes
3. Detailed Binary Evolution: Integrating complete binary evolution models
4. Multi-Element Extension: Including additional types of short-lived isotopes

In-Depth Evaluation

Strengths

1. Innovative Methodology: First three-dimensional propagation model for multi-source SLRs, providing new perspectives for understanding galactic chemical evolution
2. Mechanism Discovery: Reveals the important "shock wave surfing" isotope propagation mechanism, explaining the multi-source convergence phenomenon
3. Theory-Observation Integration: Combining theoretical models with latest deep-sea detection data, enhancing practical significance
4. Predictive Value: Provides theoretical expectations for future ¹⁸²Hf detection

Limitations

1. Parameter Dependence: The model contains multiple empirical parameters (such as binary probability, yields, etc.), whose uncertainties may affect results
2. Spatial Constraints: The 2 kpc simulation range is relatively small compared to galactic scales, potentially missing large-scale effects
3. Temporal Fitting: Requires time-shifting observational data to achieve good fits, indicating model uncertainty in absolute timescales

Impact

1. Disciplinary Contribution: Provides new research tools for galactic chemical evolution and nuclear astrophysics
2. Observational Guidance: Provides theoretical framework for future deep-sea SLR detection experiments
3. Method Generalization: Three-dimensional GCE modeling methodology can be applied to studies of other nuclides

Applicable Scenarios

1. Galactic Chemical Evolution Research: Applicable to studying distribution and evolution of various nuclides in the Milky Way
2. Nuclear Astrophysics Verification: Can be used to verify nucleosynthesis theories from different astrophysical processes
3. Earth Science Applications: Provides theoretical basis for explaining cosmogenic nuclides in Earth's geological records

References

The paper cites 35 important references covering various aspects of galactic chemical evolution, stellar evolution, nucleosynthesis theory, and observational astrophysics. Particularly noteworthy are the series of works by Wallner et al. on deep-sea SLR detection 11-14 and classical research on stellar evolution by the Geneva group 22-24.

Overall Assessment: This is a high-quality astrophysics research paper with important contributions in both methodological innovation and scientific discovery. The "shock wave surfing" mechanism revealed through three-dimensional numerical simulations provides a new physical picture for understanding short-lived isotope propagation in the Milky Way, possessing significant scientific value and application prospects.