A review of the theoretical and computational modeling of bimolecular reactions is given. The review is divided into several sections which are as follows: gas-phase thermal reactions; gas-phase state-selected reactions and product state distributions; and condensed-phase bimolecular reactions. The section on gas-phase thermal reactions covers the enthalpies and free energies of reaction, kinetics, saddle points and potential energy surfaces, rate theory for simple barrier reactions and bimolecular reactions over potential wells. The section on gas-phase state-selected reactions focuses on electronically adiabatic reactions and electronically nonadiabatic reactions. Finally, the section on condensed-phase bimolecular reactions covers reactions in liquids, reactions on surfaces and in solids and tunneling at low temperature.

Modeling the Kinetics of Bimolecular Reactions

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car–Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

/pdf/computational-electrochemistry-prediction-of-liquid-phase-1djdtuew4j.pdf

Computational electrochemistry: prediction of liquid-phase reduction potentials

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.

/pdf/symmetry-numbers-and-chemical-reaction-rates-5f8lq90m03.pdf

Symmetry numbers and chemical reaction rates

This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications. The theoretical methods reviewed here include multidimensional quantum mechanical tunneling, multistructural VTST (MS-VTST), multi-path VTST (MP-VTST), both reaction-path VTST (RP-VTST) and variable reaction coordinate VTST (VRC-VTST), system-specific quantum Rice–Ramsperger–Kassel theory (SS-QRRK) for predicting pressure-dependent rate constants, and VTST in the solid phase, liquid phase, and enzymes. We also provide some perspectives regarding the general applicability of VTST.

Variational transition state theory: theoretical framework and recent developments

Many methods for correcting harmonic partition functions for the presence of torsional motions employ some form of one-dimensional torsional treatment to replace the harmonic contribution of a specific normal mode. However, torsions are often strongly coupled to other degrees of freedom, especially other torsions and low-frequency bending motions, and this coupling can make assigning torsions to specific normal modes problematic. Here, we present a new class of methods, called multi-structural (MS) methods, that circumvents the need for such assignments by instead adjusting the harmonic results by torsional correction factors that are determined using internal coordinates. We present three versions of the MS method: (i) MS-AS based on including all structures (AS), i.e., all conformers generated by internal rotations; (ii) MS-ASCB based on all structures augmented with explicit conformational barrier (CB) information, i.e., including explicit calculations of all barrier heights for internal-rotation barriers between the conformers; and (iii) MS-RS based on including all conformers generated from a reference structure (RS) by independent torsions. In the MS-AS scheme, one has two options for obtaining the local periodicity parameters, one based on consideration of the nearly separable limit and one based on strongly coupled torsions. The latter involves assigning the local periodicities on the basis of Voronoi volumes. The methods are illustrated with calculations for ethanol, 1-butanol, and 1-pentyl radical as well as two one-dimensional torsional potentials. The MS-AS method is particularly interesting because it does not require any information about conformational barriers or about the paths that connect the various structures.

/pdf/practical-methods-for-including-torsional-anharmonicity-in-2kgx55k002.pdf

Practical methods for including torsional anharmonicity in thermochemical calculations on complex molecules: the internal-coordinate multi-structural approximation.

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.

/pdf/statistical-thermodynamics-of-bond-torsional-modes-tests-of-4avlp3l756.pdf

Statistical thermodynamics of bond torsional modes: tests of separable, almost-separable, and improved Pitzer-Gwinn approximations.

Accurate quantum mechanical partition functions and absolute free energies of H2O2 are determined using a realistic potential energy surface [J. Koput, S. Carter, and N. C. Handy, J. Phys. Chem. A 102, 6325 (1998)] for temperatures ranging from 300 to 2400 K by using Monte Carlo path integral calculations with new, efficient polyatomic importance sampling methods. The path centroids are sampled in Jacobi coordinates via a set of independent ziggurat schemes. The calculations employed enhanced-same-path extrapolation of trapezoidal Trotter Fourier path integrals, and the paths were constructed using fast Fourier sine transforms. Importance sampling was also used in Fourier coefficient space, and adaptively optimized stratified sampling was used in configuration space. The free energy values obtained from the path-integral calculations are compared to separable-mode approximations, to the Pitzer–Gwinn approximation, and to values in thermodynamic tables. Our calculations support the recently proposed revisi...

/pdf/accurate-vibrational-rotational-partition-functions-and-1giytaf7gu.pdf

Accurate vibrational-rotational partition functions and standard-state free energy values for H2O2 from Monte Carlo path-integral calculations.

Accurate quantum mechanical (QM) vibrational-rotational partition functions for HOOD, D(2)O(2), H(18)OOH, H(2)(18)O(2), D(18)OOH, and H(18)OOD are determined using a realistic potential energy surface for temperatures ranging from 300 to 2400 K by using the TT-FPI-ESPE path-integral Monte Carlo method. These data, together with our prior results for H(2)O(2), provide benchmarks for testing approximate methods of estimating isotope effects for systems with torsional motions. Harmonic approximations yield poor accuracy for these systems, and although the well-known Pitzer-Gwinn (PG) approximation provides better results for absolute partition functions, it yields the same results as the harmonic approximation for isotope effects because these are intrinsically quantal phenomena. We present QM generalizations of the PG approximation that can provide high accuracy for both isotope effects and absolute partition functions. These approximations can be systematically improved until they approach the accurate result and converge rapidly. These methods can also be used to obtain affordable estimates of zero-point energies from accurate partition functions-even those at relatively high temperatures.

High-precision quantum thermochemistry on nonquasiharmonic potentials : Converged path-integral free energies and a systematically convergent family of generalized Pitzer-Gwinn approximations

Erratum: High-precision quantum thermochemistry on nonquasiharmonic potentials: Converged path-integral free energies and a systematically convergent family of generalized pitzer-gwinn approximations (Journal of Physical Chemistry A (2005) 109A)

Vanessa Audette Lynch

Papers

Statistical thermodynamics of bond torsional modes: tests of separable, almost-separable, and improved Pitzer-Gwinn approximations.

Accurate vibrational-rotational partition functions and standard-state free energy values for H2O2 from Monte Carlo path-integral calculations.

High-precision quantum thermochemistry on nonquasiharmonic potentials : Converged path-integral free energies and a systematically convergent family of generalized Pitzer-Gwinn approximations

Erratum: High-precision quantum thermochemistry on nonquasiharmonic potentials: Converged path-integral free energies and a systematically convergent family of generalized pitzer-gwinn approximations (Journal of Physical Chemistry A (2005) 109A)