Seven-dimensional quantum dynamics study of the H+NH3-->H2+NH2 reaction

An efficient and universal parallel algorithm for high-dimensional quantum dynamics in poly-atomic reactions

This study presents a parallel algorithm for high-dimensional quantum dynamics simulations in poly atomic reactions, integrating distributed- and shared-memory models. The distributions of the wave function and potential energy matrix across message passing interface processes are based on bundled radial and angular dimensions, with implementations featuring either two- or one-sided communication schemes. Using realistic parameters for the H + NH3 reaction, performance assessment reveals linear scalability, exceeding 90% efficiency with up to 600 processors. In addition, owing to the universal and concise structure, the algorithm demonstrates remarkable extensibility to diverse reaction systems, as demonstrated by successes with six-atom and four-atom reactions. This work establishes a robust foundation for high-dimensional dynamics studies, showcasing the algorithm's efficiency, scalability, and adaptability. The algorithm's potential as a valuable tool for unraveling quantum dynamics complexities is underscored, paving the way for future advancements in the field.

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 NH₃ 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 NH₃ 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 NH₃ have comparable efficacies on driving the reaction.

Recent advances in quantum scattering calculations on polyatomic bimolecular reactions

This review surveys quantum scattering calculations on chemical reactions of polyatomic molecules in the gas phase published in the last ten years.

A nine-dimensional global potential energy surface for NH4(X2A1) and kinetics studies on the H + NH3 ↔ H2 + NH2 reaction

A nine-dimensional global potential energy surface (PES) for the NH₄ system is developed from ∼10⁵ high-level ab initio points and the hydrogen abstraction kinetics on the PES agree with experiment.

SEVEN-DIMENSIONAL QUANTUM DYNAMICS STUDY OF THE H2 + NH2 → H + NH3 REACTION ON AN INTERPOLATED POTENTIAL ENERGY SURFACE

Initial-state-selected time-dependent wave packet dynamics studies have been performed for the H 2 + NH 2 → H + NH 3 reaction with a seven-dimensional model on a new interpolated ab initio potential energy surface (PES). The PES is constructed using modified Shepard interpolation Scheme and contains 1967 data points with ab initio calculations carried out on UCCSD(T)/aug-cc-pVTZ level. In the seven-dimensional model, NH 2 group keeps C2v symmetry and two NH bonds are fixed at their equilibrium values. The total reaction probabilities are calculated when (1) the two reactants are initially at their ground states; (2) NH 2 bending mode is excited, and (3) H 2 is on its first vibrational excited state. The integral cross sections are also reported for these initial states with centrifugal-sudden approximation. Thermal rate constants are calculated for the temperature range of 200–2000 K and compared with the previous calculated values and available experimental data. Good agreements between theory and experiments for the rate constants at intermediate temperature are achieved on this PES.

Constructing Potential Energy Surfaces for Polyatomic Systems: Recent Progress and New Problems

Different methods of constructing potential energy surfaces in polyatomic systems are reviewed, with the emphasis put on fitting, interpolation, and analytical (defined by functional forms) approaches, based on quantum chemistry electronic structure calculations. The different approaches are reviewed first, followed by a comparison using the benchmark H + CH4 and the H + NH3 gas-phase hydrogen abstraction reactions. Different kinetics and dynamics properties are analyzed for these reactions and compared with the available experimental data, which permits one to estimate the advantages and disadvantages of each method. Finally, we analyze different problems with increasing difficulty in the potential energy construction: spin-orbit coupling, molecular size, and more complicated reactions with several maxima and minima, which test the soundness and general applicability of each method. We conclude that, although the field of small systems, typically atom-diatom, is mature, there still remains much work to be done in the field of polyatomic systems.

QCT and QM calculations of the Cl(2P) + NH3 reaction: influence of the reactant well on the dynamics

