Permutationally invariant potential energy surfaces in high dimensionality

Global potential energy surface of ground state singlet spin O4

A new global potential energy for the singlet spin state O₄ system is reported using CASPT2/aug-cc-pVTZ ab initio calculations. The geometries for the six-dimensional surface are constructed using a novel point generation scheme that employs randomly generated configurations based on the beta distribution. The advantage of this scheme is apparent in the reduction of the number of required geometries for a reasonably accurate potential energy surface (PES) and the consequent decrease in the overall computational effort. The reported surface matches well with the recently published singlet surface by Paukku et al. [J. Chem. Phys. 147, 034301 (2017)]. In addition to the O₄ PES, the ground state N₄ PES is also constructed using the point generation scheme and compared with the existing PES [Y. Paukku et al., J. Chem. Phys. 139, 044309 (2013)]. The singlet surface is constructed with the aim of studying high energy O₂-O₂ collisions and predicting collision induced dissociation cross section to be used in simulating non-equilibrium aerothermodynamic flows.

Full dimensional potential energy surface and low temperature dynamics of the H2CO + OH → HCO + H2O reaction

A new method is proposed to analytically represent the potential energy surface of reactions involving polyatomic molecules capable of accurately describing long-range interactions and saddle points, needed to describe low-temperature collisions. It is based on two terms, a reactive force field term and a many-body term. The reactive force field term accurately describes the fragments, long-range interactions among them and the saddle points for reactions. The many-body term increases the desired accuracy everywhere else. This method has been applied to the OH + H2CO → H2O + HCO reaction, giving a barrier of 27.4 meV. The simulated classical rate constants with this potential are in good agreement with recent experimental results [Ocaña et al., Astrophys. J., 2017, submitted], showing an important increase at temperatures below 100 K. The reaction mechanism is analyzed in detail here, and explains the observed behavior at low energy by the formation of long-lived collision complexes, with roaming trajectories, with a capture observed for very long impact parameters, >100 a.u., determined by the long-range dipole-dipole interaction.

Photodissociation of CH3CHO at 248 nm: identification of the channels of roaming, triple fragmentation and the transition state

Quasi-classical trajectory (QCT) calculations are performed on the molecular products CO + CH₄via the tight transition state (TS) and global minimum configurations. With the aid of this theoretical evidence, we have re-examined the experimental results published previously to clarify the controversial issue of photodissociation dynamics of CH₃CHO at 248 nm. For the CO (v = 0 and 1) bimodal rotational distributions obtained previously [K.-C. Hung, P.-Y. Tsai, H.-K. Li, and K.-C. Lin, J. Chem. Phys., 2014, 140, 064313], the low-rotational (J) component is re-assigned to the contribution of triple fragmentation (H + CO + CH₃), whereas the high-J component is ascribed to the CH₃-roaming pathway. The H-roaming pathway is not found in the calculations. Further, the QCT results have confirmed that the CO vibrational population especially at higher states and the low-energy component of CH₄ vibrational bimodality obtained experimentally are mainly produced following the TS pathway, which has never been identified before. While taking into account both the theoretical and experimental results, the ratio of the molecular products (CO(v = 1) + CH₄) obtained by the triple fragmentation/roaming/TS processes is evaluated to be 0.23 : 1 : 0.29.

Quantum and classical IR spectra of (HCOOH)2, (DCOOH)2 and (DCOOD)2 using ab initio potential energy and dipole moment surfaces

Full-dimensional (24 modes) quantum calculation of the IR spectrum of (DCOOD)₂, and comparison with classical MD one.

An improved two-rotor function for conformational potential energy surfaces of 20 amino acid diamides

Predicting the three-dimensional structure of a protein from its amino acid sequence requires a complete understanding of the molecular forces that influences the protein folding process. Each possible conformation has its corresponding potential energy, which characterizes its thermodynamic stability. This is needed to identify the primary intra- and inter-molecular interactions, so that we can reduce the dimensionality of the problem, and create a relatively simple representation of the system. Investigating this problem using quantum chemical methods produces accurate results; however, this also entails large computational resources. In this study, an improved two-rotor potential energy function is proposed to represent the backbone interactions in amino acids through a linear combination of a Fourier series and a mixture of Gaussian functions. This function is applied to approximate the 20 amino acid diamide Ramachandran-type PESs, and results yielded an average RMSE of 2.36 kJ mol−1, which suggest that the mathematical model precisely captures the general topology of the conformational potential energy surface. Furthermore, this paper provides insights on the conformational preferences of amino acid diamides through local minima geometries and energy ranges, using the improved mathematical model. The proposed mathematical model presents a simpler representation that attempts to provide a framework on building polypeptide models from individual amino acid functions, and consequently, a novel method for rapid but accurate evaluation of potential energies for biomolecular simulations.

