Field-Level Comparison and Robustness Analysis of Cosmological N-body Simulations

Abstract

We present the first field-level comparison of cosmological N-body simulations, considering various widely used codes: Abacus, CUBEP${}^3$M, Enzo, Gadget, Gizmo, PKDGrav, and Ramses. Unlike previous comparisons focused on summary statistics, we conduct a comprehensive field-level analysis: evaluating statistical similarity, quantifying implications for cosmological parameter inference, and identifying the regimes in which simulations are consistent. We begin with a traditional comparison using the power spectrum, cross-correlation coefficient, and visual inspection of the matter field. We follow this with a statistical out-of-distribution (OOD) analysis to quantify distributional differences between simulations, revealing insights not captured by the traditional metrics. We then perform field-level simulation-based inference (SBI) using convolutional neural networks (CNNs), training on one simulation and testing on others, including a full hydrodynamic simulation for comparison. We identify several causes of OOD behavior and biased inference, finding that resolution effects, such as those arising from adaptive mesh refinement (AMR), have a significant impact. Models trained on non-AMR simulations fail catastrophically when evaluated on AMR simulations, introducing larger biases than those from hydrodynamic effects. Differences in resolution, even when using the same N-body code, likewise lead to biased inference. We attribute these failures to a CNN's sensitivity to small-scale fluctuations, particularly in voids and filaments, and demonstrate that appropriate smoothing brings the simulations into statistical agreement. Our findings motivate the need for careful data filtering and the use of field-level OOD metrics, such as PQMass, to ensure robust inference.

15 pages, 7 figures, 2 tables

References in corpus (15)

• Simulating galaxy formation with black hole driven thermal and kinetic feedback
• A non-linear solution to the $S_8$ tension?
• The $\text{Abacus}$ Cosmological $N$-body Code
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence
• The CAMELS project: public data release
• How much information can be extracted from galaxy clustering at the field level?
• Consistency tests of field level inference with the EFT likelihood
• LtU-ILI: An All-in-One Framework for Implicit Inference in Astrophysics and Cosmology
• Impacts of the physical data model on the forward inference of initial conditions from biased tracers
• The DESI $N$-body Simulation Project I: Testing the Robustness of Simulations for the DESI Dark Time Survey
• EFTofLSS meets simulation-based inference: ${\sigma }_8$ from biased tracers
• Joint velocity and density reconstruction of the Universe with nonlinear differentiable forward modeling
• A comparative study of cosmological constraints from weak lensing using Convolutional Neural Networks
• Massive $\nu$s through the CNN lens: interpreting the field-level neutrino mass information in weak lensing

Cited by in corpus (1)

• Mitigating Model Misspecification in Simulation-Based Inference for Galaxy Clustering