Anharmonicities calculation

It is nice to have a distinctive notation for the curvilinear co-ordinates, which emphasizes their difference from and yet their one-to-one correlation with the Rt co-ordinates. Most authors reporting anharmonic calculations do not in fact make any distinction they denote the curvilinear co-ordinates by the same symbols customarily used to denote the corresponding rectilinear coordinates in harmonic calculations. For many purposes this is satisfactory, particularly since the harmonic force constants are not altered by the change from rectilinear to curvilinear co-ordinates. However, in a general discussion it is important to distinguish the two sets, and so for the remainder of this section we shall follow Hoy et al.12 and write the curvilinear co-ordinates with the symbol Hi. 

The first derivative L tensor elements, are determined as the elements of the L matrix from the preliminary harmonic calculation. The second and third derivative L tensor elements have been determined in one of two ways in the literature of anharmonic calculations. The first method involves setting up an... 

The whole of this section has been concerned with the problem of transforming the potential energy V from a representation in geometrically defined internal co-ordinates H to dimensionless normal co-ordinates q, a transformation achieved in the single equation (14) for a diatomic molecule. It will be clear to the reader that programming this transformation is a considerable part of the task of performing an anharmonic calculation on any polyatomic molecule. 

It is clear that a strong Fermi resonance of the type described above may in principle be analysed directly to give an observed value of the single anharmonic constant rrs, in Section 5 (p. 143) we give examples of such analyses, and their use in anharmonic calculations. In the absence of Fermi resonance, information on the cubic constants like rrK comes mainly from the vibrational dependence of the rotational constants a , which determine linear combinations of the cubic constants 4>m as described earlier. 

Often the original structure determination will have involved some uncorrected vibrational averaging effects it may be an r0 or an r, structure.28 However, once /2, or some approximation to /3, has been determined it is possible to correct r to re and obtain an improved equilibrium structure (in most cases this correction can be made directly from the cubic anharmonic calculation, but in some cases the calculation will allow unobserved a. values to be determined, perhaps for other isotopic species, etc.). Similarly, it is often true that the harmonic field /2 is calculated from the observed fundamentals (the v values) rather than the harmonic vibration wavenumbers (the to values), for want of information on the corrections. However, once /4, or some approximation to/4, has been determined, it may be used to calculate a complete set of x values and hence to calculate all the corrections to obtain the co values. Thus the calculation of re and may be improved from a knowledge of /3 and /4. 

Ab initio anharmonic calculations using perturbation theory... 

Perturbation theory has been applied to anharmonic calculations of spectroscopy from ab initio potentials in a large number of studies [19-25,115-121]. In nearly all cases so far, second-order perturbation theory was employed. The representation of the anharmonic potential generally used in these studies is a polynomial in the normal modes, most often a quartic force field. A code implementing this vibrational method was recently incorporated by V. Barone in gaussian [24]. Calculations were carried out for relatively large molecules, such as pyrrole and furan [25], uracil and thiouracil [118], and azabenzenes [119]. We note that in addition to spectroscopy, the ab initio perturbation theoretic algorithms were also applied to the calculation of thermodynamic properties... 

Brauer B, Gerber RB, Kabelac M, Hobza P, Bakker JM, Riziq AGA, de Vries MS (2005) Vibrational spectroscopy of the G center dot center dot center dot C base pair experiment, harmonic and anharmonic calculations, and the nature of the anharmonic couplings. J Phys Chem A 109 6974-6984... 

Figure 4. Comparison between the experimental (solid bold line) and theoretical v spectra of 1-methyluracil crystal (anharmonic calculation—Dirac 5 functions and thin solid line harmonic calculation—dashed line) for the 1-methyluracil crystal (a) and the deuterated crystal (b). Figure reprinted with permission from Ref. 58. Copyright American Institute of Physics Publishing LLC.

Hollenstein H, Marquardt R, Quack M and Suhm M A 1994 Dipole moment function and equilibrium structure of methane In an analytical, anharmonic nine-dimenslonal potential surface related to experimental rotational constants and transition moments by quantum Monte Carlo calculations J. Chem. Phys. 101 3588-602... 

After transforming to Cartesian coordinates, the position and velocities must be corrected for anharmonicities in the potential surface so that the desired energy is obtained. This procedure can be used, for example, to include the effects of zero-point energy into a classical calculation. 

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. 

