Full-dimensional Quantum Dynamics Study Research Articles

Mode specificity of water dissociating on Ni(100): An approximate full-dimensional quantum dynamics study.

The mode-specific dynamics for the dissociative chemisorption of H2O on rigid Ni(100) is investigated by approximate nine-dimensional (9D) quantum dynamics calculations. The vibrational state-specific 9D dissociation probabilities are obtained by site-averaging the site-specific seven-dimensional results based on an accurate full-dimensional potential energy surface newly developed by neural network fitting to density functional theory energy points with the revised version of the Perdew, Burke, and Ernzerhof functional. The mode specificity of H2O/Ni(100) is very different from that of H2O/Ni(111) or H2O/Cu(111) whose reactivity enhancement by vibrational excitations is quite efficient. For H2O/Ni(100), it is found that the excitation in the symmetric stretching mode is more efficacious than increasing the translational energy in promoting the reaction, while the excitations in the asymmetric stretching mode and bending mode are less efficacious than the translational energy at low collision energies. These interesting observations can be attributed to the near central-barrier reaction for H2O/Ni(100), as well as large discrepancies between the site-specific mode specificities at different impact sites. The mode-specific dynamics obtained in this study is different from that obtained with the reaction path Hamiltonian approach, indicating the importance of full-dimensional quantum dynamics for gas-surface reactions.

Toward a Comprehensive Understanding of Mode-Specific Dynamics of Polyatomic Reactions: A Full-Dimensional Quantum Dynamics Study of the H + NH3 Reaction.

Mode specificity not only sheds light on reaction dynamics but also opens the door for chemical reaction control. This work reports a state-of-the-art full-dimensional quantum dynamics study on the prototypical hydrogen abstraction reaction of hydrogen with ammonia, which serves as a benchmark for advancing our fundamental understanding of polyatomic reaction dynamics. By taking advantage of the (3 + 1) Radau-Jacobi coordinates, the bond-specific probabilities are resolved with the reactant NH3 initiated from either a non-degenerate or degenerate stretching vibrational state. The observed different atom-specific abstraction probabilities from individual states of the degenerate pair are rationalized in the local mode representation according to the different vibrational energy deposited in each N-H bond. It is verified that the three H atoms in NH3 have the same abstraction probability as that from the degenerate pair and the linear combination of the degenerate pair gives the same reaction probability. In addition, the symmetric and asymmetric stretching modes of the reactant NH3 have comparable efficacies on driving the reaction.

Quantum Effects in the Dissociative Chemisorption of N2 on Fe(111): Full-Dimensional Quantum Dynamics and Quasi-Classical Trajectory Study

Quantum Effects in the Dissociative Chemisorption of N₂ on Fe(111): Full-Dimensional Quantum Dynamics and Quasi-Classical Trajectory Study

Full-dimensional quantum dynamics study of isotope effects for the H2 + NH2/ND2/NHD and H2/D2/HD + NH2 reactions.

A full-dimensional quantum dynamical study for the bimolecular reactions of hydrogen molecules with amino radicals for different isotopologues is reported. The nonreactive amino radical is described by two Radau vectors that are very close to the valence bond coordinates. Potential-optimized discrete variable representation basis is used for the vibrational coordinates of the amino radical. Starting from the reaction H2 + NH2, we study the isotope effects for the two reagents separately, i.e., H2 + NH2/ND2/NHD and H2/D2/HD + NH2. The effects of different vibrational mode excitations of the reagents on the reactivities are studied. Physical explanations about the isotope effects are also provided thoroughly including the influence of vibrational energy differences between the different isotopologues and the impact of the tunneling effect.

Dissociative photodetachment of H3O2-: a full-dimensional quantum dynamics study.

The transition state is a central concept of chemistry. Photoelectron-photofragment coincidence (PPC) spectroscopy has been proven as an attractive method to study the transition state dynamics. Within a state-of-the-art full-dimensional quantum mechanical model, the dissociative photodetachment dynamics of H3O2- is investigated on accurate anion and neutral potential energy surfaces. The calculated PPC spectrum of H3O2- agrees well with the experimental measurement. The dissociative product OH is exclusively populated on the ground vibrational state, implying the character of the spectator bond. In contrast, the product H2O is predominantly populated in the ground and fundamental states of the symmetric and antisymmetric stretching modes, which is caused by the strong coupling between the antisymmetric motion of the transferred H atom in the transient intermediate [H3O2]* and both stretching modes of the product H2O.

A rigorous full-dimensional quantum dynamics study of tunneling splitting of rovibrational states of vinyl radical C2H3.

