A Data-Driven Census of Cislunar Orbital Stability Enabled by Volunteer Computing

The cislunar space, governed by the circular restricted three-body problem (CR3BP), presents significant challenges for mission design due to its complex stability structure. Traditional high-fidelity numerical integration is computationally prohibitive for a systematic stability census of millions of orbits. Here, we present a novel approach based on global volunteer computing via the BOINC platform to overcome this barrier. Using the public "Million Orbit" dataset from Lawrence Livermore National Laboratory, we distributed the computation of Jacobi constant time series across thousands of volunteer devices, producing over 16 billion individual values. The resulting dataset is freely available. Analysis reveals that 91.68% of orbits belong to the high-energy Region V, 8.07% to the stable Region I, and only 0.24% to Region III, with Region II completely absent. A single rare Region IV orbit (ID 754482) was identified and analyzed. This work demonstrates the transformative potential of volunteer computing for large-scale astrodynamics research, providing a detailed stability map and a benchmark for future machine-learning applications.