Data-driven construction of machine-learning-based interatomic potentials for gas-surface scattering dynamics: the case of NO on graphite

Samuel Del Fré; Gilberto A. Alou Angulo; Maurice Monnerville; Alejandro Rivero Santamaría

physics.chem-phMar 19, 2026

Data-driven construction of machine-learning-based interatomic potentials for gas-surface scattering dynamics: the case of NO on graphite

Authors: Samuel Del Fré, Gilberto A. Alou Angulo, Maurice Monnerville, Alejandro Rivero Santamaría

Organizations: Univ. Lille, CNRS, UMR 8523 – PhLAM – Physique des Lasers Atomes et Molécules, F-59000 Lille, France

Abstract

Accurate atomistic simulations of gas-surface scattering require potential energy surfaces that remain reliable over broad configurational and energetic ranges while retaining the efficiency needed for extensive trajectory sampling. Here, we develop a data-driven workflow for constructing a machine-learning interatomic potential (MLIP) tailored to gas-surface scattering dynamics, using nitric oxide (NO) scattering from highly oriented pyrolytic graphite (HOPG) as a benchmark system. Starting from an initial ab initio molecular dynamics (AIMD) dataset, local atomic environments are described by SOAP descriptors and analyzed in a reduced feature space obtained through principal component analysis. Farthest point sampling is then used to build a compact training set, and the resulting Deep Potential model is refined through a query-by-committee active-learning strategy using additional configurations extracted from molecular dynamics simulations over extended ranges of incident energies and surface temperatures. The final MLIP reproduces reference energies and forces with high fidelity and enables large-scale molecular dynamics simulations of NO scattering from graphite at a computational cost far below that of AIMD. The simulations provide detailed insight into adsorption energetics, trapping versus direct scattering probabilities, translational energy loss, angular distributions, and rotational excitation. Overall, the results reproduce the main experimental trends and demonstrate that descriptor-guided sampling combined with active learning offers an efficient and transferable strategy for constructing MLIPs for gas-surface interactions.

Figures & tables

Figure 1: PCA projection of the SOAP descriptor space for the full AIMD dataset, with FPS-selected environments highlighted as blue triangles (NO) and dark green dots (graphite). The inset shows the coverage evolution.

Figure 2: Energy (upper panel) and force (lower panel) distributions for the full AIMD dataset, the FPS-sampled dataset A, the QBC-selected configurations, and the final refined dataset B. Force magnitudes are shown on a logarithmic scale.

Panel of Figure 3 — Figure 3: Three-dimensional representation of the $N_{\mathrm{st}}^{\mathrm{QBC}}$ configurations selected during the active learning phase as a function of surface temperature $T_{\mathrm{surf}}$ and incident energy $E_{\mathrm{inc}}$ (upper panel). Distribution of the QBC model-deviation score ( $\Delta F$ ) over the extracted configurations on a logarithmic scale, with the selected uncertainty window highlighted (lower panel).

Figure 4: Parity plots between MLIP predictions of energies (upper panel) and forces (lower panel) and the corresponding DFT reference data for validation set B)

Figure 5: Scattering probability

P_{\mathrm{scat}}

obtained from molecular dynamics simulations. Top: dependence on the incident energy

E_{\mathrm{inc}}

for

T_{\mathrm{surf}}=100

K. Bottom: dependence on the surface temperature

T_{\mathrm{surf}}

for

E_{\mathrm{inc}}=0.1

eV.

Figure 6: Top: Mean kinetic energy loss of scattered NO as a function of incident energy

E_{\mathrm{inc}}

. Error bars denote one standard deviation. Bottom: Corresponding speed distributions of the scattered molecules. Solid lines show fits to a three-dimensional flux distribution [ 21 ] . The analytical fitted function and its parameters are given in table 1 . All data correspond to scattering from graphite at

T_{\mathrm{surf}}

= 100 K for incident energies between 0.05 and 2.0 eV.

$E_{\mathrm{inc}}$ (eV)	$A$	$c_{0}$	$\alpha$	$v_{\mathrm{peak}}$
		(m s -1 )	(m.s -1 )	(m.s -1 )
0.05	$2.9\!\times\!10^{-10}$	1.7	284.5	349.4
0.10	$1.6\!\times\!10^{-10}$	3.2	329.5	405.1
0.30	$1.7\!\times\!10^{-11}$	275.2	380.1	623.0
0.50	$5.1\!\times\!10^{-12}$	543.2	380.4	810.8
1.00	$1.3\!\times\!10^{-12}$	977.1	365.1	1150.8

Table 1: Parameters obtained from the three-dimensional flux distribution fits of the scattered NO velocity distributions at

