Experimental and theoretical investigations of the inelastic and reactive scattering dynamics of O(3p) + D2

This paper presents a combined experimental and theoretical study of the dynamics of O((3)P) + D(2) collisions, with emphasis on a center-of-mass (c.m.) collision energy of 25 kcal mol(-1). The experiments were conducted with a crossed-molecular-beams apparatus, employing a laser detonation source to produce hyperthermal atomic oxygen and mass spectrometric detection to measure the product angular and time-of-flight distributions. The novel beam source, which enabled these experiments to be conducted, contributed unique challenges to the experiments and to the analysis, so the experimental methods and approach to the analysis are discussed in detail. Three different levels of theory were used: (1) quasiclassical trajectories (QCT), (2) time-independent quantum scattering calculations based on high-quality potential surfaces for the two lower-energy triplet states, and (3) trajectory-surface-hopping (TSH) studies that couple the triplet surfaces with the lowest singlet surface using a spin-orbit Hamiltonian derived from ab-initio calculations. The latter calculations explore the importance of intersystem crossing in the dynamics. Both experiment and theory show that inelastically scattered O atoms scatter almost exclusively in the forward direction, with little or no loss of translational energy. For the reaction, O((3)P) + D(2) --> OD + D, the experiment shows that, on average, approximately 50% of the available energy goes into product translation and that the OD product angular distributions are largely backward-peaked. These results may be interpreted in light of the QCT and TSH calculations, leading to the conclusion that the reaction occurs mainly on triplet potential energy surfaces with, at most, minor intersystem crossing to a singlet surface. Reaction on either of the two low-lying reactive triplet surfaces proceeds through a rebound mechanism in which the angular distributions are backward-peaked and the product OD is both vibrationally and rotationally excited. The quantum scattering results are in good agreement with QCT calculations, indicating that quantum effects are relatively small for this reaction at a collision energy of 25 kcal mol(-1).     

The role of molecular rotation in activated dissociative adsorption on metal surfaces

The role of molecular rotation in dissociative adsorption of H2 on the activated NiAl(110) metal surface is systematically investigated by means of classical dynamics calculations performed on ab initio six-dimensional potential energy surfaces. The calculations show that molecules rotate abruptly when they are close to the surface and that this rotation allows the molecules to adopt the orientation that is more convenient for dissociation (i.e., nearly parallel to the surface). Also, in reactive sectors of the NiAl(110) unit cell, there is an "angular threshold" below which molecules cannot dissociate. This angular threshold goes down as the incidence energy increases, which explains the rise of the dissociation probability and the fact that it reaches a value close to 1 at incidence energies of the order of 2 eV. The fact that switching on molecular rotation favors dissociation establishes a competition between dissociation and rotational excitation of reflected molecules above the dissociation threshold. Measurements on rotational excitation might thus bring indirect evidence on the dissociation dynamics. Sample calculations for nonactivated Pd(111) and activated Cu(110) metal surfaces suggest that some of these conclusions may be of general validity.     

多原子反应体系的高精度拟合势能面

本文介绍了近几年来我们组构建多原子反应体系的高精度拟合势能面的进展。我们基于神经网络(NN)方法,成功构建了多原子气相分子体系和气相分子在金属表面解离的一系列势能面。这些势能面的拟合精度相当高,并且经过了严格的量子动力学测试,能广泛应用到动力学研究中。我们还提出了一种新的置换不变势能面的拟合方法,即基本不变量神经网络方法(FI-NN)。基本不变量的使用极大地减少了神经网络输入层多项式的个数,有效提高了势能面的计算速度。     

Generalization of classical mechanics for nuclear motions on nonadiabatically coupled potential energy surfaces in chemical reactions

Classical trajectory study of nuclear motion on the Born-Oppenheimer potential energy surfaces is now one of the standard methods of chemical dynamics. In particular, this approach is inevitable in the studies of large molecular systems. However, as soon as more than a single potential energy surface is involved due to nonadiabatic coupling, such a naive application of classical mechanics loses its theoretical foundation. This is a classic and fundamental issue in the foundation of chemistry. To cope with this problem, we propose a generalization of classical mechanics that provides a path even in cases where multiple potential energy surfaces are involved in a single event and the Born-Oppenheimer approximation breaks down. This generalization is made by diagonalization of the matrix representation of nuclear forces in nonadiabatic dynamics, which is derived from a mixed quantum-classical representation of the electron-nucleus entangled Hamiltonian [Takatsuka, K. J. Chem. Phys. 2006, 124, 064111]. A manifestation of quantum fluctuation on a classical subsystem that directly contacts with a quantum subsystem is discussed. We also show that the Hamiltonian thus represented gives a theoretical foundation to examine the validity of the so-called semiclassical Ehrenfest theory (or mean-field theory) for electron quantum wavepacket dynamics, and indeed, it is pointed out that the electronic Hamiltonian to be used in this theory should be slightly modified.     

