A Hybrid Fast Multipole Method for Cosmological N-body Simulations

Abstract

We investigate a hybrid numerical algorithm aimed at the large-scale cosmological N-body simulation for the on-going and the future high precious sky surveys. It makes use of a truncated Fast Multiple Method (FMM) for short-range gravity, incorporating with a Particle Mesh (PM) method for long-range potential, which is applied to deal with extremely large particle number. In this work, we present a specific strategy to modify a conventional FMM by a Gaussian shaped factor and provide quantitative expressions for the interaction kernels between multipole expansions. Moreover, a proper multipole acceptance criteria for the hybrid method is introduced to solve potential precision loss induced by the truncation. Such procedures reduce the mount of computation than an original FMM and decouple the global communication. A simplified version of code is introduced to verify the hybrid algorithm, accuracy and parallel implementation.

13 pages, 8 figures, Accepted for publication in Research in Astronomy and Astrophysics

Cited by in corpus (5)

• Simulating cosmic structure formation with the GADGET-4 code
• Large-scale dark matter simulations
• The Ultramarine Simulation: properties of dark matter haloes before redshift 5.5
• On the statistical theory of self-gravitating collisionless dark matter flow
• PhotoNs-GPU:A GPU accelerated cosmological simulation code