Six-dimensional classical dynamics of H2 dissociative adsorption on Pd(111)

H2 Detection Using Plasmonically Generated Surface Photoacoustic Waves in Pd Nanoparticle-Deposited GaN Microcantilevers

Detection of H₂ using plasmonic amplification of surface photoacoustic (SPA) waves generated in Pd nanoparticle-deposited GaN piezotransistive microcantilevers has been investigated using a pulsed 520 nm laser. Using 1.5 nm thickness of the Pd functionalization layer, H₂ detection down to 1.5 ppm was demonstrated with a high signal-to-noise ratio, underscoring the feasibility of sub-ppm level detection using this novel sensing method. Adsorption of H₂ in Pd nanoparticles (NPs) changes their plasmonic absorption spectra because of Pd lattice expansion, in addition to changing their work function. The high sensitivity exhibited by the SPA-based H₂ detection method is attributed to a combination of changes in the plasmonic spectrum and work function of Pd NPs and was observed to be a strong function of Pd thickness, biasing conditions, and probe laser power. A comparison of the SPA-based detection technique with traditional chemidiode and chemiresistor sensors, integrated in the functionalized piezotransistor, indicated a superior detection performance of the former.

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.

Efficient hybrid evolutionary optimization of interatomic potential models

The lack of adequately predictive atomistic empirical models precludes meaningful simulations for many materials systems. We describe advances in the development of a hybrid, population based optimization strategy intended for the automated development of material specific interatomic potentials. We compare two strategies for parallel genetic programming and show that the Hierarchical Fair Competition algorithm produces better results in terms of transferability, despite a lower training set accuracy. We evaluate the use of hybrid local search and several fitness models using system energies and/or particle forces. We demonstrate a drastic reduction in the computation time with the use of a correlation-based fitness statistic. We show that the problem difficulty increases with the number of atoms present in the systems used for model development and demonstrate that vectorization can help to address this issue. Finally, we show that with the use of this method, we are able to "rediscover" the exact model for simple known two- and three-body interatomic potentials using only the system energies and particle forces from the supplied atomic configurations.

Ab initio molecular dynamics study of hot atom dynamics after dissociative adsorption of H2 on Pd(100)

The relaxation of hot hydrogen atoms upon the dissociative adsorption of H2 on Pd(100) was studied by ab initio molecular dynamics simulations based on density functional theory, modeling the full dissociative adsorption process in a consistent manner. In spite of the nonlinear dependence of every single trajectory on the run conditions, on the average it is the energy dissipation to the substrate that determines the mean distance of the two H atoms after adsorption which amounts to three to four lattice units and provides an upper bound for heavier species such as oxygen.

Representing molecule-surface interactions with symmetry-adapted neural networks

The accurate description of molecule-surface interactions requires a detailed knowledge of the underlying potential-energy surface (PES). Recently, neural networks (NNs) have been shown to be an efficient technique to accurately interpolate the PES information provided for a set of molecular configurations, e.g., by first-principles calculations. Here, we further develop this approach by building the NN on a new type of symmetry functions, which allows to take the symmetry of the surface exactly into account. The accuracy and efficiency of such symmetry-adapted NNs is illustrated by the application to a six-dimensional PES describing the interaction of oxygen molecules with the Al(111) surface.

Low sticking probability in the nonactivated dissociation of N2 molecules on W(110)

The six-dimensional potential energy surface for the dissociation of N2 molecules on the W(110) surface has been determined by density functional calculations and interpolated using the corrugation reducing procedure. Examination of the resulting six-dimensional potential energy surface shows that nonactivated paths are available for dissociation. In spite of this, the dissociation probability goes to a very small value when the impact energy goes to zero and increases with increasing energy, a behavior usually associated with activated systems. Statistics on the dynamics indicate that this unconventional result is a consequence of the characteristics of the potential energy surface at long distances. Furthermore, two distinct channels are identified in the dissociation process, namely, a direct one and an indirect one. The former is responsible for dissociation at high energies. The latter, which includes long-lasting dynamic trapping in the vicinity of a potential well above the W top position, is the leading mechanism at low and intermediate energies.

Ab initio based tight-binding molecular dynamics simulation of the sticking and scattering of O2∕Pt(111)

The sticking and scattering of O(2)Pt(111) has been studied by tight-binding molecular dynamics simulations based on an ab initio potential energy surface. We focus, in particular, on the sticking probability as a function of the angle of incidence and the energy and angular distributions in scattering. Our simulations provide an explanation for the seemingly paradox experimental findings that adsorption experiments suggest that the O(2)Pt(111) interaction potential should be strongly corrugated while scattering experiments indicate a rather small corrugation. The potential energy surface is indeed strongly corrugated which leads to a pronounced dependence of the sticking probability on the angle of incidence. The scattered O(2) molecules, however, experience a rather flat surface due to the fact that they are predominantly scattered at the repulsive tail of the potential.