T_{\mathrm{surf}}=100

K for different incident energies

E_{\mathrm{inc}}

. The fitted function is

F(c)=A\,c^{3}\exp[-(c-c_{0})^{2}/\alpha^{2}]

[ 21 ] .

v_{\text{peak}}

corresponds to the velocity at which the flux distribution reaches its maximum, i.e., the most probable speed of the scattered molecules in the flux-weighted velocity distribution.

Figure 7: Evolution of the angular scattering exponent

n

as a function of the incident translational energy

E_{\mathrm{inc}}

for

T_{\mathrm{surf}}

= 100 K. The exponent

n

is obtained by fitting the polar angular distribution of scattered trajectories with a

\cos^{n}(\theta_{\mathrm{scatt}})

function.

Figure 8: Rotational state distributions of scattered NO for incident energies between 0.05 and 2.0 eV at a surface temperature of 100 K. Maxwell–Boltzmann fits are shown as discrete markers superimposed on the distributions.

Figure 9: Rotational temperature of scattered NO as a function of surface temperature

T_{\mathrm{surf}}

E_{\mathrm{inc}}~=~0.1

eV. The dashed line corresponds to

T_{\mathrm{rot}}=T_{\mathrm{surf}}

, indicating the limit of complete rotational accommodation.

Machine Learning Interatomic Potentials

Machine learning interatomic potentials predict atomic energies and forces for fast simulations of molecules and materials.

See all 91 papers

Universal and efficient graph neural networks with dynamic attention for machine learning interatomic potentials
Mar 24, 2026Shuyu Bi, Zhede Zhao, Qiangchao Sun +3Machine Learning Interatomic Potentials Graph Neural Networks
School of Materials Science and Engineering & State Key Laboratory of Advanced Refractories & State Key Laboratory of Advanced Special Steels, Shanghai University, Shanghai, 200444, P. R. China.
Beyond Adam: SOAP and Muon for Faster, Label-Efficient Training of Machine Learning Interatomic Potentials
Jul 2, 2026Gil Harari, Yoel Zimmermann, Ola Tangen Kulseng +4Machine Learning Interatomic Potentials Neural Network Optimization
John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA · Robert Bosch LLC Research and Technology Center, Watertown, MA, USA
Physics-Aligned Electronic Ground-State Learning Improves Generalization
Oct 7, 2026Eike S. Eberhard, Xaver Kainz, Viktor Kotsev +2Physics-Informed ML OOD Generalization
Technical University of Munich (TUM) · Munich Data Science Institute (MDSI) · Munich Center for Machine Learning (MCML) +1

Mar 24, 2026J/K move · Enter open · S save

Universal and efficient graph neural networks with dynamic attention for machine learning interatomic potentials

Shuyu Bi, Zhede Zhao, Qiangchao Sun, Tao Hu+2

School of Materials Science and Engineering & State Key Laboratory of Advanced Refractories & State Key Laboratory of Advanced Special Steels, Shanghai University, Shanghai, 200444, P. R. China.

Machine Learning Interatomic Potentials Graph Neural Networks Molecular Dynamics Simulation Graph Attention Networks

Open paper AskPDF arXiv

The core of molecular dynamics simulation fundamentally lies in the interatomic potential. Traditional empirical potentials lack accuracy, while first-principles methods are computationally prohibitive. Machine learning interatomic potentials (MLIPs) promise near-quantum accuracy at linear cost, but existing models still face challenges in efficiency and stability. We presents Machine Learning Advances Neural Network (MLANet), an efficient and robust graph neural network framework. MLANet introduces a dual-path dynamic attention mechanism for geometry-aware message passing and a multi-perspective pooling strategy to construct comprehensive system representations. This design enables highly accurate modeling of atomic environments while achieving exceptional computational efficiency, making high-fidelity simulations more accessible. Tested across a wide range of datasets spanning diverse systems, including organic molecules (e.g., QM7, MD17), periodic inorganic materials (e.g., Li-containing crystals), two-dimensional materials (e.g., bilayer graphene, black phosphorus), surface catalytic reactions (e.g., formate decomposition), and charged systems, MLANet maintains competitive prediction accuracy while its computational cost is markedly lower than mainstream equivariant models, and it enables stable long-time molecular dynamics simulations. MLANet provides an efficient and practical tool for large-scale, high-accuracy atomic simulations.

Abstract

Machine Learning Interatomic Potentials

Explore similar work

Universal and efficient graph neural networks with dynamic attention for machine learning interatomic potentials

Beyond Adam: SOAP and Muon for Faster, Label-Efficient Training of Machine Learning Interatomic Potentials

Physics-Aligned Electronic Ground-State Learning Improves Generalization