Related papers: Direct N-body Simulations
We describe an algorithm for constructing N-body realisations of equilibrium stellar systems. The algorithm complements existing orbit-based modelling techniques using linear programming or other optimization algorithms. The equilibria are…
In setting up initial conditions for ensembles of cosmological N-body simulations there are, fundamentally, two choices: either maximizing the correspondence of the initial density field to the assumed fourier-space clustering or, instead,…
Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In…
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…
Large redshift surveys of galaxies and clusters are providing the first opportunities to search for distortions in the observed pattern of large-scale structure due to such effects as gravitational redshift. We focus on non-linear scales…
In the weak field limit general relativity reduces, as is well known, to the Newtonian gravitation. Alternative theories of gravity, however, do not necessarily reduce to Newtonian gravitation; some of them, for example, reduce to…
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…
Over the course of the recent decades, $N$-body simulations have become a standard tool for quantifying the gravitational perturbations that ensue in planet-forming disks. Within the context of such simulations, massive non-central bodies…
The dynamical evolution of dense stellar systems is simulated using a two-dimensional Fokker-Planck method, with the goal of providing a model for the formation of supermassive stars which could serve as seed objects for the supermassive…
Context. Simulating stellar dynamics in a molecular cloud environment is numerically challenging due to the strong coupling between young stars and their surrounding gas, and the large range of length and time scales. Aims. This paper is…
I consider a self-gravitating, N-body system assuming that the N constituents follow regular orbits about the center of mass of the cluster, where a central massive object may be present. I calculate the average over a characteristic…
We follow the sinking of two massive black holes in a spherical stellar system where the black holes become bound under the influence of dynamical friction. Once bound, the binary hardens by three-body encounters with surrounding stars. We…
A revision of Stod\'o{\l}kiewicz's Monte Carlo code is used to simulate evolution of large star clusters. The survey on the evolution of multi-mass N-body systems influenced by the tidal field of a parent galaxy and by stellar evolution is…
The approximate computation of all gravitational forces between $N$ interacting particles via the fast multipole method (FMM) can be made as accurate as direct summation, but requires less than $\mathcal{O}(N)$ operations. FMM groups…
We study the bar instability in collisionless, rotating, anisotropic, stellar systems, using N-body simulations and also the matrix technique for calculation of modes with the perturbed collisionless Boltzmann equation. These methods are…
A wide variety of outstanding problems in astrophysics involve the motion of a large number of particles ($N\gtrsim 10^{6}$) under the force of gravity. These include the global evolution of globular clusters, tidal disruptions of stars by…
We present the first post core collapse models of initially rotating star clusters, using the numerical solution of an orbit-averaged 2D Fokker-Planck equation. Based on the code developed by Einsel & Spurzem (1999), we have improved the…
We present the results of an attempt to adapt the distribution function formalism to characterize large-scale structures like clusters of galaxies that form in a cosmological N-body simulation. While galaxy clusters are systems that are not…
We describe the construction of an N-body simulation that mimics the true velocity and mass-density fields in a box of side 256\hmpc about the Local Group , and the production of mock catalogs that mimic in detail current catalogs of…
We introduce a GPU-accelerated hybrid hydro/N-body code (Enzo-N) designed to address the challenges of concurrently simulating star clusters and their parent galaxies. This task has been exceedingly challenging, primarily due to the…