Source-linked AI summary

Adaptive Interaction Graphs for Particle Simulation

Aiden Zhou

arXiv:2609.30822v1cs.LGcs.CE

TL;DR

Learned particle simulators can be accurate over one step but accumulate error over long rollouts, while fixed interaction graphs do not account for local confidence. ADAPTGNS uses a jointly trained per-particle variance head and prior-step uncertainty to adapt graph connectivity. It improves the edge-count and MSE@200 trade-off on WATERDROP and shows a modest gain on SAND.

  • Problem

    Learned GNN particle simulators can have strong one-step accuracy, but long-horizon errors compound, and fixed geometric interaction graphs ignore local model confidence.

    Source: §1 · p. 1
  • Method

    ADAPTGNS jointly trains acceleration and per-particle variance predictions, then uses prior-step uncertainty to expand interaction neighborhoods in a single-pass adaptive graph loop.

    Source: §4 · p. 3, §7 · p. 5
  • Results

    ADAPTGNS uses 20% fewer edges and has lower MSE@200 than fixed k = 5 on WATERDROP; on SAND it has 12% lower MSE@200 than fixed k = 5 using 10% fewer edges than k = 10.

    Source: §1 · pp. 1–2
  • Takeaways & Limitations

    The authors attribute WATERDROP’s stronger results to its evolving sparse graph and more spatially heterogeneous uncertainty compared with SAND’s more uniform particle dynamics.

    Source: §5.2 · p. 4

    Runtime is represented by average edge count, which tracks the principal cost of a GNN forward pass.

    Source: §6 · p. 4

Abstract

from arXiv · show

Learned particle simulators based on graph neural networks achieve strong one-step accuracy, but errors compound over long horizons. An underexplored variable is the interaction graph: existing methods fix its topology via k-nearest neighbors or a static radius rule, regardless of local model confidence. We propose making this graph adaptive: a per-particle variance head, trained jointly with the acceleration head under a heteroscedastic Gaussian NLL loss, drives a trajectory in which high-uncertainty particles receive an expanded neighborhood. This is done at little extra inference cost by using the previous step's uncertainty estimate. A key discovery is that the variance head learns a meaningful notion of uncertainty: high-variance particles concentrate near complex regions, such as splash zones or free surfaces. When this signal drives graph topology, the resulting AdaptGNS simulator achieves a strict Pareto improvement on WaterDrop and a modest gain on Sand. Given the model's stronger performance on WaterDrop, we hypothesize that adaptive graphs are most useful when complexity is concentrated in space. Our code can be found at https://github.com/aidenzhou8/AdaptGNS.

1. Introduction

Learned particle simulators can perform well one step at a time, but long-horizon errors compound, while fixed interaction graphs ignore local confidence. ADAPTGNS adapts connectivity using predicted uncertainty and improves the accuracy–edge-count trade-off, especially on WATERDROP.

1. Introduction

  • The method selectively expands interaction edges per timestep according to predicted uncertainty, aiming to add relational capacity where confidence is low while preserving sparse connectivity elsewhere.The paper describes this as a way to improve stability while retaining the speed and efficiency of learned simulation.Source: §1 · p. 1
  • Learned simulators offer faster inference than classical solvers, but accumulated errors can erase that benefit over long simulations.The introduction cites speedups of up to 5,000× over MPM on granular flow simulations and orders-of-magnitude cost reductions for neural fluid surrogates.Source: §1 · p. 1
  • The approach is motivated by spatially uneven uncertainty, such as unpredictable fluid splash zones surrounded by quieter regions.The learning task predicts per-particle acceleration and trains a per-particle variance estimate that drives graph adaptation.Source: §1 · p. 1
  • ADAPTGNS uses 20% fewer edges and lower MSE@200 than fixed k = 5 on WATERDROP, while on SAND it has 12% lower MSE@200 than fixed k = 5 using 10% fewer edges than k = 10.Against a fixed GNS trained on NLL loss, the adaptive policy also achieves 21% fewer edges and 8% lower MSE@200.Source: §1 · pp. 1–2

2. Dataset

