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.     

Rotational spectroscopy and molecular structure of the 1,1,2-trifluoroethylene-acetylene complex

Guided by ab initio calculations, Fourier transform microwave rotational spectra in the 6-22 GHz region are obtained for the complex formed between 1,1,2-trifluoroethylene and acetylene, including the normal isotopomer, three of four singly substituted (13)C species obtained in natural abundance, and using commercially available isotopic varieties of acetylene, species containing HCCD and H(13)C(13)CH. Although the ab initio calculations suggest two possible low energy planar arrangements for the molecules in the complex, only a single, unique structure is obtained from a combined analysis of the rotational constants derived from the spectra and atomic positions determined using Kraitchman [Am. J. Phys. 21, 17 (1953)] substitution coordinates. This structure is similar to that obtained for the CF(2)CHF[Single Bond]HF complex [H. O. Leung and M. D. Marshall, J. Chem. Phys. 126, 114310 (2007)] in which both the primary and secondary interactions occur between the HCCH molecule and a F atom and a H atom bonded to the same carbon of CF(2)CHF. The 2.748(15) A hydrogen bond has acetylene as the donor and 1,1,2-trifluoroethylene as the acceptor and forms a 104.49(15) degrees C[Single Bond]Fcdots, three dots, centeredH angle. The 2.8694(9) A secondary interaction between the pi bond of acetylene and the H atom geminal to the acceptor F atom causes the hydrogen bond to deviate 69.24(67) degrees from linearity. This large deviation from linearity and the similarity of the two intermolecular bond lengths suggest that the two interactions are becoming comparable in importance.     

Steering and ro-vibrational effects on dissociative adsorption and associative desorption of H2/Pd(100)

The interaction of hydrogen with many transition metal surfaces is characterized by a coexistence of activated with non-activated paths to adsorption with a broad distribution of barrier heights. By performing six-dimensional quantum dynamical calculations using a potential energy surface derived from ab initio calculations for the system H₂/Pd(100) we show that these features of the potential energy surface lead to strong steering effects in the dissociative adsorption and associative desorption dynamics.

In particular, we focus on the coupling of the translational, rotational and vibrational degrees of freedom of the hydrogen molecule in the reaction dynamics.     

Dynamically biased RRKM model of activated gas-surface reactivity: vibrational efficacy and rotation as a spectator in the dissociative chemisorption of CH4 on Pt(111)

The reactivity of CH(4) impinging on a Pt(111) surface was examined using a precursor-mediated microcanonical trapping model of dissociative chemisorption wherein the effects of rotational and vibrational energy could be explored. Dissociative sticking coefficients for a diverse range of non-equilibrium effusive beam, supersonic beam, and eigenstate-resolved experiments were simulated and an average relative discrepancy between theory and experiment of better than 50% was achieved by treating molecular rotations and translation parallel to the surface as spectator degrees of freedom, and introducing a dynamically-biased vibrational efficacy. The model parameters are {E(0) = 57.9 kJ mol(-1), s = 2, η(v) = 0.40} where E(0) is the apparent threshold energy for reaction, s is the number of surface oscillators participating in energy exchange within each gas-surface collision complex formed, and η(v) is the mean vibrational efficacy for reaction relative to normal translational energy which figures in the assembly of the active exchangeable energy which is available to surmount the activation barrier to dissociative chemisorption. GGA-DFT electronic structure calculations provided vibrational frequencies for the transition state for dissociative chemisorption. The asymmetry of the rotational state populations in supersonic and effusive molecular beam experiments allowed kinetic analysis to establish that taking rotation as a spectator degree of freedom is a good approximation. Surface phonons, rather than the incident molecules, are calculated to play the dominant role in supplying the energy required to overcome the activation barrier for dissociative chemisorption under the thermal equilibrium conditions relevant to high pressure catalysis. Over the temperature range 300 K ≤T≤ 1000 K, the thermal dissociative sticking coefficient is predicted to be well described by S(T) = S(0) exp(-E(a)/RT) where S(0) = 0.62 and E(a) = 62.6 kJ mol(-1).     

D_2在Ni(100)表面离解吸附的量子动力学研究

用量子含时波包法研究了Ｄ２在镍表面顶－桥位上离解吸附量子动力学．计算了不同入射动能及初始振转态的离解几率．讨论了分子的同核对称性、转动取向和振动激发对离解几率的影响，并与其他理论计算结果做了比较．     

Surface bonding and dynamical behavior of the CH3SH molecule on Au(111)

