Anharmonicity quantum mechanical calculations

Since the PE surfaces are anharmonic in this limit, the Franck-Condon contributions can be difficult to treat. The approaches to this limit can be separated into two classes (a) perturbation theory corrections of the weak-coupling limit and (b) quantum-mechanical calculations of the reaction coordinate. The latter tend to be done on a reaction-by-reaction basis this makes it difficult to generalize. The perturbation theory approaches have the advantage that they make use of the parameters used in the weak-coupling limit, and this can provide usefiil insights into general trends and patterns. 

S. Erkoc, J, N. Murrell, and D. C. Clary, Quantum mechanical calculation of collinear atom-triatom transition probabilities for anharmonic triatom potentials, Chem. Phys. Lett. 72 624 (1980). 

Of course, one is not really interested in classical mechanical calculations. Thus in normal practice the partition functions used in TST, as discussed in Chapter 4, are evaluated using quantum partition functions for harmonic frequencies (extension to anharmonicity is straightforward). On the other hand rotations and translations are handled classically both in TST and in VTST, which is a standard approximation except at very low temperatures. Later, by introducing canonical partition functions one can direct the discussion towards canonical variational transition state theory (CVTST) where the statistical mechanics involves ensembles defined in terms of temperature and volume. There is also a form of variational transition state theory based on microcanonical ensembles referred to by the symbol p,. Discussion of VTST based on microcanonical ensembles pVTST is beyond the scope of the discussion here. It is only mentioned that in pVTST the dividing surface is... 

Another hmitation is inherent to the harmonic approximation on which standard quantum mechanical force-field calculations are invariably based. Due to a fortui-tious (but surpisingly systematic) cancellation of errors, the harmonic frequencies calculated by modem density functional methods often match very well with the experimental ones, in spite of the fact that the latter involve necessarily more or less anharmonic potentials. Thus one is tempted to forget that the harmonic approx-imaton can become perilous when strong anharmonicity prevails along one or another molecular deformation coordinate. 

In the estimation of Acon(t), only the first two terms are considered, neglecting the higher-order terms. (Q - Goo) and (Q m - Goo) 810 die quantum mechanical expectation values of the anharmonic oscillator. They can be calculated using perturbation theory and is given by... 

As a rule the quantum-mechanical force-fields and the corresponding normal frequencies are calculated in a harmonic approximation, while the experimentally accessible frequencies are influenced by anharmonic contributions. The Puley s scaling factors are also found to incorporate the relevant empirical corrections for the vibrational anharmonicity. 

Calculation of Matrix Elements for One-Dimensional Quantum-Mechanical Problems and the Application to Anharmonic Oscillators. 

While the classical model of an anharmonic oscillator describes the effects of non-linearity, it cannot provide information on molecular properties. Calculation of molecular properties requires a quantum mechanical model. Application of perturbation theory (Boyd, 2003) leads to the following expression ... 

Because of their importance to nucleation kinetics, there have been a number of attempts to calculate free energies of formation of clusters theoretically. The most important approaches for the current discussion are harmonic models, " Monte Carlo studies, and molecular dynamics calcula-tions. In the harmonic model the cluster is assumed to be composed of constituent atoms with harmonic intermolecular forces. The most recent calculations, which use the harmonic model, have taken the geometries of the clusters to be those determined by the minimum in the two-body additive Lennard-Jones potential surface. The oscillator frequencies have been obtained by diagonalizing the Lennard-Jones force constant matrix. In the harmonic model the translational and rotational modes of the clusters are treated classically, and the vibrational modes are treated quantum mechanically. The harmonic models work best at low temjjeratures where anharmonic-ity effects are least important and the system is dominated by a single structure. 

Ab initio constmction of such anharmonic energy surfaces for quite complicated systems has been feasible for some years, thanks to increased compnting power and veiy time-efficient quantum-mechanical software. Many years before such calculations became possible, however, Marcus (1959, 1960, 1965) showed that, by treating the internal vibrational modes of D and and A and A , as symmetrised classical harmonic... 

The shape of the potential for the proton motion, particularly when the proton is engaged in hydrogen bond formation, is one of the most fascinating problems of molecular physics and chemistry. Because of very low mass the stretching protonic vibrations are characterized by high frequencies and, if independent of additional interactions, they are anharmonic. It is commonly known that expression of quantum-chemical calculations in the harmonic approximation, on various levels of the quantum-mechanical approach, needs the application of some scaling factors [1, 2]. For stretching protonic vibrations this factor is... 

Chaban GM, Gerber RB. Anharmonic vibrational spectroscopy of the glycine-water complex calculations for ab initio, empirical and hybrid quantum mechanics/molecular mechanics potentials. J Chem Phys 2001 115 1340-1348. 

---