Permutationally invariant potential energy surfaces in high dimensionality

Trapped Abstraction in the O((1)D) + CHD3 → OH + CD3 Reaction

Despite significant progress made in past decades, it is still challenging to elucidate dynamics mechanisms for polyatomic reactions, in particular, involving complex formation. The reaction of O((1)D) with methane has long been regarded as a prototypical polyatomic system of direct insertion reaction in which the O((1)D) atom can insert into the C-H bond of methane to form a "hot" methanol intermediate before decomposition. Here, we report a combined theoretical and experimental study on the O((1)D) + CHD3 reaction, on which good agreement between theory and experiment is achieved. Our study revealed that this complex-forming reaction actually proceeds via a trapped abstraction mechanism, rather than an insertion mechanism as has long been thought. We anticipate that this reaction mechanism should also be responsible for the reaction of O((1)D) with ethane and propane, as well as many other chemical reactions with deep wells in the interaction region.

Quantum mechanical fragment methods based on partitioning atoms or partitioning coordinates

Conspectus The development of more efficient and more accurate ways to represent reactive potential energy surfaces is a requirement for extending the simulation of large systems to more complex systems, longer-time dynamical processes, and more complete statistical mechanical sampling. One way to treat large systems is by direct dynamics fragment methods. Another way is by fitting system-specific analytic potential energy functions with methods adapted to large systems. Here we consider both approaches. First we consider three fragment methods that allow a given monomer to appear in more than one fragment. The first two approaches are the electrostatically embedded many-body (EE-MB) expansion and the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE), which we have shown to yield quite accurate results even when one restricts the calculations to include only electrostatically embedded dimers. The third fragment method is the electrostatically embedded molecular tailoring approach (EE-MTA), which is more flexible than EE-MB and EE-MB-CE. We show that electrostatic embedding greatly improves the accuracy of these approaches compared with the original unembedded approaches. Quantum mechanical fragment methods share with combined quantum mechanical/molecular mechanical (QM/MM) methods the need to treat a quantum mechanical fragment in the presence of the rest of the system, which is especially challenging for those parts of the rest of the system that are close to the boundary of the quantum mechanical fragment. This is a delicate matter even for fragments that are not covalently bonded to the rest of the system, but it becomes even more difficult when the boundary of the quantum mechanical fragment cuts a bond. We have developed a suite of methods for more realistically treating interactions across such boundaries. These methods include redistributing and balancing the external partial atomic charges and the use of tuned fluorine atoms for capping dangling bonds, and we have shown that they can greatly improve the accuracy. Finally we present a new approach that goes beyond QM/MM by combining the convenience of molecular mechanics with the accuracy of fitting a potential function to electronic structure calculations on a specific system. To make the latter practical for systems with a large number of degrees of freedom, we developed a method to interpolate between local internal-coordinate fits to the potential energy. A key issue for the application to large systems is that rather than assigning the atoms or monomers to fragments, we assign the internal coordinates to reaction, secondary, and tertiary sets. Thus, we make a partition in coordinate space rather than atom space. Fits to the local dependence of the potential energy on tertiary coordinates are arrayed along a preselected reaction coordinate at a sequence of geometries called anchor points; the potential energy function is called an anchor points reactive potential. Electrostatically embedded fragment methods and the anchor points reactive potential, because they are based on treating an entire system by quantum mechanical electronic structure methods but are affordable for large and complex systems, have the potential to open new areas for accurate simulations where combined QM/MM methods are inadequate.

Experimental and theoretical investigations of energy transfer and hydrogen-bond breaking in small water and HCl clusters