Six-dimensional potential energy surface for H2at Ru(0001)

The six-dimensional (6D) potential energy surface (PES) for the H(2) molecule interacting with a clean Ru(0001) surface has been computed accurately for the first time. Density functional theory (DFT) and a pseudopotential based periodic plane-wave approach have been used to calculate the electronic interactions between the molecule and the surface. Two different generalized gradient approximation (GGA) exchange-correlation functionals, PW91 and RPBE, have been adopted. Based on the DFT/GGA calculated potential energies, an analytical 6D PES has been constructed using the corrugation reducing procedure. A very accurate representation of the DFT/GGA data has been achieved, with an average error in the interpolation of about 3 meV and a maximum error not larger than about 30 meV. The top site is found to be the most reactive site for both functionals used, but PW91 predicts a higher reactivity than RPBE, with lower-energy and earlier-located dissociation barriers. The energetic corrugation displayed by the RPBE PES is larger than the PW91 PES while the geometric corrugation is smaller. The differences between the two PESs increase as the distance of the molecular center of mass to the surface decreases. A direct comparison with experimental investigations on H(2)/Ru(0001) could shed light on the suitability of these XC potentials often used in DFT calculations.

Mechanisms of H2 dissociative adsorption on the Pt(211) stepped surface

We utilize classical trajectory calculations to study the reaction dynamics of the dissociative adsorption of H2 on the stepped Pt(211) surface. The potential-energy surface has been obtained through an accurate interpolation of density-functional theory data at the generalized gradient approximation level, using the corrugation reduction procedure. New techniques for visualizing the collective dynamics of trajectories are introduced to elucidate the reaction mechanisms involved. Reaction exhibits a nonmonotonic dependence on collision energy, first decreasing with energy, and then increasing. A strong component of direct nonactivated reaction exists at the top edge of the step over the entire range of energies. The inverse relationship between reaction and collision energy at low energies is attributed to trapping in weak chemisorption wells. These wells also influence the direct reaction at the step, leading to a strong asymmetric dependence on incidence angle. Reaction on the terrace is activated, and only contributes significantly at high energies. Agreement with experiments on Pt(533) [A. T. Gee, B. E. Hayden, C. Mormiche, and T. S. Nunney, J. Chem. Phys. 112, 7660 (2000); Surf. Sci. 512, 165 (2002)] is good, and we are able to suggest new interpretations of the experimental data.

High-dimensional quantum dynamical study of the dissociation of H2 on Pd(110)

We report the first six-dimensional quantum dynamical study of the dissociative adsorption of H(2) on a (110) surface. We have performed quantum coupled-channel calculations for the system H(2)/Pd(110) based on a potential energy surface (PES) that was derived from ab initio electronic structure calculations. In particular, we have focused on the effects of the corrugation and anisotropy of the PES on the H(2) dissociation probability. Our results agree well with the available experimental data for the sticking probability as a function of the initial kinetic energy and the angle of incidence. Because of the coupling between the anisotropy and corrugation of the potential energy surface our calculations predict an unusual rotational heating and a rather small rotational alignment in desorption.

Theoretical analysis of the relation between H2 dissociation and reflection on Pd surfaces

We study the scattering of H2 (v=0, J=0) molecules by the Pd(110) surface using classical trajectory methods. We show that the dissociative adsorption probability barely depends on incidence angle (total energy scaling) up to an impact energy of 200 meV. This is the consequence of a "loss of memory" of the initial incidence angle, mostly due to dynamic trapping, which also reflects itself in a cosinelike angular distribution of reflected molecules. Consequently, a cosinelike distribution can be the result of a subpicosecond process that involves neither energy dissipation to the surface nor transient thermal accommodation.

Surface temperature dependence of rotational excitation of H(2) scattered from Pd(111)

Classical dynamics simulations are performed to study rotational excitation of H(2) scattered from Pd(111), taking into account energy exchange with surface phonons through a 3D surface oscillator model. We show that dynamic trapping, identified recently in the study of dissociation dynamics, plays a prominent role. The corresponding long interaction time due to several recollisions allows an efficient energy exchange between H(2) molecules and surface phonons. This microscopic mechanism explains the puzzling experimental finding on the role of surface temperature in H(2)(nu = 0,J = 1-->3) excitation.