The First Star-by-star -body/Hydrodynamics Simulation of Our Galaxy Coupling with a Surrogate Model
arXiv:2510.23330 · doi:10.1145/3712285.3759866
Abstract
A major goal of computational astrophysics is to simulate the Milky Way Galaxy with sufficient resolution down to individual stars. However, the scaling fails due to some small-scale, short-timescale phenomena, such as supernova explosions. We have developed a novel integration scheme of -body/hydrodynamics simulations working with machine learning. This approach bypasses the short timesteps caused by supernova explosions using a surrogate model, thereby improving scalability. With this method, we reached 300 billion particles using 148,900 nodes, equivalent to 7,147,200 CPU cores, breaking through the billion-particle barrier currently faced by state-of-the-art simulations. This resolution allows us to perform the first star-by-star galaxy simulation, which resolves individual stars in the Milky Way Galaxy. The performance scales over CPU cores, an upper limit in the current state-of-the-art simulations using both A64FX and X86-64 processors and NVIDIA CUDA GPUs.
arXiv comment12 pages, 7 figures, 7 tables, IEEE/ACM Supercomputing Conference (SC25)
References in corpus (20)
- The cosmological simulation code GADGET-2
- FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
- The mass distribution and gravitational potential of the Milky Way
- AGAMA: Action-based galaxy modelling architecture
- The Arepo public code release
- Ultra-faint dwarfs in a Milky Way context: Introducing the Mint Condition DC Justice League Simulations
- Variable interstellar radiation fields in simulated dwarf galaxies: supernovae versus photoelectric heating
- Supernova feedback in numerical simulations of galaxy formation: separating physics from numerics
- Efficient early stellar feedback can suppress galactic outflows by reducing supernova clustering
- Determining the full satellite population of a Milky Way-mass halo in a highly resolved cosmological hydrodynamic simulation
- 24.77 Pflops on a Gravitational Tree-Code to Simulate the Milky Way Galaxy with 18600 GPUs
- Equilibrium Initialization and Stability of Three-Dimensional Gas Disks
- Modeling the Milky Way as a Dry Galaxy
- Field Level Neural Network Emulator for Cosmological N-body Simulations
- Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods
- Co-evolution of Dust and Chemistry in Galaxy Simulations with a Resolved Interstellar Medium
- The effects of local stellar radiation and dust depletion on non-equilibrium interstellar chemistry
- 4.45 Pflops Astrophysical N-Body Simulation on K computer -- The Gravitational Trillion-Body Problem
- Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning
- Surrogate Modeling for Computationally Expensive Simulations of Supernovae in High-Resolution Galaxy Simulations