MSTAR - a fast parallelized algorithmically regularized integrator with minimum spanning tree coordinates

Show full item record



Permalink

http://hdl.handle.net/10138/320660

Citation

Rantala , A , Pihajoki , P , Mannerkoski , M , Johansson , P H & Naab , T 2020 , ' MSTAR - a fast parallelized algorithmically regularized integrator with minimum spanning tree coordinates ' , Monthly Notices of the Royal Astronomical Society , vol. 492 , no. 3 , pp. 4131-4148 . https://doi.org/10.1093/mnras/staa084

Title: MSTAR - a fast parallelized algorithmically regularized integrator with minimum spanning tree coordinates
Author: Rantala, Antti; Pihajoki, Pauli; Mannerkoski, Matias; Johansson, Peter H.; Naab, Thorsten
Contributor organization: Particle Physics and Astrophysics
Department of Physics
Date: 2020-03
Language: eng
Number of pages: 18
Belongs to series: Monthly Notices of the Royal Astronomical Society
ISSN: 0035-8711
DOI: https://doi.org/10.1093/mnras/staa084
URI: http://hdl.handle.net/10138/320660
Abstract: We present the novel algorithmically regularized integration method MSTAR for high-accuracy (vertical bar Delta E/E vertical bar greater than or similar to 10(-14)) integrations of N-body systems using minimum spanning tree coordinates. The twofold parallelization of the O(N-part(2)) force loops and the substep divisions of the extrapolation method allow for a parallel scaling up to N-CPU = 0.2 x N-part. The efficient parallel scaling of MSTAR makes the accurate integration of much larger particle numbers possible compared to the traditional algorithmic regularization chain (AR-CHAIN) methods, e.g. N-part = 5000 particles on 400 CPUs for 1 Gyr in a few weeks of wall-clock time. We present applications of MSTAR on few particle systems, studying the Kozai mechanism and N-body systems like star clusters with up to N-part = 10(4) particles. Combined with a tree or fast multipole-based integrator, the high performance of MSTAR removes a major computational bottleneck in simulations with regularized subsystems. It will enable the next-generation galactic-scale simulations with up to 109 stellar particles (e.g. m(star) = 100 M-circle dot) for an M-star = 10(11) M-circle dot galaxy), including accurate collisional dynamics in the vicinity of nuclear supermassive black holes.
Subject: gravitation
methods: numerical
quasars: supermassive black holes
galaxies: star clusters: general
SUPERMASSIVE BLACK-HOLES
GALACTIC NUCLEI
GALAXY MERGERS
EVOLUTION
BINARIES
IMPLEMENTATION
PERTURBATIONS
SYSTEMS
SPACE
CODE
115 Astronomy, Space science
Peer reviewed: Yes
Rights: unspecified
Usage restriction: openAccess
Self-archived version: acceptedVersion


Files in this item

Total number of downloads: Loading...

Files Size Format View
staa084.pdf 2.490Mb PDF View/Open
2001.03180.pdf 1.941Mb PDF View/Open

This item appears in the following Collection(s)

Show full item record