An accurate full-dimensional potential energy surface and quasiclassical trajectory dynamics of the H + H2O2 two-channel reaction

Quasiclassical trajectory calculations reveal interesting dynamics features based on an accurate FI-NN PES for the H + H₂O₂ two-channel reaction.

Dynamics and kinetics of the reaction OH + H2S → H2O + SH on an accurate potential energy surface

Mode specificity and product energy disposal are unveiled in the reaction OH + H₂S → H₂O + SH.

Ab initio instanton rate theory made efficient using Gaussian process regression

Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into the calculations of chemical reaction rates based on a potential-energy surface computed on the fly from electronic-structure theory. This approach is necessary to extend conventional transition-state theory into the deep-tunnelling regime, but it is also more computationally expensive as it requires many more ab initio calculations. We propose an approach which uses Gaussian process regression to fit the potential-energy surface locally around the dominant tunnelling pathway. The method can be converged to give the same result as from an on-the-fly ab initio instanton calculation but it requires far fewer electronic-structure calculations. This makes it a practical approach for obtaining accurate rate constants based on high-level electronic-structure methods. We show fast convergence to reproduce benchmark H + CH4 results and evaluate new low-temperature rates of H + C2H6 in full dimensionality at a UCCSD(T)-F12b/cc-pVTZ-F12 level.

Diabatic potential energy surfaces of MgH2+ and dynamic studies for the Mg+(3p) + H2 → MgH+ + H reaction

The global diabatic potential energy surfaces for the Mg⁺(3p) + H₂ → MgH⁺ + H reaction are structured for the first time.

Development of a Flexible-Monomer Two-Body Carbon Dioxide Potential and Its Application to Clusters up to (CO2 )13

A flexible-monomer two-body potential energy function was developed that approaches the high level CCSD(T)/CBS potential energy surface (PES) of carbon dioxide (CO₂ ) systems. This function was generated by fitting the electronic energies of unique CO₂ monomers and dimers to permutationally invariant polynomials. More than 200,000 CO₂ configurations were used to train the potential function. Comparisons of the PESs of six orientations of flexible CO₂ dimers were evaluated to demonstrate the accuracy of the potential. Furthermore, the potential function was used to determine the minimum energy structures of CO₂ clusters containing as many as 13 molecules. For isomers of (CO₂ )₃ , the potential demonstrated energetic agreement with the M06-2X functional and structural agreement of the B2PLYP-D functional at substantially reduced computational costs. A separate function, fit to MP2/aug-cc-pVDZ reference energies, was developed to directly compare the two-body potential to the ab initio MP2 level of theory. © 2017 Wiley Periodicals, Inc.

Dynamics and Novel Mechanisms of SN2 Reactions on ab Initio Analytical Potential Energy Surfaces

We describe a novel theoretical approach to the bimolecular nucleophilic substitution (S_N2) reactions that is based on analytical potential energy surfaces (PESs) obtained by fitting a few tens of thousands high-level ab initio energy points. These PESs allow computing millions of quasi-classical trajectories thereby providing unprecedented statistical accuracy for S_N2 reactions, as well as performing high-dimensional quantum dynamics computations. We developed full-dimensional ab initio PESs for the F^- + CH₃Y [Y = F, Cl, I] systems, which describe the direct and indirect, complex-forming Walden-inversion, the frontside attack, and the new double-inversion pathways as well as the proton-transfer channels. Reaction dynamics simulations on the new PESs revealed (a) a novel double-inversion S_N2 mechanism, (b) frontside complex formation,

Accurate potential energy surfaces for hydrogen abstraction reactions: A benchmark study on the XYG3 doubly hybrid density functional