The study evaluates adaptive graphs on granular SAND and fluid WATERDROP simulations. Their graph dynamics differ over time: WATERDROP changes around impact and disperses, while SAND settles into a denser heap.

2. Dataset

  • SAND trajectories come from a Taichi-MPM solver, while WATERDROP trajectories come from an SPH solver, providing granular-flow and fluid regimes for evaluation.The datasets are sourced from Kumar & Vantassel (2023) and Sanchez-Gonzalez et al. (2020), respectively.Source: §2 · p. 2
  • At the standard interaction radius r = 0.015, both datasets produce sparse graphs with moderate average degree.The graph statistics are measured at Step 100 and averaged across trajectories.Source: §2 · p. 2
  • In WATERDROP, edge count and average degree rise at impact and taper as particles disperse; in SAND, the graph slumps around Step 50 before recovering to a stable plateau as particles settle.The paper identifies these temporal signatures as central to the performance gap between datasets.Source: §2 · p. 2

3. Related Work

Prior particle simulators commonly use fixed geometric interaction graphs, while adaptive graph methods and non-graph approaches address related challenges through different mechanisms. ADAPTGNS instead uses predicted uncertainty to adapt point-cloud particle connectivity.

3. Related Work

  • Radius-based GNS interaction graphs are fixed by geometric rules, which do not account for local dynamics; fixed sparse graphs may also restrict information exchange across bottlenecks.The paper relates this concern to over-squashing and possible interaction with compounding error in autoregressive simulation.Source: §3 · p. 2
  • NRI infers a discrete graph from particle histories but keeps it fixed at inference and incurs overhead, while AdaMeshNet rewires mesh-based GNNs using geometric heuristics.ADAPTGNS targets point-cloud particle systems, uses no separate encoder, and adapts topology with its own predicted uncertainty.Source: §3 · p. 2
  • Other long-horizon approaches use physics priors, particle-grid emulation, or global attention rather than adapting a particle interaction graph.The paper describes ADAPTGNS as orthogonal to these directions and potentially integrable with stronger message functions or physics-informed losses.Source: §3 · p. 2

4. Methodology

ADAPTGNS adds a per-particle variance head and uses lagged uncertainty to adapt graph topology in a single inference pass. Joint Gaussian NLL training encourages uncertainty estimates to concentrate where prediction error is large.

4. Methodology

  • A variance head predicts a scalar uncertainty for each particle alongside the acceleration output, adding fewer than 20,000 parameters to a standard GNS.The head is a two-layer MLP on final node embeddings with a softplus activation to enforce positivity.Source: §4 · p. 3
  • The acceleration and variance heads are trained jointly with a Gaussian NLL loss that rewards concentrating uncertainty where error is large and penalizes inflating it elsewhere.The acceleration dimension is d = 2.Source: §4 · p. 3
  • At each timestep, the adaptive loop builds the graph using the previous step’s uncertainty estimates, avoiding a second GNN forward pass.The forward pass produces the next position and a fresh uncertainty estimate for the following step.Source: §4 · p. 3
  • The rule constructs the standard graph, marks particles above the per-step τ-th percentile of uncertainty, and adds edges within a larger radius when an edge is incident to a high-uncertainty particle.The expanded radius is rlarge = rad_factor · r; the selected settings were r = 0.015, rad_factor = 1.267, and τ = 70.Source: §4 · p. 3

5. Experiments and Results

ADAPTGNS is compared with fixed-graph baselines and ablations on SAND and WATERDROP. It improves on the same NLL-trained fixed-graph model on both datasets, with stronger results on WATERDROP; the discussion links the difference to spatially concentrated complexity and notes limitations of NLL and edge count as a compute proxy.

