Anharmonicity of Coupled Torsions: The Extended Two-Dimensional Torsion Method and Its Use To Assess More Approximate Methods.

In this work we present the extended two-dimensional torsion (E2DT) method and use it to analyze the performance of several methods that incorporate torsional anharmonicity more approximately for calculating rotational-vibrational partition functions. Twenty molecules having two hindered rotors were studied for temperatures between 100 and 2500 K. These molecules present several kinds of situations; they include molecules with nearly separable rotors, molecules in which the reduced moments of inertia change substantially with the internal rotation, and molecules presenting compound rotation. Partition functions obtained by the rigid-rotor harmonic oscillator approximation, a method involving global separability of torsions and the multistructural methods without explicit potential coupling [MS-T(U)] and with explicit potential coupling [MS-T(C)] of torsions, are compared to those obtained with a quantized version - called the extended two-dimensional torsion (E2DT) method - of the extended hindered rotor approximation of Vansteenkiste et al. ( Vansteenkiste et al. J. Chem. Phys. 2006 , 124 , 044314 ). In the E2DT method, quantum effects due to the torsional modes were incorporated by the two-dimensional nonseparable method, which is a method that is based on the solution of the torsional Schrödinger equation and that includes full coupling in both the kinetic and potential energy. By comparing other methods to the E2DT method and to experimental thermochemical data, this study concludes that the harmonic approximation yields very poor results at high temperatures; the global separation of torsions from the rest of the degrees of freedom is not justified even when an accurate method to treat the torsions is employed; it is confirmed that methods based on less complete potential energy coupling of torsions, such as MS-T(U), are not accurate when dealing with rotors with different barrier heights, and more complete inclusion of torsional coupling to the method in MS-T(C) improves substantially the results in such a way that it could be used in cases where the E2DT method is unaffordable.

[1]  D. Truhlar,et al.  Computational study of the reactions of methanol with the hydroperoxyl and methyl radicals. 2. Accurate thermal rate constants. , 2010, The journal of physical chemistry. A.

[2]  Benjamin A. Ellingson,et al.  Statistical thermodynamics of bond torsional modes: tests of separable, almost-separable, and improved Pitzer-Gwinn approximations. , 2006, The Journal of chemical physics.

[3]  Antonio Fernández-Ramos,et al.  Calculation of the two-dimensional non-separable partition function for two molecular systems , 2014, Journal of Molecular Modeling.

[4]  D. Truhlar,et al.  Including Torsional Anharmonicity in Canonical and Microcanonical Reaction Path Calculations. , 2013, Journal of chemical theory and computation.

[5]  Donald G. Truhlar,et al.  Hybrid Meta Density Functional Theory Methods for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions: The MPW1B95 and MPWB1K Models and Comparative Assessments for Hydrogen Bonding and van der Waals Interactions , 2004 .

[6]  K. Pitzer,et al.  Thermodynamic Functions for Molecules with Internal Rotation , 1941 .

[7]  Rubén Meana-Pañeda,et al.  Tunneling and conformational flexibility play critical roles in the isomerization mechanism of vitamin D. , 2012, Journal of the American Chemical Society.

[8]  Jorge M. C. Marques,et al.  Symmetry numbers and chemical reaction rates , 2007 .

[9]  Kenneth S. Pitzer,et al.  Energy Levels and Thermodynamic Functions for Molecules with Internal Rotation: II. Unsymmetrical Tops Attached to a Rigid Frame , 1946 .

[10]  Robert J. Renka,et al.  Algorithm 751: TRIPACK: a constrained two-dimensional Delaunay triangulation package , 1996, TOMS.

[11]  Kenneth L. Clarkson,et al.  MSTor: A program for calculating partition functions, free energies, enthalpies, entropies, and heat capacities of complex molecules including torsional anharmonicity , 2012, Comput. Phys. Commun..

[12]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[13]  Donald G Truhlar,et al.  Multi-structural thermodynamics of C-H bond dissociation in hexane and isohexane yielding seven isomeric hexyl radicals. , 2011, Physical chemistry chemical physics : PCCP.

[14]  V Van Speybroeck,et al.  An extended hindered-rotor model with incorporation of Coriolis and vibrational-rotational coupling for calculating partition functions and derived quantities. , 2006, The Journal of chemical physics.

[15]  Donald G. Truhlar,et al.  MSTor version 2013: A new version of the computer code for the multi-structural torsional anharmonicity, now with a coupled torsional potential , 2013, Comput. Phys. Commun..

[16]  Donald G Truhlar,et al.  Computational Thermochemistry: Scale Factor Databases and Scale Factors for Vibrational Frequencies Obtained from Electronic Model Chemistries. , 2010, Journal of chemical theory and computation.

[17]  Donald G Truhlar,et al.  Practical methods for including torsional anharmonicity in thermochemical calculations on complex molecules: the internal-coordinate multi-structural approximation. , 2011, Physical chemistry chemical physics : PCCP.

[18]  Kenneth S. Pitzer,et al.  Energy Levels and Thermodynamic Functions for Molecules with Internal Rotation. III. Compound Rotation , 1949 .

[19]  Donald G Truhlar,et al.  Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions. , 2006, Journal of chemical theory and computation.

[20]  A. Fernández-Ramos,et al.  Hindered rotor tunneling splittings: an application of the two-dimensional non-separable method to benzyl alcohol and two of its fluorine derivatives. , 2016, Physical chemistry chemical physics : PCCP.

[21]  Donald G. Truhlar,et al.  A simple approximation for the vibrational partition function of a hindered internal rotation , 1991 .

[22]  D. Truhlar,et al.  Multipath variational transition state theory: rate constant of the 1,4-hydrogen shift isomerization of the 2-cyclohexylethyl radical. , 2012, The journal of physical chemistry. A.

[23]  Antonio Fernández-Ramos,et al.  Accurate treatment of two-dimensional non-separable hindered internal rotors. , 2013, The Journal of chemical physics.

[24]  D. Truhlar,et al.  Quantum Thermochemistry: Multistructural Method with Torsional Anharmonicity Based on a Coupled Torsional Potential. , 2013, Journal of chemical theory and computation.

[25]  Yvan Notay,et al.  JADAMILU: a software code for computing selected eigenvalues of large sparse symmetric matrices , 2007, Comput. Phys. Commun..

[26]  A. Fernández-Ramos,et al.  Accounting for conformational flexibility and torsional anharmonicity in the H + CH3CH2OH hydrogen abstraction reactions: a multi-path variational transition state theory study. , 2014, The Journal of chemical physics.