Related papers: On the reliability of simulations of collisional s…
The numerical simulations of massive collisional stellar systems, such as globular clusters (GCs), are very time-consuming. Until now, only a few realistic million-body simulations of GCs with a small fraction of binaries (5%) have been…
We report results of collisional N-body simulations aimed to study the N-dependance of the dynamical evolution of star clusters. Our clusters consist of equal-mass stars and are in virial equilibrium. Clusters moving in external tidal…
Stellar systems are broadly divided into collisional and non-collisional. The latter are large-N systems with long relaxation timescales and can be simulated disregarding two-body interactions, while either computationally expensive direct…
The past few years have seen dramatic improvements in the scope and realism of star cluster simulations. Accurate treatments of stellar evolution, coupled with robust descriptions of all phases of binary evolution, have been incorporated…
We report results of N-body simulations of isolated star clusters, performed up to the point where the clusters are nearly completely dissolved. Our main focus is on the post-collapse evolution of these clusters. We find that after core…
This paper, and its companion, investigate the evolution of dense stellar systems due to the influence of two-body gravitational encounters, physical collisions and stellar evolution. Our goal is the simulation of the densest centers of…
The P$^3$T scheme is a hybrid method for simulating gravitational $N$-body systems. It combines a fast particle-tree (PT) algorithm for long-range forces with a high-accuracy particle-particle (PP, direct $N$-body) solver for short-range…
A revision of Stodolkiewicz's Monte-Carlo code is used to simulate evolution of star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The first…
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…
Stellar collisions are an important formation channel for blue straggler stars in globular and old open clusters. Hydrodynamical simulations have shown that the remnants of such collisions are out of thermal equilibrium, are not strongly…
Both simulations and observations indicate that stars form in filamentary, hierarchically clustered associations, most of which disperse into their galactic field once feedback destroys their parent clouds. However, during their early…
We address the dynamical evolution of an isolated self--gravitating system with two stellar mass groups. We vary the individual ratio of the heavy to light bodies, $\mu$ from 1.25 to 50 and alter also the fraction of the total heavy mass…
The presence of protostellar disks can greatly increase the dissipation during close stellar encounters, leading to the formation of a significant population of binaries during the initial collapse and virialization of a cluster. We have…
Planet-planet collisions are a common outcome of instability in systems of transiting planets close to the star, as well as occurring during in-situ formation of such planets from embryos. Previous N-body studies of instability amongst…
Direct collisions between finite-sized particles occur commonly in many areas of astrophysics. Such collisions are typically mediated by chaotic, bound gravitational interactions involving small numbers of particles. An important…
In a Keplerian system, a large number of bodies orbit a central mass. Accretion disks, protoplanetary disks, asteroid belts, and planetary rings are examples. Simulations of these systems require algorithms that are computationally…
$N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be…
We study the dynamical evolution of idealised stellar systems by averaging results from many $N$-body simulations, each having modest numbers of stars. For isolated systems with stars of uniform mass, we discuss aspects of evolution up to…
We discuss the performance of direct summation codes used in the simulation of astrophysical stellar systems on highly distributed architectures. These codes compute the gravitational interaction among stars in an exact way and have an…
We compare the results for the dynamical evolution of star clusters derived from anisotropic gaseous models with the data from N-body simulations of isolated and one-component systems, each having modest number of stars. The statistical…