Water is one of the most pervasive molecules on earth and other planetary bodies; it is the molecule that is searched for as the presumptive precursor to extraterrestrial life. It is also the paradigm substance illustrating ubiquitous hydrogen bonding (H-bonding) in the gas phase, liquids, crystals, and amorphous solids. Moreover, H-bonding with other molecules and between different molecules is of the utmost importance in chemistry and biology. It is no wonder, then, that for nearly a century theoreticians and experimentalists have tried to understand all aspects of H-bonding and its influence on reactivity. It is somewhat surprising, therefore, that several fundamental aspects of H-bonding that are particularly important for benchmarking theoretical models have remained unexplored experimentally. For example, even the binding strength between two gas-phase water molecules has never been determined with sufficient accuracy for comparison with high-level electronic structure calculations. Likewise, the effect of cooperativity (nonadditivity) in small H-bonded networks is not known with sufficient accuracy. An even greater challenge for both theory and experiment is the description of the dissociation dynamics of H-bonded small clusters upon acquiring vibrational excitation. This is because of the long lifetimes of many clusters, which requires running classical trajectories for many nanoseconds to achieve dissociation. In this Account, we describe recent progress and ongoing research that demonstrates how the combined and complementary efforts of theory and experiment are enlisted to determine bond dissociation energies (D0) of small dimers and cyclic trimers of water and HCl with unprecedented accuracy, describe dissociation dynamics, and assess the effects of cooperativity. The experimental techniques rely on IR excitation of H-bonded X-H stretch vibrations, measuring velocity distributions of fragments in specific rovibrational states, and determining product state distributions at the pair-correlation level. The theoretical methods are based on high-level ab initio potential energy surfaces used in quantum and classical dynamical calculations. We achieve excellent agreement on D0 between theory and experiments for all of the clusters that we have compared, as well as for cooperativity in ring trimers of water and HCl. We also show that both the long-range and the repulsive parts of the potential must be involved in bond breaking. We explain why H-bonds are so resilient and hard to break, and we propose that a common motif in the breaking of cyclic trimers is the opening of the ring following transfer of one quantum of stretch excitation to form open-chain structures that are weakly bound. However, it still takes many vibrational periods to release one monomer fragment from the open-chain structures. Our success with water and HCl dimers and trimers led us to embark on a more ambitious project: studies of mixed water and HCl small clusters. These clusters eventually lead to ionization of HCl and serve as prototypes of acid dissociation in water. Measurements and calculations of such ionizations are yet to be achieved, and we are now characterizing these systems by adding monomers one at a time. We describe our completed work on the HCl-H2O dimer and mention our recent theoretical results on larger mixed clusters.

Anharmonic rovibrational calculations of singlet cyclic C4 using a new ab initio potential and a quartic force field

We report a CCSD(T)/cc-pCV5Z quartic force field (QFF) and a semi-global CCSD(T)-F12b/aug-cc-pVTZ potential energy surface (PES) for singlet, cyclic C4. Vibrational fundamentals, combinations, and overtones are obtained using vibrational second-order perturbation theory (VPT2) and the vibrational configuration-interaction (VCI) approach. Agreement is within 10 cm(-1) between the VCI calculated fundamentals on the QFF and PES using the MULTIMODE (MM) program, and VPT2 and VCI results agree for the fundamentals. The agreement between VPT2-QFF and MM-QFF results is also good for the C4 combinations and overtones. The J = 1 and J = 2 rovibrational energies are reported from both VCI (MM) on the PES and VPT2 on the QFF calculations. The spectroscopic constants of (12)C4 and two C2v-symmetry, single (13)C-substituted isotopologues are presented, which may help identification of cyclic C4 in future experimental analyses or astronomical observations.

Effects of Zero-Point Delocalization on the Vibrational Frequencies of Mixed HCl and Water Clusters

We demonstrate the significant effect that large-amplitude zero-point vibrational motion can have on the high-frequency fundamental vibrations of molecular clusters, specifically small (HCl)n-(H2O)m clusters. Calculations were conducted on a many-body potential, constructed from a mix of new and previously reported semiempirical and high-level ab initio potentials. Diffusion Monte Carlo simulations were performed to determine ground-state wave functions. Visualization of these wave functions indicates that the clusters exhibit delocalized ground states spanning multiple stationary point geometries. The ground states are best characterized by planar ring configurations, despite the clusters taking nonplanar configurations at their global minima. Vibrational calculations were performed at the global minima and the Diffusion Monte Carlo predicted configurations and also using an approach that spans multiple stationary points along a rectilinear normal-mode reaction path. Significantly better agreement was observed between the calculated vibrational frequencies and experimental peak positions when the delocalized ground state was accounted for.

High-Level, First-Principles, Full-Dimensional Quantum Calculation of the Ro-vibrational Spectrum of the Simplest Criegee Intermediate (CH2OO)

The ro-vibrational spectrum of the simplest Criegee intermediate (CH2OO) has been determined quantum mechanically based on nine-dimensional potential energy and dipole surfaces for its ground electronic state. The potential energy surface is fitted to more than 50 000 high-level ab initio points with a root-mean-square error of 25 cm(-1), using a recently proposed permutation invariant polynomial neural network method. The calculated rotational constants, vibrational frequencies, and spectral intensities of CH2OO are in excellent agreement with experiment. The potential energy surface provides a valuable platform for studying highly excited vibrational and unimolecular reaction dynamics of this important molecule.

Global analytical ab initio ground-state potential energy surface for the C((1)D)+H2 reactive system

