Here, we present the algorithm and validation for OMEGA, a systematic, knowledge-based conformer generator. The algorithm consists of three phases: assembly of an initial 3D structure from a library of fragments; exhaustive enumeration of all rotatable torsions using values drawn from a knowledge-based list of angles, thereby generating a large set of conformations; and sampling of this set by geometric and energy criteria. Validation of conformer generators like OMEGA has often been undertaken by comparing computed conformer sets to experimental molecular conformations from crystallography, usually from the Protein Databank (PDB). Such an approach is fraught with difficulty due to the systematic problems with small molecule structures in the PDB. Methods are presented to identify a diverse set of small molecule structures from cocomplexes in the PDB that has maximal reliability. A challenging set of 197 high quality, carefully selected ligand structures from well-solved models was obtained using these methods. This set will provide a sound basis for comparison and validation of conformer generators in the future. Validation results from this set are compared to the results using structures of a set of druglike molecules extracted from the Cambridge Structural Database (CSD). OMEGA is found to perform very well in reproducing the crystallographic conformations from both these data sets using two complementary metrics of success.

This article shows how to evaluate rotational symmetry numbers for different molecular configurations and how to apply them to transition state theory. In general, the symmetry number is given by the ratio of the reactant and transition state rotational symmetry numbers. However, special care is advised in the evaluation of symmetry numbers in the following situations: (i) if the reaction is symmetric, (ii) if reactants and/or transition states are chiral, (iii) if the reaction has multiple conformers for reactants and/or transition states and, (iv) if there is an internal rotation of part of the molecular system. All these four situations are treated systematically and analyzed in detail in the present article. We also include a large number of examples to clarify some complicated situations, and in the last section we discuss an example involving an achiral diasteroisomer.

No abstract

Practical approximation schemes for calculating partition functions of torsional modes are tested against accurate quantum mechanical results for H(2)O(2) and six isotopically substituted hydrogen peroxides. The schemes are classified on the basis of the type and amount of information that is required. First, approximate one-dimensional hindered-rotator partition functions are benchmarked against exact one-dimensional torsion results obtained by eigenvalue summation. The approximate one-dimensional methods tested in this stage include schemes that only require the equilibrium geometries and frequencies, schemes that also require the barrier heights of internal rotation, and schemes that require the whole one-dimensional torsional potential. Then, three classes of approximate full-dimensional vibrational-rotational partition functions are calculated and are compared with the accurate full-dimensional path integral partition functions. These three classes are (1) separable approximations combining harmonic oscillator-rigid rotator models with the one-dimensional torsion schemes, (2) almost-separable approximations in which the nonseparable zero-point energy is used to correct the separable approximations, and (3) improved nonseparable Pitzer-Gwinn-type methods in which approaches of type 1 are used as reference methods in the Pitzer-Gwinn approach. The effectiveness of these methods for the calculation of isotope effects is also studied. Based on the results of these studies, the best schemes of each type are recommended for further use on systems where a corresponding amount of information is available.

Rate constants and (12)C/(13)C kinetic isotope effects are calculated by direct dynamics for the OH + CH(4) --> H(2)O + CH(3) reaction. The electronic structure calculations required to generate the implicit potential energy surface were carried out by the high-level multicoefficient Gaussian-3/version-3 (MCG3) method and compared to two other multilevel methods, MC3BB and MC3MPW, and three density functional methods, M06-2X, BB1K, and MPW1K. The rate constants and (12)C/(13)C kinetic isotope effects are shown to depend strongly on the coordinate system used to calculate the frequencies as well as on the method used to account for the torsional anharmonicity of the lowest-frequency vibrational mode of the generalized transition states.

Rate coefficients are calculated using canonical variational transition state theory with multidimensional tunneling (CVT/SCT) for the reactions H + H2O2 --> H2O + OH (1a) and H + H2O2 --> HO2 + H2 (1b). Reaction barrier heights are determined using two theoretical approaches: (i) comparison of parametrized rate coefficient calculations employing CVT/SCT to experiment and (ii) high-level ab initio methods. The evaluated experimental data reveal considerable variations of the barrier height for the first reaction: although the zero-point-exclusive barrier for (1a) derived from the data by Klemm et al. (First Int. Chem. Kinet. Symposium 1975, 61) is 4.6 kcal/mol, other available measurements result in a higher barrier of 6.2 kcal/mol. The empirically derived zero-point-exclusive barrier for (1b) is 10.4 kcal/mol. The electronic structure of the system at transition state geometries in both reactions was found to have "multireference" character; therefore special care was taken when analyzing electronic structure calculations. Transition state geometries are optimized by multireference perturbation theory (MRMP2) with a variety of one-electron basis sets, and by a multireference coupled cluster (MR-AQCCSD) method. A variety of single-reference benchmark-level calculations have also been carried out; included among them are BMC-CCSD, G3SX(MP3), G3SX, G3, G2, MCG3, CBS-APNO, CBS-Q, CBS-QB3, and CCSD(T). Our data obtained at the MRMP2 level are the most complete; the barrier height for (1a) using MRMP2 at the infinite basis set limit is 4.8 kcal/mol. Results are also obtained with midlevel single-reference multicoefficient correlation methods, such as MC3BB, MC3MPW, MC-QCISD/3, and MC-QCISD-MPWB, and with a variety of hybrid density functional methods, which are compared with high-level theory. On the basis of the evaluated experimental values and the benchmark calculations, two possible recommended values are given for the rate coefficients.

Rate constants for the OH + H2S --> H2O + HS reaction, which is important for both atmospheric chemistry and combustion, are calculated by direct dynamics with the M06-2X density functional using the MG3S basis set. Energetics are compared to high-level MCG3/3//MC-QCISD/3 wave function theory and to results obtained by other density functionals. We employ canonical variational transition-state theory with multidimensional tunneling contributions and scaled generalized normal-mode frequencies evaluated in redundant curvilinear coordinates with anharmonicity included in the torsion. The transition state has a quantum mechanically distinguishable, nonsuperimposable mirror image that corresponds to a separate classical reaction path; the effect of the multiple paths is examined through use of a symmetry number and by torsional methods. Calculations with the reference-potential Pitzer-Gwinn treatment of the torsional mode agree with experiment, within experimental scatter, and predict a striking temperature dependence of the activation energy, increasing from -0.1 kcal/mol at 200 K to 0.2, 1.0, 3.4, and 9.8 kcal/mol at 300, 500, 1000, and 2400 K. The unusual temperature dependence arises from a dynamical bottleneck at an energy below reactants, following an addition complex on the reaction path with a classical binding energy of 4.4 kcal/mol. As a way to check the mechanism, kinetic isotope effects of the OH + D2S and OD + D2S reactions have been predicted.

We report here a theoretical study of the 13C kinetic isotope effect (KIE) and its temperature dependence for the reaction OH + CH4 --> H2O + CH3, the major sink of atmospheric methane in the troposphere. The KIE values at various atmospherically significant temperatures were determined by direct dynamics using variational transition state theory with multidimensional tunneling contributions (VTST/MT). The potential energy surfaces (PESs) were generated by hybrid density functional theory as well as by recently developed doubly hybrid density functional theory methods. Comparisons of our calculated KIEs with experimental data and theoretical values in the literature reveal the critical contributions due to multidimensional tunneling and torsion anharmonicity as well as the critical issue of the choice of internal rotational axis.