The chemisorption of the undissociated CH3SH molecule on the Au(111) surface has been studied at 5 K using scanning tunneling microscopy (STM) and density functional theory (DFT) calculations. The molecule was found to adsorb on atop Au sites on the defect-free surface. CH3SH undergoes hindered rotation about the Au-S bond on the defect-free surface which is seen in STM as a time-averaged 6-fold pattern. The pattern suggests that the potential minima directions occur for the rotating molecule at the six hollow sites surrounding the atop adsorption site. The barrier for rotation, obtained by DFT calculations, is approximately 0.1 kcal.mol(-1). At low coverages, preferential adsorption occurs at defect sites in the surface, namely, the herringbone "elbows" and random atomic step sites. Molecules adsorbed on these sites do not exhibit rotational freedom.     

A multilayer multiconfigurational time-dependent Hartree approach for quantum dynamics on general potential energy surfaces

The multiconfigurational time-dependent Hartree (MCTDH) approach facilitates multidimensional quantum dynamics calculations by representing the wavepacket in an optimal set of time-dependent basis functions, called single-particle functions. Choosing these single-particle functions to be themselves multidimensional wavefunctions which are represented using a MCTDH representation, a multilayer MCTDH scheme has been constructed and used for quantum dynamics calculations treating up to 1000 degrees of freedom rigorously [Wang and Thoss, J. Chem. Phys. 199, 1289 (2003)]. The present work gives a practical scheme which facilitates the application of the multilayer MCTDH approach, which previously has only been employed to study systems described by model-type Hamiltonians, to molecular systems described by more complicated Hamiltonians and general potential energy surfaces. A multilayer extension of the correlation discrete variable representation (CDVR) scheme employed in MCTDH calculations studying quantum dynamics on general potential energy surfaces is developed and tested in a simple numerical application. The resulting multilayer MCTDH/CDVR approach might offer a perspective to rigorously describe the quantum dynamics of larger polyatomic systems.     

Electronic properties and dissociative photoionization of thiocyanates. Part II. Valence and shallow-core (sulfur and chlorine 2p) regions of chloromethyl thiocyanate, CH2ClSCN

A combination of photoelectron spectroscopy and synchrotron based photoelectron photoion coincidence (PEPICO) spectra has been applied to investigate the electronic structure and the dissociative ionization of the CH(2)ClSCN molecule in the valence region. The PES is assigned with the electronic structure calculations at the outer-valence Green's function and symmetry adapted cluster/configuration interaction (SAC-CI) levels offer an explanation of our experimental results. Upon vacuum ultraviolet irradiation the low-lying radical cation, located at 10.39 eV is formed. The molecular ion is observed in the time-of-flight mass spectra, together with the CH(2)SCN(+) and CH(2)Cl(+) daughter ions. The total ion yield spectra have been measured in the S 2p and Cl 2p regions and several channels have been determined in dissociative photoionization events for the core-excited species. Thus, by using time-of-flight mass spectrometry and synchrotron radiation the relative abundances of the ionic fragments and their kinetic energy release values were obtained from both PEPICO and photoelectron photoion photoion coincidence spectra. Possible fragmentation processes are discussed and compared with that found for the related CH(3)SCN species.     

A 4D wave packet study of the CH3I photodissociation in the A-band. Comparison with femtosecond velocity map imaging experiments

The time-resolved photodissociation dynamics of CH(3)I in the A-band has been studied theoretically using a wave packet model including four degrees of freedom, namely the C-I dissociation coordinate, the I-CH(3) bending mode, the CH(3) umbrella mode, and the C-H symmetric stretch mode. Clocking times and final product state distributions of the different dissociation (nonadiabatic) channels yielding spin-orbit ground and excited states of the I fragment and vibrationless and vibrationally excited (symmetric stretch ν(1) and umbrella ν(2) modes) CH(3) fragments have been obtained and compared with the results of femtosecond velocity map imaging experiments. The wave packet calculations are able to reproduce with very good agreement the experimental reaction times for the CH(3)(ν(1), ν(2))+I*((2)P(1/2)) dissociation channels with ν(1) = 0 and ν(2) = 0,1,2, and also for the channel CH(3)(ν(1) = 0, ν(2) = 0)+I((2)P(3/2)). However, the model fails to predict the experimental clocking times for the CH(3)(ν(1), ν(2))+I((2)P(3/2)) channels with (ν(1), ν(2)) = (0, 1), (0, 2), and (1, 0), that is, when the CH(3) fragment produced along with spin-orbit ground state I atoms is vibrationally excited. These results are similar to those previously obtained with a three-dimensional wave packet model, whose validity is discussed in the light of the results of the four-dimensional treatment. Possible explanations for the disagreements found between theory and experiment are also discussed.     

Differential reaction cross sections from rotationally resolved quantum scattering calculations: application to gas-phase S(N)2 reactions

