Vibrations harmonic,/anharmonic

These harmonic-oscillator solutions predict evenly spaced energy levels (i.e., no anharmonicity) that persist for all v. It is, of course, known that molecular vibrations display anharmonicity (i.e., the energy levels move closer together as one moves to higher v) and that quantized vibrational motion ceases once the bond dissociation energy is reached. 

If we extend this last example to the modelling of molecular vibrations, we need to include additional terms in the differential equation to account for non-harmonic (anharmonic) forces. 

The harmonic oscillator provides equally spaced eigenenergies in molecular vibrations, additional anharmonic contributions (V = bx3 + cx4 +...), computed numerically, spread the higher-energy vibrational eigenvalues further apart if b > 0, or closer together if b < 0. 

The performance of the (zT) correction is essentially identical to that of the conventional ROHF-CCSD(T) method. Application of both to a series of diatomic molecules in ground and excited states indicates insignificant differences between the two in the prediction of bond lengths, harmonic vibrational frequencies, anharmonic constants, and so on. Unfortunately, the complicated equations associated with the (zT) correction have thus far precluded its large scale implementation and, as a result, further systematic studies involving larger basis sets have not yet been carried out. 

Note the regular variation in ry as v changes indeed, note that r shrinks as the vibrational quantum number decreases. It is the well-known non-harmonic (anharmonic), potential curve that leads to this characteristic variation. Indeed, quantum mechanics shows to a high approximation that ry is a well-deflned average... 

To calculate n E-E, the non-torsional transitional modes have been treated as vibrations as well as rotations [26]. The former approach is invalid when the transitional mode s barrier for rotation is low, while the latter is inappropriate when the transitional mode is a vibration. Harmonic frequencies for the transitional modes may be obtained from a semi-empirical model [23] or by performing an appropriate normal mode analysis as a function of the reaction path for the reaction s potential energy surface [26]. Semiclassical quantization may be used to determine anharmonic energy levels for the transitional modes [27]. 

Just as we corrected the expressions for the rigid rotor to allow for the centrifugal effect and an interaction with the vibration, we also must adjust the expression for the harmonic oscillator to account for the anharmonicity in the oscillation. The potential energy surface for the molecule is not symmetrical (Fig. 25.2). The parabola (dotted figure) represents the potential energy of the harmonic oscillator. The correct potential energy is shown by the full lines the vibration is anharmonic. The vibrational energy levels for such a system can be approximated by a series ... 

