Search

To study the scattering of CO off a movable Cu(100) surface, extensive multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) calculations are performed based on the SAP [R. Marquardt et al., J. Chem. Phys. 132, 074108 (2010)] potential energy surface in conjunction with a recently developed expansion model [Q. Meng and H.-D. Meyer, J. Chem. Phys. 143, 164310 (2015)] for including lattice motion. The surface vibration potential is constructed by a sum of Morse potentials where the parameters are determined by simulating the vibrational energies of a clean Cu(100) surface. Having constructed the total Hamiltonian, extensive dynamical calculations in both timeindependent and time-dependent schemes are performed. Two-layer MCTDH (i.e., normal MCTDH) block-improved-relaxations (time-independent scheme) show that increasing the number of included surface vibrational dimensions lets the vibrational energies of CO/Cu(100) decrease for the frustrated translation (T mode), which is of low energy but increase those of the frustrated rotation (R mode) and the CO-Cu stretch (S mode), whose vibrational energies are larger than the energies of the in-plane surface vibrations (~79 cm-1). This energy-shifting behavior was predicted and discussed by a simple model in our previous publication [Q. Meng and H.-D. Meyer, J. Chem. Phys. 143, 164310 (2015)]. By the flux analysis of the MCTDH/ML-MCTDH propagated wave packets, we calculated the sticking probabilities for the X + 0D, X + 1D, X + 3D, X + 5D, and X + 15D systems, where "X" stands for the used dimensionality of the CO/rigid-surface system and the second entry denotes the number of surface degrees of freedom included. From these sticking probabilities, the X + 5D/15D calculations predict a slower decrease of sticking with increasing energy as compared to the sticking of the X + 0D/1D/3D calculations. This is because the translational energy of CO is more easily transferred to surface vibrations, when the vibrational dimensionality of the surface is enlarged. [ABSTRACT FROM AUTHOR]

To accurately compute the rates of the abstraction channels of hydrogen plus ethane (Et), propane (Pr), and dimethyl ether (DME), ring-polymer molecular dynamics (RPMD) method is used in conjunction with the recently constructed local permutation invariant polynomial neural-networks potential energy surface of the parent H + CH4 system [Q. Meng et al., J. Chem. Phys. 144, 154312 (2016)]. For H + Et, one of the H atoms in CH4 of the parent system is replaced by a methyl group, while for the H + DME reaction, it is replaced by the methoxyl group. For the H + Pr reaction, replacing one of the H atoms in CH4 by an ethyl group, the terminal channel is built, meanwhile the middle channel is considered through replacing two H atoms in CH4 by two methyl groups. Since the potential energy barriers of the title reactions must differ from the H + CH4 barrier, the corrections have to be made by computing the ratio of free-energy barriers between H + CH4 and the title reactions at coupled cluster with a full treatment singles and doubles (where the triples contribution is calculated by perturbation theory, that is, CCSD(T)) level. Comparing the present RPMD rates with the previous theoretical and experimental results, good agreement can be found. Moreover, probable reasons for the deviation between the present RPMD rates and the previous experimental ones are discussed. [ABSTRACT FROM AUTHOR]

To fast and accurately compute rate coefficients of the H/D + CH4→H2/HD + CH3 reactions, we propose a segmented strategy for fitting suitable potential energy surface (PES), on which ring-polymer molecular dynamics (RPMD) simulations are performed. On the basis of recently developed permutation invariant polynomial neural-network approach [J. Li et al., J. Chem. Phys. 142, 204302 (2015)], PESs in local configuration spaces are constructed. In this strategy, global PES is divided into three parts, including asymptotic, intermediate, and interaction parts, along the reaction coordinate. Since less fitting parameters are involved in the local PESs, the computational efficiency for operating the PES routine is largely enhanced by a factor of ~20, comparing with that for global PES. On interaction part, the RPMD computational time for the transmission coefficient can be further efficiently reduced by cutting off the redundant part of the child trajectories. For H + CH4, good agreements among the present RPMD rates and those from previous simulations as well as experimental results are found. For D + CH4, on the other hand, qualitative agreement between present RPMD and experimental results is predicted. [ABSTRACT FROM AUTHOR]