We report a rigorous quantum mechanical study of the rovibrational energy levels of vinyl radical C2H3. The calculations are carried out using a real two-component multi-layer Lanczos algorithm in a set of orthogonal polyspherical coordinates based on a recently developed accurate ab initio potential energy surface of C2H3. All well converged 158 vibrational bands up to 3200 cm-1 are determined, together with a comparison to previous calculations and experimental results. Results show a remarkable multi-dimensional tunneling effect on the vibrational spectra of the radical. The vibrational tunneling splitting is substantially different from that of previous reduced dimensional calculations. The rotational constants of the fundamental vibrational bands of C2H3 are also given. It was found that the rovibrational states are strongly coupled, especially among those bending vibrational modes. In addition, the perturbative iteration approach of Gruebele has been extended to assign the rovibrational energy levels of C2H3 without the requirement of explicit wavefunctions.

An approximate full-dimensional quantum dynamics study of the mode specificity in the dissociative chemisorption of D2O on rigid Cu(111).

The mode specificity of the dissociative chemisorption of D2O on a rigid Cu(111) surface was investigated by calculating the approximate full-dimensional (9D) dissociation probabilities with the implementation of the seven-dimensional (7D) quantum dynamics calculations and site averaging approximations. The approximate 9D dissociation probabilities for D2O initially in various vibrational states were obtained by averaging the site-specific 7D results over 9 impact sites on an accurate 9D potential energy surface (PES). Strong mode specificity was observed for the title reaction, where vibrational excitations of the two stretching modes of D2O are more efficacious than increasing the translational energy in promoting the reaction, while the bending excitations with one and two quanta are less efficacious than the same amount of translational energy in enhancing the reactivity at low kinetic energies.

Communication: Equivalence between symmetric and antisymmetric stretching modes of NH3 in promoting H + NH3 → H2 + NH2 reaction.

Vibrational excitations of reactants sometimes promote reactions more effectively than the same amount of translational energy. Such mode specificity provides insights into the transition-state modulation of reactivity and might be used to control chemical reactions. We report here a state-of-the-art full-dimensional quantum dynamical study of the hydrogen abstraction reaction H + NH3 → H2 + NH2 on an accurate ab initio based global potential energy surface. This reaction serves as an ideal candidate to study the relative efficacies of symmetric and degenerate antisymmetric stretching modes. Strong mode specificity, particularly for the NH3 stretching modes, is demonstrated. It is further shown that nearly identical efficacies of the symmetric and antisymmetric stretching modes of NH3 in promoting the reaction can be understood in terms of local-mode stretching vibrations of the reactant molecule.

First-principles quantum dynamical theory for the dissociative chemisorption of H2O on rigid Cu(111).

Despite significant progress made in the past decades, it remains extremely challenging to investigate the dissociative chemisorption dynamics of molecular species on surfaces at a full-dimensional quantum mechanical level, in particular for polyatomic-surface reactions. Here we report, to the best of our knowledge, the first full-dimensional quantum dynamics study for the dissociative chemisorption of H2O on rigid Cu(111) with all the nine molecular degrees of freedom fully coupled, based on an accurate full-dimensional potential energy surface. The full-dimensional quantum mechanical reactivity provides the dynamics features with the highest accuracy, revealing that the excitations in vibrational modes of H2O are more efficacious than increasing the translational energy in promoting the reaction. The enhancement of the excitation in asymmetric stretch is the largest, but that of symmetric stretch becomes comparable at very low energies. The full-dimensional characterization also allows the investigation of the validity of previous reduced-dimensional and approximate dynamical models.

Full-dimensional quantum dynamics study of the H2 + C2H → H + C2H2 reaction on an ab initio potential energy surface

This work performs a time-dependent wavepacket study of the H2 + C2H → H + C2H2 reaction on a new ab initio potential energy surface (PES). The PES is constructed using neural network method based on 68 478 geometries with energies calculated at UCCSD(T)-F12a/aug-cc-pVTZ level and covers H2 + C2H↔H + C2H2, H + C2H2 → HCCH2, and HCCH2 radial isomerization reaction regions. The reaction dynamics of H2 + C2H → H + C2H2 are investigated using full-dimensional quantum dynamics method. The initial-state selected reaction probabilities are calculated for reactants in eight vibrational states. The calculated results showed that the H2 vibrational excitation predominantly enhances the reactivity while the excitation of bending mode of C2H slightly inhibits the reaction. The excitations of two stretching modes of C2H molecule have negligible effect on the reactivity. The integral cross section is calculated with J-shift approximation and the mode selectivity in this reaction is discussed. The rate constants over 200-2000 K are calculated and agree well with the experimental measured values.

