Anharmonicity quantum corrections

The approach to the evaluation of vibrational spectra described above is based on classical simulations for which quantum corrections are possible. The incorporation of quantum effects directly in simulations of large molecular systems is one of the most challenging areas in theoretical chemistry today. The development of quantum simulation methods is particularly important in the area of molecular spectroscopy for which quantum effects can be important and where the goal is to use simulations to help understand the structural and dynamical origins of changes in spectral lineshapes with environmental variables such as the temperature. The direct evaluation of quantum time- correlation functions for anharmonic systems is extremely difficult. Our initial approach to the evaluation of finite temperature anharmonic effects on vibrational lineshapes is derived from the fact that the moments of the vibrational lineshape spectrum can be expressed as functions of expectation values of positional and momentum operators. These expectation values can be evaluated using extremely efficient quantum Monte-Carlo techniques. The main points are summarized below. 

In an attempt both to include quantum corrections and to treat anharmonic effects, Etters, Kanney, Gillis, and Kaelberger have calculated the thermodynamic properties of clusters using the self-consistent phonon approximation. - Within this approximation they calculated internal energies, free energies, specific heats, and entropies for argon clusters. 

R. Ramrez, T. Lopez-Ciudad, P. Kumar, and D. Marx (2004) Quantum corrections to classical time-correlation functions Hydrogen bonding and anharmonic floppy modes. J. Chem. Phys. 121, p. 3973... 

Methods to apply quantum corrections for this inadequacy of classical mechanics when applied to crystals, as well as other phases, have been explored (124-127), but there is not yet general agreement on how best to calculate thermodynamic properties for these types of systems. The practical solution of Berens and co-workers (124) entails. . calculating the velocity spectrum S(v) from molecular dynamics and then integrating S(v) over all frequencies with a weighting fimction which is the difference between the quantum and classical harmonic weighting functions for the thermodynamic variable of interest. However, as the statement avers, anharmonic contributions to the vibrational modes will not be captured correctly this relatively minor shortcoming was addressed by Hardy and co-workers (128). 

Eqs. (D.5)-(D.7). However, when perturbations occur due to anharmonicity, the wave functions in Eqs. (D.11)-(D.13) will provide the correct zeroth-order ones. The quantum numbers V2a and 2b are therefore not physically significant, while V2 and mV2 or V2 and I2 = mV2 are. It should also be pointed out that the degeneracy in the vibrational levels will be split due to anharmonicity [28]. 

Isotope effects on anharmonic corrections to ZPE drop off rapidly with mass and are usually neglected. The ideas presented above obviously carry over to exchange equilibria involving polyatomic molecules. Unfortunately, however, there are very few polyatomics on which spectroscopic vibrational analysis has been carried in enough detail to furnish spectroscopic values for Go and o)exe. For that reason anharmonic corrections to ZPE s of polyatomics have been generally ignored, but see Section 5.6.3.2 for a discussion of an exception also theoretical (quantum package) calculations of anharmonic constants are now practical (see above), and in the future one can expect more attention to anharmonic corrections of ZPE s. 

In theory, the wave equations of quantum mechanics can be used to derive near-correct potential-energy curves for molecular vibrations. Unfortunately, the mathematical complexity of these equations precludes quantitative application to all but the very simplest of systems. Qualitatively, the curves must take the anharmonic form. Such curves depart from harmonic behavior by varying degrees, depending on the nature of the bond and the atom involved. However, the harmonic and anharmonic curves are almost identical at low potential energies, which accounts for the success of the approximate methods described. 

It should be noted that (4.28) is only an approximation for the nuclear wave function. The perturbation terms (4.36) will mix into the nuclear wave function small contributions from harmonic-oscillator functions with quantum numbers other than v. These anharmonicity corrections to the vibrational wave function will add further to the probability of transitions with At) > 1. 

Actually, many other infrared transitions occur besides those allowed by the selection rule (6.74). The neglected terms in the expansion (6.66) will give transitions with a change of 2 or more in a given vibrational quantum number and transitions in which more than one vibrational quantum number changes moreover, anharmonicity corrections to the vibrational wave function will add to the probability of such transitions. Generally, the transitions (6.74) are the strongest. 

Classical anharmonic spring models with or without damping [9], and the corresponding quantum oscillator models seem well removed from the molecular problems of interest here. The quantum systems are frequently described in terms of coulombic or muffin tin potentials that are intrinsically anharmonic. We will demonstrate their correspondence after first discussing the quantum approach to the nonlinear polarizability problem. Since we are calculating the polarization of electrons in molecules in the presence of an external electric field, we will determine the polarized molecular wave functions expanded in the basis set of unperturbed molecular orbitals and, from them, the nonlinear polarizability. At the heart of this strategy is the assumption that perturbation theory is appropriate for treating these small effects (see below). This is appropriate if the polarized states differ in minor ways from the unpolarized states. The electric dipole operator defines the interaction between the electric field and the molecule. Because the polarization operator (eq lc) is proportional to the dipole operator, there is a direct link between perturbation theory corrections (stark effects) and electronic polarizability [6,11,12]. 

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. 

---