Search

We present a new application for the Torsional Path Integral Monte Carlo (TPIMC) method in which the TPI partition functions are introduced into the calculation of Transition State Theory (TST) rate constants. In this way, an explicit treatment of torsional anharmonicity is included in the TST calculations and the magnitude of these effects can be assessed. The new method is tested on the C(2)H(6) + H hydrogen abstraction reaction and concerted hydrogen transfer in the carbonic acid dimer, for which we have developed torsional potential energy surfaces. For the C(2)H(6) + H reaction the rate constants are halved at room temperature on including a treatment of torsional anharmonicity, while the effects are found to be much smaller for the hydrogen transfer reaction in the carbonic acid dimer.

In this paper we calculate the relative conformer populations of the tryptamine molecule. Our approach combines high level electronic structure conformer energies with harmonic frequencies and an anharmonic treatment of the torsional motions using the torsional path integral Monte Carlo method. We have developed a 3-D potential energy surface as a function of the torsional coordinates at the B3LYP/6-31+G(d) level using 2535 grid points. Eight conformers of tryptamine were found to be significantly populated at 430 K as opposed to the experimental observation of seven. This, along with further comparisons with various experimental data, leads us to suppose that conformer interconversion occurs during the cooling phases of many of the experiments. The ordering of the calculated populations fits well with available experimental data. Torsional anharmonicity is found to affect conformer populations more significantly at 430 K than at 100 K (although overall the effects are small), while quantum mechanical effects are not important at either temperature.

Schemes to include a treatment of torsional anharmonicity in the conformational analysis of biological molecules are introduced. The approaches combine ab initio electronic energies and harmonic frequencies with anharmonic torsional partition functions calculated using the torsional path integral Monte Carlo method on affordable potential energy surfaces. The schemes are applied to the conformational study of the monosaccharide beta-d-galactose in the gas phase. The global minimum structure is almost exclusively populated at 100 K, but a large number of conformers are present at ambient and higher temperatures. Both quantum mechanical and anharmonic effects in the torsional modes have little effect on the populations at all temperatures considered, and it is, therefore, expected that standard harmonic treatments are satisfactory for the conformational study of monosaccharides.