Isotope Effect and Heavy-Light-Heavy Reactivity Oscillation in the Cl + CHD3/CHT3 Reaction.

Previous experiments and theories have shown the existence of heavy-light-heavy (HLH) reactivity oscillation in the Cl + CH4 reaction and anticipated that similar oscillations should exist in many HLH reactions involving polyatomic reagents. However, the total reaction probabilities for the Cl + CHD3 → HCl + CD3 reaction exhibit only a step-like feature, and the total reaction probabilities for Cl + CHT3 → HCl + CT3 do not show any structure at all. Here, we report seven-dimensional state-to-state quantum dynamics studies for this reaction on the FI-NN PES, and we demonstrate that HLH reactivity oscillations also exist in these two reactions, manifesting as peaks in the reaction probabilities for low product rotational states. These oscillations, however, are obscured in the total reaction probability because of the higher excitation of j ≥ 2 product rotational states. Furthermore, the isotope replacement of nonreactive hydrogen with deuterium and tritium significantly enhances reactivity at collision energies above 0.112 eV, indicating an inverse secondary isotope effect on the probabilities, which is proved to be also caused by HLH mass combination. We also demonstrate that the highly rotational excitation of CHD3 substantially enhances reactivity and the HLH oscillations, similar to HLH triatomic reactions. These observations are completely different from those in the H + CHD3 reaction, which is also a late-barrier reaction. Therefore, the HLH mass combination is very important, which affects not only the reactivity oscillation but also the amplitude and product rotational state distribution and makes the initial rotation excitation play a pivotal role in the reaction.

Supercollisions of NaCl + NaCl on an Accurate Full-Dimensional Potential Energy Surface.

An accurate, global, full-dimensional potential energy surface (PES) of NaCl + NaCl has been constructed by the fundamental invariant-neural network (FI-NN) fitting based on roughly 13,000 ab initio energies at the level of CCSD(T)-F12a/aug-cc-pVTZ, with the small fitting error of 0.16 meV. Extensive quasiclassical trajectory (QCT) calculations were performed on this PES to investigate the energy transfer process of the NaCl + NaCl collision at four different collision energies. Various quantities were obtained, including the cross-sections, energy transfer probability, average energy transfer, and collision lifetime. The probabilities of energy transfer (P(ΔE)) for prompt trajectories, nonreactive trajectories, and reactive trajectories deviate from a simple exponential decay pattern. Instead, a noteworthy probability is observed in the high-energy transfer region, indicative of supercollisions. The formation of the (NaCl)2 complex, coupled with a comparatively extended collision lifetime, promotes vibrational excitation in NaCl molecules. The reactive trajectories exhibit enhanced energy transfer, attributed to the longer lifetime of the NaCl dimer. This study not only provides an accurate and extensive understanding of the NaCl + NaCl collision dynamics but also reveals intriguing phenomena, such as supercollisions and enhanced energy transfer in reactive trajectories, shedding light on the complex intricacies of molecular interactions.

High vibrational excitation of the reagent transforms the late-barrier H + HOD reaction into an early-barrier reaction.

Polanyi's rules predict that a late-barrier reaction yields vibrationally cold products; however, experimental studies showed that the H2 product from the late-barrier H + H2O(|04⟩-) and H + HOD(vOH = 4) reactions is vibrationally hot. Here, we report a potential-averaged five-dimensional state-to-state quantum dynamics study for the H + HOD(vOH = 0-4) → H2 + OD reactions on a highly accurate potential energy surface with the total angular momentum J = 0. It is found that with the HOD vibration excitation increasing from vOH = 1 to 4, the product H2 becomes increasingly vibrationally excited and manifests a typical characteristic of an early barrier reaction for vOH = 3 to 4. Analysis of the scattering wave functions revealed that vibrational excitation in the breaking OH bond moves the location of dynamical saddle point from product side to reactant side, transforming the reaction into an early barrier reaction. Interestingly, pronounced oscillatory structures in the total and product vibrational-state-resolved reaction probabilities were observed for the H + HOD(vOH = 3, 4) reactions, in particular at low collision energies, which originate from the Feshbach resonance states trapped in the bending/torsion excited vibrational adiabatic potential wells in the entrance region due to van der Waals interactions.

Bond-selective effect for the dissociative chemisorption of HOD on the Ni(100) surface revealed at the full-dimensional quantum dynamical level.

