Permutationally invariant potential energy surfaces in high dimensionality

A global coupled cluster potential energy surface for HCl + OH ↔ Cl + H2O

A new and more accurate full-dimensional global potential energy surface (PES) for the ground electronic state of the ClH₂O system is developed by using the permutation invariant polynomial-neural network (PIP-NN) method to fit 15 777 points obtained using an explicitly correlated unrestricted coupled-cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b).

Calculations of the active mode and energetic barrier to electron attachment to CF3 and comparison with kinetic modeling of experimental results

To provide a deeper understanding of the kinetics of electron attachment to CF₃, the six-dimensional potential energy surfaces of both CF₃ and CF₃^- were developed by fitting ∼3000 ab initio points per surface at the AE-CCSD(T)-F12a/AVTZ level using the permutation invariant polynomial-neural network (PIP-NN) approach. The fitted potential energy surfaces for CF₃ and CF₃^- had root mean square fitting errors relative to the ab initio calculations of 1.2 and 1.8 cm^-1, respectively. The main active mode for the crossing between the two potential energy surfaces was identified as the umbrella bending mode of CF₃ in C_3v symmetry. The lowest energy crossing point is located at R_CF = 1.306 Å and θ_FCF = 113.6° with the energy of 0.051 eV above the minimum of the CF₃ electronic surface. This value is only slightly larger than the experimental data 0.026 ± 0.01 eV determined by kinetic modeling of electron attachment to CF₃. The small discrepancy between the theoretical and experimentally measured values is analyzed.

Getting the Right Answers for the Right Reasons: Toward Predictive Molecular Simulations of Water with Many-Body Potential Energy Functions

The central role played by water in fundamental processes relevant to different disciplines, including chemistry, physics, biology, materials science, geology, and climate research, cannot be overemphasized. It is thus not surprising that, since the pioneering work by Stillinger and Rahman, many theoretical and computational studies have attempted to develop a microscopic description of the unique properties of water under different thermodynamic conditions. Consequently, numerous molecular models based on either molecular mechanics or ab initio approaches have been proposed over the years. However, despite continued progress, the correct prediction of the properties of water from small gas-phase clusters to the liquid phase and ice through a single molecular model remains challenging. To large extent, this is due to the difficulties encountered in the accurate modeling of the underlying hydrogen-bond network in which both number and strength of the hydrogen bonds vary continuously as a result of a subtle interplay between energetic, entropic, and nuclear quantum effects. In the past decade, the development of efficient algorithms for correlated electronic structure calculations of small molecular complexes, accompanied by tremendous progress in the analytical representation of multidimensional potential energy surfaces, opened the doors to the design of highly accurate potential energy functions built upon rigorous representations of the many-body expansion (MBE) of the interaction energies. This Account provides a critical overview of the performance of the MB-pol many-body potential energy function through a systematic analysis of energetic, structural, thermodynamic, and dynamical properties as well as of vibrational spectra of water from the gas to the condensed phase. It is shown that MB-pol achieves unprecedented accuracy across all phases of water through a quantitative description of each individual term of the MBE, with a physically correct representation of both short- and long-range many-body contributions. Comparisons with experimental data probing different regions of the water potential energy surface from clusters to bulk demonstrate that MB-pol represents a major step toward the long-sought-after "universal model" capable of accurately describing the molecular properties of water under different conditions and in different environments. Along this path, future challenges include the extension of the many-body scheme adopted by MB-pol to the description of generic solutes as well as the integration of MB-pol in an efficient theoretical and computational framework to model acid-base reactions in aqueous environments. In this context, given the nontraditional form of the MB-pol energy and force expressions, synergistic efforts by theoretical/computational chemists/physicists and computer scientists will be critical for the development of high-performance software for many-body molecular dynamics simulations.

Two Pathways for Dissociation of Highly Energized syn-CH3CHOO to OH Plus Vinoxy

Ozonolysis of alkenes is an important nonphotolytic source of hydroxl radicals in the troposphere. The reaction proceeds through cycloaddition and subsequent decomposition to a carbonyl oxide, known as Criegee intermediates. Ozonolysis of alkene releases about 50 kcal/mol excess energy to form highly energized Criegee molecules, which can be stabilized and undergo further reaction or dissociate to OH+vinoxy products. The dissociation dynamics of partially stabilized Criegee (syn-CH3CHOO) has been thoroughly studied recently, in which the molecules dissociate by first isomerizing to vinyl hydroperoxide (VHP). Here we examine the dissociation dynamics of highly energized syn-CH3CHOO (42 kcal/mol), and a second, prompt dissociation path is discovered. The dissociation dynamics of these two paths are carefully examined through the animation of trajectories and the energy distributions of products. The new prompt path reveals a distinctly different translational energy and internal energy distributions of products compared to the known path through VHP.

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.

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.

Global diabatic potential energy surfaces and quantum dynamical studies for the Li(2p) + H2(X(1)Σ(+)g) → LiH(X(1)Σ(+)) + H reaction