A new global ab initio potential energy surface (called ZMB-a) for the 1(1)A' state of the C((1)D)+H2 reactive system has been constructed. This is based upon ab initio calculations using the internally contracted multireference configuration interaction approach with the aug-cc-pVQZ basis set, performed at about 6300 symmetry unique geometries. Accurate analytical fits are generated using many-body expansions with the permutationally invariant polynomials, except that the fit of the deep well region is taken from our previous fit. The ZMB-a surface is unique in the accurate description of the regions around conical intersections (CIs) and of van der Waals (vdW) interactions. The CIs between the 1(1)A' and 2(1)A' states cause two kinds of barriers on the ZMB-a surface: one is in the linear H-CH dissociation direction with a barrier height of 9.07 kcal/mol, which is much higher than those on the surfaces reported before; the other is in the C((1)D) collinearly attacking H2 direction with a barrier height of 12.39 kcal/mol. The ZMB-a surface basically reproduces our ab initio calculations in the vdW interaction regions, and supports a linear C-HH vdW complex in the entrance channel, and two vdW complexes in the exit channel, at linear CH-H and HC-H geometries, respectively.

Representing potential energy surfaces by high-dimensional neural network potentials

The development of interatomic potentials employing artificial neural networks has seen tremendous progress in recent years. While until recently the applicability of neural network potentials (NNPs) has been restricted to low-dimensional systems, this limitation has now been overcome and high-dimensional NNPs can be used in large-scale molecular dynamics simulations of thousands of atoms. NNPs are constructed by adjusting a set of parameters using data from electronic structure calculations, and in many cases energies and forces can be obtained with very high accuracy. Therefore, NNP-based simulation results are often very close to those gained by a direct application of first-principles methods. In this review, the basic methodology of high-dimensional NNPs will be presented with a special focus on the scope and the remaining limitations of this approach. The development of NNPs requires substantial computational effort as typically thousands of reference calculations are required. Still, if the problem to be studied involves very large systems or long simulation times this overhead is regained quickly. Further, the method is still limited to systems containing about three or four chemical elements due to the rapidly increasing complexity of the configuration space, although many atoms of each species can be present. Due to the ability of NNPs to describe even extremely complex atomic configurations with excellent accuracy irrespective of the nature of the atomic interactions, they represent a general and therefore widely applicable technique, e.g. for addressing problems in materials science, for investigating properties of interfaces, and for studying solvation processes.

Multisurface Adiabatic Reactive Molecular Dynamics

Adiabatic reactive molecular dynamics (ARMD) simulation method is a surface-crossing algorithm for modeling chemical reactions in classical molecular dynamics simulations using empirical force fields. As the ARMD Hamiltonian is time dependent during crossing, it allows only approximate energy conservation. In the current work, the range of applicability of conventional ARMD is explored, and a new multisurface ARMD (MS-ARMD) method is presented, implemented in CHARMM and applied to the vibrationally induced photodissociation of sulfuric acid (H2SO4) in the gas phase. For this, an accurate global potential energy surface (PES) involving 12 H2SO4 and 4 H2O + SO3 force fields fitted to MP2/6-311G++(2d,2p) reference energies is employed. The MS-ARMD simulations conserve total energy and feature both intramolecular H-transfer reactions and water elimination. An analytical treatment of the dynamics in the crossing region finds that conventional ARMD can approximately conserve total energy for limiting cases. In one of them, the reduced mass of the system is large, which often occurs for simulations of solvated biomolecular systems. On the other hand, MS-ARMD is a general approach for modeling chemical reactions including gas-phase, homogeneous, heterogeneous, and enzymatic catalytic reactions while conserving total energy in atomistic simulations.

Quasiclassical Trajectory Calculations of the Rate Constant of the OH + HBr → Br + H2O Reaction Using a Full-Dimensional Ab Initio Potential Energy Surface Over the Temperature Range 5 to 500 K

We report a permutationally invariant, ab initio potential energy surface (PES) for the OH + HBr → Br + H2O reaction. The PES is a fit to roughly 26 000 spin-free UCCSD(T)/cc-pVDZ-F12a energies and has no classical barrier to reaction. It is used in quasiclassical trajectory calculations with a focus on the thermal rate constant, k(T), over the temperature range 5 to 500 K. Comparisons with available experimental data over the temperature range 23 to 416 K are made using three approaches to treat the OH rotational and associated electronic partition function. All display an inverse temperature dependence of k(T) below roughly 160 K and a nearly constant temperature dependence above 160 K, in agreement with experiment. The calculated rate constant with no treatment of spin-orbit coupling is overall in the best agreement with experiment, being (probably fortuitously) within 20% of it.

Ab initio potential energy and dipole moment surfaces of the F(-)(H2O) complex

We present full-dimensional, ab initio potential energy and dipole moment surfaces for the F(-)(H2O) complex. The potential surface is a permutationally invariant fit to 16,114 coupled-cluster single double (triple)/aVTZ energies, while the dipole surface is a covariant fit to 11,395 CCSD(T)/aVTZ dipole moments. Vibrational self-consistent field/vibrational configuration interaction (VSCF/VCI) calculations of energies and the IR-spectrum are presented both for F(-)(H2O) and for the deuterated analog, F(-)(D2O). A one-dimensional calculation of the splitting of the ground state, due to equivalent double-well global minima, is also reported.

