Related papers: An adaptive algorithm for n-body field expansions
This report describes a modification of the orthogonal function Poisson solver for n-body simulations that minimizes relaxation caused by small particle number fluctuations. With the standard algorithm, the noise leading to relaxation can…
The stability of a galaxy model is most easily assessed through N-body simulation. Particle-mesh codes have been widely used for this purpose, since they enable the largest numbers of particles to be employed. We show that the functional…
We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be…
N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new…
We discuss a basis set developed to calculate perturbation coefficients in an expansion of the general N-body problem. This basis has two advantages. First, the basis is complete order-by-order for the perturbation series. Second, the…
We present the N-body simulation techniques in EXP. EXP uses empirically-chosen basis functions to expand the potential field of an ensemble of particles. Unlike other basis function expansions, the derived basis functions are adapted to an…
This work discusses the main analogies and differences between the deterministic approach underlying most cosmological N-body simulations and the probabilistic interpretation of the problem that is often considered in mathematics and…
Group theoretic and graphical techniques are used to derive the N-body wave function for a system of identical bosons with general interactions through first-order in a perturbation approach. This method is based on the maximal symmetry…
The oscillator bases expansion stands as an efficient approximation method for the time-independent Schr\"odinger equation. The method, originally formulated with one non-linear variational parameter, can be extended to incorporate two such…
Gas-poor galaxies can be modelled as composite collisionless stellar systems, with a dark matter halo and one or more stellar components, representing different stellar populations. The dynamical evolution of such composite systems is often…
We have developed an adaptive multigrid code for solving the Poisson equation in gravitational simulations. Finer rectangular subgrids are adaptively created in locations where the density exceeds a local level-dependent threshold. We…
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…
We show how to develop an expansion of nearly oblate systems in terms of a set of potential-density pairs. A harmonic (multipole) structure is imposed on the potential set at infinity, and the density can be made everywhere regular. We…
We present a new, simple, fast algorithm to numerically evolve disks of inelastically colliding particles surrounding a central star. Our algorithm adds negligible computational cost to the fastest existing collisionless N-body codes, and…
Over the last two decades, several fast, robust, and high-order accurate methods have been developed for solving the Poisson equation in complicated geometry using potential theory. In this approach, rather than discretizing the partial…
A particle method for reproducing the phase space of collisionless stellar systems is described. The key idea originates in Liouville's theorem which states that the distribution function (DF) at time t can be derived from tracing necessary…
We introduce the N-body simulation technique to follow structure formation in linear and nonlinear regimes for the extended quintessence models (scalar-tensor theories in which the scalar field has a self-interaction potential and behaves…
In the gravitational evolution of a cold infinite particle distribution, two-body interactions can be predominant at early times: we show that, by treating the simple case of a Poisson particle distribution in a static universe as an…
We present a new high-resolution N-body algorithm for cosmological simulations. The algorithm employs a traditional particle-mesh technique on a cubic grid and successive multilevel relaxations on the finer meshes, introduced recursively in…
We describe a new, adaptive solver for the two-dimensional Poisson equation in complicated geometries. Using classical potential theory, we represent the solution as the sum of a volume potential and a double layer potential. Rather than…