The potential energy surface (PES) for the H + CH₄ system has been constructed with the recently developed XYG3 doubly hybrid functional, while those with the standard B3LYP hybrid functional, and the Møller-Plesset perturbation theory up to the second order (MP2) are also presented for comparison. Quantum dynamics studies demonstrated that satisfactory results on the reaction probabilities and the rate coefficients can be obtained on top of the XYG3-PES, as compared to the results based on the highly accurate, yet expensive, CCSD(T)-PES (Li et al., J. Chem. Phys. 2015, 142, 204302). Further investigation suggested that the XYG3 functional is useful in providing accurate rate coefficients for some larger systems involving H atom abstractions. © 2017 Wiley Periodicals, Inc.

First Principles Neural Network Potentials for Reactive Simulations of Large Molecular and Condensed Systems

Modern simulation techniques have reached a level of maturity which allows a wide range of problems in chemistry and materials science to be addressed. Unfortunately, the application of first principles methods with predictive power is still limited to rather small systems, and despite the rapid evolution of computer hardware no fundamental change in this situation can be expected. Consequently, the development of more efficient but equally reliable atomistic potentials to reach an atomic level understanding of complex systems has received considerable attention in recent years. A promising new development has been the introduction of machine learning (ML) methods to describe the atomic interactions. Once trained with electronic structure data, ML potentials can accelerate computer simulations by several orders of magnitude, while preserving quantum mechanical accuracy. This Review considers the methodology of an important class of ML potentials that employs artificial neural networks.

Fitting potential energy surfaces to sum-of-products form with neural networks using exponential neurons

In this paper, the use of the neural network (NN) method with exponential neurons for directly fitting ab initio data to generate potential energy surfaces (PESs) in sum-of-product form will be discussed. The utility of the approach will be highlighted using fits of CS2, HFCO, and HONO ground state PESs based upon high-level ab initio data. Using a generic interface between the neural network PES fitting, which is performed in MATLAB, and the Heidelberg multi-configuration time-dependent Hartree (MCTDH) software package, the PESs have been tested via comparison of vibrational energies to experimental measurements. The review demonstrates the potential of the PES fitting method, combined with MCTDH, to tackle high-dimensional quantum dynamics problems.

A neural network potential energy surface for the NaH2 system and dynamics studies on the H(2S) + NaH(X1Σ+) → Na(2S) + H2(X1Σg+) reaction

In order to study the dynamics of the reaction H(²S) + NaH(X¹Σ⁺) → Na(²S) + H₂(X¹Σ_g⁺), a new potential energy surface (PES) for the ground state of the NaH₂ system is constructed based on 35 730 ab initio energy points. Using basis sets of quadruple zeta quality, multireference configuration interaction calculations with Davidson correction were carried out to obtain the ab initio energy points. The neural network method is used to fit the PES, and the root mean square error is very small (0.00639 eV). The bond lengths, dissociation energies, zero-point energies and spectroscopic constants of H₂(X¹Σ_g⁺) and NaH(X¹Σ⁺) obtained on the new NaH₂ PES are in good agreement with the experiment data. On the new PES, the reactant coordinate-based time-dependent wave packet method is applied to study the reaction dynamics of H(²S) + NaH(X¹Σ⁺) → Na(²S) + H₂(X¹Σ_g⁺), and the reaction probabilities, integral cross-sections (ICSs) and differential cross-sections (DCSs) are obtained. There is no threshold in the reaction due to the absence of an energy barrier on the minimum energy path. When the collision energy increases, the ICSs decrease from a high value at low collision energy. The DCS results show that the angular distribution of the product molecules tends to the forward direction. Compared with the LiH₂ system, the NaH₂ system has a larger mass and the PES has a larger well at the H-NaH configuration, which leads to a higher ICS value in the H(²S) + NaH(X¹Σ⁺) → Na(²S) + H₂(X¹Σ_g⁺) reaction. Because the H(²S) + NaH(X¹Σ⁺) → Na(²S) + H₂(X¹Σ_g⁺) reaction releases more energy, the product molecules can be excited to a higher vibrational state.

Accuracy of Frequencies Obtained with the Aid of Explicitly Correlated Wave Function Based Methods