On the dynamics of H2 adsorption on the Pt(111) surface

The dynamics and kinetics of the dissociation of hydrogen over the hexagonal close packed platinum (Pt(111)) surface are investigated using Car–Parrinello molecular dynamics and static density functional theory calculations of the potential energy surfaces. The calculations model the reference energy‐resolved molecular beam experiments, considering the degrees of freedom of the catalytic surface. Two‐dimensional potential energy surfaces above the main sites on Pt(111) are determined. Combined with Car–Parrinello trajectories, they confirm the dissociative adsorption of H₂ as the only adsorption pathway on this surface at H₂ incindence energies above 5 kJ/mol. A direct determination of energy‐resolved sticking coefficients from molecular dynamics is also performed, showing an excellent agreement with the experimental data at incidence energies in the 5–30 kJ/mol range. Application of dispersion corrections does not lead to an improvement in the prediction of the H₂ sticking coefficient. The adsorption reaction rate obtained from the calculated sticking coefficients is consistent with experimentally derived literature values.     

Quantum mechanical calculation of energy dependence of OCl/OH product branching ratio and product quantum state distributions for the O(1D) + HCl reaction on all three contributing electronic state potential energy surfaces

OCl/OH product branching ratios are calculated as a function of total energy for the O( (1) D) + HCl reaction using quantum wavepacket methods. The calculations take account of reaction on all the three electronic state potential energy surfaces which correlate with both reactants and products. Our results show that reaction on the excited electronic state surfaces has a large effect on the branching ratio at higher energies and that these surfaces must therefore be fully taken into account. The calculations use the potential energy surfaces of Nanbu and co-workers. Product vibrational and rotational quantum state distributions are also calculated as a function of energy for both product channels. Inclusion of the excited electronic state potential energy surfaces improves the agreement of the predicted product vibrational quantum state distributions with experiment for the OH product channel. For OCl agreement between theory and experiment is retained for the vibrational quantum state distributions when the excited electronic state potential energy surfaces are included in the analysis. For the rotational state distributions good agreement between theory and experiment is maintained for energies at which experimental results are available. At higher energies, above 0.7 eV of total energy, the OCl rotational state distributions predicted using all three electronic state potential energy surfaces shift to markedly smaller rotational quantum numbers.     

Dissociative chemisorption dynamics of small molecules on metal surfaces

Much progress has been achieved for both experimental and theoretical studies on the dissociative chemisorption of molecules on surfaces.Quantum state-resolved experimental data has provided unprecedented details for these fundamental steps in heterogeneous catalysis,while the quantitative dynamics is still not fully understood in theory.An in-depth understanding of experimental observations relies on accurate dynamical calculations,in which the potential energy surface and adequate quantum mechanical implementation are desired.This article summarizes the current methodologies on the construction of potential energy surfaces and the quantum mechanical treatments,some of which are promising for future applications.The challenges in this field are also addressed.     

A coupled states calculation of accurate quantum rate constants for H + H2

In this article we present the results of converged quantum reactivescattering calculations of thermal rate constants for H + H₂ using the Liu-Siegbahn-Truhlar-Horowitz (LSTH) potential energy surface. These calculations are based on the coupled states (CS) approximation wherein rotational states having different body fixed angular momentum projection quantum numbers are decoupled. By comparision with accurate coupled channel results on the Porter-Karplus No. 2 (PK2) potential surface, we estimate that the maximumerror in thermal rate constants arising from both this approximation and from other numerical approximations in the calculation is less than 25%. We also show that the sum over projection quantum numbers Ω associated with the CS calculation may be approximated quite accurately in terms of the Ω = 0 rate constants by assuming that the |Ω| > 0 rate constants differ from Ω = 0 by a shift in activation energy, which reflects the vibrationally adiabatic bending energy associated with each Ω. Comparison of the LSTH rate constants with experiment indicates average errors of 16% and 24% relative to the two modern measurementsof the rate constants for H + H₂. These errors are reduced to 18% and 9% if the CS rate constants are multiplied by exp(0.0065 eV/kT). The expected error (based on recent quantum structure calculations) associated with the 0.425 eV barrier of theLSTH potential surface is 0.0065 eV. Overall, the agreement of either the LSTH or modified LSTH rate constants with experiment is within the 32% maximum disagreement between the two experimental measurements at all butthe lowest temperature that has been studied. Comparison of our CS rate constants with the results from simpler theories is considered using both the LSTH and PK2 potential surfaces. The best overall agreement is with transition state theories that use accurate dynamical methods to calculate tunnelling factors. These include reduced dimensionality quantum dynamics methods and variational transition state theory using either the Marcus-Coltrin or least action ground-state tunnelling paths. Comparison with the results of quasiclassical trajectory calculations indicates substantial errors at low temperatures.     

