Physically-motivated Recursively Embedded Atom Neural Networks: Incorporating Local Completeness and Nonlocality

Recent advances in machine-learned interatomic potentials largely benefit from the atomistic representation and locally invariant many-body descriptors. It was however recently argued that including three- (or even four-) body features is incomplete to distinguish specific local structures. Utilizing an embedded density descriptor made by linear combinations of neighboring atomic orbitals and realizing that each orbital coefficient physically depends on its own local environment, we propose a recursively embedded atom neural network model. We formally prove that this model can efficiently incorporate complete many-body correlations without explicitly computing high-order terms. This model not only successfully addresses challenges regarding local completeness and nonlocality in representative systems, but also provides an easy and general way to update local many-body descriptors to have a message-passing form without changing their basic structures.

Determining the effect of hot electron dissipation on molecular scattering experiments at metal surfaces

Nonadiabatic effects that arise from the concerted motion of electrons and atoms at comparable energy and time scales are omnipresent in thermal and light-driven chemistry at metal surfaces. Excited (hot) electrons can measurably affect molecule-metal reactions by contributing to state-dependent reaction probabilities. Vibrational state-to-state scattering of NO on Au(111) has been one of the most studied examples in this regard, providing a testing ground for developing various nonadiabatic theories. This system is often cited as the prime example for the failure of electronic friction theory, a very efficient model accounting for dissipative forces on metal-adsorbed molecules due to the creation of hot electrons in the metal. However, the exact failings compared to experiment and their origin from theory are not established for any system, because dynamic properties are affected by many compounding simulation errors of which the quality of nonadiabatic treatment is just one. We use a high-dimensional machine learning representation of electronic structure theory to minimize errors that arise from quantum chemistry. This allows us to perform a comprehensive quantitative analysis of the performance of nonadiabatic molecular dynamics in describing vibrational state-to-state scattering of NO on Au(111) and compare directly to adiabatic results. We find that electronic friction theory accurately predicts elastic and single-quantum energy loss, but underestimates multi-quantum energy loss and overestimates molecular trapping at high vibrational excitation. Our analysis reveals that multi-quantum energy loss can potentially be remedied within friction theory, whereas the overestimation of trapping constitutes a genuine breakdown of electronic friction theory. Addressing this overestimation for dynamic processes in catalysis and surface chemistry will require more sophisticated theories.

Accelerating Atomistic Simulations with Piecewise Machine Learned Ab Initio Potentials at Classical Force Field-like Cost

Machine learning methods have nowadays become easy-to-use tools for constructing high-dimensional interatomic potentials with ab initio accuracy. Although machine learned interatomic potentials are generally orders of magnitude faster than first-principles calculations, they remain much slower than classical force fields, at the price of using more complex structural descriptors. To bridge this efficiency gap, we propose an embedded atom neural network approach with simple piecewise switching function based descriptors, resulting in a favorable linear scaling with the number of neighbor atoms. Numerical examples validate that this piecewise machine learning model can be over an order of magnitude faster than various popular machine learned potentials with comparable accuracy for both metallic and covalent materials, approaching the speed of the fastest embedded atom method (i.e. several {mu}s/atom per CPU core). The extreme efficiency of this approach promises its potential in first-principles atomistic simulations of very large systems and/or in long timescale.

Efficient and Accurate Simulations of Vibrational and Electronic Spectra with Symmetry-Preserving Neural Network Models for Tensorial Properties

Machine learning has revolutionized the high-dimensional representations for molecular properties such as potential energy. However, there are scarce machine learning models targeting tensorial properties, which are rotationally covariant. Here, we propose tensorial neural network (NN) models to learn both tensorial response and transition properties, in which atomic coordinate vectors are multiplied with scalar NN outputs or their derivatives to preserve the rotationally covariant symmetry. This strategy keeps structural descriptors symmetry invariant so that the resulting tensorial NN models are as efficient as their scalar counterparts. We validate the performance and universality of this approach by learning response properties of water oligomers and liquid water, and transition dipole moment of a model structural unit of proteins. Machine learned tensorial models have enabled efficient simulations of vibrational spectra of liquid water and ultraviolet spectra of realistic proteins, promising feasible and accurate spectroscopic simulations for biomolecules and materials.

Symmetry-Adapted High Dimensional Neural Network Representation of Electronic Friction Tensor of Adsorbates on Metals

Nonadiabatic effects in chemical reaction at metal surfaces, due to excitation of electron-hole pairs, stand at the frontier of the studies of gas-surface reaction dynamics. However, the first principles description of electronic excitation remains challenging. In an efficient molecular dynamics with electronic friction (MDEF) method, the nonadiabatic couplings are effectively included in a so-called electronic friction tensor (EFT), which can be computed from first-order time-dependent perturbation theory (TDPT) in terms of density functional theory (DFT) orbitals. This second-rank tensor depends on adsorbate position and features a complicated transformation with regard to the intrinsic symmetry operations of the system. In this work, we develop a new symmetry-adapted neural network representation of EFT, based on our recently proposed embedded atom neural network (EANN) framework. Inspired by the derivation of the nonadiabatic coupling matrix, we represent the tensorial friction by the first and second derivatives of multiple outputs of NNs with respect to atomic Cartesian coordinates. This rigorously preserves the positive semidefiniteness, directional property, and correct symmetry-equivariance of EFT. Unlike previous methods, our new approach can readily include both molecular and surface degrees of freedom, regardless of the type of surface. Tests on the H2+Ag(111) system show that this approach yields an accurate, efficient, and continuous representation of EFT, making it possible to perform large scale TDPT-based MDEF simulations to study both adiabatic and nonadiabatic energy dissipation in a unified framework.

Embedded Atom Neural Network Potentials: Efficient and Accurate Machine Learning with a Physically Inspired Representation

We propose a simple, but efficient and accurate machine learning (ML) model for developing high-dimensional potential energy surface. This so-called embedded atom neural network (EANN) approach is inspired by the well-known empirical embedded atom method (EAM) model used in condensed phase. It simply replaces the scalar embedded atom density in EAM with a Gaussian-type orbital based density vector, and represents the complex relationship between the embedded density vector and atomic energy by neural networks. We demonstrate that the EANN approach is equally accurate as several established ML models in representing both big molecular and extended periodic systems, yet with much fewer parameters and configurations. It is highly efficient as it implicitly contains the three-body information without an explicit sum of the conventional costly angular descriptors. With high accuracy and efficiency, EANN potentials can vastly accelerate molecular dynamics and spectroscopic simulations in complex systems at ab initio level.

Adiabatic and Nonadiabatic Energy Dissipation during Scattering of Vibrationally Excited CO from Au(111)

A high-dimensional potential energy surface (PES) for CO interaction with the Au(111) surface is developed using a machine-learning algorithm. Including both molecular and surface coordinates, this PES enables the simulation of the recent experiment on scattering of vibrationally excited CO from Au(111). Trapping in a physisorption well is observed to increase with decreasing incidence energy. While energy dissipation of physisorbed CO is slow, due to weak coupling with both the phonons and electron-hole pairs, its access to the chemisorption well facilitates fast vibrational relaxation of CO through nonadiabatic coupling with surface electron-hole pairs.