We asses the basis set convergence of harmonic frequencies using different explicitly correlated wave function based methods. All commonly available CCSD(T) variants as well as MP2-F12 and MP4(F12*) are considered, and a hierarchy of the different approaches is established. As for reaction and atomization energies, CCSD(F12*)(T*) is a close approximation to CCSD(F12)(T*) and clearly superior to the other tested approximations. The used scaling for the triples correction enhances the accuracy relative to CCSD(F12*)(T) especially for small basis sets and is very attractive since no additional computational costs are added. However, this scaling slightly breaks size consistency, and therefore we additionally study the accuracy of CCSD(F12*)(T*) and CCSD(F12*)(T) in the context of calculating anharmonic frequencies to check if this causes problems in the generation of the potential energy surface (PES). We find a fast basis set convergence for harmonic and anharmonic frequencies. Already in the cc-pVDZ-F12 basis, the RMSD to the CBS limit is only around 4-5 cm^-1.

Toolkit for the Construction of Reproducing Kernel-Based Representations of Data: Application to Multidimensional Potential Energy Surfaces

In the early days of computation, slow processor speeds limited the amount of data that could be generated and used for scientific purposes. In the age of big data, the limiting factor usually is the method with which large amounts of data are analyzed and useful information is extracted. A typical example from chemistry are high-level ab initio calculations for small systems, which have nowadays become feasible even if energies at many different geometries are required. Molecular dynamics simulations often require several thousand distinct trajectories to be run. Under such circumstances suitable analytical representations of potential energy surfaces (PESs) based on ab initio calculations are required to propagate the dynamics at an acceptable cost. In this work we introduce a toolkit which allows the automatic construction of multidimensional PESs from gridded ab initio data based on reproducing kernel Hilbert space (RKHS) theory. The resulting representations require no tuning of parameters and allow energy and force evaluations at ab initio quality at the same cost as empirical force fields. Although the toolkit is primarily intended for constructing multidimensional potential energy surfaces for molecular systems, it can also be used for general machine learning purposes. The software is published under the MIT license and can be downloaded, modified, and used in other projects for free.

Semiclassical "Divide-and-Conquer" Method for Spectroscopic Calculations of High Dimensional Molecular Systems

A new semiclassical "divide-and-conquer" method is presented with the aim of demonstrating that quantum dynamics simulations of high dimensional molecular systems are doable. The method is first tested by calculating the quantum vibrational power spectra of water, methane, and benzene-three molecules of increasing dimensionality for which benchmark quantum results are available-and then applied to C_{60}, a system characterized by 174 vibrational degrees of freedom. Results show that the approach can accurately account for quantum anharmonicities, purely quantum features like overtones, and the removal of degeneracy when the molecular symmetry is broken.

On-the-Fly ab Initio Semiclassical Calculation of Glycine Vibrational Spectrum

We present an on-the-fly ab initio semiclassical study of vibrational energy levels of glycine, calculated by Fourier transform of the wavepacket correlation function. It is based on a multiple coherent states approach integrated with monodromy matrix regularization for chaotic dynamics. All four lowest-energy glycine conformers are investigated by means of single-trajectory semiclassical spectra obtained upon classical evolution of on-the-fly trajectories with harmonic zero-point energy. For the most stable conformer I, direct dynamics trajectories are also run for each vibrational mode with energy equal to the first harmonic excitation. An analysis of trajectories evolved up to 50 000 atomic time units demonstrates that, in this time span, conformers II and III can be considered as isolated species, while conformers I and IV show a pretty facile interconversion. Therefore, previous perturbative studies based on the assumption of isolated conformers are often reliable but might be not completely appropriate in the case of conformer IV and conformer I for which interconversion occurs promptly.

A new (multi-reference configuration interaction) potential energy surface for H2CO and preliminary studies of roaming

We report a new global potential energy surface (PES) for H2CO, based on precise fitting of roughly 67 000 MRCI/cc-pVTZ energies. This PES describes the global minimum, the cis- and trans-HCOH isomers, and barriers relevant to isomerization, formation of the molecular (H2+CO) and radical (H+HCO) products, and the loose so-called roaming transition-state saddle point. The key features of the PES are reviewed and compared with a previous PES, denoted by PES04, based on five local fits that are 'stitched' together by switching functions (Zhang et al. 2004 J. Phys. Chem. A108, 8980-8986 (doi:10.1021/jp048339l)). Preliminary quasi-classical trajectory calculations are performed at the total energy of 36 233 cm-1 (103 kcal mol-1), relative to the H2CO global minimum, using the new PES, with a particular focus on roaming dynamics. When compared with the results from PES04, the new PES findings show similar rotational distributions, somewhat more roaming and substantially higher H2 vibrational excitation.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

