Calculation of rotational partition functions by an efficient Monte Carlo importance sampling technique

The evaluation of the classical rotational partition function represented by a configuration integral over all external and internal rotational degrees of freedom of nonrigid chain polyatomic molecules is described. The method of Pitzer and Gwinn is used to correct the classical partition function for quantum mechanical effects at low temperatures. The internal rotor hindrance and all coupling arising from the external and internal rotational degrees of freedom are explicitly taken into account. Importance sampling Monte Carlo based on the adaptive VEGAS algorithm to perform multidimensional integration is implemented within the TINKER program package. A multidimensional potential energy hypersurface is calculated with the MM3(2000) molecular mechanics force field. Numerical tests are performed on a number of small n-alkanes (from ethane to octane), for which the absolute entropies calculated at three different temperatures are compared both with the experimental values and with the previous theoretical results. The application of a more efficient importance sampling technique developed here results in a substantial reduction of statistical errors in the evaluation of the configuration integral for a given number of Monte Carlo steps. Error estimates for the calculated entropies are given, and possible sources of systematic errors, and their importance for a reliable prediction of the absolute entropy, are discussed.