H + H2 quantum dynamics using potential energy surfaces based on the XYG3 type of doubly hybrid density functionals: Validation of the density functionals

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.

Bayesian machine learning for quantum molecular dynamics

This article discusses applications of Bayesian machine learning for quantum molecular dynamics.

The harpooning mechanism as evidenced in the oxidation reaction of the Al atom

The harpooning model is firstly evidenced by the maximum impact parameter derived from AlO(N_max) products with speed v′ ≈ 0 from the Al + O₂ reaction.

The harpooning mechanism has long been proposed for elementary reaction dynamics involving metals. It is characterized by an initial electron transfer (ET) process from the metal to the oxidant molecule. For the titled reaction Al + O₂, the ET distance can be predicted to be 2.6 Å by simply calculating the energy difference between the ionization energy of the Al atom and the electron affinity of the O₂ molecule. Hereby we experimentally derived the maximum impact parameter b_max of 2.5 ± 0.2 Å for the titled reaction, in consistency with the predicted ET distance. This derivation of b_max was achieved by using the crossed molecular beam experiment at a collision energy of 507 cm^–1 (i.e. 1.45 kcal mol^–1) with a high resolution time-sliced ion velocity imaging detection of the state-selective AlO products based on the (1 + 1) resonance-enhanced multiphoton ionization. The small rotational constant of the AlO(X²Σ⁺) radical (B_e = 0.6413 cm^–1) facilitated the formation of the AlO(v = 0) products in high rotational levels up to the energetically limited state, N_max = 52, with an almost zero velocity mapping. Hence, in this extreme angular momentum disposal case, the collisional orbital angular momentum l was nearly completely channeled into the product rotational angular momentum as a consequence of the conservations of energy and angular momentum, offering a reaction system that breaks the restriction of kinematically favored mass combination in order to obtain information on the impact parameters. The present study yields the first direct derivation of b_max from the maximum rotational level of products under the experimental condition with the recoil energy E′_T ≈ 0. This, in turn, provides solid evidence in supporting the harpooning mechanism.

Accurate potential energy surfaces for hydrogen abstraction reactions: A benchmark study on the XYG3 doubly hybrid density functional

The potential energy surface (PES) for the H + CH₄ system has been constructed with the recently developed XYG3 doubly hybrid functional, while those with the standard B3LYP hybrid functional, and the Møller-Plesset perturbation theory up to the second order (MP2) are also presented for comparison. Quantum dynamics studies demonstrated that satisfactory results on the reaction probabilities and the rate coefficients can be obtained on top of the XYG3-PES, as compared to the results based on the highly accurate, yet expensive, CCSD(T)-PES (Li et al., J. Chem. Phys. 2015, 142, 204302). Further investigation suggested that the XYG3 functional is useful in providing accurate rate coefficients for some larger systems involving H atom abstractions. © 2017 Wiley Periodicals, Inc.

When does a functional correctly describe both the structure and the energy of the transition state?

Requiring that several properties are well reproduced is a severe test on density functional approximations. This can be assessed through the estimation of joint and conditional success probabilities. An example is provided for a small set of molecules, for properties characterizing the transition states (geometries and energies).

Beyond energies: geometry predictions with the XYG3 type of doubly hybrid density functionals

The scaled mean absolute deviations (s-MADs) of the optimized geometric parameters for covalent bondings (the CCse set), nonbonded interactions (the S22G30 set) and the transition state structures (the TSG36 set), with Tot referring to the averaged s-MAD for general performances.