We present a comprehensive investigation into the dissociative chemisorption of HOD on a rigid Ni(100) surface using an approximate full-dimensional (9D) quantum dynamics approach, which was based on the time-dependent wave-packet calculations on a full-dimensional potential energy surface obtained through neural network fitting to density functional theory energy points. The approximate-9D probabilities were computed by averaging the seven-dimensional (7D) site-specific dissociation probabilities across six impact sites with appropriate relative weights. Our results uncover a distinctive bond-selective effect, demonstrating that the vibrational excitation of a specific bond substantially enhances the cleavage of that excited bond. The product branching ratios are substantially influenced by which bond undergoes excitation, exhibiting a clear preference for the product formed through the cleavage of the excited bond over the alternative product.

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.

Six-dimensional quantum dynamics study for the dissociative chemisorption of H2 on pure and alloyed AgAu surfaces.

The 6D time-dependent wave packet calculations were performed to explore H2 dissociation on Ag, Au, and two AgAu alloy surfaces, using four newly fitted potential energy surfaces based on the neural network fitting to density functional theory energy points. The ligand effect resulting from the Ag-Au interaction causes a reduction in the barrier height for H2+Ag/Au(111) compared to H2+Ag(111). However, the scenario is reversed for H2+Au/Ag(111) and H2+Au(111). The 6D dissociation probabilities of H2 on Ag/Au(111) surfaces are significantly higher than those on the pure Ag(111) surface, but the corresponding results for H2 on Au/Ag(111) surfaces are substantially lower than those on the pure Au(111) surface. The reactivity of H2 on Au(111) is larger than that on Ag(111), despite Ag(111) having a slightly lower static barrier height. This can be attributed to the exceptionally small dissociation probabilities at the hcp and fcc regions, which are at least 100 times smaller compared to those at the bridge or top site for H2+Ag(111). Due to the late barrier being more pronounced, the vibrational excitation of H2 on Ag(111) is more effective in promoting the reaction than on Au(111). Moreover, a high degree of alignment dependence is detected for the four reactions, where the H2 dissociation has the highest probability at the helicopter alignment, as opposed to the cartwheel alignment.

Quantum Tunneling in Peroxide O-O Bond Breaking Reaction.

The importance of quantum-mechanical tunneling becomes increasingly recognized in chemical reactions involving hydrogen as well as heavier atoms. Here we report concerted heavy-atom tunneling in an oxygen-oxygen bond breaking reaction from cyclic beryllium peroxide to linear dioxide in cryogenic Ne matrix, as evidenced by subtle temperature-dependent reaction kinetics and unusually large kinetic isotope effects. Furthermore, we demonstrate that the tunneling rate can be tuned through noble gas atom coordination on the electrophilic beryllium center of Be(O2), as the half-life dramatically increased from 0.1 h for NeBe(O2) at 3 K to 12.8 h for ArBe(O2). Quantum chemistry and instanton theory calculations reveal that noble gas coordination notably stabilizes the reactants and transition states, increases the barrier heights and widths, and consequently reduces the reaction rate drastically. The calculated rates and in particular kinetic isotope effects are in good agreement with experiment.

Neural network potential energy surfaces and dipole moment surfaces for SO2(H2O) and SO2(H2O)2 complexes.

Full-dimensional, ab initio-based many-body potential energy surfaces and dipole moment surfaces constructed using the neural network method for SO2(H2O)n (n = 1,2) complexes are reported. The database of the SO2 1-body PES, SO2(H2O) 2-body PES and SO2(H2O)2 3-body PES consists of 11 952, 79 882 and 84 159 ab initio energies, respectively. All 1-body energies were calculated at the CCSD(T)/CBS(AVTZ:AVQZ) level and all 2,3-body energies were calculated at the DSD-PBEP86/AVTZ level. The database of DMSs is the same as that of PESs and all dipole moments were calculated at the MP2/AVTZ level. Harmonic frequencies and dissociation energies of SO2(H2O) and SO2(H2O)2 were calculated on these PESs and compared with ab initio results to examine the fidelity of these PESs.

A highly accurate full-dimensional ab initio potential surface for the rearrangement of methylhydroxycarbene (H3C-C-OH).

We report here a full-dimensional machine learning global potential surface (PES) for the rearrangement of methylhydroxycarbene (H3C-C-OH, 1t). The PES is trained with the fundamental invariant neural network (FI-NN) method on 91 564 ab initio energies calculated at the UCCSD(T)-F12a/cc-pVTZ level of theory, covering three possible product channels. FI-NN PES has the correct symmetry properties with respect to permutation of four identical hydrogen atoms and is suitable for dynamics studies of the 1t rearrangement. The averaged root mean square error (RMSE) is 11.4 meV. Six important reaction pathways, as well as the energies and vibrational frequencies at the stationary geometries on these pathways are accurately preproduced by our FI-NN PES. To demonstrate the capacity of the PES, we calculated the rate coefficient of hydrogen migration in -CH3 (path A) and hydrogen migration of -OH (path B) with instanton theory on this PES. Our calculations predicted the half-life of 1t to be 95 min, which is excellent in agreement with experimental observations.

