A novel, fast, and accurate code for cosmological loop calculations
Abstract
We introduce neckar, a public code for the one-loop matter and galaxy power spectrum, as well as correlation function, in the Eulerian Effective Field Theory of Large-Scale Structure (EFTofLSS). It is based on the FFTLog expansion, and is designed for fast and accurate parameter inference from large-scale structure data, making it a powerful tool for high-dimensional Bayesian analyses of current and upcoming galaxy surveys. neckar combines several numerical and analytical improvements, including the separation between pre-computable and cosmology-dependent contributions, a semi-analytic treatment of the Alcock-Paczynski effect for power spectrum multipoles, and an implementation of the symmetry-based bootstrap kernels, allowing straightforward applications to a broader class of cosmological models. The code also incorporates an improved infrared resummation scheme that corrects previous implementation inaccuracies, and which led to errors of the order of 1% for the quadrupole, and of 10% for the hexadecapole at scales k ≳ 0.1 h Mpc−1, as well as an asymptotic treatment that extends the range of accurate predictions into mildly non-linear scales k ∼ 0.4 h Mpc−1. This leads to relevant corrections in the correlation function around the BAO peak. Taken altogether, these improvements render neckar up to two orders of magnitude faster than similar existing codes without relying on emulators, and at the same time improving numerical accuracy. The code is written in C and is distributed with a Python wrapper for convenience.
Related
- Magnetising the quokka code with fast, second-order accurate, error-correcting schemes for magnetohydrodynamics on GPUs
- Fast and accurate differentiable code of the galaxy power spectrum based on Eulerian and Lagrangian one-loop perturbation theories
- FastSpec : Fast and versatile stellar spectrum synthesis code
- Neutrino Direct Simulation Monte Carlo: Accurate modeling of out-of-equilibrium decaying cosmological relics
- Bridging perturbation theory and simulations: initial conditions and fast integrators for cosmological simulations
BibTeX
@misc{harder2026novel,
title = {{A novel, fast, and accurate code for cosmological loop calculations}},
author = {Tjark Harder and Luca Amendola and Thiago S. Pereira and Miguel Quartin},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.03466},
eprint = {2610.03466},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.03466},
}
RIS
TY - UNPB
TI - A novel, fast, and accurate code for cosmological loop calculations
AU - Tjark Harder
AU - Luca Amendola
AU - Thiago S. Pereira
AU - Miguel Quartin
PY - 2026
DA - 2026/10/02
PB - arXiv
DO - 10.48550/arxiv.2610.03466
UR - https://arxiv.org/abs/2610.03466
AB - We introduce neckar, a public code for the one-loop matter and galaxy power spectrum, as well as correlation function, in the Eulerian Effective Field Theory of Large-Scale Structure (EFTofLSS). It is based on the FFTLog expansion, and is designed for fast and accurate parameter inference from large-scale structure data, making it a powerful tool for high-dimensional Bayesian analyses of current and upcoming galaxy surveys. neckar combines several numerical and analytical improvements, including the separation between pre-computable and cosmology-dependent contributions, a semi-analytic treatment of the Alcock-Paczynski effect for power spectrum multipoles, and an implementation of the symmetry-based bootstrap kernels, allowing straightforward applications to a broader class of cosmological models. The code also incorporates an improved infrared resummation scheme that corrects previous implementation inaccuracies, and which led to errors of the order of 1% for the quadrupole, and of 10% for the hexadecapole at scales k ≳ 0.1 h Mpc^−1, as well as an asymptotic treatment that extends the range of accurate predictions into mildly non-linear scales k ∼ 0.4 h Mpc^−1. This leads to relevant corrections in the correlation function around the BAO peak. Taken altogether, these improvements render neckar up to two orders of magnitude faster than similar existing codes without relying on emulators, and at the same time improving numerical accuracy. The code is written in C and is distributed with a Python wrapper for convenience.
ER -