Dissociative Chemisorption of O2 on Al(111): Dynamics on a Correlated Wave-Function-Based Potential Energy Surface

Dissociative chemisorption of O2 on Al(111): dynamics on a potential energy surface computed with a non-self-consistent screened hybrid density functional approach

Density functional theory (DFT) at the generalized gradient approximation (GGA) level is often considered the best compromise between feasibility and accuracy for reactions of molecules on metal surfaces. Recent work, however, strongly suggests that density functionals (DFs) based on GGA exchange are not able to describe molecule-metal surface reactions for which the work function of the metal surface minus the electron affinity of the molecule is less than 7 eV. Systems for which this is true exhibit an increased charge transfer from the metal to the molecule at the transition state, increasing the delocalisation of the electron density. This enlarged delocalisation can cause GGA-DFT to underestimate energy values relative to the gas-phase and thus underestimate the barrier height, similar to what has been observed for several gas-phase reactions. An example of such a molecule-metal surface system is O2 + Al(111). Following a similar strategy as for gas-phase reactions, previous work showed results of increased accuracy when using a screened hybrid DF for O2 + Al(111). However, even screened hybrid DFs are computationally expensive to use for metal surfaces. To resolve this, we test a non-self-consistent field (NSCF) screened hybrid DF approach. This approach computes screened hybrid DFT energies based on self-consistent-field (SCF) GGA electronic densities. Here, we explore the accuracy of the NSCF screened hybrid DF approach by implementing the NSCF HSE03-1/3x@RPBE DF for O2 + Al(111). We compute and analyse molecular beam sticking probabilities as well as a set of sticking probabilities for rotationally aligned O2. Our results show that the NSCF approach results in reaction probability curves that reproduce SCF results with near-chemical accuracy, suggesting that the NSCF approach can be used advantageously for exploratory purposes. An analysis of the potential energy surface and the barriers gives insight into the cause of the disagreement between the SCF and NSCF reaction probabilities and into the changes needed in theoretical modelling to further improve the description of the O2 + Al(111) system. Finally, the hole model yields fair agreement with dynamics results for the reaction probability curve, but results in an increased slope of the reaction probability curve compared to the molecular dynamics, with a shift to lower or higher energies depending on whether the vibrational energy of the molecule is included in the initial energy of the molecule or not.

On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K

We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.

Temperature-Dependent Kinetics for the Reactions of Fen- (n = 2-17) and FexNiy- (x + y = 3-9) with O2: Comparison of Pure and Mixed Metal Clusters with Relevance to Meteor Radio Afterglows and Surface Oxidation

Rate constants and product branching fractions were measured from 300-600 K for Fen- + O2 (n = 2-17) and for 300-500 K for FexNiy- + O2 (x + y = 3-9) using a selected-ion flow tube (SIFT) apparatus. Rate constants for 46 species are reported. All rate constants increased with increasing temperature, and several were in excess of the Langevin-Gioumousis-Stevenson (LGS) capture rate at elevated temperatures. As with previously studied transition metal anion oxidation reactions, the collision limit is treated as the sum of the LGS limit along with a hard-sphere contribution, allowing for determination of activation energies. These values are compared to each other along with previous results for Nin-. Measured rate constants for all three series (Fen-, Nin-, and FexNy-) vary over a relatively narrow range (1-5 × 10-10 cm3 s-1 at 300 K) being at least 15% of the collision rate constant. All reaction rate constants increase with temperature, described by small activation energies of 0.5-4 kJ mol-1. The data are consistent with an anticorrelation between the electron binding energy and rate constant, previously noted in other systems. The Fen- reaction produces a larger population of higher energy electrons than do the Nin- reactions, with FexNiy- producing an intermediate amount. The results suggest that the overall rate constant is limited by a small energetic barrier located at a large internuclear distance where electrostatic forces dominate, causing the potentials to be similar across systems, while the product formation is determined by the shorter-range, valence portion of the potential, which varies widely between systems.

Best-of-Both-Worlds Predictive Approach to Dissociative Chemisorption on Metals