A theoretical study on excited state double proton transfer reaction of a 7-azaindole dimer: an ab initio potential energy surface and its empirical valence bond model

Double proton transfer (DPT) reaction of a 7-azaindole dimer in the first ππ* electronically excited state was studied theoretically. We investigated the reaction mechanism through constructing a full dimensional empirical valence bond potential energy function (PEF) based on potential energies evaluated by ab initio molecular orbital methods, and carrying out quantum dynamics calculations with the PEF. Potential energy surfaces of the DPT obtained at the multi-reference perturbation level of theory favors a concerted DPT mechanism, although a stepwise channel is suggested to open for an excited initial vibrational state. Reduced two dimensional quantum dynamics calculations for a reaction surface Hamiltonian of DPT coordinates were performed. Time constants of the reaction were evaluated to be on the order of picoseconds, which is consistent with experiments. On the other hand, the computed kinetic isotope effect deviates from experimental evidence, suggesting the importance of intermolecular stretching motion, which is not explicit in the present calculations for the quantum effect.     

Local physical quantities for spin based on the relativistic quantum field theory in molecular systems

Local physical quantities for spin are investigated on the basis of the four‐ and two‐component relativistic quantum theory. In the quantum field theory, local physical quantities for spin such as the spin angular momentum density, spin torque density, zeta force density, and zeta potential play important roles in spin dynamics. We discuss how to calculate these local physical quantities based on the two‐component relativistic quantum theory. Some different types of relativistic numerical calculations of local physical quantities in Li atom and C₆H₆ are demonstrated and compared. Local physical quantities for each orbital are also discussed, and it is seen that a total local zeta potential is given as a result of some cancellation of large contributions from each orbital. © 2016 Wiley Periodicals, Inc.     

Wave packet dynamics of photon- and electron-stimulated desorption of ammonia from surfaces

Time-dependent quantum mechanical calculations have been carried out for the photon- and electron-stimulated desorption of ammonia from metal or semiconductor surfaces. The desorption is facilitated by a short-lived complex which excites the N---H₃ inversion mode. The desorption yield and its isotope effect have been obtained from a wave packet on two-dimensional empirical potential energy surfaces. The translational and vibrational distributions of the desorbed ammonia are also calculated. The desorption mechanism includes both direct and predesorption, but the latter predominates. It is shown that the quantum desorption dynamics is much more complex than the simple MGR model.     

Spin-orbit effect in the energy pooling reaction O2(a 1Delta)+O2(a 1Delta)-->O2(b 1Sigma)+O2(X 3Sigma)

