Related papers: Multi-mass schemes for collisionless N-body simula…
We review the field of collisionless numerical simulations for the large-scale structure of the Universe. We start by providing the main set of equations solved by these simulations and their connection with General Relativity. We then…
Next generation galaxy surveys demand the development of massive ensembles of galaxy mocks to model the observables and their covariances, what is computationally prohibitive using $N$-body simulations. COLA is a novel method designed to…
Special high-accuracy direct force summation N-body algorithms and their relevance for the simulation of the dynamical evolution of star clusters and other gravitating N-body systems in astrophysics are presented, explained and compared…
We describe an algorithm for constructing N-body realizations of equilibrium spherical systems. A general form for the mass density rho(r) is used, making it possible to represent most of the popular density profiles found in the…
We report on a series of tests of agreement between three types of N-body simulations: PM, P$^3$M, and Tree codes. We find good agreement in both the individual and the statistical properties only on scales larger than the mean…
We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against $N$-body simulations…
We perform a detailed comparison of the phase-space density traced by the particle distribution in Gadget simulations to the result obtained with a spherical Vlasov solver using the splitting algorithm. The systems considered are apodized…
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…
This paper describes a method for quantitatively comparing an N-body model with a sample of discrete kinematic data. The comparison has two stages: (i) finding the optimum scaling and orientation of the model relative to the data; and (ii)…
Galaxy observations and N-body cosmological simulations produce conflicting dark matter halo density profiles for galaxy central regions. While simulations suggest a cuspy and universal density profile (UDP) of this region, the majority of…
We present a new method for constructing equilibrium phase models for stellar systems, which we call the iterative method. It relies on constrained, or guided evolution, so that the equilibrium solution has a number of desired parameters…
Numerical simulations of self-gravitating systems are generally based on N-body codes, which solve the equations of motion of a large number of interacting particles. This approach suffers from poor statistical sampling in regions of low…
Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with…
(Abridged) We examine the theoretical relationship between Omega_0 and substructure in galaxy clusters which are formed by the collapse of high density peaks in a gaussian random field. The radial mass distributions of the clusters are…
We use high resolution dissipationless N-body simulations to examine the spatial distribution of galaxy clusters on large scales. The Standard CDM model and two of its main competitors, Low density CDM and Mixed Dark Matter are compared.The…
Purpose: To develop and apply multiscale Monte Carlo (MC) simulations to assess variations in nucleus and cytoplasm dose enhancement factors (n,cDEFs) over tumor-scale volumes. Methods: The intrinsic variation of n,cDEFs (due to variations…
Gravitational N-body simulations, that is numerical solutions of the equations of motions for N particles interacting gravitationally, are widely used tools in astrophysics, with applications from few body or solar system like systems all…
N-body dynamical simulations are used to analyze the conditions for the gravitational stability of a three-dimensional stellar disk in the gravitational field of two rigid spherical components--a bulge and a halo whose central…
Performing a shell model calculation for heavy nuclei has been a long-standing problem in nuclear physics. Here we propose one possible solution. The central idea of this proposal is to take the advantages of two existing models, the…
We use N-body simulation to study the structure formation in the Cubic Galileon Gravity model where along with the usual kinetic and potential term we also have a higher derivative self-interaction term. We find that the large scale…