Predictive capability, accuracy, and affordability are essential features of a theory that is capable of describing dissociative chemisorption on a metal surface. This type of reaction is important for heterogeneous catalysis. Here we present an approach in which we use diffusion Monte Carlo (DMC) to pin the minimum barrier height and construct a density functional that reproduces this value. This predictive approach allows the construction of a potential energy surface at the cost of density functional theory while retaining near DMC accuracy. Scrutinizing effects of energy dissipation and quantum tunneling, dynamics calculations suggest the approach to be of near chemical accuracy, reproducing molecular beam sticking experiments for the showcase H2 + Al(110) system to ∼1.4 kcal/mol.

Constructing Mixed Density Functionals for Describing Dissociative Chemisorption on Metal Surfaces: Basic Principles

The production of a majority of chemicals involves heterogeneous catalysis at some stage, and the rates of many heterogeneously catalyzed processes are governed by transition states for dissociative chemisorption on metals. Accurate values of barrier heights for dissociative chemisorption on metals are therefore important to benchmarking electronic structure theory in general and density functionals in particular. Such accurate barriers can be obtained using the semiempirical specific reaction parameter (SRP) approach to density functional theory. However, this approach has thus far been rather ad hoc in its choice of the generic expression of the SRP functional to be used, and there is a need for better heuristic approaches to determining the mixing parameters contained in such expressions. Here we address these two issues. We investigate the ability of several mixed, parametrized density functional expressions combining exchange at the generalized gradient approximation (GGA) level with either GGA or nonlocal correlation to reproduce barrier heights for dissociative chemisorption on metal surfaces. For this, seven expressions of such mixed density functionals are tested on a database consisting of results for 16 systems taken from a recently published slightly larger database called SBH17. Three expressions are derived that exhibit high tunability and use correlation functionals that are either of the PBE GGA form or of one of two limiting nonlocal forms also describing the attractive van der Waals interaction in an approximate way. We also find that, for mixed density functionals incorporating GGA correlation, the optimum fraction of repulsive RPBE GGA exchange obtained with a specific GGA density functional is correlated with the charge-transfer parameter, which is equal to the difference in the work function of the metal surface and the electron affinity of the molecule. However, the correlation is generally not large and not large enough to obtain accurate guesses of the mixing parameter for the systems considered, suggesting that it does not give rise to a very effective search strategy.

Spin-dependent reactivity and spin-flipping dynamics in oxygen atom scattering from graphite

The formation of two-electron chemical bonds requires the alignment of spins. Hence, it is well established for gas-phase reactions that changing a molecule's electronic spin state can dramatically alter its reactivity. For reactions occurring at surfaces, which are of great interest during, among other processes, heterogeneous catalysis, there is an absence of definitive state-to-state experiments capable of observing spin conservation and therefore the role of electronic spin in surface chemistry remains controversial. Here we use an incoming/outgoing correlation ion imaging technique to perform scattering experiments for O(3P) and O(1D) atoms colliding with a graphite surface, in which the initial spin-state distribution is controlled and the final spin states determined. We demonstrate that O(1D) is more reactive with graphite than O(3P). We also identify electronically nonadiabatic pathways whereby incident O(1D) is quenched to O(3P), which departs from the surface. With the help of molecular dynamics simulations carried out on high-dimensional machine-learning-assisted first-principles potential energy surfaces, we obtain a mechanistic understanding for this system: spin-forbidden transitions do occur, but with low probabilities.

Nonadiabatic Potential Energy Surfaces for a Molecule on a Surface as Found by Constrained Complete Active Space Theory

In order to study electron-transfer mediated chemical processes on a metal surface, one requires not one but two potential energy surfaces (one ground state and one excited state) as in Marcus theory. In this letter, we report that a novel, dynamically weighted, state-averaged constrained CASSCF(2,2) (DW-SA-cCASSCF(2,2)) can produce such surfaces for the Anderson impurity model. Both ground and excited state potentials are smooth, they incorporate states with a charge transfer character, and the accuracy of the ground state surface can be verified for some model problems by renormalization group theory. Future development of gradients and nonadiabatic derivative couplings should allow for the study of nonadiabatic dynamics for molecules near metal surfaces.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