Molecular-surface studies are often done by assuming a corrugated, static (i.e., rigid) surface. To be able to investigate the effects that vibrations of surface atoms may have on spectra and cross sections, an expansion Hamiltonian model is proposed on the basis of the recently reported [R. Marquardt et al., J. Chem. Phys. 132, 074108 (2010)] SAP potential energy surface (PES), which was built for the CO/Cu(100) system with a rigid surface. In contrast to other molecule-surface coupling models, such as the modified surface oscillator model, the coupling between the adsorbed molecule and the surface atoms is already included in the present expansion SAP-PES model, in which a Taylor expansion around the equilibrium positions of the surface atoms is performed. To test the quality of the Taylor expansion, a direct model, that is avoiding the expansion, is also studied. The latter, however, requests that there is only one movable surface atom included. On the basis of the present expansion and direct models, the effects of a moving top copper atom (the one to which CO is bound) on the energy levels of a bound CO/Cu(100) system are studied. For this purpose, the multiconfiguration time-dependent Hartree calculations are carried out to obtain the vibrational fundamentals and overtones of the CO/Cu(100) system including a movable top copper atom. In order to interpret the results, a simple model consisting of two coupled harmonic oscillators is introduced. From these calculations, the vibrational levels of the CO/Cu(100) system as function of the frequency of the top copper atom are discussed. [ABSTRACT FROM AUTHOR]

The ring polymer molecular dynamics (RPMD) calculations are performed to calculate rate constants for the title reaction on the recently constructed potential energy surface based on permutation invariant polynomial (PIP) neural-network (NN) fitting [J. Li et al., J. Chem. Phys. 142, 204302 (2015)]. By inspecting convergence, 16 beads are used in computing free-energy barriers at 300 K ≤ T ≤ 1000 K, while different numbers of beads are used for transmission coefficients. The present RPMD rates are in excellent agreement with quantum rates computed on the same potential energy surface, as well as with the experimental measurements, demonstrating further that the RPMDis capable of producing accurate rates for polyatomic chemical reactions even at rather low temperatures. [ABSTRACT FROM AUTHOR]

Employing the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method in conjunction with the multistate multimode vibronic coupling Hamiltonian (MMVCH) model, we perform a full dimensional (9D) quantum dynamical study on the simplest Criegee intermediate, formaldehyde oxide, in five lower-lying singlet electronic states. The ultraviolet (UV) spectrum is then simulated by a Fourier transform of the auto-correlation function. The MMVCH model is built based on extensive MRCI(8e,8o)/aug-cc-pVTZ calculations. To ensure a fast convergence of the final calculations, a large number of ML-MCTDH test calculations is performed to find an appropriate multilayer separations (ML-trees) of the ML-MCTDH nuclear wave functions, and the dynamical calculations are carefully checked to ensure that the calculations are well converged. To compare the computational efficiency, standard MCTDH simulations using the same Hamiltonian are also performed. A comparison of the MCTDH and ML-MCTDH calculations shows that even for the present not-too-large system (9D here) the ML-MCTDH calculations can save a considerable amount of computational resources while producing identical spectra as the MCTDH calculations. Furthermore, the present theoretical B¹A' ← X¹A' UV spectral band and the corresponding experimental measurements [J. M. Beames, F. Liu, L. Lu, and M. I. Lester, J. Am. Chem. Soc.134, 20045-20048 (2012); L. Sheps, J. Phys. Chem. Lett.4, 4201-4205 (2013); W.-L. Ting, Y.-H. Chen, W. Chao, M. C. Smith, and J. J.-M. Lin, Phys. Chem. Chem. Phys.16, 10438-10443 (2014)] are discussed. To the best of our knowledge, this is the first theoretical UV spectrum simulated for this molecule including nuclear motion beyond an adiabatic harmonic approximation. [ABSTRACT FROM AUTHOR]