COmoving Computer Acceleration (COCA): N-body simulations in an emulated frame of reference
Astronomy & Astrophysics EDP Sciences 694 (2025) ARTN A287
Abstract:
<jats:p><jats:italic>Context.N</jats:italic>-body simulations are computationally expensive and machine learning (ML) based emulation techniques have thus emerged as a way to increase their speed. Surrogate models are indeed fast, however, they are limited in terms of their trustworthiness due to potentially substantial emulation errors that current approaches are not equipped to correct.</jats:p> <jats:p><jats:italic>Aims.</jats:italic> To alleviate this problem, we have introduced COmoving Computer Acceleration (COCA), a hybrid framework interfacing ML algorithm with an <jats:italic>N</jats:italic>-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. Thus, we are able to find a solution for the perturbation of particle trajectories around the ML solution. This approach is computationally cheaper than obtaining the full solution and it is guaranteed to converge to the truth as the number of force evaluations is increased.</jats:p> <jats:p><jats:italic>Methods.</jats:italic> Even though it is applicable to any ML algorithm and <jats:italic>N</jats:italic>-body simulator, we assessed this approach in the particular case of particle-mesh (PM) cosmological simulations in a frame of reference predicted by a convolutional neural network. In such cases, the time dependence is encoded as an additional input parameter to the network.</jats:p> <jats:p><jats:italic>Results.</jats:italic> We find that COCA efficiently reduces emulation errors in particle trajectories, requiring far fewer force evaluations than running the corresponding simulation without ML. As a consequence, we were able to obtain accurate final density and velocity fields for a reduced computational budget. We demonstrate that this method exhibits 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.</jats:p> <jats:p><jats:italic>Conclusions.</jats:italic> Therefore, COCA makes <jats:italic>N</jats:italic>-body simulations cheaper by skipping unnecessary force evaluations, while still solving the correct equations of motion and correcting for emulation errors made by ML.</jats:p>The Velocity Field Olympics: Assessing velocity field reconstructions with direct distance tracers
(2025)
Bye-bye, Local-in-matter-density Bias: The Statistics of the Halo Field Are Poorly Determined by the Local Mass Density
The Astrophysical Journal Letters American Astronomical Society 977:2 (2024) ARTN L44
Abstract:
<jats:title>Abstract</jats:title> <jats:p>Bias models relating the dark matter field to the spatial distribution of halos are widely used in current cosmological analyses. Many models predict halos purely from the local Eulerian matter density, yet bias models in perturbation theory require other local properties. We assess the validity of assuming that only the local dark matter density can be used to predict the number density of halos in a model-independent way and in the nonperturbative regime. Utilizing <jats:italic>N</jats:italic>-body simulations, we study the properties of the halo counts field after spatial voxels with near-equal dark matter density have been permuted. If local-in-matter-density (LIMD) biasing were valid, the statistical properties of the permuted and unpermuted fields would be indistinguishable since both represent equally fair draws of the stochastic biasing model. If the Lagrangian radius is greater than approximately half the voxel size and for halos less massive than ∼10<jats:sup>15</jats:sup> <jats:italic>h</jats:italic> <jats:sup>−1</jats:sup> <jats:italic>M</jats:italic> <jats:sub>☉</jats:sub>, we find the permuted halo field has a scale-dependent bias with greater than 25% more power on scales relevant for current surveys. These bias models remove small-scale power by not modeling correlations between neighboring voxels, which substantially boosts large-scale power to conserve the field’s total variance. This conclusion is robust to the choice of initial conditions and cosmology. Assuming LIMD halo biasing cannot, therefore, reproduce the distribution of halos across a large range of scales and halo masses, no matter how complex the model. One must either allow the biasing to be a function of other quantities and/or remove the assumption that neighboring voxels are statistically independent.</jats:p>Scant evidence for thawing quintessence
Physical Review D American Physical Society (APS) 110:8 (2024) 83528