Five-dimensional nonadiabatic quantum dynamics studies have been carried out on two new potential energy surfaces of S(2)((1)A(')) and T(7)((3)A(")) states for the title oxygen molecules collision with coplanar configurations, along with the spin-orbit coupling between them. The ab initio calculations are based on complete active state second-order perturbation theory with the 6-31+G(d) basis set. The calculated spin-orbit induced transition probability as a function of collision energy is found to be very small for this energy pooling reaction. The rate constant obtained from a uniform J-shifting approach is compared with the existing theoretical and experimental data, and the spin-orbit effect is also discussed in this electronic energy-transfer process.     

Application of integrated computer simulation approach to solid surfaces and interfaces

The atomistic understanding of the structure, reactivity, and electronic properties of solid surfaces and interfaces are essential for the design of novel catalysts and electronics/photonics devices which have high-performance and unexplored properties. Computational chemistry is expected not only to rationalize the experimental results but also to predict new features. We have applied integrated computer simulation methods including quantum chemistry, periodic density functional theory, molecular dynamics, embedded atom method, and atomic force microscopy simulation to various topics related to solid surfaces and interfaces. In the present paper, we reviewed our recent activities on supported metal catalysts, metal clusters, atomic force microscopy simulation, high-temperature superconductors, tribology, Si semiconductor and V₂O₅ catalysts. Our activities also involve the generation of a lot of new computer simulation codes. We emphasize that the integrated computer simulation system provides not only methods for scientific studies but also a key technology for industrial innovations in research and development.     

Cover Image

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe–O distance, Fe–O–O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin–orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin–orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau–Zener model underestimated the nonadiabatic transition probability.     

Efficient evaluation of the accuracy of molecular quantum dynamics on an approximate analytical or interpolated ab initio potential energy surface

Ab initio electronic structure methods have reached a satisfactory accuracy for the calculation of static properties, but remain too expensive for quantum dynamical calculations. Recently, an efficient semiclassical method was proposed to evaluate the accuracy of quantum dynamics on an approximate potential without having to perform the expensive quantum dynamics on the accurate potential. Here, this method is applied for the first time to evaluate the accuracy of quantum dynamics on an approximate analytical or interpolated potential in comparison to the quantum dynamics on an accurate potential obtained by an ab initio electronic structure method. Specifically, the vibrational dynamics of H₂ on a Morse potential is compared with that on the full CI potential, and the photodissociation dynamics of CO₂ on a LEPS potential with that on the excited ¹Π surface computed at the EOM‐CCSD/aug‐cc‐pVDZ level of theory. Finally, the effect of discretization of a potential energy surface on the quantum dynamics is evaluated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2426–2435, 2010     

The importance of tunneling in the first hydrogenation step in ammonia synthesis over a Ru(0001) surface

The hydrogenation of nitrogen (N(ads)+H(ads)-->NH(ads)) on metal surfaces is an important step in ammonia catalysis. We investigate the reaction dynamics of this hydrogenation step by time independent scattering theory and variational transition state theory (VTST) including tunneling corrections. The potential energy surface is derived by hybrid density functional theory on a model cluster composed of 12 ruthenium atoms resembling a Ru(0001) surface. The scattering calculations are performed on a reduced dimensionality potential energy hypersurface, where two dimensions are treated explicitly and all others are included implicitly by the zero-point correction. The VTST calculations include quantum effects along the reaction coordinate by applying the small curvature tunneling scheme. Even at room temperature (where ruthenium already shows catalytic activity) we find rate enhancement by tunneling by a factor of approximately 70. Inspection of the reaction probabilities shows that the major contribution to reactivity comes from the vibrational ground state of the reactants into vibrationally excited product states. The reaction rates are higher than determined in previous studies, and are compatible with experimental overall rates for ammonia synthesis.     

Classical and quantum studies of the photodissociation of a HX (X=Cl,F) molecule adsorbed on ice

The photodissociation dynamics of a HX (X = Cl,F) molecule adsorbed on a hexagonal ice surface at T = 0 K is studied using time-dependent quantum wave packets and quasiclassical trajectories. The relevant potential energy surfaces are calculated using high-level ab initio methods. We present here two dimensional calculations for the dynamics of the hydrogen photofragment for both HCl and HF molecules. The purpose of this paper is to compare the photodissociation dynamics of the two molecules which are adsorbed on the ice surface with different equilibrium geometries. The total photodissociation cross section and the angular distribution are calculated. The comparison with classical trajectory calculations provides evidence for typical quantum effects and reveals rainbow structures.     

Investigation of ligand exchange reactions in aqueous uranyl carbonate complexes using computational approaches

Carbonate anion exchange reactions with water in the uranyl-carbonate and calcium-uranyl-carbonate aqueous systems have been investigated using computational methods. Classical molecular dynamics (MD) simulations with the umbrella sampling technique were employed to determine potentials of mean force for the exchange reactions of water and carbonate. The presence of calcium counter-ions is predicted to increase the stability of the uranyl-carbonate species in accordance with previous experimental observations. However, the free energy barrier to carbonate exchange with water is found to be comparable both in the presence and absence of calcium cations. Possible implications of these results for uranyl adsorption on mineral surfaces are discussed. Density functional theory (DFT) calculations were also used to confirm the trends observed in classical molecular dynamics simulations and to corroborate the validity of the potential parameters employed in the MD scheme.     

Atomistic modeling of enantioselection in chromatography

A review of atomistic molecular modeling studies related to chromatographic separations of enantiomers is presented. Only those types of calculations where direct interactions between a selector and a selectand are involved are described in this review; omitted are regression models. An emphasis is placed on comparing methods used for sampling potential energy surfaces implementing different methodologies like quantum and molecular mechanics for energy calculations, and molecular dynamics and Monte Carlo sampling strategies for simulations. Type I-V chiral stationary phases and additives for capillary electrophoresis and ion-pair chromatography are covered in this review.