5.1. Quantitative Results

  • The predicted uncertainty acts as a rank-based signal for physical structure, although it is not perfectly calibrated in absolute value.On SAND, bins 0–7 have |∆| ≤0.06; total ECE is large because error is severely underestimated for the hardest particles, bins 8 and 9.Source: §5.1 · p. 3, Figure 2 · p. 4
  • On SAND, ADAPTGNS reduces edges by 21% and lowers MSE@200 by 8% (0.0417 vs. 0.0382) against its fixed-graph NLL counterpart.At matched edge counts (∼23,000 each), the MSE-trained model has MSE@200 of 0.0355 vs. 0.0382.Source: §5.1 · p. 4
  • On WATERDROP, ADAPTGNS uses 20% fewer edges than k = 5 (7,821 vs. 9,741) and achieves lower MSE@200 (0.0412 vs. 0.0414).Compared with k = 10, it uses 36% fewer edges for a 6% relative MSE@200 increase.Source: §5.1 · p. 4
  • ADAPTGNS strictly beats the same NLL model without adaptive graphs on both datasets, but does not defeat the best MSE-trained models on SAND.The comparisons include fixed-graph GNS models with k = 5 and k = 10, plus two ablations; the matched-edge ablation was not conducted for WATERDROP.Source: §5.1 · pp. 3–4, §5 · p. 3

5.2. Analysis and Discussion

  • The authors attribute WATERDROP’s stronger results to dispersing particles and more spatially heterogeneous uncertainty, whereas SAND particles undergo more similar dynamics.They note that NLL can be harder than MSE for raw acceleration accuracy, and that edge count is a noisy compute proxy because worse models can drift into denser configurations.Source: §5.2 · p. 4

6. Limitations and Future Work

The paper identifies runtime measurement, reaction lag, and limited evaluation scope as restrictions on interpreting its results. It calls for broader evaluation and future methods that address these constraints.

6. Limitations and Future Work

  • Average edge count tracks the principal cost of a GNN forward pass, but it is not a complete runtime metric.Adaptive graph construction also searches an expanded radius for selected particles, and its overhead remains to be quantified.Source: §6 · p. 4, WATERDROP · p. 5

WATERDROP

  • Figure 3 tracks per-step edge count and MSE on SAND and WATERDROP.The table reports results over 320 SAND steps and 1000 WATERDROP steps.Source: Figure 3 · p. 5, Table 4 · p. 5
  • Using the previous step’s uncertainty to construct the current graph avoids a second forward pass but may react too late during abrupt transformations.A next-step predictor or a physical prior based on acceleration or velocity divergence could reduce this lag.Source: WATERDROP · p. 5
  • The experiments cover one SPH fluid dataset and one MPM granular dataset; broader evaluation across materials, 3D contexts, and hybrid rigid–fluid scenes is needed.The results table covers SAND at 320 steps and WATERDROP at 1000 steps.Source: WATERDROP · p. 5, Table 4 · p. 5
  • Potential extensions include molecular dynamics, mesh-based PDE solvers, and robotics, where local regions may benefit from larger neighborhoods or adaptive refinement.The paper presents these as possible extensions, not evaluated settings.Source: WATERDROP · p. 5

7. Conclusion

The paper introduces a single-pass adaptive interaction graph driven by a jointly trained per-particle variance head. It reports that the signal corresponds to physical uncertainty and that AdaptGNS achieves modest but notable results on WATERDROP and SAND.

7. Conclusion

  • AdaptGNS uses a per-particle variance head, trained jointly with the acceleration head under NLL loss, to adapt the interaction graph.The adaptive mechanism avoids the cost of two GNN computations per step.Source: §7 · p. 5
  • The predicted variance corresponds to physical uncertainty, and AdaptGNS achieves modest but notable results on WATERDROP and SAND.Source: §7 · p. 5
  • Future work should examine variance-aware learning objectives, next-step uncertainty prediction, learnable expansion radii, broader physical regimes, longer horizons, and multi-seed evaluations.The paper says these are needed to examine the general usefulness of adaptive interaction graphs more carefully.Source: §7 · p. 5

Impact Statement

The paper identifies downstream uses of faster and more accurate simulation systems as potential societal consequences and reports no additional risks beyond those typically associated with machine learning research.

Impact Statement

  • Potential societal consequences are tied to downstream uses of faster and more accurate simulation systems; no additional risks beyond typical machine-learning research risks are identified.Source: Impact Statement · p. 5
Loading 2609.30822v1…