Anharmonicity perturbation theory

Starting from the normal mode approximation, one can introduce anharmonicity in different ways. Anharmonic perturbation theory [206] and local mode models [204] may be useful in some cases, where anharmonic effects are small or mostly diagonal. Vibrational self-consistent-field and configuration-interaction treatments [207, 208] can also be powerful and offer a hierarchy of approximation levels. Even more rigorous multidimensional treatments include variational calculations [209], diffusion quantum Monte Carlo, and time-dependent Hartree approaches [210]. 

The anharmonic terms, i.e. the cubic and higher terms in the displacement expansion of the intermolecular potential and the rotational kinetic energy terms, which are neglected in the harmonic Hamiltonian, can be considered as perturbations. They affect the vibrational excitations of the crystal in two ways they shift the excitation frequencies and they lead to finite lifetimes of the excited states, which are visible as spectral line broadening. By means of anharmonic perturbation theory based on a Green s function approach [64, 65] it is possible to calculate the frequency shifts, as well as the line widths. 

Figure 7. The striking property of these curves is that they are not of the form prescribed by leading-order anharmonic perturbation theory. From leading-order anharmonic perturbation theory, in the high temperature region, both and should go as, with volume-dependent coefficients. Our results for and cannot be fitted to T curves, within the estimated statistical errors of our MD averages.

The quasiharmonic part of S contains the temperature-independent term 3Nk(l - in0). To determine, we assume that contains no temperature-independent term. Such a term would follow from a contribution to of the form Aj T, where Aj may be volume-dependent. From our MD calculations of and P, we would not be able to determine A. However, anharmonic perturbation theory suggests that Fy in the high-temperature region should be a power series in T, of lowest order T. Therefore, to analyze the sodium calculations, we write... 

The interatomic potential function for the diatomic molecule was described in Section 6 5. In the Taylpr-series development of this function (6-72)3 cubic and higher terms were neglected in the harmonic approximation. It is now of interest to evaluate the importance of these so-called anharmonic terms with the aid of the perturbation theory outlined above. If cubic and quartic... 

To invoke the perturbation theory for a small anharmonic coupling coefficient, we use the Wick theorem for the coupling of the creation and annihilation operators of low-frequency modes in expression (A3.19). Retaining the terms of the orders y and y2, we are led to the following expressions for the shift AQ and the width 2T of the high-frequency vibration spectral line 184... 

Although the harmonic ZPVE must always be taken into account in the calculation of AEs, the anharmonic contribution is much smaller (but oppositely directed) and may sometimes be neglected. However, for molecules such as H2O, NH3, and CH4, the anharmonic corrections to the AEs amount to 0.9, 1.5, and 2.3 kJ/mol and thus cannot be neglected in high-precision calculations of thermochemical data. Comparing the harmonic and anharmonic contributions, it is clear that a treatment that goes beyond second order in perturbation theory is not necessary as it would give contributions that are small compared with the errors in the electronic-structure calculations. 

M. Quack In answer to the question by Prof. Kellman on exact analytical treatments of anharmonic resonance Hamiltonians, I might point out that to the best of my knowledge no fully satisfactory result beyond perturbation theory is known. Interesting efforts concern very recent perturbation theories by Sibert and co-workers and by Duncan and co-workers as well as by ourselves using as starting point internal coordinate Hamiltonians, normal coordinate Hamiltonians, and perhaps best, Fermi modes [1]. Of course, Michael Kellman himself has contributed substantial work on this question. Although all the available analytical results are still rather rough approximations, one can always... 

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... 

The emphasis of the theoretical discussion is (1) derivation and interpretation of the sum on states perturbation theory for charge polarization (2) development of physical models for the hyperpolarizability to assist molecular design (e.g., reduction of molecular orbital representations to the corresponding anharmonic oscillator description for hyperpolarizability). 

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]. 