The experimental vibrational harmonic frequencies cue for Au Cl (and Au Cl equal to 382.8 cm (and 373.9 cm , respectively) and anharmonicity constants coupled cluster procedure theory QCISD(T) which gives cue values of 369.5 cm (and 360.9 cm ) and weXe values of 1.32 cm (and 1.26 cm ), respectively, for Au Cl (and Au Cl) isotopomers. The estimated " dissociation energy of 3.0 0.7 eV and the value of 2.85 estimated at the QCISD level indicate that the AuCl dissociation energy should be below 3.5 eV and that the experimentally obtained value" of 3.5 0.1 eV is probably overestimated by about 0.5 eV. The vibrational-state dependencies of the molecular properties for Au Cl have been established (equations 77-79) ... 

In our discussion so far, we have assumed that the motions of atoms in a vibrating molecule are harmonic. Although making this assumption made the mathematics easier, it is not a realistic view of the motion of atoms in a real vibrating molecule. Anharmonic motion is the type of motion that really takes place in vibrating molecules. The energy levels of such an anharmonic oscillator are approximately given by... 

Figure 1 illustrates the simplest theoretical treatment for the energy absorption behavior of covalent bonds. Hooke s law curve corresponds to bonds vibrating in simple harmonic motion. In practice, covalent bond vibrations are anharmonic and follow the potential energies depicted in the Morse curve, also illustrated in Figure 1. It is this phenomenon of anharmonicity that allows the overtones or harmonics to occur in the NIR region. Thus, the frequencies of the overtones are slightly less than whole number...

In some cases, it is possible to see the transitions to higher vibrational states (overtones). In the harmonic approximation, these transitions are forbidden. Hence, the appearance of overtones is a sign that the vibration is anharmonic. 

In the discussion which follows, the harmonic oscillator functions will be employed as vibrational wave functions. The results obtained, however, are equally applicable to anharmonic oscillators for the folloAving reasons. The wave funcitions for a molecule whose vibrations are anharmonic can always be expressed as linear combinations of the harmonic o,scillator functions. 

In this chapter, our goal is to present theoretical methods applied to gas phase vibrational spectroscopy. This is reviewed in Sect. 3 where we present harmonic and anharmonic spectra calculations, with special emphasis on dynamical approaches to anharmonic spectroscopy. In particular, we present the many reasons and advantages of dynamical anharmonic theoretical spectroscopy over harmonic/ anharmonic non-dynamical methods in Sect. 3.1. Illustrations taken from our work on dynamical theoretical spectroscopy are then presented in Sects. 4-6 in relation to action spectroscopy experiments. All examples presented here are conducted either in relation to finite temperature IR-MPD and IR-PD experiments, or to cold IR-MPD experiments. Beyond the conformational dynamics provided by the finite temperature trajectories, the chosen examples illustrate how the dynamical spectra manage to capture vibrational anharmonicities of different origins, of different strengths, in various domains of the vibrations from 100 to 4,000 cm and on various molecular systems. Other comprehensive reviews on theoretical anharmonic spectroscopy can be found in [7-9]. 

Here we pause for a more general discussion on vibrational anharmonicities. As already pointed out in the introduction, static quantum chemistry calculations can be performed beyond the vibrational harmonic approximation, though they are more rarely applied in the communities of IR-MPD and IR-PD action spectroscopy. There are exceptions however see works by the group of A.B. McCoy in collaboration with IR-PD experiments from the group of M.A. Johnson at Yale University... 

The calculated barrier to the transition from a pyramidal to a planar configuration and in the opposite direction is very low and varies depending on the calculation procedure. In particular, according to Vetere et al. (2000), the DFT/PBEO procedure yields barrier values of 0.71, 0.33, 0.08, and 0.04 kj/mol for lanthanum trihalides from trifluoride to triiodide, respectively. Of greater importance is another conclusion made by Kovacs et al. (1997) and substantiated by theoretical calculations in more recent works. The low out-of-plane bending vibrational frequencies V2 obtained in calculations are caused by the use of the harmonic approximation. In reality, these vibrations are anharmonic, and the V2 values are much higher. Unfortunately, taking into account anharmonic effects is currently impracticable for lanthanide trichlorides. 

Up to now we have been discussing only harmonic vibrations. Mechanical anharmonicity results if the restoring force is not linearly proportional to the nuclear displacement coordinate. Electrical anharmonicity results if the change in dipole moment is not linearly proportional to the nuclear displacement coordinate. If a vibration is mechanically harmonic, the classical picture of a plot of nuclear displacement versus time is a sine or cosine wave (see Fig. 1.8). If mechanical anharmonicity is present this plot will be periodic but not a simple sine or cosine wave. One result of mechanical anharmonicity... 

Each of the partition functions is now regarded as a product of independent translational, rotational, and vibrational partition functions—the implication being that vibration-rotation interaction is negligible, and it is then assumed that the rotations are classical and the vibrations harmonic. If the structure of each molecular species is known, the moments of inertia can be calculated, and if necessary— as may well be the case for hydrogen isotopes—a correction can be applied to account for the fact that the rotational partition function has not reached its classical value. If complete vibrational analjrses of all the molecules are also available, the vibrational partition functions can be set up, and an approximate correction for neglect of anharmonicity can also be made. Having done all this, we can calculate the isotope effect. 

R. M. Levy, O. de la Luz Rojas, and R. A. Friesner. Quasi-harmonic method for calculating vibrational spectra from classical simulations on multidimensional anharmonic potential surfaces. J. Phys. Chem., 88 4233-4238, 1984. 

Among the main theoretical methods of investigation of the dynamic properties of macromolecules are molecular dynamics (MD) simulations and harmonic analysis. MD simulation is a technique in which the classical equation of motion for all atoms of a molecule is integrated over a finite period of time. Harmonic analysis is a direct way of analyzing vibrational motions. Harmonicity of the potential function is a basic assumption in the normal mode approximation used in harmonic analysis. This is known to be inadequate in the case of biological macromolecules, such as proteins, because anharmonic effects, which MD has shown to be important in protein motion, are neglected [1, 2, 3]. 

This is a check on the reasonableness of the method chosen. For example, it would not be reasonable to select a method to investigate vibrational motions that are very anharmonic with a calculation that uses a harmonic oscillator approximation. To avoid such mistakes, it is important the researcher understand the method s underlying theory. 

The reason that does not change with isotopic substitution is that it refers to the bond length at the minimum of the potential energy curve (see Figure 1.13), and this curve, whether it refers to the harmonic oscillator approximation (Section 1.3.6) or an anharmonic oscillator (to be discussed in Section 6.1.3.2), does not change with isotopic substitution. Flowever, the vibrational energy levels within the potential energy curve, and therefore tq, are affected by isotopic substitution this is illustrated by the mass-dependence of the vibration frequency demonstrated by Equation (1.68). 

Equations (6.5) and (6.12) contain terms in x to the second and higher powers. If the expressions for the dipole moment /i and the polarizability a were linear in x, then /i and ot would be said to vary harmonically with x. The effect of higher terms is known as anharmonicity and, because this particular kind of anharmonicity is concerned with electrical properties of a molecule, it is referred to as electrical anharmonicity. One effect of it is to cause the vibrational selection mle Au = 1 in infrared and Raman spectroscopy to be modified to Au = 1, 2, 3,. However, since electrical anharmonicity is usually small, the effect is to make only a very small contribution to the intensities of Av = 2, 3,. .. transitions, which are known as vibrational overtones. 

However, unlike electrical anharmonicity, mechanical anharmonicity modifies the vibrational term values and wave functions. The harmonic oscillator term values of Equation (6.3) are modified to a power series in (u + ) ... 

Owing to the effects of mechanical anharmonicity - to which we shall refer in future simply as anharmonicity since we encounter electrical anharmonicity much less frequently -the vibrational wave functions are also modified compared wifh fhose of a harmonic oscillator. Figure 6.6 shows some wave functions and probabilify densify functions (iA A ) for an anharmonic oscillator. The asymmefry in and (iA A ) 5 compared wifh fhe harmonic oscillator wave functions in Figure f.i3, increases fheir magnitude on the shallow side of the potential curve compared with the steep side. 

---