Towards optimal softening in threedimensional Nbody codes  I. Minimizing the force error
Abstract
In Nbody simulations of collisionless stellar systems, the forces are softened to reduce the large fluctuations owing to shot noise. Softening essentially modifies the law of gravity at r=x_{i}x_{j} smaller than some softening length ɛ. Therefore, the softened forces are increasingly biased for ever larger ɛ, and there is some optimal ɛ which yields the best compromise between reducing the fluctuations and introducing a bias. Here, analytical relations are derived for the amplitudes of the bias and the fluctuations in the limit of ɛ being much smaller than any physical scale of the underlying stellar system. In particular, it is shown that the fluctuations of the force are generated locally, in contrast to the variations of the potential, which originate from noise in the whole system. Based on these asymptotic relations and using numerical simulations, I study the dependence of the resulting force error on the number of bodies, the softening length, and on the functional form by which Newtonian gravity is replaced. The widely used Plummer softening, where each body is replaced by a Plummer sphere of scale radius ɛ, yields significantly larger force errors than do methods in which the bodies are replaced by density kernels of finite extent. I also give special kernels, which reduce the errors even further. These kernels largely compensate for the errors made with too small interparticle forces at r<ɛ by exceeding Newtonian forces at r~ɛ. Additionally, the possibilities of locally adapting ɛ and of using unequal weights for the bodies are investigated. By using these various techniques without increasing N, the rms force error can be reduced by a factor ~2 when compared with the standard Plummer softening with constant ɛ. The results of this study are directly relevant to Nbody simulations using direct summation techniques or the tree code for force evaluation.
 Publication:

Monthly Notices of the Royal Astronomical Society
 Pub Date:
 June 2001
 DOI:
 10.1046/j.13658711.2001.04237.x
 arXiv:
 arXiv:astroph/0011568
 Bibcode:
 2001MNRAS.324..273D
 Keywords:

 STELLAR DYNAMICS;
 METHODS: NBODY SIMULATIONS;
 METHODS: NUMERICAL;
 Astrophysics
 EPrint:
 20 pages, 8 figs, accepted for publication in MNRAS