Numerical irreversibility in self-gravitating small N-body systems, II: Influence of instability affected by softening parameters

Nobuyoshi Komatsu, Takahiro Kiwata and Shigeo Kimura

Physica A: Statistical Mechanics and its Applications, 2009, vol. 388, issue 5, 639-650

Abstract: We investigate the fundamental characteristics of numerical irreversibility appearing in self-gravitating small N-body systems by means of a molecular dynamics method from the viewpoint of time-reversible dynamics. We reconsider a closed spherical system consisting of 250 point-particles interacting through the Plummer softened potential. To investigate the characteristics of numerical irreversibility, we examine the influence of the instability affected by the softening parameter for the softened potential (the instability considered here is the instability of a dynamical system in chaos theory, e.g., a separation rate of the distance between two nearby trajectories in phase space or speed space). To this end, under the restriction of constant initial energy, the softening parameter for the Plummer softened potential is varied from 0.005R to 0.050R, where R is the radius of the spherical container. We first confirm that the size of the softening parameter, i.e., the deviation of the potential from a pure gravitational potential, influences the virial ratio, the density, the pressure on the spherical container, etc., during an early stage of the relaxation process. Through a time-reversible simulation based on a velocity inversion technique, we demonstrate that numerical irreversibility due to round-off errors appears more rapidly with decreasing softening parameter. This means that the higher the instability of the system or the higher the separation rate of the distance between two nearby trajectories, the earlier the memory of the initial conditions is lost. We show that the memory loss time tm, when the simulated trajectory completely forgets its initial conditions, increases approximately linearly with the timescale of the chaotic system, i.e., the Lyapunov time tλ. In a small self-gravitating system, propagation of numerical irreversibility or loss of reversibility depends on both the energy state of the system and the instability affected by the softening parameter.

Keywords: Self-gravitating system; Irreversibility; Instability; Round-off errors; N-body simulation (search for similar items in EconPapers)
Date: 2009
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S037843710800945X
Full text for ScienceDirect subscribers only. Journal offers the option of making the article available online on Science direct for a fee of $3,000

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:phsmap:v:388:y:2009:i:5:p:639-650

DOI: 10.1016/j.physa.2008.11.012

Access Statistics for this article

Physica A: Statistical Mechanics and its Applications is currently edited by K. A. Dawson, J. O. Indekeu, H.E. Stanley and C. Tsallis

More articles in Physica A: Statistical Mechanics and its Applications from Elsevier
Bibliographic data for series maintained by Catherine Liu ().