Chaos in self-gravitating many-body systems: Lyapunov time dependence of $N$ and the influence of general relativity

In self-gravitating $N$-body systems, small perturbations introduced at the start, or infinitesimal errors produced by the numerical integrator or due to limited precision in the computer, grow exponentially with time. For Newton's gravity, we confirm earlier results by \cite{1992ApJ...386..635K} and \cite{1993ApJ...415..715G}, that for relatively homogeneous systems, this rate of growth per crossing time increases with $N$ up to $N \sim 30$, but that for larger systems, the growth rate has a weaker dependency with $N$. For concentrated systems, however, the rate of exponential growth continues to scale with $N$. In relativistic self-gravitating systems, the rate of growth is almost independent of $N$. This effect, however, is only noticeable when the system's mean velocity approaches the speed of light to within three orders of magnitude. The chaotic behavior of systems with $\apgt 10$ bodies for the usually adopted approximation of only solving the pairwise interactions in the Einstein-Infeld-Hoffmann equation of motion, is qualitatively different than when the interaction terms (or cross terms) are taken into account. This result provides a strong motivation for follow-up studies on the microscopic effect of general relativity on orbital chaos, and the influence of higher-order cross-terms in the Taylor-series expansion of the EIH equations of motion.