Theoretical calculation of vibrational frequencies

The procedure for calculating harmonic vibrational frequencies and force constants by GF matrix method has been described in Sec. 1.12. In this method, both G (kinetic energy) and F (potential energy) matrices are expressed in terms of internal coordinates (R) such as increments of bond distances and bond angles. Then, the kinetic (7) and potential (V) energies are written as  

Since the total electronic enegy ( ) of the H2O molecule is a function of the nine Cartesian coordinates (x i where i = 1, 2,. .., 9), the equilibrium configuration that renders all the nine first derivatives (dEldxj) zero can be determined (energy optimization). Accuracy of the results depends on the level of the approximation method 

The next step is to calculate the second derivatives near the equilibrium positions [(dE /dXf dxj)o4J= 1 9], namely, the force constants in terms of Cartesian coordinates (fij), using analytical methods [116a]. The force constants thus obtained are converted to mass-weighted Cartesian coordinates (Xmj) by 

These/m,j terms are the elements of the force constant (Hessian) matrix. 

Using matrix expression, the relationship between mass-weighted Cartesian (Xm) and Cartesian (X) coordinates is written as follows  

---

Chapter 3 is devoted to dipole dispersion laws for collective excitations on various planar lattices. For several orientationally inequivalent molecules in the unit cell of a two-dimensional lattice, a corresponding number of colective excitation bands arise and hence Davydov-split spectral lines are observed. Constructing the theory for these phenomena, we exemplify it by simple chain-like orientational structures on planar lattices and by the system CO2/NaCl(100). The latter is characterized by Davydov-split asymmetric stretching vibrations and two bending modes. An analytical theoretical analysis of vibrational frequencies and integrated absorptions for six spectral lines observed in the spectrum of this system provides an excellent agreement between calculated and measured data. 

Theoretical Calculations. The vibrational frequencies and infrared and Raman intensities can be calculated for benzene using programs such as Gaussian/Gaussview and HyperChem. These programs also permit animated visualization of the normal modes such as the symmetric CC and CH stretches studied in this experiment. 

This review has summarize the applications of neutron inelastic scattering to the study of pol3uners. The technique has proven useful for measuring and characterizii low-frequency intramolecular and inter-molecular vibrations, particularly for three systems, such as polyethylene and the n-paraffins, for which theoretical calculations of phase-frequency relations are available. More calculations of this type, and extension of them to include the effects of departures of chain conformations from their ideal transplanar or helical configurations, are needed for an optimum application of the method. 

Theoretical calculation of the frequencies of the normal modes of vibration is possible and has been demonstrated for moderately complex molecules such as acetazolamide [34]. This facilitates assignment of the observed spectral bands beyond what may be possible by comparison with correlation tables and manual interpretation. The ability to predict the visual appearance of IR and Raman spectra is, however, more challenging since the inclusion of selection rules to ascertain which bands will be IR or Raman active and the calculation of their intensity must be considered. 

Now that we see that we can combine partition functions for all the quantized energy systems into a total partition function, we can think of other ways to use the quantized energy formulas. There is a curious history for this approach. We can see above that gvib is an important part of the total partition function and yet for many years low-resolution infrared spectra blurred many of the 3N — 6 vibrational modes of molecules typically larger than benzene. Thus the equations for quantum thermodynamics were known before 1940 but could only be applied to cases of small molecules in the gas phase using experimental vibrational frequencies. Since about 1985, quantum chemistry programs have included the calculation of vibrational frequencies with some correction factors that now make it possible to write down the full partition function by including theoretical... 

Theoretical prediction of vibrational frequencies serves as an aid to experiment for the characterization of chemical species. When experimental observation is difficult or impossible, theoretical approaches are able to offer useful information. The calculated vibrational frequencies can be used to verify whether a molecular structure is a genuine minimum, or is a transition state structure. The vibrational frequencies of a molecule are all real for a true minimum, and should have a single imaginary frequency for a saddle point. The characterization of transition state structures can only be achieved accurately by theoretical methods. 

The entropy difference A5tot between the HS and the LS states of an iron(II) SCO complex is the driving force for thermally induced spin transition [97], About one quarter of AStot is due to the multiplicity of the HS state, whereas the remaining three quarters are due to a shift of vibrational frequencies upon SCO. The part that arises from the spin multiplicity can easily be calculated. However, the vibrational contribution AS ib is less readily accessible, either experimentally or theoretically, because the vibrational spectrum of a SCO complex, such as [Fe(phen)2(NCS)2] (with 147 normal modes for the free molecule) is rather complex. Therefore, a reasonably complete assignment of modes can be achieved only by a combination of complementary spectroscopic techniques in conjunction with appropriate calculations. 

Vibrational spectroscopy is of utmost importance in many areas of chemical research and the application of electronic structure methods for the calculation of harmonic frequencies has been of great value for the interpretation of complex experimental spectra. Numerous unusual molecules have been identified by comparison of computed and observed frequencies. Another standard use of harmonic frequencies in first principles computations is the derivation of thermochemical and kinetic data by statistical thermodynamics for which the frequencies are an important ingredient (see, e. g., Hehre et al. 1986). The theoretical evaluation of harmonic vibrational frequencies is efficiently done in modem programs by evaluation of analytic second derivatives of the total energy with respect to cartesian coordinates (see, e. g., Johnson and Frisch, 1994, for the corresponding DFT implementation and Stratman etal., 1997, for further developments). Alternatively, if the second derivatives are not available analytically, they are obtained by numerical differentiation of analytic first derivatives (i. e., by evaluating gradient differences obtained after finite displacements of atomic coordinates). In the past two decades, most of these calculations have been carried... 

Theoretical calculations of R2Ge give bond parameters and vibrational frequencies. For example, Me2Ge has calculated GeC = 202 pm and CGeC = 98° with the Ge-C stretches at 560 cm-1 (Ai) and 497 cm-1 (Bi). As the latter agree reasonably with experimental values from matrix-isolated species (527 and 541 cm-1), the structural values are probably good indications. 

No characteristic IR data were reported for pyrrolizines or dihydropyrrolizines. An almost complete set of vibrational frequencies was deduced by combining an infrared (IR) and a Raman spectrum of pyrrolizinone 2 <2001J(P2)2195>. The experimental values thus obtained were used to scale the theoretical complete set of vibrational frequencies of 2. Using the same scaling constant, the authors proposed a set of calculated vibrational frequencies for dihydropyrrolizinone 3. 

Crystals lack some of the dynamic complexity of solutions, but are still a challenging subject for theoretical modeling. Long-range order and forces in crystals cause their spectrum of vibrational frequencies to appear more like a continuum than a series of discrete modes. Reduced partition function ratios for a continuous vibrational spectrum can be calculated using an integral, rather than the hnite product used in Equation (3) (Kieffer 1982),... 

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

There exists an extensive literature on theoretical calculations of the vibrational damping of an excited molecule on a metal surface. The two fundamental excitations that can be made in the metal are creation of phonons and electron-hole pairs. The damping of a high frequency mode via the creation of phonons is a process with small probability, because from pure energy conservation, it requires about 6-8 phonons to be created almost simultaneously. 

The theoretical value of the frequency of vibration, depending on the curvature of the cmrve at its minimum, is naturally more uncertain. Calculation shows that the curve gives a frequency of vibration of 5300 cm. S about 20% higher than the value 4360 cm. from experiment. As for the moment of inertia, while it is larger than most of the values from specific heat theories, it is in accord with the larger values which have been found by Richardson and Tanaka from analysis of the hydrogen bands. 

Vibrational frequencies measured in IR experiments can be used as a probe of the metal—ligand bond strength and hence for the variation of the electronic structure due to metal—radical interactions. Theoretical estimations of the frequencies are obtained from the molecular Hessian, which can be straightforwardly calculated after a successful geometry optimization. Pure density functionals usually give accurate vibrational frequencies due to an error cancellation resulting from the neglect of... 

From the theoretical point of view the work of Leszczynski and collaborators159, should be mentioned, where the calculated harmonic vibrational frequencies of thioformaldehyde were compared with those of formaldehyde and selenoformaldehyde. The largest discrepancy (120 cm-1) between calculated and experimental vibrational frequencies was found for the aforementioned CH2 in-plane deformation. 

Fig. 5.10. The theoretical calculation of the time trace of transient absorption (TRABS) for a one-mode system. The energy gap is 20 cm-1 and the vibrational mode is 420 cm-1. The dark curve is the reactant TRABS and the light-gray curve is the product TRABS. The probing frequency is set at respective peak positions of the induced absorption spectra of both reactant state and product state. For discussion see text.

Figure 15 Tj (p, T) vs. temperature for the solvent carbon dioxide at the critical density and the theoretically calculated curve. The frequency u> and the hard sphere diameters are the same as those used in the fit of the 33°C data. The theory is scaled to match the data at 33°C and the critical density, 10.6 mol/L. Unlike ethane at the critical density, there is no inverted region, and the vibrational lifetime decreases nearly linearly with temperature. The theory does not quantitatively fit the data, but it does show the correct general behavior. Most importantly, the hydrodynamic/thermodynamic theory shows the existence of the inverted region in ethane and the lack of one in carbon dioxide.

The H local vibrational mode at 3096 cm has been assigned to N-H centers. This assignment is consistent with a recent theoretical investigation by Van de Walle who calculated a vibrational frequency of approximately 3100 cm for the N-H center. It is interesting to note that this frequency deviates by about 10 % from the value observed for the N-H vibration in ammonia molecules. 

Similar sorts of conclusions apply to the frequencies. A systematic study " found that a DZP basis set yields vibrational frequencies within about 9% of experimental (harmonic) values. The discrepancy diminishes to 4% when correlation is included via CISD and to 2% with a coupled cluster treatment. Another set of calculations confirmed the eost-effec-tiveness of the MP2 treatment of vibrational frequencies, indicating better agreement with experiment than MP3 on some oceasions. Certain types of modes can be more sensitive to the level of theoretical treatment than others. For example, out-of-plane bending motions for it-bonded systems can require triple- plus two sets of polarization functions, as well as a set of/-functions in the basis set . 

The theoretical basis for a full valence force field treatment of extended lattices lies with the work of Kleinman and Spitzer published in 1962. ° These authors, who developed their force field to calculate the vibrational frequencies of quartz, felt that the most accurate way to represent the vibrational motion in a quartz crystal was to include the relative motion of oxygen and silicon atoms. The valence force field was the most effective method for treating this localized picture. 

---