A quasi-classical study in a quantum spirit of mode specificity of the H + HOD abstraction reaction

Quantum benchmark calculations of the H + HOD abstraction reaction [B. Zhao et al., Phys. Chem. Chem. Phys., 2018, 20, 17029-17037] provide an opportunity to test approximate methods, such as quasi-classical trajectories (QCTs) with Gaussian binning. However, the large mode-specific enhancements of this reaction lead to special challenges and unphysical QCT cross-section artifacts. We propose and apply a general backward-forward-backward (BFB) trajectory procedure to avoid cross-section artifacts arising from vibrational motion delocalization of the initially prepared reactants. We also develop and apply a general hybrid weighting scheme, in which vibrationally adiabatic products receive unit weights, while vibrationally non-adiabatic products are Gaussian weighted. Motivated by examination of the product vibrational actions, the hybrid weighting is necessary to avoid large cross-section errors. These QCTs with Gaussian weighting extensions are general and are recommended for future polyatomic reaction simulations.

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.

Femtochemistry of bimolecular reactions from weakly bound complexes: computational study of the H + H'OD → H'OH + D or HOD + H' exchange reactions

A full-dimensional wavepacket propagation describing the bimolecular exchange reactions H + H'OD → H'OH + D or HOD + H' initiated by photolysis of HCl in the hydrogen-bound complex (HCl)⋯(HOD) is reported. The dynamics of this reaction is carried out with the MCTDH method on an ab initio potential energy surface (PES) of H₃O and the initial state is derived from the ground state wavefunction of the complex obtained by relaxation on its own electronic ground state ab initio PES. The description of the system makes use of polyspherical coordinates parametrizing a set of Radau and Jacobi vectors. The calculated energy- and time-resolved reaction probabilities show, owing to the large collision energies at play stemming from the (almost full) photolysis of HCl, that the repulsion between oxygen in the H'OD molecule and the incoming hydrogen atom is the main feature of the collision and leads to non-reactive scattering. No abstraction reaction products are observed. However, both exchange processes are still observable, with a preference in O-H' bond dissociation over that of O-D. The selectivity is reversed upon vibrational pre-excitation of the O-D stretching mode in the H'OD molecule. It is shown that, after the collision, the hydrogen atom of HCl does most likely not encounter the almost stationary chlorine atom again but we also consider the limit case where the H atom is forced to collide multiple times against H'OD as a result of being pushed back by the Cl atom.

Feshbach resonances in D + HD(v = 1, j = 0) reaction at low collision energies

Feshbach resonances in D + HD(v = 1, j = 0) reaction are studied by using the time-independent quantum method. The integral cross section (ICS) results present three Feshbach resonance peaks, which are different from H + HD(v = 1, j = 0) reaction dominated by only one peak. These resonances are attributed to coupling with adiabatic effective potentials of D + HD(v = 1, j = 1) reaction, and the most obvious peak is contributed by J = 1 at 83.16 cm^-1 collision energy. For J = 0 and 2, the resonances are related with the same L partial wave and present a double-peak structure in total ICS. The characteristics of product angular distribution show that the resonance of J = 1 is long-lived, while the lifetimes are relatively shorter for the resonance of J = 0 and 2.

Mode- and Bond-Selected Reaction of H with Local Mode Molecule HDS

Understanding mode- and bond-selected dynamics of elementary chemical reactions is of central importance in molecular reaction dynamics. The initial state-selected time-dependent wave packet method is employed to study the mode and bond selectivity, isotopic branching ratio, and temperature dependence of rate constants of the two-channel reaction of H with local mode molecule HDS. For the abstraction channel, fundamental excitation of the HS (DS) bond of the reactant HDS significantly enhances the H-abstraction (D-abstraction) reaction, whose efficacy is higher than the same amount of translational energy except at low energies just above the energy threshold. This is in sharp contrast to the prediction of Polanyi rules: translational energy is more efficient than vibrational energy in enhancing a reaction with an early barrier. The recent sudden vector projection model is then applied to rationalize the observed mode specificity, which, however, shows that the translational mode vector has a larger coupling with the reaction coordinate than the stretching vector of the active bond, implying a reversed relative efficacy on promoting the reaction as well. In contrast, the mode and bond specificity for the exchange channel is not as strong as for the abstraction channel due to the regulation of the shallow well along the reaction path.

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.

Recent Advances in Quantum Dynamics of Bimolecular Reactions

In this review, we survey the latest advances in theoretical understanding of bimolecular reaction dynamics in the past decade. The remarkable recent progress in this field has been driven by more accurate and efficient ab initio electronic structure theory, effective potential-energy surface fitting techniques, and novel quantum scattering algorithms. Quantum mechanical characterization of bimolecular reactions continues to uncover interesting dynamical phenomena in atom-diatom reactions and beyond, reaching an unprecedented level of sophistication. In tandem with experimental explorations, these theoretical developments have greatly advanced our understanding of key issues in reaction dynamics, such as microscopic reaction mechanisms, mode specificity, product energy disposal, influence of reactive resonances, and nonadiabatic effects.