Related papers: Optimal softening for force calculations in collis…
In N-body simulations of collisionless stellar systems, the forces are softened to reduce the shot noise. Softening modifies gravity at r=|x-y| smaller than softening length epsilon and the softened forces are increasingly biased for ever…
Pairwise forces between particles in cosmological N-body simulations are generally softened to avoid hard collisions. Physically, this softening corresponds to treating the particles as diffuse clouds rather than point masses. For particles…
A criterion for the choice of optimal softening length $\epsilon$ and time-step $dt$ for $N$-body simulations of a collisionless stellar system is analyzed. Plummer and Hernquist spheres are used as models to follow how changes in various…
In any collisionless N-body code, there is an optimal choice for the smoothing parameter that minimizes the average error in the force evaluations. We show how to compute the optimal softening length in a direct-summation code and…
In this paper we describe an adaptive softening length formalism for collisionless N-body and self-gravitating Smoothed Particle Hydrodynamics (SPH) calculations which conserves momentum and energy exactly. This means that spatially…
We discuss capability of Smooth Particle Hydrodynamics to represent adequately the dynamics of self-gravitating systems, in particular for what regards the quality of approximation of force fields in the motion equations. When cubic spline…
When starting an N-body simulation of an isolated galaxy, it is desirable to select particles from a distribution function to ensure that the model is in equilibrium. Random sampling from a DF is widely used, but results in a set of…
Dynamical simulations are a fundamental tool for studying the secular evolution of disc galaxies. Even at their maximum resolution, they still follow a limited number of particles and typically resolve scales of the order of a few tens of…
We propose to determine the optimal softening length in N-body halo simulations by minimizing the ensemble-average acceleration error at a small radius r0. This strategy ensures that the error never exceeds the optimal value beyond r0.…
Using the Ricci and scalar curvatures of the configuration manifold of gravitational N-body systems, we study the exponential instability in their trajectories. It is found that the exponentiation time-scale for isotropic Plummer spheres…
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…
Analysis of self-similarity in scale-free $N$-body simulations reveals the spatial and temporal scales for which statistics measured in cosmological simulations are converged to the physical continuum limit. We examine how the range of…
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve…
In self-consistent N-body simulations of collisionless systems, gravitational interactions are modified on small scales to remove singularities and simplify the task of numerically integrating the equations of motion. This `gravitational…
We present a systematic investigation of the distribution of normal forces at the boundaries of static packings of spheres. A new method for the efficient construction of large hexagonal-close-packed crystals is introduced and used to study…
Modelling gravity is a fundamental problem that must be tackled in N-body simulations of stellar systems, and satisfactory solutions require a deep understanding of the dynamical effects of softening. In a previous paper (Romeo 1997), we…
The paper considers functional linear regression, where scalar responses $Y_1,...,Y_n$ are modeled in dependence of random functions $X_1,...,X_n$. We propose a smoothing splines estimator for the functional slope parameter based on a…
Gravitational softening length is one of the key parameters to properly set up a cosmological $N$-body simulation. In this paper, we perform a large suit of high-resolution $N$-body simulations to revise the optimal softening scheme…
We present a two-dimensional (2D) Particle-Particle-Particle-Mesh (P$^3$M) algorithm with an optimized Green function and adaptive softening length for gravitational lensing studies in N-Body simulations. The analytical form of the…
Effects of gravitational softening on the global structure of self-gravitating disks in centrifugal equillibrium are examined in relation to hydrodynamical/gravitational simulations. The one-parameter spline softening proposed by Hernquist…