Mode specificity of the dissociative chemisorption of HOD on rigid Cu(111): an approximate full-dimensional quantum dynamics study.

The validity of the site averaging approximation with exact potential (SAEP) has been confirmed in a recent work on the H2O/Cu(111) system [Z. Zhang, T. Liu, B. Fu, X. Yang, D. H. Zhang, Nat. Commun. 2016, 7, 11953]. Here, the mode specificity of the dissociative chemisorption of HOD on a rigid Cu(111) surface is investigated by carrying out the seven-dimensional (7D) quantum dynamics calculations on an accurate nine-dimensional (9D) potential energy surface together with the implementation of the SAEP. The approximate 9D dissociation probabilities for HOD initially in various vibrational states are obtained by averaging the site-specific 7D results over 9 impact sites. A strong bond-selective effect for the title reaction is observed, where vibrational excitation of a particular bond leads to a large enhancement only in the reaction in which the excited bond is broken. The product branching ratios strongly depend on which bond is excited, and the product from the cleavage of the excited bond is much more favored than the other product. The implementation of the SAEP allows us to investigate the mode-specific dynamics at a level of accuracy that can only be achieved in full-dimensional quantum dynamics calculations.

A full-dimensional quantum dynamics study of the mode specificity in the H + HOD abstraction reaction.

We employ the initial state-selected time-dependent wave packet approach to an atom-triatom reaction to study the H + HOD → OH + HD/OD + H2 reaction without the centrifugal sudden approximation, based on an accurate potential energy surface which was recently developed by neural network fitting to high level ab initio energy points. The total reaction probabilities and integral cross sections, which are the exact coupled-channel results, are calculated for the HOD reactant initially in the ground and several vibrationally excited states, including the bending excited state, OD stretching excited states, OH stretching excited states, and combined excitations of them. The reactivity enhancements from different initial states of HOD are presented, which feature strong bond-selective effects of the reaction dynamics. The current results for the product branching ratios, reactivity enhancements, and relative cross sections are largely improved over the previous calculations, in quantitatively good agreement with experiment. The thermal rate constant for the title reaction and the contributions from individual vibrational states of HOD are also obtained.

A full-dimensional multilayer multiconfiguration time-dependent Hartree study on the ultraviolet absorption spectrum of formaldehyde oxide.

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̃(1)A'←X̃(1)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.

A full-dimensional quantum dynamical study of H2+H2 collisions: Coupled-states versus close-coupling formulation

Collision-induced energy transfer involving H2 molecules plays an important role in many areas of physics. Kinetic models often require a complete set of state-to-state rate coefficients for H2+H2 collisions in order to interpret results from spectroscopic observations or to make quantitative predictions. Recent progress in full-dimensional quantum dynamics using the numerically exact close-coupling (CC) formulation has provided good agreement with existing experimental data for low-lying states of H2 and increased the number of state-to-state cross sections that may be reliably determined over a broad range of energies. Nevertheless, there exist many possible initial states (e.g., states with high rotational excitation) that still remain elusive from a computational standpoint even at relatively low collision energies. In these cases, the coupled-states (CS) approximation offers an alternative full-dimensional formulation. We assess the accuracy of the CS approximation for H2+H2 collisions by comparison with benchmark results obtained using the CC formulation. The results are used to provide insight into the orientation effects of the various internal energy transfer mechanisms. A statistical CS approximation is also investigated and cross sections are reported for transitions which would otherwise be impractical to compute.

Mode specificity in the H + H2O → H2 + OH reaction: A full-dimensional quantum dynamics study

The initial state-selected time-dependent wave packet approach to an atom-triatom reaction is employed to study the H + H2O → H2 + OH reaction without the centrifugal sudden approximation. The total reaction probabilities and integral cross sections, which are the exact coupled-channel results, are calculated for the H2O reactant initially in the ground and several vibrationally excited states, including bending excited states, first and second stretching excited states, and simultaneous excitations of both bending and stretching modes. The reactivity enhancements from different initial states of the H2O reagent are presented and discussed in detail. The thermal rate constant for the title reaction and the contributions to this coefficient from individual vibrational states of H2O are also obtained and compared with the previous theoretical and experimental data.

A multilayer MCTDH study on the full dimensional vibronic dynamics of naphthalene and anthracene cations

