COmoving Computer Acceleration (COCA): $N$-body simulations in an emulated frame of reference
Abstract: $N$-body simulations are computationally expensive, so machine-learning (ML)-based emulation techniques have emerged as a way to increase their speed. Although fast, surrogate models have limited trustworthiness due to potentially substantial emulation errors that current approaches cannot correct for. To alleviate this problem, we introduce COmoving Computer Acceleration (COCA), a hybrid framework interfacing ML with an $N$-body simulator. The correct physical equations of motion are solved in an emulated frame of reference, so that any emulation error is corrected by design. This approach corresponds to solving for the perturbation of particle trajectories around the machine-learnt solution, which is computationally cheaper than obtaining the full solution, yet is guaranteed to converge to the truth as one increases the number of force evaluations. Although applicable to any ML algorithm and $N$-body simulator, this approach is assessed in the particular case of particle-mesh cosmological simulations in a frame of reference predicted by a convolutional neural network, where the time dependence is encoded as an additional input parameter to the network. COCA efficiently reduces emulation errors in particle trajectories, requiring far fewer force evaluations than running the corresponding simulation without ML. We obtain accurate final density and velocity fields for a reduced computational budget. We demonstrate that this method shows robustness when applied to examples outside the range of the training data. When compared to the direct emulation of the Lagrangian displacement field using the same training resources, COCA's ability to correct emulation errors results in more accurate predictions. COCA makes $N$-body simulations cheaper by skipping unnecessary force evaluations, while still solving the correct equations of motion and correcting for emulation errors made by ML.
- (Buchert, 1989) T. Buchert, A class of solutions in Newtonian cosmology and the pancake theory, Astron. & Astrophys. 223, 9 (1989).
- (Leclercq, 2015) F. Leclercq, Bayesian large-scale structure inference and cosmic web analysis, Ph.D. thesis, Institut d’Astrophysique de Paris (2015), arXiv:1512.04985 [astro-ph.CO] .
- (Leclercq, 2018) F. Leclercq, Bayesian optimization for likelihood-free cosmological inference, Phys. Rev. D 98, 063511 (2018), arXiv:1805.07152 .
- (Leclercq, 2022) F. Leclercq, Simulation-based inference of Bayesian hierarchical models while checking for model misspecification, Physical Sciences Forum 5, arXiv:2209.11057 (2022), arXiv:2209.11057 [stat.ME] .
- (Planck Collaboration, 2020) Planck Collaboration, Planck 2018 results. VI. Cosmological parameters, Astron. & Astrophys. 641, A6 (2020), arXiv:1807.06209 [astro-ph.CO] .
- (Zel’dovich, 1970) Y. B. Zel’dovich, Gravitational instability: An approximate theory for large density perturbations., Astron. & Astrophys. 5, 84 (1970).
Paper Prompts
Sign up for free to create and run prompts on this paper using GPT-5.
Top Community Prompts
Collections
Sign up for free to add this paper to one or more collections.