The abacus cosmological N-body code

Lehman H. Garrison, Daniel J. Eisenstein, Douglas Ferrer, Nina A. Maksimova, Philip A. Pinto

Research output: Contribution to journalArticlepeer-review

7 Scopus citations


We present abacus, a fast and accurate cosmological N-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near-and far-field forces, reducing the former to direct 1/r2 summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of 10-5. We express the simulation time-step as an event-driven 'pipeline', incorporating asynchronous events such as completion of co-processor work, input/output, and network communication. abacus has been used to produce the largest suite of N-body simulations to date, the abacussummit suite of 60 trillion particles, incorporating on-The-fly halo finding. abacus enables the production of mock catalogues of the volume and resolution required by the coming generation of cosmological surveys.

Original languageEnglish (US)
Pages (from-to)575-596
Number of pages22
JournalMonthly Notices of the Royal Astronomical Society
Issue number1
StatePublished - Nov 1 2021


  • cosmology: Theory
  • methods: numerical

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science


Dive into the research topics of 'The abacus cosmological N-body code'. Together they form a unique fingerprint.

Cite this