The effect of corrugation on the quantum dynamics of dissociative and diffractive scattering of H2 from Pt(111)

Ring Polymer Molecular Dynamics Approach to Quantum Dissociative Chemisorption Rates

A ring polymer molecular dynamics (RPMD) method is proposed for the calculation of the dissociative chemisorption rate coefficient on surfaces. The RPMD rate theory is capable of handling quantum effects such as the zero-point energy and tunneling in dissociative chemisorption, while it relies on classical trajectories for the simulation. Applications to H2 dissociative chemisorption are demonstrated. For the highly activated process on Ag(111), strong deviations from Arrhenius behavior are found at low temperatures and attributed to tunneling. On Pt(111), where the dissociation has a barrierless pathway, the RPMD rate coefficient is found to agree with the experimentally derived thermal sticking coefficient within a factor of 2 over a large temperature range. Significant quantum effects are also found.

Test of the Transferability of the Specific Reaction Parameter Functional for H2 + Cu(111) to D2 + Ag(111)

The accurate description of the dissociative chemisorption of a molecule on a metal surface requires a chemically accurate description of the molecule-surface interaction. Previously, it was shown that the specific reaction parameter approach to density functional theory (SRP-DFT) enables accurate descriptions of the reaction of dihydrogen with metal surfaces in, for instance, H₂ + Pt(111), H₂ + Cu(111), and H₂ + Cu(100). SRP-DFT likewise allowed a chemically accurate description of dissociation of methane on Ni(111) and Pt(111), and the SRP functional for CH₄ + Ni(111) was transferable to CH₄ + Pt(111), where Ni and Pt belong to the same group. Here, we investigate whether the SRP density functional derived for H₂ + Cu(111) also gives chemically accurate results for H₂ + Ag(111), where Ag belongs to the same group as Cu. To do this, we have performed quasi-classical trajectory calculations using the six-dimensional potential energy surface of H₂ + Ag(111) within the Born-Oppenheimer static surface approximation. The computed reaction probabilities are compared with both state-resolved associative desorption and molecular beam sticking experiments. Our results do not yet show transferability, as the computed sticking probabilities and initial-state selected reaction probabilities are shifted relative to experiment to higher energies by about 2-3 kcal/mol. The lack of transferability may be due to the different character of the SRP functionals for H₂ + Cu and CH₄ + group 10 metals, the latter containing a van der Waals correlation functional and the former not.

SIX-DIMENSIONAL DYNAMICS OF DISSOCIATIVE CHEMISORPTION OF H2 ON METAL SURFACES

The theory of time-dependent quantum dynamics of dissociative chemisorption of hydrogen on metal surfaces is reviewed, in the framework of electronically adiabatic scattering from static surfaces. Four implementations of the time-dependent wave packet (TDWP) method are discussed. In the direct product pseudo-spectral and the spherical harmonics pseudo-spectral methods, no use is made of the symmetry associated with the surface unit cell. This symmetry is exploited by the symmetry adapted wave packet and the symmetry adapted pseudo-spectral (SAPS) method, which are efficient for scattering at normal incidence. The SAPS method can be employed for potential energy surfaces of general form.

Comparison to experiment shows that the TDWP method yields good, but not yet excellent, quantitative accuracy for dissociation of (ν = 0, j = 0) H 2 if the calculations are based on accurately fitted density functional theory calculations that are performed using the generalized gradient approximation. The influence of the molecule's vibration (rotation) is well (reasonably well) described. The theory does not yet yield quantitatively accurate results for rovibrationally inelastic scattering. The theory has helped with the interpretation of existing experimental results, for instance, by solving a parodox regarding the corrugation of Pt(111) as seen by reacting and scattering H 2. The theory has also provided some exciting new predictions, for instance, concerning where on the surface of Cu(100) H2 reacts depending on its vibrational state. Future theoretical studies of H 2 reacting on metal surfaces will likely be aimed at validating GGAs for molecule-surface interactions, and understanding trends in collisions of H 2 with complex metal surfaces.

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

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.

Off-normal incidence dissociative sticking of H2 on Cu(100) studied using six-dimensional quantum calculations

Six-dimensional quantum calculations of the sticking probability for H2 hitting a Cu(100) surface with off-normal incidence are presented. The multiconfiguration time-dependent Hartree approach is employed for an efficient wave-packet propagation. The sticking probability is calculated for different initial momenta parallel to the surface. In contrast with the picture described in the literature, the sticking probability was found to depend on the parallel momentum. The results are explained by the topology of the potential-energy surface, which shows significant corrugation with a moderate variation of the barrier height with the surface site.

Six-dimensional quantum dynamics of (v=0,j=0)D2 and of (v=1,j=0)H2 scattering from Cu(111)

We report six-dimensional quantum dynamics calculations of the dissociative scattering of molecular hydrogen from the copper111 surface. Two potential energy surfaces are investigated and the results are compared with experiment. Our study completes the preliminary work of Somers et al. [Chem. Phys. Lett. 360, 390 (2002)] and focuses on the role of initial vibrational excitation and on isotopic effects. None of the two investigated potential energy surfaces is found satisfactory: the use of neither potential yields reaction and vibrational excitation probabilities and vibrational efficacies that are in close agreement with experiment. In addition to showing the shortcomings of existing potential energy surfaces we point out an inconsistency in the experimental fits for D2.

Multiconfigurational time-dependent Hartree calculations for dissociative adsorption of H2 on Cu(100)

The efficiency of the multiconfigurational time-dependent Hartree (MCTDH) method for calculating the initial-state selected dissociation probability of H(2)(v=0,j=0) on Cu(100) is investigated. The MCTDH method is shown to be significantly more efficient than standard wave packet methods. A large number of single-particle functions is required to converge the initial-state selected reaction probability for dissociative adsorption. Employing multidimensional coordinates in the MCTDH ansatz (mode combination) is found to be crucial for the efficiency of these MCTDH calculations. Perspectives towards the application of the MCTDH approach to study dissociative adsorption of polyatomic molecules on surfaces are discussed.

Theoretical study of the H2 reaction with a Pt4 (111) cluster

The C(s) symmetry reaction of the H(2) molecule on a Pt(4) (111) clusters, has been studied using ab initio multiconfiguration self-consistent field plus extensive multireference configuration interaction variational and perturbative calculations. The H(2) interaction by the vertex and by the base of a tetrahedral Pt(4) cluster were studied in ground and excited triplet and singlet states (closed and open shells), where the reaction curves are obtained through many avoided crossings. The Pt(4) cluster captures and activates the hydrogen molecule; it shows a similar behavior compared with other Pt(n) (n=1,2,3) systems. The Pt(4) cluster in their lowest five open and closed shell electronic states: (3)B(2), (1)B(2), (1)A(1) (3)A(1), (1)A(1), respectively, may capture and dissociate the H(2) molecule without activation barriers for the hydrogen molecule vertex approach. For the threefolded site reaction, i.e., by the base, the situation is different, the hydrogen adsorption presents some barriers. The potential energy minima occur outside and inside the cluster, with strong activation of the H-H bond. In all cases studied, the Pt(4) cluster does not absorb the hydrogen molecule.