Equivariant optimisation for the gravitational n-body problem: A computational factory of symmetric orbits

In this paper we present SymOrb.jl (2024), a software which combines group representation theory and variational methods to provide numerical solutions of singular dynamical systems of paramount relevance in Celestial Mechanics and other interacting particles models. Among all, it prepares for large-scale search of symmetric periodic orbits for the classical n-body problem and their classification, paving the way towards a computational validation of Poincaré conjecture about the density of periodic orbits. Through the accessible language of Julia, Symorb.jl offers a unified implementation of an earlier version (dlfer – symorb, 2017). This paper provides theoretical and practical guidelines for the specific approach we adopt, complemented with examples.