Incomplete relaxation in a two-mass one-dimensional self-gravitating system.

Abstract

Due to the apparent ease with which they can be numerically simulated, one-dimensional gravitational systems were first introduced by astronomers to explore different modes of gravitational evolution. These include violent relaxation and the approach to thermal equilibrium. Careful work by dynamicists and statistical physicists has shown that several claims made by astronomers regarding these models were incorrect. Unusual features of the evolution include the development of long lasting structures on large scales, which can be thought of as one-dimensional analogs of Jupiter's red spot or a galactic spiral density wave or bar. The existence of these structures demonstrates that in gravitational systems evolution is not entirely dominated by the second law of thermodynamics and also appears to contradict the Arnold diffusion ansatz. Thus it is correct to assert that the one-dimensional planar sheet gravitational system is the nonextensive analog of the Fermi-Pasta-Ulam model of dynamical systems. This paper is an extension of a preliminary study where we conclusively showed mass segregation and equipartition of kinetic energy in a two-mass planar sheet system for the first time. Here we employ both mean-field theory and dynamical simulation to more thoroughly probe the statistical and ergodic properties of these systems. Valuable information is obtained from local and global time averaging, and temporal and spatial correlation functions. Using these tools we show that the system appears to approach the equilibrium distribution on very long time scales, but the relaxation is incomplete.