Abstract:Transition state theory (TST) based on activation parameters computed using quantum mechanics calculations combined with the polarizable continuum model (QM/PCM) is a fundamental tool for investigating reaction rates in the liquid phase. In conventional QM/PCM methods, thermodynamic data and partition functions for a solute are often derived from a quasi-ideal gas treatment (IGT) widely implemented in commercially available computation packages. This approach tends to overestimate entropy because calculations of thermodynamic parameters in the liquid phase ignore hindered translational and rotational modes in real solutions. The present work formulated partition functions for more realistic solutes hindered by surrounding solvent molecules in conjunction with the basic QM/PCM concept. In addition, a configuration partition function for solute molecules at a standard concentration of 1 mol dm -3 was incorporated using a simple lattice model. The canonical partition function and thermodynamic functions were derived based on statistical thermodynamics for non-localized systems. Expressions for rate coefficients within TST were also derived with a consistent formalism based on the standard state selected in partition function calculations. The performance of the proposed method was assessed by predicting rate coefficients for three different Diels-Alder reactions and comparing these with experimental results. QM/PCM calculations at the G4//ωB97X-D/6-311++G(d,p)/IEF-PCM level of theory were performed to obtain the electronic structures of stationary points on potential energy surfaces as a means of finding activation enthalpy, entropy and Gibbs energy values based on revised partition functions as well as predicting rate coefficients. This new technique demonstrated superior rate coefficient predictions compared with the IGT technique. The activation Gibbs energies obtained from our proposed method were lower than those obtained from the IGT method due to reasonable entropy computations. However, the obtained rate coefficients were one to two orders of magnitude lower than those obtained from the IGT method. This apparent paradox is attributed to the choice of standard states in the two approaches, which has a significant impact on the computation of the activation Gibbs energy.

Two-component transformation inclusive contraction scheme in the relativistic molecular orbital theory

A new computational framework for spinor-based relativistic exact two-component calculations using contracted basis functions

Improved Version of a Local Contracted Configuration Interaction of Singles and Doubles with Partial Inclusion of Triples and Quadruples

Unified Implementation of Relativistic Wave Function Methods: 4C-iCIPT2 as a Showcase

Open-shell Tensor Hypercontraction

Recovering four-component solutions by the inverse transformation of the infinite-order two-component wave functions

4-component Relativistic Calculations in a Multiwavelet Basis with Improved Convergence

Computation of Rate Coefficients in Solutions Based on Transition State Theory Combined with a Heuristically Corrected Polarizable Continuum Model: Intermolecular Diels-Alder Reactions as Case Studies

Exact Two-Component Relativistic Theory for NMR Parameters: General Formulation and Pilot Application.

Metric-Based Assessment Method for MS-T Formalism with Small Subsets of Torsional Conformers

Renormalized Internally-Contracted Multireference Coupled Cluster with Perturbative Triples

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

Squeezing in the Triatomic Linear Molecule Model Revealed by Virtue of the Iwop Technique

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

Time-dependent restricted-active-space self-consistent-field theory for laser-driven many-electron dynamics. II. Extended formulation and numerical analysis

Exact Two-Component Relativistic Theory for Nuclear Magnetic Resonance Parameters

A (T) Correction for Multicomponent Coupled-Cluster Theory

Double Exponential Transformation For Computing Three-Centre Nuclear Attraction Integrals

Unitary Basis Transformations in Mixed Quantum-Classical Dynamics

Fully self-consistent optimization of the Jastrow-Slater-type wave function using a similarity-transformed Hamiltonian

Double Configuration Interaction Singles: Scalable and size-intensive approach for orbital relaxation in excited states and bond-dissociation