“…The evolution under Newtonian self-gravity of an isolated distribution of particles that start from initial conditions (ICs) in which the particles are initially at rest and are randomly distributed in a sphere with uniform mean density has been studied by numerous authors since the earliest numerical simulations of such systems. Numerous variants have also been explored, notably with initial velocity dispersion (using radially dependent density profiles) and for spheroidal boundaries (Henon, 1964(Henon, , 1973van Albada, 1982;Merritt & Aguilar, 1985;Aarseth et al, 1988;Aguilar & Merritt, 1990;Theis & Spurzem, 1999;Boily et al, 2002;Roy & Perez, 2004;Boily & Athanassoula, 2006;Barnes et al, 2009;Joyce et al, 2009;Sylos Labini, 2012, 2013Worrakitpoonpon, 2015;Sylos Labini et al, 2015;Benhaiem et al, 2016;Benhaiem & Sylos Labini, 2017;Spera & Capuzzo-Dolcetta, 2017;Worrakitpoonpon, 2020). Such systems provide a natural and simple setting to explore the rich complexity of the nonlinear phase of gravitational dynamics relevant to many different problems in astrophysics and cosmology.…”