Q2DTor: A program to treat torsional anharmonicity through coupled pair torsions in flexible molecules

Abstract The Q 2DTor program (Quantum 2-Dimensional Torsions) is designed to calculate accurate rotational–vibrational partition functions (also called rovibrational partition functions) and thermodynamic functions for molecular systems having two [1] or more torsions. Systems with more than two torsions can also be studied by treating the torsions by pairs. The program searches for all the torsional conformers and evaluates the rovibrational partition function using the multi-structural harmonic oscillator (MS-HO) approximation and the extended two-dimensional torsion (E2DT) approximation. The latter incorporates full coupling of the two torsions by means of the two-dimensional non-separable (2D-NS) approximation [2], and it also includes their influence on the remaining degrees of freedom. The program also calculates the ideal gas-phase standard-state thermodynamic functions at the requested temperatures. Twenty molecules have been used to test Q 2DTor. Program summary Program Title: Q2DTor Program Files doi: http://dx.doi.org/10.17632/wbgchgc2kp.1 Licensing provisions: GNU GPL v3 Programming language: Python 2.7 Nature of problem: Calculation of accurate partition functions and thermodynamic functions in molecular systems involving two torsional modes. Torsional anharmonicity is treated quantically and includes full coupling in the kinetic and potential energies between the torsions and between the torsions and the rest of the degrees of freedom. Solution method: The program uses the variational method to solve the Schrodinger equation of a two-dimensional torsional potential using Fourier series. All of the remaining degrees of freedom (non-torsional) are incorporated through a projected (the torsional modes are removed) rigid-rotator harmonic-oscillator partition function which is calculated at every torsional stationary point and that is allowed to vary with the torsional motion. The integration of the rovibrational partition function over the torsional space leads to a mixed quantum-classical vibrational partition function, which is transformed into a full quantum partition function by including the quantum contribution due to the torsions. For the evaluation of the integral, the rovibrational partition function at nonstationary points is carried out through a Delaunay triangulation procedure using the calculated rovibrational partition functions at the stationary points as nodes. Additional comments including Restrictions and Unusual features: The program is limited to two coupled torsional modes. References: [1] L. Simon-Carballido et al., J. Chem. Theory Comput. 13 (2017) 3478. [2] A. Fernandez-Ramos, J. Chem. Phys. 138 (2013) 134112.

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

[2]  Rubén Meana-Pañeda,et al.  Anharmonicity of Coupled Torsions: The Extended Two-Dimensional Torsion Method and Its Use To Assess More Approximate Methods. , 2017, Journal of chemical theory and computation.

[3]  María Luisa Senent,et al.  Ab Initio Study of the Rotational-Torsional Spectrum of Methyl Formate , 2005 .

[4]  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.

[5]  Kenneth S. Pitzer,et al.  Energy Levels and Thermodynamic Functions for Molecules with Internal Rotation I. Rigid Frame with Attached Tops , 1942 .

[6]  William H Green,et al.  Intramolecular hydrogen migration in alkylperoxy and hydroperoxyalkylperoxy radicals: accurate treatment of hindered rotors. , 2010, The journal of physical chemistry. A.

[7]  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.

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

[9]  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..

[10]  Sean C. Smith.  Angular‐momentum resolution in transitional‐mode state counting for loose transition states , 1992 .

[11]  Ekaterina I Izgorodina,et al.  How Accurate Are Approximate Methods for Evaluating Partition Functions for Hindered Internal Rotations? , 2008, The journal of physical chemistry. A.

[12]  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.

[13]  Jim Pfaendtner,et al.  The 1-D hindered rotor approximation , 2007 .

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

[15]  Michel Waroquier,et al.  Why does the uncoupled hindered rotor model work well for the thermodynamics of n-alkanes? , 2005 .

[16]  Peter Groner,et al.  Analysis of torsional spectra of molecules with two internal C3v rotors. II - Far infrared and low frequency Raman spectra of dimethylether isotopes , 1977 .

[17]  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..

[18]  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.

[19]  Donald G. Truhlar,et al.  Statistical thermodynamics of bond torsional modes , 2000 .

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

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

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