Full-dimensional quantum calculations of the dissociation energy, zero-point, and 10 K properties of H7+/D7+ clusters using an ab initio potential energy surface

Diffusion Monte Carlo (DMC) and path-integral Monte Carlo computations of the vibrational ground state and 10 K equilibrium state properties of the H7 (+)/D7 (+) cations are presented, using an ab initio full-dimensional potential energy surface. The DMC zero-point energies of dissociated fragments H5 (+)(D5 (+))+H2(D2) are also calculated and from these results and the electronic dissociation energy, dissociation energies, D0, of 752 ± 15 and 980 ± 14 cm(-1) are reported for H7 (+) and D7 (+), respectively. Due to the known error in the electronic dissociation energy of the potential surface, these quantities are underestimated by roughly 65 cm(-1). These values are rigorously determined for first time, and compared with previous theoretical estimates from electronic structure calculations using standard harmonic analysis, and available experimental measurements. Probability density distributions are also computed for the ground vibrational and 10 K state of H7 (+) and D7 (+). These are qualitatively described as a central H3 (+)/D3 (+) core surrounded by "solvent" H2/D2 molecules that nearly freely rotate.

Full-dimensional potentials and state couplings and multidimensional tunneling calculations for the photodissociation of phenol

Full-dimensional potentials and state couplings were developed for the photodissociation of phenol. We also present multidimensional tunneling calculations at the transition state on the shoulder of the first conical intersection.

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.

Six-dimensional quantum dynamics for dissociative chemisorption of H2 and D2 on Ag(111) on a permutation invariant potential energy surface

Quantum dynamics on a permutation invariant potential energy surface for H₂ dissociation on Ag(111) yield satisfactory agreement with experiment.

Mode selective dynamics and kinetics of the H2 + F2 → H + HF + F reaction

The reactivity is significantly enhanced by vibrational excitation of F₂ whereas excitation of H₂ vibration has a moderate effect.

About the variational property of generalized discrete variable representation

It is explained why the generalized discrete variable representation is variational neither with respect to the size (truncation) of the basis set nor with respect to the selection of grid points.

Quantum mechanical characterization of the He4ICl weakly bound complex

Vibrational calculations are performed for the 12-dimensional He4ICl van der Waals complex using the multiconfiguration time-dependent Hartree (MCTDH) method. The potential energy surface of the cluster is represented as a sum of the triatomic He-ICl ab initio parametrized terms plus the He-He interactions. The topology of the potential presents higher anisotropy compared to the one with a homonuclear dopant, and this is clearly reflected in the structure and energetics of the low-lying conformers of the system. In order to take advantage of the MCTDH method, natural potential fits are employed for the potential energy operator, and also, a mode combination scheme is introduced in order to speed up the computations. Zero-point energy, binding energies, and vibrationally averaged structures of different isomers of the He4ICl cluster are obtained. The present results predict that the (3,1,0) structure, involving three He atoms in the near T-shaped and one He atom in the linear configurations, to be the most stable one in accord with recent experimental findings. Comparisons with previous theoretical and experimental data are presented, and the stability of the high-order conformers is discussed in connection with the multiple minima (global and local) of the underlying potential surface.

Ab initio based potential energy surface and kinetics study of the OH + NH3 hydrogen abstraction reaction

A full-dimensional analytical potential energy surface (PES) for the OH + NH3 → H2O + NH2 gas-phase reaction was developed based exclusively on high-level ab initio calculations. This reaction presents a very complicated shape with wells along the reaction path. Using a wide spectrum of properties of the reactive system (equilibrium geometries, vibrational frequencies, and relative energies of the stationary points, topology of the reaction path, and points on the reaction swath) as reference, the resulting analytical PES reproduces reasonably well the input ab initio information obtained at the coupled-cluster single double triple (CCSD(T)) = FULL/aug-cc-pVTZ//CCSD(T) = FC/cc-pVTZ single point level, which represents a severe test of the new surface. As a first application, on this analytical PES we perform an extensive kinetics study using variational transition-state theory with semiclassical transmission coefficients over a wide temperature range, 200-2000 K. The forward rate constants reproduce the experimental measurements, while the reverse ones are slightly underestimated. However, the detailed analysis of the experimental equilibrium constants (from which the reverse rate constants are obtained) permits us to conclude that the experimental reverse rate constants must be re-evaluated. Another severe test of the new surface is the analysis of the kinetic isotope effects (KIEs), which were not included in the fitting procedure. The KIEs reproduce the values obtained from ab initio calculations in the common temperature range, although unfortunately no experimental information is available for comparison.