🔭 astrophysics

The DORCHA suite: nature, nurture, and the phase space distribution of the Milky Way's high redshift progenitors today

Using the DORCHA suite of high-resolution cosmological simulations, this study demonstrates that the remnants of the Milky Way's earliest high-redshift progenitors are robustly concentrated in the galaxy's innermost regions with distinctive kinematic signatures, guiding future searches for ancient stellar populations like Population III stars toward the inner bulge and halo.

Original authors: Sreedhar Balu, Chris Power, Kris Walker, J. Stuart B. Wyithe

Published 2026-03-19

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sreedhar Balu, Chris Power, Kris Walker, J. Stuart B. Wyithe

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Cosmic Detective Story: Finding the Milky Way's Oldest Ancestors

Imagine the Milky Way not just as a galaxy, but as a giant, ancient family tree. For decades, astronomers have been trying to find the "great-great-grandparents" of our galaxy—the very first stars that formed billions of years ago. The big question is: Where in our galaxy today should we look to find the fossils of these ancient stars?

This paper, titled The Dorcha suite, acts like a cosmic time machine and a detective's map to answer that question.

1. The Setup: Building 25 Digital Universes

To solve this mystery, the researchers didn't just look at the real sky; they built 25 digital universes in a computer. They called this collection the "Dorcha suite" (pronounced DUR-uh-khuh, which means "Dark" in Gaelic).

The Analogy: Think of these simulations as 25 different "what-if" scenarios. In each one, they created a galaxy that looks and acts just like our Milky Way.
The Variety: Some of these digital galaxies were lonely (isolated), while others had a massive neighbor nearby (like our Milky Way and the Andromeda galaxy). Some formed very early in the universe's history, and others formed a bit later. This allowed the scientists to see if the galaxy's history or its neighbors changed where the oldest stars ended up.

2. The Method: Tagging the "VIPs"

In the real universe, we can't see the very first stars directly because they are too faint or have already died. So, the scientists used a clever trick in their simulations.

The Analogy: Imagine a massive, chaotic dance party (the early universe) where thousands of people are moving around. The scientists wanted to track the "VIPs"—the people who were in the very center of the most crowded, energetic rooms at the start of the party.
The Technique: They went back in time in their simulation to when the universe was very young (high redshift, $z \ge 5$ ). They identified the densest, most tightly packed clumps of matter (the "VIPs") that were likely to form the first stars. They put a digital tag on these specific particles.
The Journey: Then, they let the simulation run forward for 13 billion years, watching how the universe evolved. They watched to see where these "tagged" particles ended up in the present day.

3. The Discovery: The "Inner Sanctum"

The results were surprisingly consistent, no matter which of the 25 digital galaxies they looked at.

The Finding: The ancient, tagged material didn't scatter randomly. Instead, it sank deep into the center of the galaxy, like a heavy stone sinking to the bottom of a pond.
The Numbers: About 90% to 100% of these ancient remnants are now living within the innermost 15,000 light-years of the galaxy's center.
The Behavior: These ancient particles are "calm." They aren't zooming around wildly; they are moving in smooth, predictable, radial paths (like spokes on a wheel) and are tightly packed together.

4. The "Nature vs. Nurture" Conclusion

The paper asked a crucial question: Does it matter if a galaxy grew up in a quiet neighborhood or a chaotic one with a giant neighbor?

The Answer: No. Whether the digital galaxy was lonely or had a massive neighbor, the ancient stars still ended up in the same place: the deep center.
The Metaphor: It's like dropping a handful of marbles into a bowl. Whether you drop them gently (a quiet galaxy) or shake the bowl violently (a galaxy with a massive neighbor), the heavy marbles (the ancient stars) always end up at the very bottom of the bowl. The "nurture" (environment) didn't change the "nature" (where the old stars live).

5. What This Means for Real Astronomers

This is a game-changer for real-world astronomy.