Differential reaction cross sections have been computed based on previous rotationally resolved time-independent quantum-mechanical scattering calculations for the complex-forming S(N)2 reaction Cl(-) + CH(3)Br → ClCH(3) + Br(-). The results show almost isotropic cross sections for reactant molecules with high rotational quantum numbers. Backward scattering is disfavoured for reaction out of states with small rotational excitation, in particular the rovibrational ground state. This is a quantum-mechanical effect (interference of partial waves) that can partly be rationalized by simple classical arguments. In particular for higher vibrational excitations, an umbrella effect can be observed that favours the backward direction. It can be explained by the strong enhancement of the reactivity by opening a direct mechanism. The ion-dipole interaction exerts a torque onto the molecule which carries out a rotation by about 90° and then completes the reaction.     

Nonequilibrium activated dissociative chemisorption: SiH4 on Si(100)

A three-parameter local hot spot model of gas-surface reactivity is employed to analyze and predict dissociative sticking coefficients for SiH4 incident on Si(100) under varied nonequilibrium conditions. Two Si surface oscillators and the molecular vibrations, rotations, and translational energy directed along the local surface normal are active degrees of freedom in the 15 dimensional microcanonical kinetics. The threshold energy for SiH4 dissociative chemisorption is found to be 19 kJ/mol, in quantitative agreement with recent GGA-DFT calculations. A simple scheme for increasing the rate of chemical vapor deposition of silicon from SiH4 at low surface temperatures is modeled.     

Dissociative chemisorption and energy transfer for methane on Ir(111)

A 3-parameter local hot spot model of gas-surface reactivity is employed to analyze and predict dissociative sticking coefficients for CH(4) incident on Ir(111) under varied nonequilibrium and equilibrium conditions. One Ir surface oscillator and the molecular vibrations, rotations, and translational energy directed along the surface normal are treated as active degrees of freedom in the 14 dimensional microcanonical kinetics. The threshold energy for CH(4) dissociative chemisorption on Ir(111) derived from modeling molecular beam experiments is E(0) = 39 kJ/mol. Over more than 4 orders of magnitude of variation in sticking, the average relative discrepancy between the beam and theoretically derived sticking coefficients is 88%. The experimentally observed enhancement in dissociative sticking as beam translational energies decrease below approximately 10 kJ/mol is consistent with a parallel dynamical trapping/energy transfer channel that likely fails to completely thermalize the molecules to the surface temperature. This trapping-mediated sticking, indicative of specific energy transfer pathways from the surface under nonequilibrium conditions, should be a minor contributor to the overall dissociative sticking at thermal equilibrium. Surprisingly, the CH(4) dissociative sticking coefficient predicted for Ir(111) surfaces at thermal equilibrium, based on the molecular beam experiments, is roughly 4 orders of magnitude higher than recent measurements on supported nanoscale Ir catalysts at 1 bar pressure, which suggests that substantial improvements in catalyst turnover rates may be possible.     

Multidimensional density operator propagations in open systems: model studies on vibrational relaxations and surface sticking processes

An efficient method for the numerical treatment of multidimensional dynamics of open systems is presented: the multiconfiguration time-dependent Hartree (MCTDH) method extended to the propagation of density operators. With this method we investigate the relaxation process of a CO molecule adsorbed on a copper surface, i.e., CO/Cu(100), after the excitation with an infrared (IR) pulse. The interaction potential was taken from the literature. Lifetime estimations and thermalization studies were performed on the IR excited CO molecule. We were able to treat this system with all six degrees of freedom (DOF) and thus 12 dynamical variables, but most of our studies used a two or four DOF model. Finally, we demonstrate the applicability of MCTDH to the analysis of scattering processes in an open environment. We calculate sticking coefficients of a scattered particle to a model surface, the latter acting as heat bath. The surface corrugation and the initial particle energy have been varied, and six different relaxation strengths have been studied. These calculations were done under the inclusion of three DOFs: the two surface coordinates and the distance between the particle and the surface.     

Reduced dimensionality quantum scattering calculations on the F + CH(4) --> FH + CH(3) reaction

The existence of recently observed scattering resonances in the hydrogen abstraction reaction F + CH4 --> FH + CH3 was investigated using the reduced dimensionality rotating line umbrella (RLU) quantum scattering model and employing an analytical potential energy surface, PES-2006, recently developed by our group. The calculations were performed in hyperspherical coordinates. The wells found in the hyperspherical adiabats, the oscillatory pattern in the cumulative and state-to-state reaction probabilities, the forward/backward predominance in the differential cross section at a collision energy of 1.8 kcal mol(-1), and the dramatic change of the scattering angle with energy are related to scattering resonances, and they are assigned to a quasi-bound complex on the vibrationally adiabatic ground-state potential.     

A multi-configurational time-dependent Hartree approach to the eigenstates of multi-well systems

