Where are the Population III star relics in the simulated Milky Way analogues?

Basic Information

• Paper ID: 2502.17088
• Title: Where are the Population III star relics in the simulated Milky Way analogues?
• Authors: Hang Yang, Liang Gao, Qi Guo, Haining Li, Shi Shao, Gang Zhao
• Classification: astro-ph.GA (Astrophysics - Galaxies)
• Publication Date: October 17, 2025 (Draft version)
• Paper Link: https://arxiv.org/abs/2502.17088

Abstract

This paper investigates the total mass, spatial distribution, and kinematic properties of Population III star relics at z=0 using six Milky Way analogues from the Auriga simulations with two different numerical resolutions. These relics (primarily second-generation stars) formed over a wide redshift range from z≈22 to z≈4, with an average formation redshift of z∼10.0, comprising approximately 2×10⁻⁵ of the total stellar population in the galaxy. The disk and bulge components contain only a small fraction of these relics, totaling less than 12%. The stellar halo, particularly the outer stellar halo at galactic radius r>30 kpc, contains the largest proportion (averaging ~46%), with approximately one relic star per 4,000 to 10,000 stars, making it a promising region for observational searches.

Research Background and Motivation

Problem Definition

1. Core Question: The spatial distribution and kinematic characteristics of Population III star relics in the present-day Milky Way
2. Scientific Significance: Population III stars mark the end of the cosmic dark ages and initiate chemical enrichment and cosmic reionization processes
3. Observational Challenges: The initial mass function of Population III stars is highly uncertain, and most Population III stars evolve rapidly, ending their lives through supernovae or direct collapse to black holes

Research Motivation

1. Stellar Archaeology Needs: By studying the chemical composition of surviving second-generation stars within the Milky Way, we can infer the properties of Population III stars
2. Limitations of Existing Methods: Most cosmological hydrodynamic simulations cannot explicitly model Population III star formation, typically using extremely metal-poor stars (Fe/H < -3) as proxies
3. Theoretical Prediction Verification: Provides valuable guidance for ongoing and future stellar archaeology surveys

Core Contributions

1. Improved Tracer Method: Employs stellar particles with zero metallicity as a more reasonable tracer for Population III star relics
2. Systematic Spatial Distribution Analysis: First detailed investigation of Population III star relic distribution across different galactic components
3. Kinematic Decomposition: Decomposes relic stars into disk, bulge, inner halo, and outer halo components based on orbital characteristics
4. Satellite Galaxy Analysis: Systematically investigates the distribution patterns of Population III star relics in satellite galaxies
5. Numerical Convergence Verification: Validates result convergence through two different resolutions

Methodology Details

Simulation Setup

Auriga Project:

• Uses magnetohydrodynamic moving mesh code AREPO
• Cosmological parameters: Ωₘ = 0.307, Ωᵦ = 0.048, ΩΗ = 0.693, h = 0.6777
• Two resolutions: L3 (dark matter mass resolution ~5×10⁴M⊙) and L4 (~4×10⁵M⊙)

Population III Star Relic Sample Selection

Key Innovation: Selects stellar particles with metallicity Fe/H = 0 as tracers for Population III star relics (PSR), justified by:

1. These stars formed in pristine, uncontaminated gas
2. Represent the most direct stellar carriers of Population III information in the first galaxies
3. May include both true Population III stars and second-generation stars

Orbital Decomposition Method

Classifies stars into four categories based on galactic radius r and circularity ε = jz/rVc(r):

1. Disk Stars: r < 30 kpc and ε > 0.7
2. Bulge Stars: r < 3 kpc and ε < 0.7
3. Inner Halo Stars: r < 30 kpc and ε < 0.7
4. Outer Halo Stars: r > 30 kpc

Experimental Setup

Dataset

• Six Milky Way Analogues: Au6, Au16, Au21, Au23, Au24, Au27
• Two Resolutions: L3 and L4, totaling 12 simulations
• Galactic Mass Range: log₁₀M₂₀₀ ≈ 12.0-12.24 M⊙
• Stellar Mass Range: log₁₀M* ≈ 10.7-11.0 M⊙

Analysis Tools

• Friends-of-Friends and SUBFIND algorithms for constructing halo and subhalo catalogs
• Galactic coordinate system established based on total angular momentum direction of all stars within 10 kpc at z=0

Experimental Results

Formation Time and Total Mass

1. Formation Redshift Range: z≈22 to z≈4, spanning a wide redshift range
2. Average Formation Redshift: z≈10.0 for Auriga-L3, z≈9.3 for L4
3. Mass Fraction: Approximately 2×10⁻⁵ (L3) to 5×10⁻⁵ (L4) of total galactic stellar population
4. Total Mass: Median approximately 10⁶·²M⊙ (L3) and 10⁶·⁶M⊙ (L4)

Spatial Distribution Characteristics