It is possible to use computational techniques to gain insight into the vibrational motion of molecules. There are a number of computational methods available that have varying degrees of accuracy. These methods can be powerful tools if the user is aware of their strengths and weaknesses. The user is advised to use ah initio or DFT calculations with an appropriate scale factor if at all possible. Anharmonic corrections should be considered only if very-high-accuracy results are necessary. Semiempirical and molecular mechanics methods should be tried cautiously when the molecular system prevents using the other methods mentioned. 

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. 

Accurate intensity measurements have been made in many cases and calculations of r — r" made, including the effects of anharmonicity and even allowing for breakdown of the Bom-Oppenheimer approximation. 

The MP2 treatment recovers the majority of the correlation effect, and the CCSD(T) results with the cc-pVQZ basis sets are in good agreement with the experimental values. The remaining discrepancies of 9cm , 13cm and lOcm are mainly due to basis set inadequacies, as indicated by the MP2/cc-pV5Z results. The MP2 values are in respectable agreement with the experimental harmonic frequencies, but of course still overestimate the experimental fundamental ones by the anharmonicity. For this reason, calculated MP2 harmonic frequencies are often scaled by 0.97 for comparison with experimental results. ... 

For comparison with experimental frequencies (which necessarily are anharmonic), there is normally little point in improving the theoretical level beyond MP2 with a TZ(2df,2pd) type basis set unless anharmonicity constants are calculated explicitly. Although anharmonicity can be approximately accounted for by scaling the harmonic frequencies by 0.97, the remaining errors in the harmonic force constants at this level are normally smaller than the corresponding errors due to variations in anharmonicity. 

C) The error in AE" /AEq is 0.1 kcal/mol. Corrections from vibrations, rotations and translation are clearly necessary. Explicit calculation of the partition functions for anharmonic vibrations and internal rotations may be considered. However, at this point other factors also become important for the activation energy. These include for example ... 

Vibrational Spectra Many of the papers quoted below deal with the determination of vibrational spectra. The method of choice is B3-LYP density functional theory. In most cases, MP2 vibrational spectra are less accurate. In order to allow for a comparison between computed frequencies within the harmonic approximation and anharmonic experimental fundamentals, calculated frequencies should be scaled by an empirical factor. This procedure accounts for systematic errors and improves the results considerably. The easiest procedure is to scale all frequencies by the same factor, e.g., 0.963 for B3-LYP/6-31G computed frequencies [95JPC3093]. A more sophisticated but still pragmatic approach is the SQM method [83JA7073], in which the underlying force constants (in internal coordinates) are scaled by different scaling factors. 

For diatomic molecules, B0 is the rotational constant to use with equation (10.125), while Be applies to equation (10.124). They are related by Bq = Be 2 - The moment of inertia 70(kg-m2) is related to 50(cm ) through the relationship /0 = h/ 8 x 10 27r22 oc), with h and c expressed in SI units. For polyatomic molecules, /a, /b, and Iq are the moments of inertia to use with Table 10.4 where the rigid rotator approximation is assumed. For diatomic molecules, /0 is used with Table 10.4 to calculate values to which we add the anharmonicity and nonrigid rotator corrections. 

For diatomic molecules, lj0 is the vibrational constant to use with equation (10.125) for calculating anharmonicity and nonrigid rotator corrections, while J)e and tDe-Ve... 

Under most circumstances the equations given in Table 10.4 accurately calculate the thermodynamic properties of the ideal gas. The most serious approximations involve the replacement of the summation with an integral [equations (10.94) and (10.95)] in calculating the partition function for the rigid rotator, and the approximation that the rotational and vibrational partition functions for a gas can be represented by those for a rigid rotator and harmonic oscillator. In general, the errors introduced by these approximations are most serious for the diatomic molecule." Fortunately, it is for the diatomic molecule that corrections are most easily calculated. It is also for these molecules that spectroscopic information is often available to make the corrections for anharmonicity and nonrigid rotator effects. We will summarize the relationships... 

One of the more interesting results of these calculations is the contribution to the heat capacity. Figure 10.10 shows the temperature dependence of this contribution to the heat capacity for CH3-CCU as calculated from Pitzer s tabulation with 7r = 5.25 x 10-47 kg m2 and VQ/R — 1493 K. The heat capacity increases initially, reaches a maximum near the value expected for an anharmonic oscillator, but then decreases asymptotically to the value of / expected for a free rotator as kT increases above Vo. The total entropy calculated for this molecule at 286.53 K is 318.86 J K l-mol l, which compares very favorably with the value of 318.94T 0.6 TK-1-mol 1 calculated from Third Law measurements.7... 

---