SBH17: Benchmark Database of Barrier Heights for Dissociative Chemisorption on Transition Metal Surfaces

Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.

Stereodynamics of adiabatic and non-adiabatic energy transfer in a molecule surface encounter

Molecular energy transfer and reactions at solid surfaces depend on the molecular orientation relative to the surface. While such steric effects have been largely understood in electronically adiabatic processes, the orientation-dependent energy transfer in NO scattering from Au(111) was complicated by electron-mediated nonadiabatic effects, thus lacking a clear interpretation and posing a great challenge for theories. Herein, we investigate the stereodynamics of adiabatic and nonadiabatic energy transfer via molecular dynamics simulations of NO(v = 3) scattering from Au(111) using realistic initial orientation distributions based on accurate neural network fitted adiabatic potential energy surface and electronic friction tensor. Our results reproduce the observed stronger vibrational relaxation for N-first orientation and enhanced rotational rainbow for O-first orientation, and demonstrate how adiabatic anisotropic interactions steer molecules into the more attractive N-first orientation to experience more significant energy transfer. Remaining disagreements with experiment suggest the direction for further developments of nonadiabatic theories for gas-surface scattering.

Capture Theory Models: An overview of their development, experimental verification, and applications to ion-molecule reactions

Since Arrhenius first proposed an equation to account for the behaviour of thermally activated reactions in 1889, significant progress has been made in our understanding of chemical reactivity. A number of capture theory models have been developed over the past several decades to predict the rate coefficients for reactions between ions and molecules-ranging from the Langevin equation (for reactions between ions and non-polar molecules) to more recent fully quantum theories (for reactions at ultra-cold temperatures). A number of different capture theory methods are discussed, with the key assumptions underpinning each approach clearly set out. The strengths and limitations of these capture theory methods are examined through detailed comparisons between low-temperature experimental measurements and capture theory predictions. Guidance is provided on the selection of an appropriate capture theory method for a given class of ion-molecule reaction and set of experimental conditions-identifying when a capture-based model is likely to provide an accurate prediction. Finally, the impact of capture theories on fields such as astrochemical modelling is noted, with some potential future directions of capture-based approaches outlined.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Structures and Electron Affinities of Aluminum Hydride Clusters AlnH (n = 3-13)

Low-energy structures and electron affinities (EAs) for aluminum hydride clusters Al_nH (n = 3-13) have been calculated using ab initio and density functional calculations. Geometries were optimized at the PBE0/def-2-TZVPP level of theory, which has been shown to match the currently accepted lowest-energy structures for the all-aluminum clusters Al_n and their anions. Neutral hydride clusters with n = 4, 7, and 9-12 are predicted to adopt terminal structures with the hydrogen atom bound to only one aluminum atom and with only minor alterations of the aluminum atom arrangement from that of the all-aluminum cluster. Clusters with n = 3 and 13 are predicted to adopt "face-centered" geometries, and the n = 6 cluster is predicted to prefer an isomer with the hydrogen atom bridging two aluminum atoms, also with little or no distortion to the aluminum atom arrangement from the all-aluminum cluster. Addition of a hydrogen atom to clusters with n = 5 and 8 is predicted to distort the aluminum atom arrangement significantly from that of the corresponding all-aluminum cluster. In the anionic clusters, terminal clusters are preferred for all cluster sizes except for n = 6 that prefers a face-centered arrangement. Minor distortions in the aluminum scaffolding for Al₁₁ and Al₁₂ were found, while all other anionic clusters adopt structures with little or no deviation in the aluminum atom arrangement from the corresponding all-aluminum cluster. Raw adiabatic electron affinities were computed using CCSD(T)/aug-cc-pVTZ single-point energies for the anionic and neutral hydride clusters at their respective DFT geometries. Isodesmic electron affinities for the hydride clusters were computed relative to their all-aluminum counterparts and show an even-odd alternation with cluster size. Derived EAs alternate in magnitude between even- and odd-numbered clusters, with the even-numbered clusters having relatively larger EAs.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Computational approaches to dissociative chemisorption on metals: towards chemical accuracy