1. Radial Distribution: Population III star relics are more dispersed than the overall stellar distribution
2. Concentration Comparison: Approximately half of relic stars are located beyond 40 kpc from the galactic center, while half of all stars are distributed within 10 kpc
3. Numerical Convergence: Good qualitative consistency between different resolutions

Kinematic Decomposition Results

Average Results at L3 Resolution:

• Bulge: ~4% of relic stars, one relic per 10⁵·²⁵ bulge stars
• Disk: ~8% of relic stars, one relic per 10⁵·⁵ disk stars
• Inner Halo: ~24% of relic stars, one relic per 10⁴·⁴ inner halo stars
• Outer Halo: ~46% of relic stars, one relic per 10³·⁷ outer halo stars (i.e., one per 4,000-10,000 stars)

Distribution in Satellite Galaxies

1. Overall Proportion: ~18% of Population III star relics are located in satellite galaxies
2. Mass Correlation: Low-mass satellite galaxies contain higher proportions of relic stars
3. Extreme Cases:
   • 16% of satellite galaxies have 100% relic star fraction (experienced only single star formation episode)
   • 16% of satellite galaxies contain no relic stars (formed later, already chemically enriched)

Comparison with Related Work

Comparison with Previous Studies

1. Scannapieco et al. (2006): Using particle tagging methods, predicted second-generation star fraction ~10⁻³·⁵, higher than this study
2. Gao et al. (2010): Semi-analytic modeling predicted ~half of relic stars within 40 kpc, consistent with this study
3. Brook et al. (2007): SPH simulations show metal-poor stars dispersed distribution, qualitatively consistent but larger total mass predictions
4. Starkenburg et al. (2017): APOSTLE simulations show 75% metal-poor stars beyond 8 kpc, slightly lower than this study's 85%

Methodological Advantages

1. More Reasonable Tracer: Uses truly zero-metallicity stars rather than extremely metal-poor stars (Fe/H < -3)
2. Higher Resolution: Provides comparative verification between L3 and L4 resolutions
3. Systematic Analysis: Includes kinematic decomposition and detailed satellite galaxy analysis

Conclusions and Discussion

Main Conclusions

1. Optimal Search Region: The outer stellar halo (r > 30 kpc) is the most promising region for finding Population III star relics
2. Satellite Galaxy Targets: Low-mass, early-formed satellite galaxies are ideal search targets
3. Rarity: Relic stars are extremely rare, with approximately one relic per 4,000-10,000 stars in the outer halo
4. Numerical Reliability: Results show good convergence across different resolutions

Observational Implications

1. Survey Strategy: Should focus on the outer stellar halo and low-mass satellite galaxies
2. Future Missions: Deep survey projects such as the Chinese Space Station Telescope (CSST) will greatly expand sample sizes
3. Statistical Expectations: Provides theoretical foundation for interpreting results from ongoing and future stellar archaeology surveys

Limitations

1. Resolution Constraints: Cannot directly simulate Population III star formation processes
2. Physical Models: Gas physics in first galaxy formation is not sufficiently realistic
3. Chemical Evolution: Lacks detailed modeling of early metal production and mixing processes
4. Sample Definition: Cannot determine precise proportions of true Population III versus second-generation stars

In-Depth Evaluation

Strengths

1. Methodological Innovation: Employs more reasonable tracer method for Population III star relics
2. Strong Systematicity: Comprehensively analyzes spatial distribution, kinematics, and satellite galaxy properties
3. Numerical Verification: Validates reliability through dual-resolution verification
4. Practical Value: Provides clear guidance for observational searches
5. Clear Writing: Well-structured and accurately expressed

Weaknesses

1. Physical Limitations: Constrained by simulation resolution, cannot directly model Population III star formation
2. Sample Purity: Selected zero-metallicity stars may include different stellar types
3. Statistical Sample: Uses only six Milky Way analogues, relatively small sample size
4. Chemical Information: Lacks detailed chemical abundance pattern predictions

Impact

1. Theoretical Contribution: Provides systematic theoretical predictions for spatial distribution of Population III star relics
2. Observational Guidance: Directly guides current and future stellar archaeology observation strategies
3. Methodological Value: Tracer method can be adopted by other similar studies
4. Reproducibility: Based on publicly available Auriga simulation data, results are verifiable

Applicable Scenarios

1. Stellar Archaeology: Provides theoretical guidance for Population III star relic searches in large-scale stellar surveys
2. Galaxy Formation Research: Understands the impact of Population III stars on galaxy evolution in the early universe
3. Observation Strategy Development: Formulates observation plans for next-generation survey telescopes such as CSST
4. Theoretical Model Verification: Provides observational constraints for Population III star theoretical models

References

The paper cites 67 important references, covering multiple aspects including Population III star formation theory, numerical simulation methods, and observational research, reflecting the comprehensiveness and rigor of the research. Key references include the original Auriga project paper (Grand et al. 2017) and Population III star formation theory reviews (Bromm & Yoshida 2011; Klessen & Glover 2023).