The global diabatic potential energy surfaces which are correlated with the ground state 1A' and the excited state 2A' of the Li(2p) + H2 reaction are presented in this study. The multi-reference configuration interaction method and large basis sets (aug-cc-pVQZ for H atom and cc-pwCVQZ for Li atom) were employed in the ab initio single-point energy calculations. The diabatic potential energies were generated by the diabatization scheme based on transition dipole moment operators. The neural network method was utilized to fit the matrix elements of the diabatic energy surfaces, and the root mean square errors were extremely small (3.69 meV for , 5.34 meV for and 5.06 meV for ). The topographical features of the diabatic potential energy surfaces were characterized and the surfaces were found to be sufficiently smooth for the dynamical calculation. The crossing seam of the conical intersections between the and surfaces were pinpointed. Based on this new analytical diabatic potential energy surfaces, time-dependent wave packet calculation were conducted to investigate the mechanism of the title reaction. At low collision energies, the product LiH molecule tends to forward scattering, while at high collision energies, the forward and backward scatterings exist simultaneously.

Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems

A new flavor of the frozen natural orbital complete active space second-order perturbation theory method (FNO-CASPT2, Aquilante et al., J. Chem. Phys. 131, 034113) is proposed herein. In this new implementation, the virtual space in Cholesky decomposition-based CASPT2 computations (CD-CASPT2) is truncated by excluding those orbitals that contribute the least toward preserving a predefined value of the trace of an approximate density matrix, as that represents a measure of the amount of dynamic correlation retained in the model. In this way, the amount of correlation included is practically constant at all nuclear arrangements, thus allowing for the computation of smooth electronic states surfaces and energy gradients-essential requirements for theoretical studies in photochemistry. The method has been benchmarked for a series of relevant biochromophores for which large speed-ups have been recorded while retaining the accuracy achieved in the corresponding CD-CASPT2 calculations. Both vertical excitation energies and gradient calculations have been carried out to establish general guidelines as to how much correlation needs to be retained in the calculation for the results to be consistent with the CD-CASPT2 findings. Our results feature errors within a tenth of an eV for the most difficult cases and have been validated to be used for gradient computations where an up to 3-fold speed-up is observed depending on the size of the system and the basis set employed.

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.

Calculations of the IR spectra of bend fundamentals of (H2O)n=3,4,5 using the WHBB_2 potential and dipole moment surfaces

Stimulated by new experiments from the Havenith group, we report IR spectra of the bend fundamentals of (H₂O)_n=3,4,5, using anharmonic, coupled-mode VSCF/VCI calculations, done in a subspace of modes consisting of all the monomer bends plus the hydrogen-bonded OH stretches.

An ab initio potential energy surface for the formic acid dimer: zero-point energy, selected anharmonic fundamental energies, and ground-state tunneling splitting calculated in relaxed 1–4-mode subspaces

We report a full-dimensional, permutationally invariant potential energy surface (PES) for the cyclic formic acid dimer.

An accurate multi-channel multi-reference full-dimensional global potential energy surface for the lowest triplet state of H2O2

A full dimensional potential energy surface for the lowest triplet state of H₂O₂ was developed at the MRCI-F12 level.

Quantum dynamics of polyatomic dissociative chemisorption on transition metal surfaces: mode specificity and bond selectivity

Recent advances in quantum dynamical characterization of polyatomic dissociative chemisorption on accurate global potential energy surfaces are critically reviewed.

Accurate nonadiabatic dynamics

This Perspective addresses the use of coupled diabatic potential energy surfaces (PESs) together with rigorous quantum dynamics in full or reduced dimensional coordinate spaces to obtain accurate solutions to problems in nonadiabatic dynamics.

Rovibrational transitions of the methane–water dimer from intermolecular quantum dynamical computations

Quantum dynamical computations of the methane–water dimer explain the far-infrared spectrum of this important prototype of a water–hydrocarbon interaction.

Pruning the Hamiltonian Matrix in MULTIMODE: Test for C2H4 and Application to CH3NO2 Using a New Ab Initio Potential Energy Surface

We report vibrational self-consistent field/virtual state configuration interaction energies of nitromethane using the code MULTIMODE and a new full-dimensional potential energy surface (PES). The PES is a precise, permutationally invariant linear least-squares fit to 17,049 electronic energies, using the CCSD(T)-F12b method with HaDZ basis (cc-pVDZ basis for H atoms, and aug-cc-pVDZ basis for C, O, N atoms). Nitromethane has 15 vibrational degrees of freedom, including one that is a nearly free internal methyl torsion, which is accurately described by the PES. This torsional mode makes vibrational calculations very challenging and here we present the results of calculations without it. Nevertheless, 14-mode calculations are still challenging and can lead to very large Hamiltonian matrices. To address this issue, we apply a pruning scheme, suggested previously by Handy and Carter, that reduces the size of matrix without sacrificing accuracy in the eigenvalues. The method is briefly described here in the context of partitioning theory. A new and more efficient implementation of it, coded in the latest version of MULTIMODE program, is described. The accuracy and efficiency are demonstrated for 12-mode C2H4 and then applied to CH3NO2. Agreement with available experimental values of the CH3NO2 14 fundamentals is very good. Diffusion Monte Carlo calculations in full dimensionality are done for the zero-point energy and wavefuction. These indicate that the torisonal motion is nearly a free-rotor in this state.