Differential Cross Sections for the H + H2S → H2 + SH Reaction: A Full-Dimensional State-to-State Quantum Dynamics Study.

We utilized the time-dependent wave packet approach to compute the first full-dimensional (6D) state-to-state differential cross sections (DCSs) for the title reaction with the initial nonrotating H2S in the ground and the (100) and (001) vibrational excited states. It is found that the fundamental symmetric and asymmetric stretching excitations of H2S exhibit almost the same influence on the DCS, but unlike the H + H2O → H2 + OH reaction, they greatly increase the vibrational excitation of both the H2 and SH products. The hot vibrational state distributions of H2 are consistent with the prediction of product energy disposal by Polanyi's rules for an early barrier reaction. Because the incident H atom reacts strongly with both the ground and excited S-H states, the large populations of product SH(v2 = 1), which are very close to the relative reactivity of the initial S-H(v = 0) state, can still be explained by the local mode picture for H2S and the nonreacting SH bond's spectator nature.

Roaming Dynamics in Hydroxymethyl Hydroperoxide Decomposition Revealed by the Full-Dimensional Potential Energy Surface of the CH2OO + H2O Reaction.

The CH2OO + H2O reaction is an important atmospheric process that leads to the formation of formic acid (HCOOH) and water via the intermediate hydroxymethyl hydroperoxide (HOCH2OOH, HMHP). We investigated the intricacies of this process by employing quasiclassical trajectory calculations on an accurate, full-dimensional ab initio potential energy surface (PES). In addition to the direct mechanism via the transition state (TS), an interesting roaming mechanism was found to play the predominant role in producing H2O and HCOOH. This roaming pathway is featured as the near direct dissociation of HMHP into OH and hydroxymethoxy radical, followed by the retraction of OH and abstraction of the H atom, culminating in the formation of H2O. Due to the longer interaction time of the roaming mechanism, less product translational energy was released, but more internal energies of HCOOH were obtained, as compared with the direct TS mechanism. The enhanced yield of H2O and formic acid achieved through roaming dynamics underscores the significance of dynamics simulations based on an accurate full-dimensional PES. This work provides new insights into the dynamics of the CH2OO + H2O reaction and its implications for atmospheric chemistry.

Semiclassical Vibrational Spectroscopy of Real Molecular Systems by Means of Cross-Correlation Filter Diagonalization.

We applied the harmonic inversion technique to extract vibrational eigenvalues from the semiclassical initial value representation (SC-IVR) propagator of molecular systems described by explicit potential surfaces. The cross-correlation filter-diagonalization (CCFD) method is used for the inversion problem instead of the Fourier transformation, which allows much shorter propagation time and is thus capable of avoiding numerical divergence issues while getting rid of approximations like the separable one to the pre-exponential factor. We also used the "Divide-and-Conquer" technique to control the total dimensions under consideration, which helps to further enhance the numerical behavior of SC-IVR calculations and the stability of harmonic inversion methods. The technique is tested on small molecules and water trimer to justify its applicability and reliability. Results show that the CCFD method can effectively extract the vibrational eigenvalues from short trajectories and reproduce the original spectra conventionally obtained from long-time ones, with no loss on accuracy while the numerical behavior is much better. This work demonstrates the possibility to apply the combined method of CCFD and SC-IVR to real molecular potential surfaces, which might be a new way to overcome the numerical instabilities caused by the increase of dimensions.

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.

Accurate calculation of tunneling splittings in water clusters using path-integral based methods.

Tunneling splittings observed in molecular rovibrational spectra are significant evidence for tunneling motion of hydrogen nuclei in water clusters. Accurate calculations of the splitting sizes from first principles require a combination of high-quality inter-atomic interactions and rigorous methods to treat the nuclei with quantum mechanics. Many theoretical efforts have been made in recent decades. This Perspective focuses on two path-integral based tunneling splitting methods whose computational cost scales well with the system size, namely, the ring-polymer instanton method and the path-integral molecular dynamics (PIMD) method. From a simple derivation, we show that the former is a semiclassical approximation to the latter, despite that the two methods are derived very differently. Currently, the PIMD method is considered to be an ideal route to rigorously compute the ground-state tunneling splitting, while the instanton method sacrifices some accuracy for a significantly smaller computational cost. An application scenario of such a quantitatively rigorous calculation is to test and calibrate the potential energy surfaces of molecular systems by spectroscopic accuracy. Recent progress in water clusters is reviewed, and the current challenges are discussed.

Reactivity oscillation in the heavy-light-heavy Cl + CH4 reaction.

