Deterministic and probabilistic neural surrogates of global hybrid-Vlasov simulations

TL;DR

This paper demonstrates that graph neural network emulators can accurately and rapidly predict plasma states in hybrid-Vlasov simulations, significantly reducing computational costs while maintaining high accuracy.

Contribution

The authors develop both deterministic and probabilistic GNN-based emulators that learn from 5D hybrid-Vlasov simulation data, enabling fast and accurate plasma state predictions.

Findings

Emulators achieve over 100x speedup compared to traditional simulations.

High correlation (above 0.95) maintained at 50 seconds lead time.

Probabilistic model improves ensemble forecast calibration.

Abstract

Hybrid-Vlasov simulations resolve ion-kinetic effects in the solar wind-magnetosphere interaction, but even 5D (2D + 3V) configurations are computationally expensive. We show that graph-based machine learning emulators can learn the spatiotemporal evolution of electromagnetic fields and lower order moments of ion velocity distribution in the near-Earth space environment from four 5D Vlasiator runs performed with identical steady solar wind conditions. The initial ion number density is systematically varied, while the grid spacing is held constant, to scan the ratio of the characteristic ion skin depth to the numerical grid size. Using a graph neural network (GNN) operating on the 2D spatial simulation grid comprising 670k cells, we demonstrate that both a deterministic forecasting model (Graph-FM) and a probabilistic ensemble forecasting model (Graph-EFM) based on a latent variable formulation are capable of producing accurate predictions of future plasma states. A divergence penalty is incorporated to encourage divergence-freeness in the magnetic fields. For the probabilistic model, a continuous ranked probability score objective is added to improve the calibration of the ensemble forecasts. The trained emulators achieve over two orders of magnitude speedup per time step on a single GPU compared to 100 CPU Vlasiator simulations. Most forecasted fields have Pearson correlations above 0.95 at 50 seconds lead time. However, we find that fields that exhibit near-zero degenerate distributions in the 5D setting are more challenging for the emulator to maintain high correlations for. Overall, these results demonstrate that GNNs provide a viable framework for rapid ensemble generation in hybrid-Vlasov modeling and highlight promising directions for future work.