We will of course be rather more focused here. We shall be concerned with the generic computational strategies needed to address the problems of phase behavior. The physical context we shall explore will not extend beyond the structural organization of the elementary phases (liquid, vapor, crystalline) of matter, although the strategies are much more widely applicable than this. We shall have nothing to say about a wide spectrum of techniques (density functional theory [1], integral equation theories [2], anharmonic perturbation... 

In the lowest order of perturbation theory, the energy levels of the three-dimensional anharmonic oscillator are... 

A non-perturbative theory of the multiphonon relaxation of a localized vibrational mode, caused by a high-order anharmonic interaction with the nearest atoms of the crystal lattice, is proposed. It relates the rate of the process to the time-dependent non-stationary displacement correlation function of atoms. A non-linear integral equation for this function is derived and solved numerically for 3- and 4-phonon processes. We have found that the rate exhibits a critical behavior it sharply increases near a specific (critical) value(s) of the interaction. 

To sum up, we have developed a general non-perturbative method that allows one to calculate the rate of relaxation processes in conditions when perturbation theory is not applicable. Theories describing non-radiative electronic transitions and multiphonon relaxation of a local mode, caused by a high-order anharmonic interaction have been developed on the basis of this method. In the weak coupling limit the obtained results agree with the predictions of the standard perturbation theory. 

The second term in aB arises from second-order perturbation theory it may be described as due to the fact that cubic anharmonicity results in a mean first-power displacement in q given by... 

In practice, the result of the perturbation treatment may be expressed as a series of formulae for the spectroscopic constants, i.e. the coefficients in the transformed or effective hamiltonian, in terms of the parameters appearing in the original hamiltonian, i.e. the wavenumbers tor, the anharmonic force constants , the moments of inertia Ia, their derivatives eft , and the zeta constants These formulae are analogous to equations (23)—(27) for a diatomic molecule. They are too numerous and too complicated to quote all of them here, but the various spectroscopic constants are listed in Table 3, with their approximate relative orders of magnitude, an indication of which parameters occur in the formula for each spectroscopic constant, and a reference to an appropriate source for the perturbation theory formula for that constant. 

A similar interaction would be observed between all Fermi polyads containing sets of vibrational levels related by the selection rule A tv = 2, A tv = +1, and the hamiltonian matrix should be diagonalized for each Fermi polyad without the use of perturbation theory. If, on the other hand, the interaction (63) were smaller, or the separation between the unperturbed levels were larger, the interaction could be treated by perturbation theory it can be shown that, in second-order perturbation theory, equation (63) would contribute a term to the vibrational anharmonic constants... 

All this applies to weak and medium strong H-bonds like those encountered for alcohols and many other systems up to carboxylic acid dimers or about 32-42 kJ/mol. (8 or 10 kcal/mol.) Unfortunately vibrational spectra of systems with very strong H-bonds could, with a few exceptions, only be measured in condensed phases. Factors that come in when such systems are examined are potential surfaces with two minima with, in certain cases, the possibility of tunnelling, or flat single minima Most of these systems are likely to be so anharmonic that second order perturbation theory breaks down and the concept of normal vibrations becomes itself question-nable. Many such systems are highly polarizable and are strongly influenced by the environment yielding extremely broad bands 92). Bratos and Ratajczak 93) has shown that even such systems can be handled by relaxation theories. 

An alternative approach widely used in polyatomic molecule studies is based on the Golden Rule and a perturbative treatment of the anharmonic coupling (57,62). This approach is not much used for diatomic molecules. In the liquid O2 example cited above, the Hamiltonian must be expanded to 30th order or so to calculate the multiphonon emission rate. But for vibrations of polyatomic molecules, which can always find relatively low-order VER pathways for each VER step, perturbation theory is very useful. In the perturbation approach, the molecule s entire ladder of vibrational excitations is the system and the phonons are the bath. Only lower-order processes are ordinarily needed (57) because polyatomic molecules have many vibrations ranging from higher to lower frequencies and only a small number of phonons, usually one or two, are excited in each VER step. The usual practice is to expand the interaction Hamiltonian (qn, Q) in Equation (2) in powers of normal coordinates (57,62) ... 

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 ... 

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. 

Later work evaluated the two-dimensional potential energy surface using various correlation treatments including many-body perturbation theory and coupled cluster techniques Evaluation of the vibrational spectrum was explicitly anharmonic in nature, mak-... 

Of course, this corresponds to an adiabatic potential-energy surface with two potential-energy minima separated by a well-defined barrier. Note that y A contains anharmonicities induced into Xr by D/A interactions as well as energy corrections that originate from D/A coupling. To some extent (i.e., as in first-order perturbation theory or in a Taylor s series expansion around the diabatic values of... 

---