KRIOS: A new basis-expansion $N$-body code for collisional stellar dynamics

Abstract

The gravitational $N$-body problem is a nearly universal problem in astrophysics which, despite its deceptive simplicity, still presents a significant computational challenge. For collisional systems such as dense star clusters, the need to resolve individual encounters between $N$ stars makes the direct summation of forces - with quadratic complexity - almost infeasible for systems with $N\gtrsim 10^6$ particles over many relaxation times. At the same time, the most common Monte Carlo $N$-body algorithm - that of Hénon - assumes the cluster to be spherically symmetric. This greatly limits the study of many important features of star clusters, including triaxiality, rotation, and the production of tidal debris. In this paper, we present a new hybrid code, KRIOS, that combines 3D collisionless relaxation using an adaptive self-consistent field method with collisional dynamics handled via Hénon's method. We demonstrate that KRIOS can accurately model the long-term evolution of clusters and provide its complete phase-space information over many relaxation times. As a test of our new code, we present detailed comparisons to well-known results from stellar dynamics: (i) the collisional evolution of a family of Plummer spheres with varying anisotropy and rotation to core collapse, and (ii) the emergence of the radial-orbit instability in radially anisotropic star clusters, including its non-spherical effects.

23 pages, 15 figures, submitted to ApJ

References in corpus (49)

• First results from the IllustrisTNG simulations: matter and galaxy clustering
• Resolved Massive Star Clusters in the Milky Way and its Satellites: Brightness Profiles and a Catalogue of Fundamental Parameters
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• NBODY6++GPU: Ready for the gravitational million-body problem
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• A stellar census in globular clusters with MUSE: The contribution of rotation to cluster dynamics studied with 200 000 stars
• Massive Black Holes in Star Clusters. II. Realistic Cluster Models
• Tidal Tails of Star Clusters
• Monte-Carlo Simulations of Globular Cluster Evolution - I. Method and Test Calculations
• Hubble Space Telescope proper motion (HSTPROMO) catalogs of Galactic globular clusters. II. Kinematic profiles and maps
• MOCCA code for star cluster simulations - I. Blue Stragglers, first results
• The State-of-the-Art HST Astro-Photometric Analysis of the core of $\omega$ Centauri. I. The Catalog
• The internal rotation of globular clusters revealed by Gaia DR2
• Evolution of star clusters in arbitrary tidal fields
• Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code
• Rotating Globular Clusters
• The eye of Gaia on globular clusters kinematics: internal rotation
• Mergers and ejections of black holes in globular clusters
• Central rotations of Milky Way Globular Clusters
• N-body Models of Rotating Globular Clusters
• An adaptive algorithm for n-body field expansions
• MIKiS: the Multi-Instrument Kinematic Survey of Galactic Globular Clusters. I. Velocity dispersion profiles and rotation signals of 11 globular clusters
• Kinetic theory of spatially inhomogeneous stellar systems without collective effects
• A fully parallel, high precision, N-body code running on hybrid computing platforms
• A new Fokker-Planck approach for relaxation-driven evolution of galactic nuclei
• Comparative study between N-body and Fokker-Planck simulations for rotating star clusters: I. Equal-mass system
• Comparative study between N-body and Fokker-Planck simulations for rotating star clusters - II. 2-component models
• A new Monte Carlo method for dynamical evolution of non-spherical stellar systems
• A Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces
• A Halo Expansion (HEX) Technique for Approximating Simulated Dark Matter Haloes
• Monte-Carlo Simulations of Star Clusters I. First Results
• The impact of stellar evolution on rotating star clusters: the gravothermal-gravogyro catastrophe and the formation of a bar of black holes
• Great Balls of FIRE II: The evolution and destruction of star clusters across cosmic time in a Milky Way-mass galaxy
• An efficient parallel algorithm for O(N^2) direct summation method and its variations on distributed-memory parallel machines
• The orbital anisotropy profiles of nearby globular clusters from Gaia Data Release 2
• EXP: N-body integration using basis function expansions
• The Kinematic Richness of Star Clusters I. Isolated Spherical Models with Primordial Anisotropy
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• On the nature of the radial orbit instability in spherically symmetric collisionless stellar systems
• Galaxy halo expansions: a new biorthogonal family of potential-density pairs
• Towards time symmetric N-body integration
• Should $N$-body integrators be symplectic everywhere in phase space?
• New biorthogonal potential--density basis functions
• Radial orbit instability: review and perspectives
• The hunt for self-similar core collapse
• Euclid: Early Release Observations -- Unveiling the morphology of two Milky Way globular clusters out to their periphery
• A general basis set algorithm for galactic haloes and discs
• Non-resonant relaxation of anisotropic globular clusters
• Non-resonant relaxation of rotating globular clusters