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.
ABSTRACT
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.
ABSTRACT
Accurate quantum mechanical partition functions and absolute free energies of H(2)O(2) 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 2,400 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 revisions to the JANAF tables.