Employing the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method in conjunction with the multistate multimode vibronic coupling Hamiltonian model, we perform a full dimensional quantum dynamical study on the naphthalene (48D) and anthracene (66D) radical cations in their six lowest-lying doublet electronic states. For easily comparing results of full and reduced dimensionalities, MCTDH simulations based on larger sizes of primitive basis functions and single-particle functions than the previous ones [S. Ghanta, V. S. Reddy, and S. Mahapatra, Phys. Chem. Chem. Phys. 13, 14531 (2011)], are also performed. Extensive ML-MCTDH test calculations are performed to find appropriate ML separations of the wave functions (so-called ML-trees), and the convergence of the dynamical calculations are carefully checked. The ML-MCTDH method was developed for efficiently simulating quantum dynamics of large systems, and in fact the full dimensional ML-MCTDH calculations save a considerable amount of CPU-time in comparison with corresponding reduced dimensional MCTDH simulations. On basis of the present full and reduced dimensional simulations, the photoelectron (PE) spectra of these two cations are simulated and compared with corresponding experimental spectra. The agreement between theoretical and experimental PE spectra is good. Both full and reduced dimensional simulations give all main bands in the PE spectra. The vibronic energy-level positions from both ML-MCTDH and MCTDH calculations agree with corresponding experimental results. These quantum dynamical studies also complement the observations on diffuse interstellar bands with the wavelength of ~7088, ~6707, ~6490, ~6120, and ~5959 Å measured by astronomers as well as laboratory experimentalists.

Full-dimensional quantum dynamics study of exchange processes for the D + H2O and D + HOD reactions

The exchange processes of D + H(2)O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H(2)O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H(2)O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D(') + H(2)O → D(')OH + H and D(') + HOD → D(')OD + H reactions are substantially larger than those for the D(') + HOD → HOD(') + D reaction, indicating that the D(')/H exchange reactions are much more favored than the D(')/D exchange.

Quantum-dynamics study of the ${\rm H}_5^+$H5+ cluster: Full dimensional benchmark results on its vibrational states

A full-dimensional quantum dynamics study is carried out for the highly fluxional H(5)(+) cation on a recent reference potential energy surface by using the multi configuration time-dependent Hartree method. With five equivalent light atoms and shallow barriers between various low-lying stationary points on the surface, the spectroscopic characterization of H(5)(+) represents a huge challenge for accurate quantum dynamics simulations. The present calculation is the first such a study on this cation, which together with its isotope analogies are of primary importance in the interstellar chemistry. The vibrational ground state properties and several vibrationally excited states corresponding to low vibrational frequency motions, not yet directly observable by the experiment, are presented and analyzed.

Shape resonance in the H + D2O → D + HOD reaction: a full-dimensional quantum dynamics study

We present exact coupled-channel (CC) results of full-dimensional quantum mechanical calculations for the H + D2O → D + HOD exchange reaction which provide the first evidence for a shape resonance in this reaction, employing the initial state-selected time-dependent wave packet method with both OD bonds in the D2O reactant treated as reactive bonds. Due to the shallow C3v minimum along the reaction path, the signature of shape resonance is manifested as the existence of a distinct step-like feature in the CC integral cross sections of the exchange reaction just above the threshold. Our study also shows that the current CC integral cross sections are in very good agreement with the experimental results, and a bit larger than those calculated in the previous quantum dynamics study, where the one reactive bond and centrifugal sudden (CS) approximations were invoked, revealing that the one reactive bond approximation overestimated, while the CS approximation underestimated the integral cross sections.

Photodissociation of methyl iodide embedded in a host-guest complex: A full dimensional (189D) quantum dynamics study of CH3I@resorc[4]arene

Accurate full dimensional quantum dynamics calculations studying the photodissociation of CH(3)I@resorc[4]arene on an ab initio based potential energy surface (PES) model are reported. The converged 189D quantum dynamics calculations are facilitated by the multilayer multi-configurational time-dependent Hartree (ML-MCTDH) approach combined with the correlation discrete variable representation (CDVR) for the evaluation of potential energy matrix elements. The potential employed combines an established ab initio PES describing the photodissociation of methyl iodide in the A band with a harmonic description of the resorc[4]arene host and a bilinear modeling of the host-guest interaction. All potential parameters required in the description of the vibrations of the host molecule and the host-guest interaction are derived from ab initio calculations on the host-guest complex. Absorption spectra at 0 K and 300 K are calculated and the electronic population dynamics during the bond breaking process occurring in the first 20-30 fs after the photoexcitation is investigated. Weak but significant effects resulting from the host-guest interaction on this time scale are found and interpreted. The present study demonstrates that accurate fully quantum mechanical dynamics calculations can be preformed for systems consisting of more than 50 atoms using the ML-MCTDH/CDVR approach. Utilizing an efficient statistical approach for the construction of the ensemble of initial wavepackets, these calculations are not restricted to zero temperature but can also study the dynamics at 300 K.