We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (E_b, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate E_b can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol^-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for E_b for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH₄. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of E_b for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate E_b for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O₂, H₂O, NH₃, CO₂, and CH₃OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.

Normal and off-normal incidence dissociative dynamics of O2(v,J) on ultrathin Cu films grown on Ru(0001)

The dissociative adsorption of molecular oxygen on metal surfaces has long been controversial, mostly due to the spin-triplet nature of its ground state, to possible non-adiabatic effects, such as an abrupt charge transfer from the metal to the molecule, or even to the role played by the surface electronic state. Here, we have studied the dissociative adsorption of O₂ on Cu_ML/Ru(0001) at normal and off-normal incidence, from thermal to super-thermal energies, using quasi-classical dynamics, in the framework of the generalized Langevin oscillator model, and density functional theory based on a multidimensional potential energy surface. Our simulations reveal a rather intriguing behavior of dissociative adsorption probabilities, which exhibit normal energy scaling at incidence energies below the reaction barriers and total energy scaling above, irrespective of the reaction channel, either direct dissociation, trapping dissociation, or molecular adsorption. We directly compare our results with existing scanning tunneling spectroscopy and microscopy measurements. From this comparison, we infer that the observed experimental behavior at thermal energies may be due to ligand and strain effects, as already found for super-thermal incidence energies.

Gas-Phase Anionic Metal Clusters are Model Systems for Surface Oxidation: Kinetics of the Reactions of Mn- with O2 (M = V, Cr, Co, Ni; n = 1-15)

The reactions of anionic metal clusters M_n^- with O₂ (M = V (n = 1-15), Cr (n = 1-15), Co (n = 1-12), and Ni (n = 1-14)) are investigated from 300 to 600 K using a selected-ion flow tube. All rate constants show a positive temperature dependence, well described by an Arrhenius equation. Rate constants exceed (or are extrapolated to exceed at higher temperatures) the Langevin-Gioumousis-Stevenson capture rate constant. Application of a capture model accounting for the finite size of the clusters reproduces the size-dependent trends in reactivity. The assumption that reactivity is further controlled by an energetic barrier early in the reaction coordinate is consistent with the experimental observations. An observed correlation of the derived barrier heights on the electron binding energy of M_n^- suggests the barrier may be formed at an avoided crossing between electronic states correlating to M_n^- + O₂ and M_n + O₂^- reactants, analogous to that previously proposed for Al_n^- + O₂ systems. The mechanism is analogous to that for reactions of O₂ with neutral metal surfaces, indicating that gas-phase reactions of anionic metal clusters can be an appropriate model systems for surface oxidation.

Nonadiabatic Dynamics at Metal Surfaces: Fewest Switches Surface Hopping with Electronic Relaxation

A new scheme is proposed for modeling molecular nonadiabatic dynamics near metal surfaces. The charge-transfer character of such dynamics is exploited to construct an efficient reduced representation for the electronic structure. In this representation, the fewest switches surface hopping (FSSH) approach can be naturally modified to include electronic relaxation (ER). The resulting FSSH-ER method is valid across a wide range of coupling strengths as supported by tests applied to the Anderson-Holstein model for electron transfer. Future work will combine this scheme with ab initio electronic structure calculations.

Density Functional Theory for Molecule-Metal Surface Reactions: When Does the Generalized Gradient Approximation Get It Right, and What to Do If It Does Not

While density functional theory (DFT) is perhaps the most used electronic structure theory in chemistry, many of its practical aspects remain poorly understood. For instance, DFT at the generalized gradient approximation (GGA) tends to fail miserably at describing gas-phase reaction barriers, while it performs surprisingly well for many molecule-metal surface reactions. GGA-DFT also fails for many systems in the latter category, and up to now it has not been clear when one may expect it to work. We show that GGA-DFT tends to work if the difference between the work function of the metal and the molecule's electron affinity is greater than ∼7 eV and to fail if this difference is smaller, with sticking of O₂ on Al(111) being a spectacular example. Using dynamics calculations we show that, for this system, the DFT problem may be solved as done for gas-phase reactions, i.e., by resorting to hybrid functionals, but using screening at long-range to obtain a correct description of the metal. Our results suggest the GGA error in the O₂ + Al(111) barrier height to be functional driven. Our results also suggest the possibility to compute potential energy surfaces for the difficult-to-treat systems with computationally cheap nonself-consistent calculations in which a hybrid functional is applied to a GGA density.

High-Fidelity Potential Energy Surfaces for Gas-Phase and Gas-Surface Scattering Processes from Machine Learning

In this Perspective, we review recent advances in constructing high-fidelity potential energy surfaces (PESs) from discrete ab initio points, using machine learning tools. Such PESs, albeit with substantial initial investments, provide significantly higher efficiency than direct dynamics methods and/or high accuracy at a level that is not affordable by on-the-fly approaches. These PESs not only are a necessity for quantum dynamical studies because of delocalization of wave packets but also enable the study of low-probability and long-time events in (quasi-)classical treatments. Our focus here is on inelastic and reactive scattering processes, which are more challenging than bound systems because of the involvement of continua. Relevant applications and developments for dynamical processes in both the gas phase and at gas-surface interfaces are discussed.

Redefining the Mechanism of O2 Etching of Aln- Superatoms: An Early Barrier Controls Reactivity, Analogous to Surface Oxidation

New insights into aluminum anion cluster reactivity with O₂ were obtained through temperature-dependent kinetics measurements. Overall reactivity is controlled by a barrier at an avoided crossing where charge is transferred from the cluster to the O₂, mechanistically similar to what occurs as O₂ approaches a bulk Al surface. Contrary to prior interpretations, spin conservation does not inhibit the reaction of clusters with an odd number of Al atoms. In fact, the only spin constraint in these systems is on the reactivity of even clusters due to repulsive surfaces, not previously recognized. Although the superatom nature of Al₁₃^- is manifest in its high electron binding energy (EBE), the mechanism of its reactivity is not special; it reacts with O₂ as if it were a small piece of bulk Al. These experiments highlight the sensitivity of Al cluster reactivity with O₂ to temperature and EBE, uncovering routes to industrial scale use of aluminum superatom-based materials.

Unified and transferable description of dynamics of H2 dissociative adsorption on multiple copper surfaces via machine learning

Schematic of the developed neural network potential energy surface enabling a unified and transferable description of dynamics of H₂ dissociative adsorption on multiple copper surfaces.

Dynamics in reactions on metal surfaces: A theoretical perspective

Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.

Reactive and Nonreactive Scattering of HCl from Au(111): An Ab Initio Molecular Dynamics Study

The HCl + Au(111) system has recently become a benchmark for highly activated dissociative chemisorption, which presumably is strongly affected by electron-hole pair excitation. Previous dynamics calculations, which were based on density functional theory at the generalized gradient approximation level (GGA-DFT) for the molecule-surface interaction, have all overestimated measured reaction probabilities by at least an order of magnitude. Here, we perform ab initio molecular dynamics (AIMD) and AIMD with electronic friction (AIMDEF) calculations employing a density functional that includes the attractive van der Waals interaction. Our calculations model the simultaneous and possibly synergistic effects of surface temperature, surface atom motion, electron-hole pair excitation, the molecular beam conditions of the experiments, and the van der Waals interaction on the reactivity. We find that reaction probabilities computed with AIMDEF and the SRP32-vdW functional still overestimate the measured reaction probabilities, by a factor 18 for the highest incidence energy at which measurements were performed (≈2.5 eV). Even granting that the experiment could have underestimated the sticking probability by about a factor three, this still translates into a considerable overestimation of the reactivity by the current theory. Likewise, scaled transition probabilities for vibrational excitation from ν = 1, j = 1 to ν = 2 are overestimated by the AIMDEF theory, by factors 3-8 depending on the initial conditions modeled. Energy losses to the surface and translational energy losses are, however, in good agreement with experimental values.