Ab initio molecular dynamics of hydrogen dissociation on metal surfaces using neural networks and novelty sampling

Ludwig, Jeffery; Vlachos, Dionisios G.
October 2007
Journal of Chemical Physics;10/21/2007, Vol. 127 Issue 15, p154716
Academic Journal
We outline a hybrid multiscale approach for the construction of ab initio potential energy surfaces (PESs) useful for performing six-dimensional (6D) classical or quantum mechanical molecular dynamics (MD) simulations of diatomic molecules reacting at single crystal surfaces. The algorithm implements concepts from the corrugation reduction procedure, which reduces energetic variation in the PES, and uses neural networks for interpolation of smoothed ab initio data. A novelty sampling scheme is implemented and used to identify configurations that are most likely to be predicted inaccurately by the neural network. This hybrid multiscale approach, which couples PES construction at the electronic structure level to MD simulations at the atomistic scale, reduces the number of density functional theory (DFT) calculations needed to specify an accurate PES. Due to the iterative nature of the novelty sampling algorithm, it is possible to obtain a quantitative measure of the convergence of the PES with respect to the number of ab initio calculations used to train the neural network. We demonstrate the algorithm by first applying it to two analytic potentials, which model the H2/Pt(111) and H2/Cu(111) systems. These potentials are of the corrugated London-Eyring-Polanyi-Sato form, which are based on DFT calculations, but are not globally accurate. After demonstrating the convergence of the PES using these simple potentials, we use DFT calculations directly and obtain converged semiclassical trajectories for the H2/Pt(111) system at the PW91/generalized gradient approximation level. We obtain a converged PES for a 6D hydrogen-surface dissociation reaction using novelty sampling coupled directly to DFT. These results, in excellent agreement with experiments and previous theoretical work, are compared to previous simulations in order to explore the sensitivity of the PES (and therefore MD) to the choice of exchange and correlation functional. Despite having a lower energetic corrugation in our PES, we obtain a broader reaction probability curve than previous simulations, which is attributed to increased geometric corrugation in the PES and the effect of nonparallel dissociation pathways.


Related Articles

  • Towards a specific reaction parameter density functional for reactive scattering of H2 from Pd(111). Boereboom, J. M.; Wijzenbroek, M.; Somers, M. F.; Kroes, G. J. // Journal of Chemical Physics;Dec2013, Vol. 139 Issue 24, p244707 

    Recently, an implementation of the specific reaction parameter (SRP) approach to density functional theory (DFT) was used to study several reactive scattering experiments of H2 on Cu(111). It was possible to obtain chemical accuracy (1 kcal/mol ≈ 4.2 kJ/mol), and therefore, accurately...

  • Molecular dynamics investigations of the dissociation of SiO2 on an ab initio potential energy surface obtained using neural network methods. Agrawal, Paras M.; Raff, Lionel M.; Hagan, Martin T.; Komanduri, Ranga // Journal of Chemical Physics;4/7/2006, Vol. 124 Issue 13, p134306 

    The neural network (NN) procedure to interpolate ab initio data for the purpose of molecular dynamics (MD) simulations has been tested on the SiO2 system. Unlike other similar NN studies, here, we studied the dissociation of SiO2 without the initial use of any empirical potential. During the...

  • Molecular dynamics investigations of chlorine peroxide dissociation on a neural network ab initio potential energy surface. Le, Anh; Vu, Nam; Dinh, Thach; Cao, Thi; Le, Hung // Theoretical Chemistry Accounts: Theory, Computation, & Modeling;Mar2012, Vol. 131 Issue 3, p1 

    Molecular dissociation of chlorine peroxide (ClOOCl), which consists of two elementary dissociation channels (of Cl-O and O-O), is investigated using molecular dynamics simulations on a neural network-fitted potential energy surface constructed by MP2 calculations with the 6-311G(d,p) basis set....

  • The role of exchange-correlation functionals in the potential energy surface and dynamics of N2 dissociation on W surfaces. Bocan, G. A.; Díez Muiño, R.; Alducin, M.; Busnengo, H. F.; Salin, A. // Journal of Chemical Physics;4/21/2008, Vol. 128 Issue 15, p154704 

    We study the dissociative adsorption of N2 on W(100) and W(110) by means of density functional theory and classical dynamics. Working with a full six-dimensional adiabatic potential energy surface (PES), we find that the theoretical results of the dynamical problem strongly depend on the choice...

  • The A 1Σu+ state of K2 up to the dissociation limit. Falke, St.; Sherstov, I.; Tiemann, E.; Lisdat, Ch. // Journal of Chemical Physics;12/14/2006, Vol. 125 Issue 22, p224303 

    We report an experimental study of the K2 A 1Σu+ state. Long-range levels up to the dissociation limit were observed in a two laser spectroscopic experiment using a highly collimated molecular beam. We derive an analytical potential energy curve for the complete A state including long-range...

  • Molecular dynamics of hydrogen dissociation on an oxygen covered Pt(111) surface. Ludwig, Jeffery; Vlachos, Dionisios G. // Journal of Chemical Physics;4/21/2008, Vol. 128 Issue 15, p154708 

    The dissociation of hydrogen on a Pt(111) surface covered with a p(2×2) oxygen phase was investigated using quasiclassical, six-dimensional molecular dynamics. The potential energy surface (PES) used in these simulations was obtained by an iterative novelty sampling algorithm. Compared to...

  • Highly excited vibrational states of HCN around 30 000 cm-1. Martínez, R. Z.; Lehmann, Kevin K.; Carter, Stuart // Journal of Chemical Physics;11/7/2006, Vol. 125 Issue 17, p174306 

    This article describes the analysis and interpretation of rovibrational spectra involving highly excited vibrational states in the molecule of HCN. The spectra were obtained by means of the vibrationally mediated photodissociation technique. Analysis of the spectra revealed four bands with...

  • Reactive force fields for surface chemical reactions: A case study with hydrogen dissociation on Pd surfaces. Xiao, Y.; Dong, W.; Busnengo, H. F. // Journal of Chemical Physics;1/7/2010, Vol. 132 Issue 1, p014704 

    An approach based on reactive force fields is applied to the parametrization of potential energy surface (PES) for chemical reactions on surfaces with a benchmark system, H2/Pd(111). We show that a simple reactive force field based on the second moment approximation does not allow for obtaining...

  • Photodissociation dynamics of indole in a molecular beam. Ming-Fu Lin; Chien-Ming Tseng; Lee, Yuan T.; Chi-Kung Ni // Journal of Chemical Physics;9/22/2005, Vol. 123 Issue 12, p124303 

    Photodissociation of indole at 193 and 248 nm under collision-free conditions has been studied in separate experiments using multimass ion imaging techniques. H atom elimination was found to be the only dissociation channel at both wavelengths. The photofragment translational energy distribution...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics