Anharmonicity scaling factors

Our quantum-chemical simulation of benzene oxidation reaction based on pseudospinel iron center (see Fig. 20.36, bottom) reveals the same structure. The characteristic feature of such intermediate is the presence of C(sp )-H bond. The presence of the C(5/7 )-H bond intermediate was confirmed by in-situ IR experiment of Panov et al. [84]. The IR band at 2874 cm appeared immediately after benzene was fed to the FeO catalyst. At the same time no phenol signals were detected. Heating of the sample resulted in complete disappearance of this band. According to our quantem-chemical simulation only the a-complex structure has the characteristic of this IR band. For benzene oxide, which also has two C(sp )-H bonds, this band is not present, since all of the vibrational frequencies are within narrow range of 3182-3218 cm . In the case of the benzene o-complex the calculated IR frequency for the C(sp )-H vibration is 2930 cm , while the other C-H vibrations are within 3178-3215 cm . Applying anharmonic scaling factor/= 0.96 one may obtain quite reasonable agreement 2813 em and 3050-3086 cm (theory estimation) versus 3037-3090 cm and 2874 em (experimental data). 

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. 

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. 

In order to obtain better agreement between theory and experiment, computed frequencies are usually scaled. Scale factors can be obtained through multiparameter fitting towards experimental frequencies. In addition to limitations on the level of calculation, the discrepancy between computed and experimental frequencies is also due to the fact that experimental frequencies include anharmonicity effects, while theoretical frequencies are computed within the harmonic approximation. These anharmonicity effects are implicitly considered through the scaling procedure. 

In order to assign more IR signals of 4a, ab initio calculations on Hbdmpza (3b) and 4a were performed. It is well known for the chosen HF/6-31G basis set that calculated harmonical vibrational frequencies are typically overestimated compared to experimental data. These errors arise from the neglecting anharmonicity effects, incomplete incorporation of electron correlation and the use of finite basis sets in the theoretical treatment (89). In order to achieve a correlation with observed spectra a scaling factor (approximately 0.84-0.90) has to be applied (90). The calculations were calibrated on the asymmetric carboxylate Vasym at 1653 cm. We were especially interested in... 

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. 

H-bonded systems may require additional diffuse or polarization functions. For example, the 6-311++G(d,p) basis set had been found to be suitable for H-bonded systems [78-81], It may be necessary to include Basis Set Superposition Errors (BSSE) [82] and Zero-Point-Energy (ZPE) corrections in evaluating the relative stabilities. Such corrections are often of the same magnitude as the energy differences among the dominant conformers. Moreover, the relative conformer energies may also differ noticeably with the basis sets used. All these factors will affect the Boltzmann factors predicted for different conformers and therefore the appearance of the population weighted VA and VCD spectra. Thus, an appropriate selection of DFT functionals and basis sets is very important for VCD simulations. A scale factor of 0.97-0.98 is usually applied to the calculated harmonic frequencies to account for the fact that the observed frequencies arise from an anharmonic force field instead of a harmonic one. A Lorentzian line shape is typically used in simulations of VA and VCD spectra. The full-width at half maximum (FWHM) used in the spectral simulation is usually based on the experimental VA line widths. 

NOTE Empirical scaling factors have been developed for each model chemistry to help correct theoretical frequencies for anharmonic effects [Scott, A. P., and L. Radom,/. Phys. Chem., 100(1996) 16502]. 

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

We will now comment briefly on the force fields that have been computed for benzene. The first one to incorporate both harmonic and anharmonic force constants (some cubic force constants) was reported by Pulay et al. (146) in 1981. This force field was determined with the 4-21P basis set at the SCF level. The computed harmonic force constants, as expected, were larger than the empirical values, so Pulay et al. developed a scaled quadratic force field. The scale factors were adjusted so that the computed frequencies agreed well with the experimental fundamental frequencies. 

Although the harmonic approximation is satisfactory for small displacements from the equifibrium position, ab initio harmonic force constants and vibrational frequencies are known to be typically overestimated as compared with those experimentally found [86]. Sources of this disagreement are the omission or incomplete incorporation of electron correlation, basis set deficiencies, and the neglect of anharmonicity effects. However, as the overestimation is fairly uniform, the appHcation of appropriate scahng procedures becomes feasible. Due to its simplicity, global scafing (using one uniform scale factor determined by a least-squares fit of the calculated to the experimental vibrational frequencies) has widely been used at different levels of theory [87]. However, for most spectro-... 

Vibrational Frequencies. As noted in Section 15.13, theoretically calculated vibrational frequencies are often multiplied by a scale factor to improve agreement with experiment. For a set of 122 molecules and 1066 vibrational frequencies, scaled theoretical harmonic vibrational frequencies showed the following rms deviations from experimental fundamental (anharmonic) vibrational frequencies also given are the optimum scale factors and the percentages of scaled frequencies with less than 6% error and more than 20% error [A. P. Scott and L. Radom,/ Phys. Chem., 100,16502 (19%)] ... 

Xu are the diagonal anharmonicity constants and Go is the polyatomic counterpart of the small Too Dunham constant [82] in diatomics. Consequently [50, 84, 90], the optimal scaling factor for ZPVEs is almost exactly midway between a 2(co) suitable for harmonic frequencies (as an approximate correction for systematic bias in the calculated frequencies) and a 2(v) suitable for fundamental frequencies (which additionally seeks to approximately corrects for anharmonicity). In fact, Alecu et al. [86] found for a large variety of basis sets and ab initio and DFT methods that 2((o)/2(ZPVE) = 1.014 0.002, which is almost exactly the ratio of 1.0143 found by Perdew and coworkers [87] between harmonic frequencies and ZPVEs derived from experimental anharmonic force fields. Note that the small uncertainty of 0.002 on a ZPVE of 140 kcal/mol still would translate to about 0.3 kcaFmol, and even that is probably optimistic for the uncertainty in an individual... 

---