Permutationally invariant potential energy surfaces in high dimensionality

Toward spin-orbit coupled diabatic potential energy surfaces for methyl iodide using effective relativistic coupling by asymptotic representation

The theoretical treatment of state-state interactions and the development of coupled multidimensional potential energy surfaces (PESs) is of fundamental importance for the theoretical investigation of nonadiabatic processes. Usually, only derivative or vibronic coupling is considered, but the presence of heavy atoms in a system can render spin-orbit (SO) coupling important as well. In the present study, we apply a new method recently developed by us (J. Chem. Phys. 2012, 136, 034103, and J. Chem. Phys. 2012, 137, 064101) to generate SO coupled diabatic PESs along the C-I dissociation coordinate for methyl iodide (CH3I). This is the first and mandatory step toward the development of fully coupled full-dimensional PESs to describe the multistate photodynamics of this benchmark system. The method we use here is based on the diabatic asymptotic representation of the molecular fine structure states and an effective relativistic coupling operator. It therefore is called effective relativistic coupling by asymptotic representation (ERCAR). This approach allows the efficient and accurate generation of fully coupled PESs including derivative and SO coupling based on high-level ab initio calculations. In this study we develop a specific ERCAR model for CH3I that so far accounts only for the C-I bond cleavage. Details of the diabatization and the accuracy of the results are investigated in comparison to reference ab initio calculations and experiments. The energies of the adiabatic fine structure states are reproduced in excellent agreement with ab initio SO-CI data. The model is also compared to available literature data, and its performance is evaluated critically. This shows that the new method is very promising for the construction of fully coupled full-dimensional PESs for CH3I to be used in future quantum dynamics studies.

A density-functional theory-based neural network potential for water clusters including van der Waals corrections

The fundamental importance of water for many chemical processes has motivated the development of countless efficient but approximate water potentials for large-scale molecular dynamics simulations, from simple empirical force fields to very sophisticated flexible water models. Accurate and generally applicable water potentials should fulfill a number of requirements. They should have a quality close to quantum chemical methods, they should explicitly depend on all degrees of freedom including all relevant many-body interactions, and they should be able to describe molecular dissociation and recombination. In this work, we present a high-dimensional neural network (NN) potential for water clusters based on density-functional theory (DFT) calculations, which is constructed using clusters containing up to 10 monomers and is in principle able to meet all these requirements. We investigate the reliability of specific parametrizations employing two frequently used generalized gradient approximation (GGA) exchange-correlation functionals, PBE and RPBE, as reference methods. We find that the binding energy errors of the NN potentials with respect to DFT are significantly lower than the typical uncertainties of DFT calculations arising from the choice of the exchange-correlation functional. Further, we examine the role of van der Waals interactions, which are not properly described by GGA functionals. Specifically, we incorporate the D3 scheme suggested by Grimme (J. Chem. Phys. 2010, 132, 154104) in our potentials and demonstrate that it can be applied to GGA-based NN potentials in the same way as to DFT calculations without modification. Our results show that the description of small water clusters provided by the RPBE functional is significantly improved if van der Waals interactions are included, while in case of the PBE functional, which is well-known to yield stronger binding than RPBE, van der Waals corrections lead to overestimated binding energies.

Quasiclassical trajectory study of CH3NO2 decomposition via roaming mediated isomerization using a global potential energy surface

We report a global potential energy surface (PES) for CH3NO2 based on fitting roughly 114,000 density functional theory (UB3LYP/6-311+g(d,p)) electronic energies. The PES is a precise, permutationally invariant fit to these energies. Properties of the PES are described, as are some preliminary quasiclassical trajectory calculations. An isomerization-roaming pathway to the CH3ONO isomer and then to the CH3O + NO products is found. Although the pathway occurs at larger distances than a related loose saddle-point on the PES, the pathway supports the supposition of such a pathway based on locating a loose first-order saddle point and associated IRC, reported previously by Zhu and Lin [Zhu, R. S. and Lin, M. C. Chem. Phys. Lett. 2009, 478, 11].

Improved Potential Energy Surface of Ozone Constructed Using the Fitting by Permutationally Invariant Polynomial Function

New global potential energy surface for the ground electronic state of ozone is constructed at the complete basis set level of the multireference configuration interaction theory. A method of fitting the data points by analytical permutationally invariant polynomial function is adopted. A small set of 500 points is preoptimized using the old surface of ozone. In this procedure the positions of points in the configuration space are chosen such that the RMS deviation of the fit is minimized. New ab initio calculations are carried out at these points and are used to build new surface. Additional points are added to the vicinity of the minimum energy path in order to improve accuracy of the fit, particularly in the region where the surface of ozone exhibits a shallow van der Waals well. New surface can be used to study formation of ozone at thermal energies and its spectroscopy near the dissociation threshold.

Full-dimensional (15-dimensional) ab initio analytical potential energy surface for the H7+ cluster