An improved spin-orbit coupling model for use within the effective relativistic coupling by asymptotic representation (ERCAR) method

An improved atomic spin-orbit model is presented, which is designed to be used within the framework of the effective relativistic coupling by asymptotic representation method. This method is used for the generation of highly accurate coupled potential energy surfaces (PESs) to represent the fine structure energies of appropriate systems. The approach is demonstrated using CH₃I and its photodissociation as a typical example. The method is based on a specific diabatization of electronic spin-space ("spin-free") states with respect to the asymptote at which a single relativistic atom is separated from a molecular non-relativistic fragment. Thus, the relativistic coupling effects can be treated entirely within the atomic framework. So far, an effective spin-orbit coupling operator which only accounts for intra-state coupling within each atomic spin-space state was used. In the present work, this approach is extended to account for inter-state couplings among different atomic spin-space states as well. It is shown that this extended approach improves the accuracy of the PESs significantly for higher excited states and also enhances the accuracy of low energy states. In particular, it improves the representation of the spin-orbit induced conical intersection among the ³Q₀ and ¹Q₁ states of CH₃I, which is of high relevance for the nonadiabatic quantum dynamics of the photodissociation.

Dynamical barrier and isotope effects in the simplest substitution reaction via Walden inversion mechanism

Reactions occurring at a carbon atom through the Walden inversion mechanism are one of the most important and useful classes of reactions in chemistry. Here we report an accurate theoretical study of the simplest reaction of that type: the H+CH4 substitution reaction and its isotope analogues. It is found that the reaction threshold versus collision energy is considerably higher than the barrier height. The reaction exhibits a strong normal secondary isotope effect on the cross-sections measured above the reaction threshold, and a small but reverse secondary kinetic isotope effect at room temperature. Detailed analysis reveals that the reaction proceeds along a path with a higher barrier height instead of the minimum-energy path because the umbrella angle of the non-reacting methyl group cannot change synchronously with the other reaction coordinates during the reaction due to insufficient energy transfer from the translational motion to the umbrella mode.

High-level ab initio potential energy surface and dynamics of the F- + CH3I SN2 and proton-transfer reactions

Bimolecular nucleophilic substitution (S_N2) and proton transfer are fundamental processes in chemistry and F^- + CH₃I is an important prototype of these reactions. Here we develop the first full-dimensional ab initio analytical potential energy surface (PES) for the F^- + CH₃I system using a permutationally invariant fit of high-level composite energies obtained with the combination of the explicitly-correlated CCSD(T)-F12b method, the aug-cc-pVTZ basis, core electron correlation effects, and a relativistic effective core potential for iodine. The PES accurately describes the S_N2 channel producing I^- + CH₃F via Walden-inversion, front-side attack, and double-inversion pathways as well as the proton-transfer channel leading to HF + CH₂I^-. The relative energies of the stationary points on the PES agree well with the new explicitly-correlated all-electron CCSD(T)-F12b/QZ-quality benchmark values. Quasiclassical trajectory computations on the PES show that the proton transfer becomes significant at high collision energies and double-inversion as well as front-side attack trajectories can occur. The computed broad angular distributions and hot internal energy distributions indicate the dominance of indirect mechanisms at lower collision energies, which is confirmed by analyzing the integration time and leaving group velocity distributions. Comparison with available crossed-beam experiments shows usually good agreement.

Perspective: Machine learning potentials for atomistic simulations

Nowadays, computer simulations have become a standard tool in essentially all fields of chemistry, condensed matter physics, and materials science. In order to keep up with state-of-the-art experiments and the ever growing complexity of the investigated problems, there is a constantly increasing need for simulations of more realistic, i.e., larger, model systems with improved accuracy. In many cases, the availability of sufficiently efficient interatomic potentials providing reliable energies and forces has become a serious bottleneck for performing these simulations. To address this problem, currently a paradigm change is taking place in the development of interatomic potentials. Since the early days of computer simulations simplified potentials have been derived using physical approximations whenever the direct application of electronic structure methods has been too demanding. Recent advances in machine learning (ML) now offer an alternative approach for the representation of potential-energy surfaces by fitting large data sets from electronic structure calculations. In this perspective, the central ideas underlying these ML potentials, solved problems and remaining challenges are reviewed along with a discussion of their current applicability and limitations.