It has long been predicted that oscillatory behavior exists in reactivity as a function of collision energy for heavy-light-heavy (HLH) chemical reactions in which a light atom is transferred between two heavy atoms or groups of atoms, but direct observation of such a behavior in bimolecular reactions remains a challenge. Here we report a joint theoretical and crossed-molecular-beam study on the Cl + CH4 → HCl + CH3 reaction. A distinctive peak at a collision energy of 0.15 eV for the CH3(v = 0) product was experimentally detected in the backward scattering direction. Detailed quantum-dynamics calculations on a highly accurate potential energy surface revealed that this feature originates from the reactivity oscillation in this HLH polyatomic reaction. We anticipate that such reactivity oscillations exist in many HLH reactions involving polyatomic reagents.

Infrared spectroscopy of neutral water clusters at finite temperature: Evidence for a noncyclic pentamer.

Infrared spectroscopic study of neutral water clusters is crucial to understanding of the hydrogen-bonding networks in liquid water and ice. Here we report infrared spectra of size-selected neutral water clusters, (H2O) n (n = 3-6), in the OH stretching vibration region, based on threshold photoionization using a tunable vacuum ultraviolet free-electron laser. Distinct OH stretch vibrational fundamentals observed in the 3,500-3,600-cm-1 region of (H2O)5 provide unique spectral signatures for the formation of a noncyclic pentamer, which coexists with the global-minimum cyclic structure previously identified in the gas phase. The main features of infrared spectra of the pentamer and hexamer, (H2O) n (n = 5 and 6), span the entire OH stretching band of liquid water, suggesting that they start to exhibit the richness and diversity of hydrogen-bonding networks in bulk water.

Torsional Tunneling Splitting in a Water Trimer.

Using a full-dimensional quantum method for nuclei and a new first-principles water potential, we show that the torsional tunneling splitting in a water trimer can be reproduced with accuracy up to ∼1 cm-1. We quantify the coupling constants of the nuclear quantum states between nonadjacent wells and show that they are the main reason for shifting the quartet-split levels in spectra from a 1:2:1 spacing. This demonstrates the limitation of treatments using simplified models such as the Hückel model and emphasizes the nonlocal nature of the quantum interactions in this system. With such an ab initio endeavor, we examine the quality of the water potential developed and provide a rigorous scheme to decipher the experimental spectra with unprecedented accuracy, which is applicable to more general systems.

Permutation invariant polynomial neural network based diabatic ansatz for the (E + A) × (e + a) Jahn-Teller and Pseudo-Jahn-Teller systems.

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of H3 + is studied, which is an (E + A) × (e + a) Jahn-Teller and Pseudo-Jahn-Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. The adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries.

Quantum dynamics reveal different ligand effects by vibrational excitation in the dissociative chemisorption of HCl on the Au/Ag(111) surface.

The reactivity and selectivity of bimetallic surfaces are of fundamental importance in industrial applications. Here, we report the first six-dimensional (6D) quantum dynamics study for the role of surface strain and ligand effects on the reactivity of HCl on a strained pseudomorphic monolayer of Au deposited onto a Ag(111) substrate, with the aid of accurate machine learning-based potential energy surfaces. The substitute of Au into Ag changes the location of the transition state; however, the static barrier height remains roughly the same as pure Au(111). The 6D quantum dynamics calculations reveal that the surface strain due to lattice expansion slightly enhances the reactivity. The ligand effect due to electronic structure interactions between Au and Ag substantially suppresses the reactivity of HCl in the ground vibrational state but promotes the reactivity via vibrational excitation at high kinetic energies. This finding can be attributed to more close interaction with Ag atoms at the transition state close to the fcc site, as well as the tight transition-state region, making the vibrational excitation highly efficient in enhancing the reactivity. Our study quantitatively unravels the dynamical origin of reactivity control by two metals, which will ultimately provide valuable insight into the selectivity of the catalyst.

A fundamental invariant-neural network representation of quasi-diabatic Hamiltonians for the two lowest states of H3.

The fundamental invariant neural network (FI-NN) approach is developed to represent coupled potential energy surfaces in quasidiabatic representations with two-dimensional irreducible representations of the complete nuclear permutation and inversion (CNPI) group. The particular symmetry properties of the diabatic potential energy matrix of H3 for the 1A' and 2A' electronic states were resolved arising from the E symmetry in the D3h point group. This FI-NN framework with symmetry adaption is used to construct a new quasidiabatic representation of H3, which reproduces accurately the ab initio energies and derivative information with perfect symmetry behaviors and extremely small fitting errors. The quantum dynamics results on the new FI-NN diabatic PESs give rise to accurate oscillation patterns in the product state-resolved differential cross sections. These results strongly support the accuracy and efficiency of the FI-NN approach to construct reliable diabatic representations with complicated symmetry problems.