Abstract:The accurate prediction of thermochemistry and kinetic parameters is an important task for reaction modeling. Unfortunately, the commonly used harmonic oscillator model is often not accurate enough due to the absence of anharmonic effects. In this work, we improve the representation of an anharmonic potential energy surface (PES) using uncoupled mode (UM) approximations, which model the full-dimensional PES as a sum of one-dimensional potentials of each mode. We systematically analyze different PES sampling schemes and coordinate systems for constructing the one-dimensional potentials, and benchmark the performance of UM methods on data sets of molecular thermochemistry and kinetic properties. The results show that the accuracy of the UM approach strongly depends on how the one-dimensional potentials are defined. If one-dimensional potentials are constructed by sampling along normal mode directions (UM-N) or along the directions that minimize intermode coupling (E- and E'-optimized), the accuracies of the predicted properties are not significantly improved compared to the harmonic oscillator model. However, significant improvements can be achieved by sampling the torsional modes separately from the vibrational modes (UM-T and UM-VT). In this work, three types of coordinate systems are examined, including redundant internal coordinates (RIC), hybrid internal coordinates (HIC), and translation-rotation-internal coordinates (TRIC). The HIC and TRIC coordinate systems can outperform RIC since transition state species may contain large-amplitude interfragmentary motions that regular internal coordinates can not describe adequately. Among all the methods we examined, the activation energies and pre-exponential factors calculated using UM-VT with either TRIC or HIC best agree with the reference values. Since UM-VT requires only a number of additional single point energy calculations for each independent mode, the scaling of computational costs of UM-VT is the same as that of the standard harmonic oscillator model, making UM-VT an appealing way of calculating the thermochemistry and kinetic properties for large-size systems.

Efficient Generation of Torsional Energy Profiles by Multifidelity Gaussian Processes for Hindered-Rotor Corrections

Efficient Generation of Torsional Energy Profiles by Multifidelity Gaussian Processes for Hindered Rotor Corrections

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

Practical Methods for Including Torsional Anharmonicity in Thermochemical Calculations on Complex Molecules: the Internal-Coordinate Multi-Structural Approximation.

Approximate Reconstruction of Torsional Potential Energy Surface based on Voronoi Tessellation

Identification of Torsional Modes in Complex Molecules Using Redundant Internal Coordinates: the Multistructural Method with Torsional Anharmonicity with a Coupled Torsional Potential and Delocalized Torsions.

Efficient vibrationally correlated calculations using n -mode expansion-based kinetic energy operators

Comparative Analysis of Uncoupled Mode Approximations for Molecular Thermochemistry and Kinetics

Gaussian process regression models for the properties of micro-tearing modes in spherical tokamaks

Aerodynamic interference analysis and modeling of heterogeneous rotors

Multiobjective optimization for flapping foil hydrodynamics with a multitask and multifidelity approach

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations

A path integral approach to molecular thermochemistry

Accurate Entropy Calculation for Large Flexible Hydrocarbons Using a Multi-Structural 2-Dimensional Torsion Method

Force-free identification of minimum-energy pathways and transition states for stochastic electronic structure theories

Improved Prediction of Aerodynamic Loss Propagation as Entropy Rise in Wind Turbines Using Multifidelity Analysis

Multi-Fidelity Simulation of Gas Turbine Overall Performance by Directly Coupling High-Fidelity Models of Multiple Rotating Components

Gaussian Process Regression models for the properties of micro-tearing modes in spherical tokamak

Quantum Thermochemistry: Multistructural Method with Torsional Anharmonicity Based on a Coupled Torsional Potential.

Accurate Vibrational and Ro-Vibrational Contributions to the Properties of Large Molecules by a New Engine Employing Curvilinear Internal Coordinates and Vibrational Perturbation Theory to Second Order

Efficient time-dependent vibrational coupled cluster computations with time-dependent basis sets at the two-mode coupling level: full and hybrid TDMVCC[2]