Advanced Megatron cosmological radiation-hydrodynamics simulations offer a physical bridge between ancient stars and modern observations, explaining why the universe’s smallest dwarf galaxies maintain a surprisingly uniform iron-to-hydrogen ratio despite containing vastly different total numbers of stars, according to research published in The Open Journal of Astrophysics 2026.
Astronomers Combine Telescope Observations with Stellar Archaeology
Astronomers investigating the early universe typically divide their efforts into two distinct approaches. The James Webb Space Telescope captures distant galactic light emitted when the cosmos was young. Meanwhile, stellar archaeology examines nearby ancient stars preserving chemical clues from their primordial gas clouds.
Dr Martin Rey and other researchers have approached the early universe through two directions. The James Webb Space Telescope observes distant galaxies whose light began traveling when the cosmos was young, while astronomers also study how elements and conditions allow later generations of stars to emerge.
Simulating Gas, Radiation, and Chemical Networks
Megatron models track gravity, flowing gas, traveling radiation, and shifting chemistry simultaneously. The simulation’s chemical network tracks more than 80 distinct atoms, ions, and molecules. This granular approach maps precisely where primordial gas cools, gathers, and achieves the density necessary to spawn new stars.
For the iron-abundance investigation, researchers simulated a region destined to evolve into a Milky-Way-like galaxy by redshift zero. Four distinct simulation runs began with identical initial conditions while varying assumptions about later stellar formation and feedback mechanisms. The calculations resolved physical structures down to a scale of just a few parsecs.
These high-resolution runs cover an interval within the first billion years of cosmic history. Rather than tracing every individual gas cloud to the present day, the team relied on auxiliary calculations to estimate which early dwarf galaxies likely survived as satellites around a larger host galaxy.
Explaining the Faint-Galaxy Iron Plateau
Larger galaxies generally exhibit higher iron-to-hydrogen ratios compared to smaller systems. Among the faintest dwarf galaxies, however, that correlation flattens. Many display remarkably similar average iron abundances despite holding declining stellar populations.
All four Megatron simulation runs successfully reproduced this low-mass plateau. The flattening appeared in dwarf galaxies containing no more than roughly 100,000 solar masses of stars. Their characteristic iron-to-hydrogen ratio measured approximately 300 times lower than that of the Sun.
The simulation data traced this chemical signature back to pair-instability supernovae. These exceptionally energetic blasts utterly destroy massive first-generation stars, scattering freshly synthesized iron into adjacent gas clouds. A single explosion from a sufficiently massive primordial star supplies the necessary iron to establish the observed plateau.
Ultraviolet starlight regulated whether that iron remained accessible for subsequent star formation. Specific energetic photons dissociated molecular hydrogen, which served as a primary coolant for early gas clouds.
Radiation thus dictates the final iron abundance by modulating the timeline and spatial distribution of early star formation.
Simulations Show Gaps in Carbon and Nitrogen Abundances
Despite matching the iron plateau, the simulations exposed notable limitations in broader chemical tracking. Only three of the four model variations successfully replicated the general mass-metallicity relation across larger dwarf galaxies. One variant retained excessive iron within massive dwarfs due to insufficiently powerful galactic outflows.
Additional chemical discrepancies emerged. Predicted carbon and magnesium abundances diverged from empirical observations, while nitrogen production fell significantly short of measured levels. Matching the iron record alone does not confirm that the simulations capture the complete chemical pedigree of early galaxies.
All four simulation runs utilized identical assumptions regarding the primordial stellar population. Consequently, the models test variations in subsequent stellar physics rather than resolving the uncertain masses and explosion yields of the absolute first generation. Different initial star populations or radiation environments could alter the projected outcomes.
The simulations also yield a distinct, testable forecast: a smaller population of ultra-iron-poor galaxies. Roughly one-fifth of the modeled faint dwarfs exhibited average iron-to-hydrogen ratios at least 1,000 times below the solar baseline. Empirical verification remains difficult due to sparse chemical measurements in extremely faint systems, leaving the existence of this predicted iron-poor population as an open question for future astronomical censuses.
