Looking for cracks in the standard cosmological model
New computer simulations follow the formation of galaxies and the cosmic large-scale structure with unprecedented statistical precision
TORONTO, July 19, 2023 – York University and an international team of astrophysicists have made an ambitious attempt to simulate the formation of galaxies and cosmic large-scale structure throughout staggeringly large swaths of space. First results of their “MillenniumTNG” project are published in a series of 10 articles in the journal Monthly Notices of the Royal Astronomical Society. The new calculations help to subject the standard cosmological model to precision tests and to unravel the full power of upcoming new cosmological observations, say the researchers including York Assistant Professor Rahul Kannan.
Over the past decades, cosmologists have gotten used to the perplexing conjecture that the universe’s matter content is dominated by enigmatic dark matter and that an even stranger dark energy field that acts as some kind of anti-gravity to accelerate the expansion of today’s cosmos. Ordinary baryonic matter makes up less than five per cent of the cosmic mix, but this source material forms the basis for the stars and planets of galaxies like our own Milky Way.
This seemingly strange cosmological model is known under the name LCDM. It provides a stubbornly successful description of a large number of observational data, ranging from the cosmic microwave radiation – the rest-heat left behind by the hot Big Bang – to the “cosmic web”, where galaxies are arranged along an intricate network of dark matter filaments. However, the real physical nature of dark matter and dark energy is still not understood, prompting astrophysicists to search for cracks in the LCDM theory. Identifying tensions to observational data could lead to a better understanding of these fundamental puzzles about our Universe. Sensitive tests are required that need both: powerful new observational data as well as more detailed predictions about what the LCDM model actually implies.
An international team of researchers led by the Max Planck Institute for Astrophysics (MPA) in Germany, Harvard University in the US, Durham University in the UK, and the Donostia International Physics Center in Spain, along with York University, have now managed to take a decisive step forward on the latter challenge. Building up on their previous successes with the “Millennium” and “IllustrisTNG” projects, they developed a new suite of simulation models dubbed “MillenniumTNG”, which trace the physics of cosmic structure formation with considerably higher statistical accuracy than what was possible with previous calculations.
Large simulations including new physical details
The team utilized the advanced cosmological code GADGET-4, custom-built for this purpose, to compute the largest high-resolution dark matter simulations to date, covering a region nearly 10 billion light-years across. In addition, they employed the moving-mesh hydrodynamical code AREPO to follow the processes of galaxy formation directly, throughout volumes still so large that they can be considered representative for the universe as a whole. Comparing both types of simulations allows a precise assessment of the impact of baryonic processes related to supernova explosions and supermassive black holes on the total matter distribution. An accurate knowledge of this distribution is key for interpreting upcoming observations correctly, such as so-called weak gravitational lensing effects, which respond to matter irrespective of whether it is of dark or baryonic type.
Furthermore, the team included massive neutrinos in their simulations, for the first time in simulations big enough to allow meaningful cosmological mock observations. Previous cosmological simulations had usually omitted them for simplicity, because they make up at most one to two per cent of the dark matter mass, and since their nearly relativistic velocities mostly prevent them from clumping together. Now, however, upcoming cosmological surveys (such as those of the recently launched Euclid satellite of the European Space Agency) will reach a precision allowing a detection of the associated percent-level effects. This raises the tantalizing prospect to constrain the neutrino mass itself, a profound open question in particle physics, so the stakes are high.
For their ground-breaking MillenniumTNG simulations, the researchers made efficient use of two extremely powerful supercomputers, the SuperMUC-NG machine at the Leibniz Supercomputing Center in Garching, and the Cosma8 machine at Durham Universe. More than 120 000 compute cores toiled away for nearly two months at SuperMUC-NG, using computing time awarded by the German Gauss Centre for Supercomputing, to produce the most comprehensive hydrodynamical simulation model to date. MillenniumTNG is tracking the formation of about one hundred million galaxies in a region of the universe around 2400 million light-years across (see Figure 1). This calculation is about 15 times bigger than the previously best is this category, the TNG300 model of the IllustrisTNG project.
Using Cosma8, the team computed an even bigger volume of the universe, filled with more than a trillion dark matter particles and more than 10 billion particles for tracking massive neutrinos (see Figure 2). Even though this simulation did not follow the baryonic matter directly, its galaxy content can be accurately predicted in MillenniumTNG with a semi-analytic model that is calibrated against the full physical calculation of the project. This procedure leads to a detailed distribution of galaxies and matter in a volume that for the first time is large enough to be representative for the universe as a whole, putting comparisons to upcoming observational surveys on a sound statistical basis.
Theoretical predictions for cosmology
The first results of the MillenniumTNG project show a wealth of new theoretical predictions that reinforce the importance of computer simulations in modern cosmology. The team has written and submitted ten introductory scientific papers for the project. Eight of them have just appeared simultaneously in the journal MNRAS, the remaining two are about to follow shortly.
One timely study examines the discovery of a population of very massive galaxies in the young universe with the James Webb Space Telescope. The masses of these galaxies are unexpectedly large just a brief time after the Big Bang, seemingly defying theoretical expectations. Dr. Kannan analyzed the predictions of MillenniumTNG for this early epoch. While the simulations agree with the observations out to redshifts of z=10 (when the universe was less than 500 million years old), he confirmed that, if they hold up, the new results by JWST at even higher redshift conflicts with the simulation predictions.
“Perhaps star formation is much more efficient shortly after the Big Bang than at later times, or maybe massive stars are formed in higher proportions back then, making these galaxies unusually bright”, says Kannan of York University’s Faculty of Science.
Another study looked at the shapes of galaxies. Nearby galaxies have the subtle tendency to orient their shapes in similar directions instead of pointing randomly, an effect called “intrinsic galaxy alignments”. This poorly understood effect distorts inferences based on weak gravitational lensing, which creates its own statistical alignment signal. The MillenniumTNG project could for the first-time measure intrinsic alignments with very high signal-to-noise directly from the shapes of the simulated galaxies, out to distances of several hundred million light-years. “Perhaps our determination of the intrinsic alignment of galaxy orientations can help to resolve the current discrepancy between the amplitude of matter clustering inferred from weak lensing and from the cosmic microwave background,” says PhD-student Ana Maria Delgado of Harvard University, first author of this study of the MillenniumTNG team. Using these results, astronomers will be able to correct for this important systematic effect much better.
Other works of the team’s initial analysis focus on the clustering signals of galaxies. For example, MPA PhD student Monica Barrera produced extremely large and highly realistic mock catalogues of galaxies on the past backwards “lightcone” of a fiducial observer (see Figure 3). In this case, galaxies that are more distant are also automatically younger, reflecting the travel time of the light that is reaching our telescopes. Using these virtual observations, she looked at the so-called baryonic acoustic oscillation (BAO) feature (which provides a cosmologically important standard ruler) in the projected two-point correlation function of galaxies. Her results showed, that measuring these BAOs is a fairly tricky endeavour that can be significantly influenced by so-called cosmic variance effects – even when extremely large volumes are studied in observational surveys. While in simulations one can observe the modelled universe from different vantage points to recover the correct statistical ensemble average, this is unfortunately not readily possible for the real Universe. “The MillenniumTNG simulations are so big and contain so many galaxies, more than 1 billion in the biggest calculation, that it was really hard to study them”, says Monica Barrera. “Analysis scripts that work just fine for smaller simulations tend to take forever for MillenniumTNG.”
Analyzing cosmological data
The flurry of first results from the MillenniumTNG simulations make it clear that they will be of great help to design better strategies for the analysis of upcoming cosmological data. The team’s principal investigator, Professor Volker Springel from MPA argues that “MillenniumTNG combines recent advances in simulating galaxy formation with the field of cosmic large-scale structure, allowing an improved theoretical modelling of the connection of galaxies to the dark matter backbone of the Universe. This may well prove instrumental for progress on key questions in cosmology, such as how the mass of neutrinos can be best constrained with large-scale structure data.” The MillenniumTNG simulations produced more than three Petabytes of simulation data, forming a rich asset for further research that will keep the participating scientists busy for many years to come.
Original scientific publications:
- The MillenniumTNG Project: The galaxy population at z ≥ 8
R. Kannan, V. Springel, L. Hernquist, R. Pakmor, A. M. Delgado, B. Hadzhiyska, C. Hernández-Aguayo, M. Barrera, F. Ferlito, S. Bose, S. D. M. White, C. Frenk, A. Smith, E. Garaldi
MNRAS, July 2023 (preprint: https://arxiv.org/abs/2210.10066)
- The MillenniumTNG Project: High-precision predictions for matter clustering and halo statistics
C. Hernández-Aguayo, V. Springel, R. Pakmor, M. Barrera, F. Ferlito, S. D. M. White, L. Hernquist, B. Hadzhiyska, A. M. Delgado, R. Kannan, S. Bose, C. Frenk
MNRAS, July 2023 (preprint: https://arxiv.org/abs/2210.10059)
- The MillenniumTNG Project: The hydrodynamical full physics simulation and a first look at its galaxy clusters
R. Pakmor, V. Springel, J. P. Coles, T. Guillet, C. Pfrommer, S. Bose, M. Barrera, A. M. Delgado, F. Ferlito, C. Frenk, B. Hadzhiyska, C. Hernández-Aguayo, L. Hernquist, R. Kannan, S. D. M. White
MNRAS, July 2023 (preprint: https://arxiv.org/abs/2210.10060)
- The MillenniumTNG Project: Semi-analytic galaxy formation models on the past lightcone
M. Barrera, V. Springel, S. White, C. Hernández-Aguayo, L. Hernquist, C. Frenk, R. Pakmor, F. Ferlito, B. Hadzhiyska, A. M. Delgado, R. Kannan, S. Bose
MNRAS, submitted (preprint: https://arxiv.org/abs/2210.10419)
- The MillenniumTNG Project: Refining the one-halo model of red and blue galaxies at different redshifts
B. Hadzhiyska, L. Hernquist, D. Eisenstein, A. M. Delgado, S. Bose, R. Kannan, R. Pakmor, V. Springel, S. Contreras, M. Barrera, F. Ferlito, C. Hernández-Aguayo, S. D. M. White, C. Frenk
MNRAS, July 2023 (preprint: https://arxiv.org/abs/2210.10068)
- The MillenniumTNG Project: An improved two-halo model for the galaxy-halo connection of red and blue galaxies
B. Hadzhiyska, D. Eisenstein, L. Hernquist, R. Pakmor, S. Bose, A. M. Delgado, S. Contreras, R. Kannan, S. D. M. White, V. Springel, C. Frenk, C. Hernández-Aguayo, F. Ferlito, M. Barrera
MNRAS, July 2023 (preprint: https://arxiv.org/abs/2210.10072)
- The MillenniumTNG Project: The large-scale clustering of galaxies
S. Bose, B. Hadzhiyska, M. Barrera, A. M. Delgado, F. Ferlito, C. Frenk, C. Hernández-Aguayo, L. Hernquist, R. Kannan, R. Pakmor, V. Springel, S. D. M. White
MNRAS, July 2023 (preprint: https://arxiv.org/abs/2210.10065)
- The MillenniumTNG Project: Inferring cosmology from galaxy clustering with accelerated N-body scaling and subhalo abundance matching
S. Contreras, R. E. Angulo, V. Springel, S. D. M. White, B. Hadzhiyska, L. Hernquist, R. Pakmor, R. Kannan, C. Hernández-Aguayo, M. Barrera, F. Ferlito, A. M. Delgado, S. Bose, C. Frenk
MNRAS, July 2023 (preprint: https://arxiv.org/abs/2210.10075)
- The MillenniumTNG Project: Intrinsic alignments of galaxies and halos
A. M. Delgado, B. Hadzhiyska, S. Bose, V. Springel, L. Hernquist, M. Barrera, R. Pakmor, F. Ferlito, R. Kannan, C. Hernández-Aguayo, S. D. M. White, C. Frenk
MNRAS, July 2023 (preprint: https://arxiv.org/abs/2304.12346)
- The MillenniumTNG Project: The impact of baryons and massive neutrinos on high-resolution weak gravitational lensing convergence maps
F. Ferlito, V. Springel, C. T. Davies, C. Hernández-Aguayo, R. Pakmor, M. Barrera, S. D. M. White, A. M. Delgado, B. Hadzhiyska, L. Hernquist, R. Kannan, S. Bose, C. Frenk
MNRAS, submitted (preprint: https://arxiv.org/abs/2304.12338)
Web site of the MillenniumTNG project
Gauss Centre for Supercomputing
SuperMUC-NG at the Leibniz Supercomputing Centre
York University is a modern, multi-campus, urban university located in Toronto, Ontario. Backed by a diverse group of students, faculty, staff, alumni and partners, we bring a uniquely global perspective to help solve societal challenges, drive positive change and prepare our students for success. York's fully bilingual Glendon Campus is home to Southern Ontario's Centre of Excellence for French Language and Bilingual Postsecondary Education. York’s campuses in Costa Rica and India offer students exceptional transnational learning opportunities and innovative programs. Together, we can make things right for our communities, our planet, and our future.
York University: Sandra McLean, York University Media Relations, 416-272-6317, firstname.lastname@example.org
MPA: Hannelore Hämmerle, +49-89-30000-3980, email@example.com
Harvard & Smithsonian Center for Astrophysics: +1- 617-721-7371, pr@, https://cfa.harvard.edu
Prof. Dr. Volker Springel, Max-Planck Institute for Astrophysics (MPA), +49-89-30000-2195, firstname.lastname@example.org
Figures / Captions:
Projections of gas (top left), dark matter (top right), and stellar light (bottom center) for a slice in the largest hydrodynamical simulation of MillenniumTNG at the present epoch. The slice is about 35 million light-years thick. The projections show the vast physical scales in the simulation from size, about 2400 million light-years across, to an individual spiral galaxy (final round inset) with a radius of ~150 000 light-years. The underlying calculation is presently the largest high-resolution hydrodynamical simulation of galaxy formation, containing more than 160 billion resolution elements. © MPA
Comparison of the neutrino (top) and dark matter (bottom) distributions on the past backwards lightcone of a fiducial observer positioned at the centre of the two horizontal stripes. As cosmic expansion slows down the neutrinos at late times (small redshift/distance), they start to weakly cluster around the biggest concentrations of dark matter as shown by a comparison of the zoomed insets. This slightly increases the mass and further growth rate of these largest structures. © MPA
Galaxy distribution on the past backwards lightcone in MillenniumTNG, where the galaxies are predicted with a sophisticated semi-analytic model on top of the dark matter backbone. Galaxies are shown down to Johnson apparent magnitude 𝑅 < 23, in a 180 degrees wide, thin wedge with opening angle 0.24 degrees, out to redshift 𝑧 = 2. The galaxy positions are drawn as circles with comoving coordinates in real space, using red for galaxies with rest frame color index 𝐵−𝑅 > 0.7, and blue otherwise. Real observations of the galaxy positions would additionally be perturbed by small shifts along the line of sight due to the Doppler effects from the galaxies’ motions, an effect that can also be easily included in the models. The two circular insets show nested zooms with diameters of around 1.25 billion light-years and 125 million light-years, and fainter apparent magnitude limits of 𝑅 < 25 and 𝑅 < 28, respectively. © MPA