A detailed dynamics study, using both quasi-classical trajectory (QCT) and reduced-dimensional quantum mechanical (QM) calculations, was carried out to understand the reactivity and mechanism of the Cl((2)P) + NH(3)→ HCl + NH(2) gas-phase reaction, which evolves through deep wells in the entry and exit channels. The calculations were performed on an analytical potential energy surface recently developed by our group, PES-2010 [M. Monge-Palacios, C. Rangel, J. C. Corchado and J. Espinosa-Garcia, Int. J. Quantum. Chem., 2011], together with a simplified model surface, mod-PES, in which the reactant well is removed to analyze its influence. The main finding was that the QCT and QM methods show a change of the reaction probability with collision energy, suggesting a change of the atomic-level mechanism of reaction with energy. This change disappeared when the mod-PES was used, showing that the behaviour at low energies is a direct consequence of the existence of the reactant well. Analysis of the trajectories showed that different mechanisms operate depending on the collision energy. Thus, while at high energies (E(coll) > 5 kcal mol(-1)) practically all trajectories are direct, at low energies (E(coll) < 3 kcal mol(-1)) the trajectories are indirect, i.e., with the mediation of a trapping complex in the entry and/or the exit wells. The reactant complex allows repeated encounters between the reactants, increasing the reaction probability at low energies. The differential cross section results reinforce this change of mechanism, showing also the influence of the reactant well on this reaction. Thus, the PES-2010 surface yields a forward-backward symmetry in the scattering, while when the reactant well is removed with the mod-PES the shape is more isotropic.

COMPARATIVE STUDY OF REACTION RATE CONSTANTS FOR THE NH3 + H → NH2 + H2 REACTION WITH GLOBE DYNAMICS AND TRANSITION STATE THEORIES

The nine-dimension quasi-classical trajectory (QCT) calculations have been carried out for the title reaction with a global potential energy surface (PES) constructed by Corchado and Espinosa-García (J Chem Phys106:4013, 1997). The detailed dynamics calculations cover the specific collision energies falling in the range of 0.62–3.04 eV, which are sufficient to fit the calculated reactive cross-sections into a barrier-type excitation function and to obtain the thermal rate constants. The present QCT rate constants are in good agreement with the recent quantum dynamics (QD) results, both of which are much lower than that of the previous variational transition state theory (VTST).

Kinetics and dynamics of the NH3 + H → NH2 + H2 reaction using transition state methods, quasi-classical trajectories, and quantum-mechanical scattering

On a recent analytical potential energy surface developed by two of the authors, an exhaustive kinetics study, using variational transition state theory with multidimensional tunneling effect, and dynamics study, using both quasi-classical trajectory and full-dimensional quantum scattering methods, was carried out to understand the reactivity of the NH(3) + H → NH(2) + H(2) gas-phase reaction. Initial state-selected time-dependent wave packet calculations using a full-dimensional model were performed, where the total reaction probabilities were calculated for the initial ground vibrational state and for four excited vibrational states of ammonia. Thermal rate constants were calculated for the temperature range 200-2000 K using the three methods and compared with available experimental data. We found that (a) the total reaction probabilities are very small, (b) the symmetric and asymmetric N-H stretch excitations enhance the reactivity, (c) the quantum-mechanical calculated thermal rate constants are about one order of magnitude smaller than the transition state theory results, which reproduce the experimental evidence, and (d) quasi-classical trajectory calculations, which were performed with the main goal of analyzing the influence of the zero-point energy problem on the final dynamics results, reproduce the quantum scattering calculations on the same surface.

A quantum reaction dynamics study of the translational, vibrational, and rotational motion effects on the HD + H3+ reaction

Time-dependent, quantum reaction dynamics wavepacket approach is employed to investigate the impacts of the translational, vibrational, and rotational motion on the HD+H(3)(+) → H(2)D(+) + H(2) reaction using the Xie-Braams-Bowman potential energy surface [Z. Xie, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 122, 224307 (2005)]. We treat this five atom reaction with a seven-degree-of-freedom model by fixing one Jacobi and one torsion angle related to H(3) (+) at the lowest saddle point geometry of the potential energy surface. The initial state selected reaction probabilities show that the rotational excitations of H(+)-H(2) greatly enhance the reactivity with the reaction probabilities increased double at high rotational states compared to the ground state. However, the vibrational excitations of H(3) (+) hinder the reactivity. The ground state reaction probability shows no reaction threshold for this exoergic reaction, and as the translational energy increases, the reaction probability decreases. Furthermore, reactive resonances and zero point energy play very important roles on the reaction dynamics. The obtained integral cross section has the character of an exoergic reaction without a threshold: it decreases with the translational energy increasing. The calculated thermal rate constants using this seven-degree-of-freedom model are in agreement with a later experiment measurement.

Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.