English

Gravitational Softening and Adaptive Mass Resolution

Astrophysics 2008-04-10 v3

Abstract

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 of unequal mass (and hence unequal softening length), computing the softened force involves a nontrivial double integral over the volumes of the two particles. We show that Plummer force softening is consistent with this interpretation of softening while spline softening is not. We provide closed-form expressions and numerical implementation for pairwise gravitational force laws for pairs of particles of general softening scales ϵ1\epsilon_1 and ϵ2\epsilon_2 assuming the commonly used cloud profiles: NGP, CIC, TSC, and PQS. Similarly, we generalize Plummer force law into pairs of particles of general softenings. We relate our expressions to the gaussian, Plummer and spline force softenings known from literature. Our expressions allow possible inclusions of pointlike particles such as stars or supermassive black holes.

Keywords

Cite

@article{arxiv.0711.2989,
  title  = {Gravitational Softening and Adaptive Mass Resolution},
  author = {Alexander Shirokov},
  journal= {arXiv preprint arXiv:0711.2989},
  year   = {2008}
}

Comments

27 pages, 6 figures; Minor additions; WSHAPE package available at http://www.gracos.org/wshape; Submitted to the Journal of Computational Physics

R2 v1 2026-06-21T09:44:58.636Z