The Old Mistake: For a long time, astronomers thought the oldest, most metal-poor stars might be scattered all over the galaxy, or that they were too hard to find in the crowded, dusty center of the Milky Way. They often looked at the outer edges.
The New Map: This paper says, "Stop looking at the edges! Look at the center."
- If you want to find the fossilized remains of the first stars (Population III stars), you need to look in the innermost bulge of the Milky Way.
- These stars are likely hiding there, moving in a specific, calm way that distinguishes them from the younger, chaotic stars around them.

6. A Word of Caution

The authors add a small but important warning: Age and "metal-poor" aren't always the same thing.

The Analogy: Think of "metal-poor" like a vintage car that has been repainted. Just because a car looks old (low metallicity) doesn't mean it's from the 1920s; it could be a newer car painted to look old. Conversely, a very old car might have been repainted many times.
The Takeaway: While the oldest stars are likely in the center, simply finding a star with low metals doesn't guarantee it's ancient. Astronomers need to look at the star's movement (kinematics) as well as its chemical makeup to be sure.

Summary

The Dorcha suite simulations tell us that the Milky Way's earliest ancestors didn't wander off to the suburbs. They stayed right in the city center, settling into a calm, dense neighborhood deep in the galactic core. If we want to find the fossils of the universe's first stars, we need to turn our telescopes toward the heart of our galaxy, not the outskirts.

1. Problem Statement

The primary scientific question addressed is: Where in the present-day Milky Way (MW) should astronomers search for the remnants of its earliest stars?

While Galactic Archaeology surveys (e.g., Gaia, APOGEE, GALAH) aim to reconstruct the assembly history of the MW by studying old, metal-poor stars, there is significant uncertainty regarding their spatial distribution. Observational campaigns have historically overlooked the inner regions of the MW due to dust obscuration and the dominance of metal-rich populations. Furthermore, theoretical models suggest that the oldest stars (formed at $z \gtrsim 5$ ) may constitute a negligible fraction of the inner stellar population, making them difficult to distinguish from younger, metal-poor stars using metallicity alone. The authors seek to determine the robustness of the dynamical signature of these high-redshift progenitors against variations in assembly history (e.g., isolated vs. paired environments) and formation time.

2. Methodology

The authors utilize the Dorcha suite (Gaelic for "Dark"), a set of 25 high-resolution, dark-matter-only (DMO) cosmological zoom-in simulations of Milky Way Analogue (MWA) haloes evolved to $z=0$ .

Simulation Setup:
- Parent Simulation: A $100 \, h^{-1}\text{Mpc}$ periodic volume ( $N=512^3$ particles) using the AREPO code with Planck 2016 cosmological parameters.
- Zoom Simulations: 25 MWA haloes were selected based on mass ( $1.04 \times 10^{12} \le M_h \le 1.95 \times 10^{12} \, h^{-1}M_\odot$ $1.04 \times 1 0^{12} \leq M_{h} \leq 1.95 \times 1 0^{12} h^{- 1} M_{⊙}$ ) and environment.
  - 15 Isolated Systems: No massive neighbours within $1.4 \, h^{-1}\text{Mpc}$ .
  - 5 Paired Systems: Analogous to the MW-M31 pair, with a massive neighbour within $1 \, h^{-1}\text{Mpc}$ .
- Resolution: The zoom regions achieve a particle mass of $1.57 \times 10^5 \, h^{-1}M_\odot$ and a softening length of $50 \, h^{-1}\text{pc}$ .
- Snapshots: 123 snapshots from $z=30$ to $z=0$ , with high cadence at high redshift ( $z>5$ ) to track merger trees accurately.
Tagging Strategy (Tracing High- $z$ Progenitors):
- The authors identify progenitor haloes at redshifts $z \in [25, 20, 15, 10, 8, 5]$ .
- Tagging Criterion: The most bound 10% of particles in each progenitor halo are tagged. This selection targets the deepest gravitational potential wells, which are the most likely sites for early star formation (atomic hydrogen cooling).
- Classification: Particles are categorized as "old" (in place at $z \ge 10$ ) and "new" (accreted at $z \le 5$ ). Haloes are also split into "early-forming" ( $z_{form} \gtrsim 1.7$ ) and "late-forming" ( $z_{form} \lesssim 1.0$ ) based on when 50% of their $z=0$ mass was assembled.
Analysis Tools:
- Spherically averaged radial profiles (density, velocity, anisotropy).
- Phase-space analysis ( $r$ vs. $v_r$ ).
- Orbital tracking to determine minimum pericentric distances ( $r_{peri}$ ).