A rigorous and efficient approach for the calculation of eigenstates in polyatomic molecular systems with potentials displaying multiple wells is introduced. The scheme is based on the multi-configurational time-dependent Hartree (MCTDH) approach and uses multiple MCTDH wavefunctions with different single-particle function bases to describe the quantum dynamics in the different potential wells. More specifically, an iterative block Lanczos-type diagonalization scheme utilizing state-averaged MCTDH wavefunctions localized in different wells is employed to obtain the energy eigenvalues and eigenstates. The approach does not impose any formal restriction on the symmetry of the potential or the number of wells. A seven-dimensional model system of tetrahedral symmetry, which is inspired by A·CH(4) type complexes and displays four equivalent potential minima, is used to study the numerical performance of the new approach. It is found that the number of configurations in the MCTDH wavefunctions required to obtain converged results is decreased by roughly one order of magnitude compared to standard MCTDH calculations employing a block-relaxation scheme.     

Reaction dynamics of Ca(4s3d 1D2)+CH4-->CaH(X 2Sigma+)+CH3: reaction pathway and energy disposal for the CaH product

The reaction pathway for Ca(4s3d 1D2)+CH4-->CaH(X 2Sigma+)+CH3 has been investigated by using a pump-probe technique in combination with potential-energy surface (PES) calculations. The nascent product distributions of CaH have been characterized with Boltzmann rotational temperatures of 1013+/-102 and 834+/-70 K for the v=0 and 1 levels, respectively, and a Boltzmann vibrational temperature of 1313+/-173 K. The rotational and vibrational energy partitions in CaH have been estimated to be 461+/-45 and 252+/-15 cm(-1), respectively. According to the PES calculations, the pathway favors an insertion mechanism. Ca(3 1D2) approaches CH4 along an attractive potential surface in a C2v (or Cs) symmetry and then the collision complex undergoes nonadiabatic transition to the reactive ground-state surface. An Arrhenius plot shows a potential-energy requirement of 2695+/-149 cm(-1), which accounts for the endothermicity of 2930 cm(-1) for the reaction scheme. The Ca-C bond distance in the transition state structure is short enough to allow for tight orbital overlap between CaH and CH3. The strong coupling between the moieties renders the energy transfer sufficient from CaH into the CH3 radical. As compared to the Ca(4 1P1) reaction, the dissociation lifetime of the intermediate complex with less excess energy is prolonged so as to cause much less vibrational energy disposal into CaH.     

The dissociative chemisorption of methane on Ni(100): reaction path description of mode-selective chemistry

We derive a model for the dissociative chemisorption of methane on a Ni(100) surface, based on the reaction path Hamiltonian, that includes all 15 molecular degrees of freedom within the harmonic approximation. The total wavefunction is expanded in the adiabatic vibrational states of the molecule, and close-coupled equations are derived for wave packets propagating on vibrationally adiabatic potential energy surfaces, with non-adiabatic couplings linking these states to each other. Vibrational excitation of an incident molecule is shown to significantly enhance the reactivity, if the molecule can undergo transitions to states of lower vibrational energy, with the excess energy converted into motion along the reaction path. Sudden models are used to average over surface impact site and lattice vibrations. Computed dissociative sticking probabilities are in good agreement with experiment, with respect to both magnitude and variation with energy. The ν(1) vibration is shown to have the largest efficacy for promoting reaction, due to its strong non-adiabatic coupling to the ground state, and a significant softening of the vibration at the transition state. Most of the reactivity at 475 K is shown to result from thermally assisted over-the-barrier processes, and not tunneling.     

Hydrogen molecule in the small dodecahedral cage of a clathrate hydrate: quantum five-dimensional calculations of the coupled translation-rotation eigenstates

We report quantum five-dimensional (5D) calculations of the energy levels and wave functions of the hydrogen molecule, para-H2 and ortho-H2, confined inside the small dodecahedral (H2O)20 cage of the sII clathrate hydrate. All three translational and the two rotational degrees of freedom of H2 are included explicitly, as fully coupled, while the cage is treated as rigid. The 5D potential energy surface (PES) of the H2-cage system is pairwise additive, based on the high-quality ab initio 5D (rigid monomer) PES for the H2-H2O complex. The bound state calculations involve no dynamical approximations and provide an accurate picture of the quantum 5D translation-rotation dynamics of H2 inside the cage. The energy levels are assigned with translational (Cartesian) and rotational quantum numbers, based on calculated root-mean-square displacements and probability density plots. The translational modes exhibit negative anharmonicity. It is found that j is a good rotational quantum number, while the threefold degeneracy of the j = 1 level is lifted completely. There is considerable translation-rotation coupling, particularly for excited translational states.