Optimal Choice of the Softening Length and Time-Step in N-body Simulations
S. A. Rodionov, N. Ya. Sotnikova
Abstract
A criterion for the choice of optimal softening length ε 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 parameters of an initially equilibrium stable model depend on ε and dt. These dependences are used to derive a criterion for choosing ε and dt. The resulting criterion is compared to Merritt's criterion for choosing the softening length, which is based on minimizing the mean irregular force acting on a particle with unit mass. Our criterion for choosing ε and dt indicate that ε must be a factor of 1.5-2 smaller than the mean distance between particles in the densest regions to be resolved. The time-step must always be adjusted to the chosen ε (the particle must, on average, travel a distance smaller than 0.5ε during one time-step). An algorithm for solving N-body problems with adaptive variations of the softening length is discussed in connection with the task of choosing ε, but is found not to be promising.
Create a lesson
Related papers
On binary pulsars and the force of gravity
Davor Palle
Tidal torques. A critical review of some techniques
Michael Efroimsky, James G. Williams
Dynamics of a Spherical Accretion Shock with Neutrino Heating and Alpha-Particle Recombination
Rodrigo Fernández, Christopher Thompson
Asymptotically FRW black holes
J. T. Firouzjaee, Reza Mansouri
Reaction of Accretion Disks to Abrupt Mass Loss During Binary Black Hole Merger
Sean M. O'Neill, M. Coleman Miller, Tamara Bogdanovic et al.
A Gamma-Ray Burst/Pulsar for Cosmic-Ray Positrons with a Dark Matter-like Spectrum
Kunihito Ioka