3. Key Contributions

The Dorcha Suite: Provides a statistically significant sample (25 MWAs) with diverse formation histories and environments (isolated vs. paired) to test the universality of dynamical signatures.
Particle Tagging Protocol: Establishes a robust method for tracing the "fossil record" of high-redshift star formation sites using the most-bound dark matter particles as proxies, independent of specific baryonic feedback models.
Decoupling Metallicity from Age: Explicitly demonstrates that while metallicity is a useful proxy, it is not a unique identifier for ancient stars; instead, kinematic and spatial signatures are more robust discriminators.

4. Key Results

A. Spatial Distribution

Central Concentration: The material associated with high-redshift progenitors ( $z \ge 5$ $z \geq 5$ ) is highly centrally concentrated at $z=0$ $z = 0$ .
- 90–100% of the material tagged at $z \ge 15$ resides within $r \lesssim 15 \, h^{-1}\text{kpc}$ (approx. $0.1 R_{200}$ ).
- The density profile of this ancient component is extremely steep, following $\rho \propto r^{-4}$ (Hernquist-like), significantly steeper than the overall halo profile.
Environment Independence: This central concentration holds true for both isolated and paired systems, as well as early- and late-forming haloes. The presence of a massive neighbour (like Andromeda) does not significantly alter the inner phase-space distribution of the oldest components.

B. Kinematic Properties

Dynamical Equilibrium: The median radial velocity $\bar{v}_r \approx 0$ within the inner regions, indicating the material is in dynamical equilibrium and has been embedded in the potential for many dynamical times.
Low Velocity Dispersion: The ancient component is "dynamically cold" compared to the bulk halo, with low radial velocity dispersion ( $\sigma_r / \sigma_{max} \lesssim 0.6$ ).
Radial Orbits: The velocity anisotropy parameter $\beta$ is high ( $\beta \gtrsim 0.5$ for $r \gtrsim 0.1 R_{200}$ ), indicating that these particles are on radially biased orbits.
Pericentric Passage: Orbital tracking reveals that the minimum pericentric distances for these particles are extremely small ( $r_{peri} \sim 10^{-3} - 10^{-2} R_{200}$ ), confirming they have sunk deep into the central potential well.

C. Phase Space Structure

The phase-space distribution ( $r$ vs. $v_r$ ) shows that the high- $z$ progenitor material is confined to the innermost potential, forming a relaxed, cold population.
While caustics and shells exist (reflecting the accretion history), the bulk of the ancient material is well-mixed in the center, distinct from the outer halo which retains more coherent substructure.

5. Significance and Implications

Observational Strategy: The results strongly suggest that searches for the first generations of stars (Population III and extremely metal-poor stars) should focus on the innermost regions of the Milky Way (the bulge and inner halo), specifically within $r \lesssim 15 \, \text{kpc}$ .
Kinematic Selection: Since metallicity alone is insufficient to identify ancient stars (due to overlap with metal-poor stars formed later), surveys should prioritize kinematic selection: looking for stars with low velocity dispersion and radially biased orbits in the inner Galaxy.
Robustness: The dynamical signature of early progenitors is robust to the assembly history of the host halo. Whether the MW formed early or late, or exists in isolation or a pair, the oldest material ends up in the same dynamically cold, central reservoir.
Future Work: The Dorcha suite serves as a baseline for the upcoming Solas suite (hydrodynamical simulations), which will incorporate baryonic physics to refine these predictions regarding star formation efficiency and feedback effects.

In conclusion, the paper establishes that the "fossil record" of the Milky Way's high-redshift assembly is preserved in a centrally concentrated, dynamically cold, and radially biased population in the inner halo, providing a clear target for future Galactic Archaeology surveys.