Towards optimal softening in 3D N-body codes: I. Minimizing the force error

Details Towards optimal softening in 3D N-body codes: I. Minimizing the force error Towards optimal softening in 3D N-body codes: I. Minimizing the force error

2000-11-30

arxiv.org/abs/astro-ph/0011568v1

Mon.Not.Roy.Astron.Soc. 324 (2001) 273

Astronomy and Astrophysics

Astrophysics

20 pages, 8 figs, accepted for publication in MNRAS

Scientific paper

10.1046/j.1365-8711.2001.04237.x

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 larger epsilon. There is, thus, some optimum 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 small epsilon and large N. 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 the asymptotic relations and using numerical experiments, I study the dependence of the resulting force error on N, epsilon, and on the functional form by which Newtonian gravity is replaced. The Plummer softening, where each body is replaced by a Plummer sphere of scale radius epsilon, 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 the errors made with too small inter-particle forces at r

Affiliated with

Also associated with

No associations

Rating

Towards optimal softening in 3D N-body codes: I. Minimizing the force error does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Towards optimal softening in 3D N-body codes: I. Minimizing the force error, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Towards optimal softening in 3D N-body codes: I. Minimizing the force error will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-249242