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

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Tep, Kerwann, Cook, Brian T., Rodriguez, Carl L., Jolly, Jiya, Sawin, Eddie, Petersen, Michael S., Gaffud, Christoph
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866912741990596608
author Tep, Kerwann
Cook, Brian T.
Rodriguez, Carl L.
Jolly, Jiya
Sawin, Eddie
Petersen, Michael S.
Gaffud, Christoph
author_facet Tep, Kerwann
Cook, Brian T.
Rodriguez, Carl L.
Jolly, Jiya
Sawin, Eddie
Petersen, Michael S.
Gaffud, Christoph
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2506_13636
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle KRIOS: A new basis-expansion $N$-body code for collisional stellar dynamics
Tep, Kerwann
Cook, Brian T.
Rodriguez, Carl L.
Jolly, Jiya
Sawin, Eddie
Petersen, Michael S.
Gaffud, Christoph
Astrophysics of Galaxies
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.
title KRIOS: A new basis-expansion $N$-body code for collisional stellar dynamics
topic Astrophysics of Galaxies
url https://arxiv.org/abs/2506.13636