Full-dimensional ab initio potential energy surface is constructed for the H(7)(+) cluster. The surface is a fit to roughly 160,000 interaction energies obtained with second-order MöllerPlesset perturbation theory and the cc-pVQZ basis set, using the invariant polynomial method [B. J. Braams and J. M. Bowman, Int. Rev. Phys. Chem. 28, 577 (2009)]. We employ permutationally invariant basis functions in Morse-type variables for all the internuclear distances to incorporate permutational symmetry with respect to interchange of H atoms into the representation of the surface. We describe how different configurations are selected in order to create the database of the interaction energies for the linear least squares fitting procedure. The root-mean-square error of the fit is 170 cm(-1) for the entire data set. The surface dissociates correctly to the H(5)(+) + H(2) fragments. A detailed analysis of its topology, as well as comparison with additional ab initio calculations, including harmonic frequencies, verify the quality and accuracy of the parameterized potential. This is the first attempt to present an analytical representation of the 15-dimensional surface of the H(7)(+) cluster for carrying out dynamics studies.

Constructing Potential Energy Surfaces for Polyatomic Systems: Recent Progress and New Problems

Different methods of constructing potential energy surfaces in polyatomic systems are reviewed, with the emphasis put on fitting, interpolation, and analytical (defined by functional forms) approaches, based on quantum chemistry electronic structure calculations. The different approaches are reviewed first, followed by a comparison using the benchmark H + CH4 and the H + NH3 gas-phase hydrogen abstraction reactions. Different kinetics and dynamics properties are analyzed for these reactions and compared with the available experimental data, which permits one to estimate the advantages and disadvantages of each method. Finally, we analyze different problems with increasing difficulty in the potential energy construction: spin-orbit coupling, molecular size, and more complicated reactions with several maxima and minima, which test the soundness and general applicability of each method. We conclude that, although the field of small systems, typically atom-diatom, is mature, there still remains much work to be done in the field of polyatomic systems.

Intersystem crossing and dynamics in O(3P) + C2H4 multichannel reaction: experiment validates theory

The O((3)P) + C(2)H(4) reaction, of importance in combustion and atmospheric chemistry, stands out as a paradigm reaction involving triplet- and singlet-state potential energy surfaces (PESs) interconnected by intersystem crossing (ISC). This reaction poses challenges for theory and experiments owing to the ruggedness and high dimensionality of these potentials, as well as the long lifetimes of the collision complexes. Primary products from five competing channels (H + CH(2)CHO, H + CH(3)CO, H(2) + CH(2)CO, CH(3) + HCO, CH(2) + CH(2)O) and branching ratios (BRs) are determined in crossed molecular beam experiments with soft electron-ionization mass-spectrometric detection at a collision energy of 8.4 kcal/mol. As some of the observed products can only be formed via ISC from triplet to singlet PESs, from the product BRs the extent of ISC is inferred. A new full-dimensional PES for the triplet state as well as spin-orbit coupling to the singlet PES are reported, and roughly half a million surface hopping trajectories are run on the coupled singlet-triplet PESs to compare with the experimental BRs and differential cross-sections. Both theory and experiment find almost equal contributions from the two PESs to the reaction, posing the question of how important is it to consider the ISC as one of the nonadiabatic effects for this and similar systems involved in combustion chemistry. Detailed comparisons at the level of angular and translational energy distributions between theory and experiment are presented for the two primary channel products, CH(3) + HCO and H + CH(2)CHO. The agreement between experimental and theoretical functions is excellent, implying that theory has reached the capability of describing complex multichannel nonadiabatic reactions.

Potential energy and dipole moment surfaces of HCO- for the search of H- in the interstellar medium

Potential energy and permanent dipole moment surfaces of the electronic ground state of formyl negative ion HCO(-) are determined for a large number of geometries using the coupled-cluster theory with single and double and perturbative treatment of triple excitations ab initio method with a large basis set. The obtained data are used to construct interpolated surfaces, which are extended analytically to the region of large separations between CO and H(-) with the multipole expansion approach. We have calculated the energy of the lowest rovibrational levels of HCO(-) that should guide the spectroscopic characterization of HCO(-) in laboratory experiments. The study can also help to detect HCO(-) in the cold and dense regions of the interstellar medium where the anion could be formed through the association of abundant CO with still unobserved H(-).

Dynamics of the O(3P) + CHD3(vCH = 0,1) reactions on an accurate ab initio potential energy surface

Recent experimental and theoretical studies on the dynamics of the reactions of methane with F and Cl atoms have modified our understanding of mode-selective chemical reactivity. The O + methane reaction is also an important candidate to extend our knowledge on the rules of reactivity. Here, we report a unique full-dimensional ab initio potential energy surface for the O((3)P) + methane reaction, which opens the door for accurate dynamics calculations using this surface. Quasiclassical trajectory calculations of the angular and vibrational distributions for the ground state and CH stretching excited O + CHD(3)(v(1) = 0,1) → OH + CD(3) reactions are in excellent agreement with the experiment. Our theory confirms what was proposed experimentally: The mechanistic origin of the vibrational enhancement is that the CH-stretching excitation enlarges the reactive cone of acceptance.