EXP: N-body integration using basis function expansions

Abstract

We present the N-body simulation techniques in EXP. EXP uses empirically-chosen basis functions to expand the potential field of an ensemble of particles. Unlike other basis function expansions, the derived basis functions are adapted to an input mass distribution, enabling accurate expansion of highly non-spherical objects, such as galactic discs. We measure the force accuracy in three models, one based on a spherical or aspherical halo, one based on an exponential disc, and one based on a bar-based disc model. We find that EXP is as accurate as a direct-summation or tree-based calculation, and in some ways is better, while being considerably less computationally intensive. We discuss optimising the computation of the basis function representation. We also detail numerical improvements for performing orbit integrations, including timesteps.

Submitted to MNRAS, comments welcome

References in corpus (5)

• The NumPy array: a structure for efficient numerical computation
• The dynamics of tidal tails from massive satellites
• Towards time symmetric N-body integration
• Should $N$-body integrators be symplectic everywhere in phase space?
• The Dynamical Response of Dark Matter to Galaxy Evolution Affects Direct-Detection Experiments

Cited by in corpus (19)

• Galactic Archaeology with Gaia
• The effect of the deforming dark matter haloes of the Milky Way and the Large Magellanic Cloud on the Orphan-Chenab stream
• A Comprehensive Perturbative Formalism for Phase-Mixing in Perturbed Disks. I. Phase spirals in an Infinite, Isothermal Slab
• Constraining the Milky Way halo kinematics via its Linear Response to the Large Magellanic Cloud
• Tidally stripped halo stars from the Large Magellanic Cloud in the Galactic North
• Effects on the local dark matter distribution due to the Large Magellanic Cloud
• Shaping the Milky Way: The interplay of mergers and cosmic filaments
• Efficient and accurate force replay in cosmological-baryonic simulations
• Dynamical Data Mining Captures Disc-Halo Couplings that Structure Galaxies
• The Non-Axisymmetric Influence: Radius and Angle-Dependent Trends in a Barred Galaxy
• LMC-driven anisotropic boosts in stream--subhalo interactions
• A general basis set algorithm for galactic haloes and discs
• EXP: a Python/C++ package for basis function expansion methods in galactic dynamics
• Collisionless relaxation from near equilibrium configurations: Linear theory and application to tidal stripping
• How do the dynamics of the Milky Way -- Large Magellanic Cloud system affect gamma-ray constraints on particle dark matter?
• The dark matter wake of a galactic bar revealed by multichannel Singular Spectral Analysis
• KRIOS: A new basis-expansion $N$-body code for collisional stellar dynamics
• lintsampler: Easy random sampling via linear interpolation
• Mapping the Distorted Dark Matter Distribution of the LMC-SMC System Prior to Milky Way Infall with Basis Function Expansions