Density functional and experimental studies on the FT-IR and FT-Raman spectra and structure of 2,6-diamino purine and 6-methoxy purine

FT-IR and FT-Raman spectra of 2,6-diamino purine (DAP) and 6-methoxy purine (MP) have been recorded in the regions of 4000-400cm(-1) and 3500-100cm(-1), respectively. The spectra were interpreted with the aid of normal coordinate analysis following full structure optimizations and force field calculations based on density functional theory (DFT) using standard B3LYP/6-31G* and B3LYP/6-311+G** method and basis set combinations. Normal coordinate calculations were performed with the DFT force field corrected by a recommended set of scaling factors yielding fairly good agreement between observed and calculated frequencies.     

An evaluation of scaling factors for multiparameter scaling procedures based on DFT force fields

An extended database of scaling factors for calculating fundamental frequencies within multiparameter scaled quantum mechanical (SQM) force field, and effective scaling frequency factor (ESFF) methods, based on various DFT force fields is reported. Twenty-six density functionals have been examined in conjunction with various Pople's and Dunning's basis sets of VDZ and VTZ quality. The calculations were based on a standard training set of 30 molecules proposed by Baker et al., for which 660 vibrational modes were assigned. Six functionals turned out to be particularly well-suited to the calculations oriented toward determination of scaled frequencies. They are B3LYP, B3PW91, B97, B97-1, B97-2, and O3LYP; they are all capable of providing reasonable scaled frequencies even for the small, 6-31G* basis set (rms <15 cm(-1)). Correlations between the quality of the scaled frequencies and the basis set quality as well as the accuracy of the predicted molecular geometry were investigated. The 6-311+G** basis set turned out to be preferable. In addition, correlation between the rms values in the scaled frequencies and the errors in the predicted geometric parameters was found. Both 11- and 9-parameter scaling frames are compared. It is shown that 9-parameter scaling is preferable in the middle range of the vibrational spectrum (500-2500 cm(-1)), provided it was based on high-quality force fields. Finally, statistical uncertainties of the calculated scaling factors are reported.     

Calculation of vibrational spectra of linear tetrapyrroles. 4. Methine bridge C-H out-of-plane modes

Quantum chemical force fields obtained by density functional theory (DFT) calculations systematically overestimate the frequencies of normal modes including ethylenic C-H out-of-plane (HOOP) coordinates. Compensation of this deviation requires a specific scaling factor for this type of coordinate that is distinctly lower than those applicable to out-of-plane coordinates in general. Such a specific scaling factor (0.900) has been optimized for the DFT(B3LYP) level of theory on the basis of vibrational analyses of training molecules including the HOOP coordinate. Thus, the root-mean-square deviation for the calculated frequencies of these modes is reduced from 16 to 8 cm(-1). Although Raman intensities are yet not reproduced in a satisfactory manner, implementation of the HOOP scaling factor into the set of global scaling factors determined previously (Magdo et al. J. Phys. Chem. A 1999, 103, 289-303) allows for a substantially improved reproduction of the experimental (resonance) Raman spectra of test molecules including linear methine-bridged tetrapyrroles. A very good agreement between calculated and experimental spectra is noted for the phycocyanobilin dimethylester dimer as well as for the protein-bound phycocyanobilin in the antenna pigment alpha-CPC. However, for the phycocyanobilin chromophore in the P(r) state of the plant photoreceptor phytochrome phyA, considerable deviations remain in the spectral range between 800 and 500 cm(-1), which are attributed to the effect of specific protein-chromophore interactions. The influence of the protein environment is not considered in the present calculations that refer to the molecule in vacuo.     

Molecular structure and vibrational investigation of benzenesulfonic acid methyl ester using DFT (LSDA, B3LYP, B3PW91 and MPW1PW91) theory calculations

The FT-Raman and FT-IR spectra for benzenesulfonic acid methyl ester (BSAME) have been recorded in the region 4000-100 cm(-1) and compared with the harmonic vibrational frequencies calculated using DFT (LSDA, B3LYP, B3PW91 and MPW1PW91) method by employing 6-311G (d, p) basis set with appropriate scale factors. IR intensities and Raman activities are also calculated by DFT (LSDA, B3LYP, B3PW91 and MPW1PW91) methods. Optimized geometries of the molecule have been interpreted and compared with the reported experimental values for sulfonic acid and some substituted sulfonic acids. The experimental geometrical parameters show satisfactory agreement with the theoretical prediction from DFT. The scaled vibrational frequencies at LSDA/B3LYP/6-311G (d, p) seem to coincide with the experimentally observed values with acceptable deviations. The theoretical spectrograms (IR and Raman) have been constructed and compared with the experimental FT-IR and FT-Raman spectra. Some of the vibrational frequencies of the sulfonic acid are effected upon profusely with the methyl substitution in comparison to benzene sulfonamide and these differences are interpreted.     

Molecular structure and vibrational analysis of 3-Ethylpyridine using ab initio HF and density functional theory (B3LYP) calculations

The FT-Raman and FT-IR spectra for 3-Ethylpyridine (3-EP) have been recorded in the region 4000-100 cm(-1) and compared with the harmonic vibrational frequencies calculated using HF/DFT (B3LYP) method by employing 6-31G(d,p) and 6-311++G(d,p) basis set with appropriate scale factors. IR intensities and Raman activities are also calculated by HF and DFT (B3LYP) methods. Optimized geometries of the molecule have been interpreted and compared with the reported experimental values of some substituted benzene. The experimental geometrical parameters show satisfactory agreement with the theoretical prediction from HF and DFT. The scaled vibrational frequencies at B3LYP/6-311++G(d,p) seem to coincide with the experimentally observed values with acceptable deviations. The theoretical spectrograms (IR and Raman) have been constructed and compared with the experimental FT-IR and FT-Raman spectra. Some of the vibrational frequencies of the pyridine are effected upon profusely with the C2H5 substitutions in comparison to pyridine and these differences are interpreted.     

Vibrational spectroscopic (FTIR and FTRaman) investigation using ab initio (HF) and DFT (LSDA and B3LYP) analysis on the structure of Toluic acid

The FTRaman and FTIR spectra for Toluic acid (TA) have been recorded in the region 4000-100 cm(-1) and compared with the harmonic vibrational frequencies calculated using HF/DFT (LSDA and B3LYP) method BY employing 6-311G (d, p) basis set with appropriate scale factors. IR intensities and Raman activities are also calculated by HF and DFT (LSDA/B3LYP) methods. Optimized geometries of the molecule have been interpreted and compared with the reported experimental values for benzoic acid and some substituted benzoic acids. The experimental geometrical parameters show satisfactory agreement with the theoretical prediction from HF and DFT. The scaled vibrational frequencies at B3LYP/6-311G (d, p) seem to coincide with the experimentally observed values with acceptable deviations. The theoretical spectrograms (IR and Raman) have been constructed and compared with the experimental FT-IR and FT-Raman spectra. Some of the vibrational frequencies of the TA are effected upon profusely with the methyl substitutions in comparison to benzoic acid and these differences are interpreted.     

Density functional theory study of vibrational spectra, and assignment of fundamental vibrational modes of 1-bromo 4-fluoronaphthalene

The FTIR and FT-Raman spectra of 1-bromo 4-fluoronaphthalene have been recorded in the regions 4000-100cm(-1) and 3500-100cm(-1), respectively. The spectra were interpreted with the aid of normal coordinate analysis based on DFT (density functional theory) using standard B3LYP/6-311+G** basis set combination for the most optimized geometry. Normal coordinate calculations performed with the DFT force field and subsequently corrected by a recommended set of scale factors, yielded fairly good agreement between observed and calculated frequencies. On the basis of the comparison between calculated and experimental results, assignments of fundamental modes were examined.     

Vibrational spectra and potential energy distributions for 4,5-dichloro-3-hydroxypyridazine by density functional theory and normal coordinate calculations

The solid phase FT-IR and FT-Raman spectra of 4,5-dichloro-3-hydroxypyridazine have been recorded in the regions 4000-400 cm(-1) and 3500-100 cm(-1), respectively. The spectra were interpreted with the aid of normal coordinate analysis following a full structure optimization and force field calculations based on the density functional theory (DFT) using the standard B3LYP/6-31G* and B3LYP/6-311+G** method and basis set combinations. The DFT force field transformed to natural internal coordinates was corrected by a well-established set of scale factors that were found to be transferable to the title compound. The IR and Raman spectra were predicted theoretically and compared with the experimental spectra.     

Simulation of IR and Raman spectral based on scaled DFT force fields: a case study of 2-amino 4-hydroxy 6-trifluoromethylpyrimidine, with emphasis on band assignment

This work deals with the vibrational spectroscopy of 2-amino 4-hydroxy 6-triflouromethylpyrimidine (AHFMP) by means of quantum chemical calculations. The mid and far FTIR and FT-Raman spectra were measured in the condensed state. The fundamental vibrational frequencies and intensity of vibrational bands were evaluated using density functional theory (DFT) with the standard B3LYP/6-31G* and B3LYP/6-311+G** method and basic set combinations. Normal co-ordinate calculations were performed with the DFT force field corrected by a recommended set of scaling factors yielding fairly good agreement between observed and calculated frequencies. Simulation of infrared and Raman spectra utilizing the results of these calculations led to excellent overall agreement with the observed spectral patterns. The SQM approach applying selective scaling of the DFT force field was shown to be superior to the uniform scaling method in its ability to allow for making modifications in the band assignment, resulting in more accurate simulation of IR and Raman Spectra.     

Molecular structure and vibrational spectra of o-chlorotoluene, m-chlorotoluene, and p-chlorotoluene by ab initio HF and DFT calculations

In this work, experimental and theoretical study on the molecular structure and the vibrational spectra of o-chlorotoluene (OCT), m-chlorotoluene (MCT) and p-chlorotoluene (PCT) are presented. The vibrational frequencies of these compounds were obtained theoretically by ab initio HF and DFT/B3LYP calculations employing the standard 6-311++G(d,p) basis set for optimized geometries and were compared with Fourier transform infrared (FTIR) in the region of 400-4000 cm(-1) and with Raman spectra in the region of 100-4000 cm(-1). Complete vibrational assignment, analysis and correlation of the fundamental modes for these compounds have been carried out. The vibrational harmonic frequencies were scaled using scale factors, yielding a good agreement between the experimentally recorded and the theoretically calculated values.     

FT-IR, NIR-FT-Raman and gas phase infrared spectra of 3-aminoacetophenone by density functional theory and ab initio Hartree-Fock calculations

The gas phase infrared spectrum of 3-aminoacetophenone (3AAP) was measured in the range 5000-500cm(-1) and with a resolution of 0.5cm(-1). The Fourier transform Raman (FT-Raman) and Fourier transform infrared (FT-IR) spectra of 3AAP were recorded in the solid phase. Geometry optimizations were done without any constraint and several thermodynamic parameters were calculated for the minimum energy conformer at ab initio and density functional theory (DFT) levels invoking 6-311G(2df 2p) basis set and the results are compared with the experimental values. Harmonic-vibrational wavenumber was also calculated for the minimum energy conformer at ab initio and DFT levels using 6-31G(d,p) basis set and the results are compared with related molecules. With the help of specific scaling procedures, the observed vibrational wavenumbers in gas phase, FT-IR and FT-Raman spectra were analyzed and assigned to different normal modes of the molecule. Most of the modes have wavenumbers in the expected range, the error obtained was in general very low. The appropriate theoretical spectrogram for the FT-IR spectra of the title molecule is also constructed.     

Calculated gas-phase infrared spectra of 2,2,5,5,8,8-hexamethylhexahydroimidazo [1,2-a]-pyrazine-3,6-dione: relative performance of Hartree-Fock and hybrid density functional theory methods.

Geometric parameters and IR spectra of 2,2,5,5,8,8-hexamethylhexahydroimidazo[1,2-a]pyrazine-3,6-dione were computed by the HF, B3LYP, B3P86 and B3PW91 methods with the 6-31G(d) basis set. All the computation methods produce a boat conformation of the six-membered ring with the orientation opposite to the one experimentally observed, that is obviously due to incapability of such calculations to take into account mutual influence of the molecules in a crystal package. Quality of predicting the interatomic distances observed experimentally decreases in the series HF, B3LYP, B3PW91 and B3P86; for bond angles, in the series B3LYP, P3PW91, B3P86 and HF. At the same time, the IR spectra computed are in a good agreement with the observed FTIR spectrum. The correlation between the calculated and experimental vibration frequencies is characterized by the coefficients higher than 0.9999 for all three DFT methods; for HF it is slightly lower than 0.9995. The calculated absolute band intensities satisfactory match the observed relative intensities as well. Optimal uniform scaling factors calculated are 0.8973, 0.9593, 0.9540 and 0.9552 for HF, B3LYP, B3P86 and B3PW91, respectively.     

Ab initio and DFT studies of the vibrational spectra of hydrogen-bonded PhOH...(H2O)4 complexes

The vibrational characteristics (vibrational frequencies and infrared intensities) for the hydrogen-bonded complex of phenol with four water molecules PhOH...(H2O)4 (structure 4A) have been predicted using ab initio and DFT (B3LYP) calculations with 6-31G(d,p) basis set. The changes in the vibrational characteristics from free monomers to a complex have been calculated. The ab initio and B3LYP calculations show that the observed four intense bands at 3299, 3341, 3386 and 3430 cm(-1) can be assigned to the hydrogen-bonded OH stretching vibrations in the complex PhOH...(H2O)4 (4A). The complexation leads to very large red shifts of these vibrations and very strong increase in their IR intensity. The predicted red shifts for these vibrations with B3LYP/6-31G(d,p) calculations are in very good agreement with the experimentally observed. It was established that the phenolic OH stretching vibration is the most sensitive to the hydrogen bonding. The predicted red-shift with the B3LYP/6-31G(d,p) calculations for the most stable ring structure 4A (-590 cm(-1)) is in better agreement with the experimentally observed than the red-shift, predicted with SCF/6-31G(d,p) calculations. The magnitude of the wavenumber shift is indicative of relatively strong OH...H hydrogen-bonded interaction. The complexation between phenol and four water molecules leads to strong increase of the IR intensity of the phenolic OH stretching vibration (up to 38 times).     

Scaling dynamical correlation energy from density functional theory correlation functionals

The scaling of dynamical correlation energy in molecules obtained by the correlation functionals of density functional theory (DFT) is examined. The approach taken is very similar to the scaled external correlation method of Brown and Truhlar but is based on the observation that DFT correlation functionals, especially the LYP, appear to represent the dynamical portion of the correlation energy in molecules. We examine whether higher accuracy in atomization energies can be gained by scaling without significant deterioration of the structural and spectroscopic properties of the molecules using four DFT functionals (BLYP, OLYP, B3LYP, and O3LYP) on 19 molecules including the six molecule AE6 database, the latter being representative of a much larger, 109 molecule training set. We show that, with molecule specific scale factors, nearly perfect agreement with experiment can be achieved in atomization energies without increasing the average errors in other molecular properties relative to the DFT calculation. We further show that it is possible to find optimal scale factors which reduce the mean unsigned error per bond to levels comparable to those of some multilevel multicoefficient methods.     

Harmonic frequency scaling factors for Hartree-Fock, S-VWN, B-LYP, B3-LYP, B3-PW91 and MP2 with the Sadlej pVTZ electric property basis set

Scaling factors for obtaining fundamental vibrational frequencies from harmonic frequencies calculated at six of the most commonly used levels of theory have been determined from regression analysis for the polarized-valence triple-zeta (pVTZ) Sadlej electric property basis set. The Sadlej harmonic frequency scaling factors for first- and second-row molecules were derived from a comparison of a total of 900 individual vibrations for 111 molecules with available experimental frequencies. Overall, the best performers were the hybrid density functional theory (DFT) methods, Becke's three-parameter exchange functional with the Lee–Yang–Parr fit for the correlation functional (B3-LYP) and Becke's three-parameter exchange functional with Perdew and Wang's gradient-corrected correlation functional (B3-PW91). The uniform scaling factors for use with the Sadlej pVTZ basis set are 0.9066, 0.9946, 1.0047, 0.9726, 0.9674 and 0.9649 for Hartree–Fock, the Slater–Dirac exchange functional with the Vosko–Wilk–Nusair fit for the correlation functional (S-VWN), Becke's gradient-corrected exchange functional with the Lee–Yang–Parr fit for the correlation functional (B-LYP), B3-LYP, B3-PW91 and second-order M?ller–Plesset theory with frozen core (MP2(fc)), respectively. In addition to uniform frequency scaling factors, dual scaling factors were determined to improve the agreement between computed and observed frequencies. The scaling factors for the wavenumber regions below 1800 cm⁻¹ and above 1800 cm⁻¹ are 0.8981 and 0.9097, 1.0216 and 0.9857, 1.0352 and 0.9948, 0.9927 and 0.9659, 0.9873 and 0.9607, 0.9844 and 0.9584 for Hartree–Fock, S-VWN, B-LYP, B3-LYP, B3-PW91 and MP2(fc), respectively. Hybrid DFT methods along with the Sadlej pVTZ basis set provides reliable theoretical vibrational spectra in a cost-effective manner. Received: 22 May 2000 / Accepted: 30 August 2000 / Published online: 28 February 2001     

Calculated gas-phase infrared spectra of imidazo[1,2-a]pyrazinediones derived from alanine

IR spectra of bicyclic and tricyclic amidine derivatives of alanine (2,5,8-trimethylhexahydroimidazo[1,2-a]pyrazine-3,6-dione, 2,5,8-trimethyltetrahydroimidazo[1,2-a]pyrazine-3,6-dione and 2,5,7,10-tetramethylhexahydroimidazo[1,2-a]imidazo[1,2-d]pyrazine-3,8-dione) were computed by the HF, B3LYP, B3P86 and B3PW91 methods in conjunction with the 6-31G(d) basis set. The IR spectra calculated are in a good agreement with the observed FTIR spectra. The correlation between the calculated and experimental vibration frequencies is characterized by the coefficients of 0.9997 for all three DFT methods; for HF it is about 0.9992. The calculated absolute band intensities satisfactory match the observed relative intensities. Optimal uniform scaling factors calculated for this series of compounds are 0.8967, 0.9598, 0.9544 and 0.9555 for HF, B3LYP, B3P86 and B3PW91, respectively. Taking into account small variations of the scaling factors for the derivatives of different amino acids, for future IR spectral predictions for unknown compounds of this class, one can recommend scaling factors of 0.897, 0.959, 0.954 and 0.955 for HF, B3LYP, B3P86 and B3PW91, respectively.     

FT-IR and FT-Raman spectroscopic investigation, computed vibrational frequency analysis and IR intensity and Raman activity peak resemblance analysis on 2-nitroanisole using HF and DFT (B3LYP and B3PW91) calculations

Fourier-transform Raman and infrared spectra of 2-nitroanisole are recorded (4000-100 cm(-1)) and interpreted by comparison with respective theoretical spectra calculated using HF and DFT method. The geometrical parameters with C(S) symmetry, harmonic vibrational frequencies, infrared and Raman scattering intensities are determined using HF/6-311++G (d, p), B3LYP/6-311+G (d, p), B3LYP/6-311++G (d, p) and B3PW91/6-311++G (d, p) level of theories. A detailed vibrational spectral analysis has been carried out and assignments of the observed fundamental bands have been proposed on the basis of peak positions and relative intensities. The results of the calculations have been used to simulate IR and Raman spectra for the molecule that showed good agreement with the observed spectra. The SQM method, which implies multiple scaling of the DFT force fields has been shown superior to the uniform scaling approach. The vibrational frequencies and the infrared intensities of the C-H modes involved in back-donation and conjugation are also investigated.     

A theoretical study of P4O10: vibrational analysis, infrared and Raman spectra

The normal mode frequencies and the corresponding vibrational assignments of tetraphosphorus decaoxide (P4O10) in tetrahedral (Td) symmetry are examined theoretically and experimentally. The Gaussian 98 set of quantum chemistry codes at the HF/6-311G*, MP2/6-311G*, and DFT/B3LYP/6-311G* levels of theory are used. By comparison to experimental normal mode frequencies deduced by Gilliam et al. [J. Phys. Chem. B 107 (2003) 2892], Chapman [Spectrochim. Acta A, 24 (1968) 1687], Beattie et al. [J. Chem. Soc. A (1970) 449], Konings et al. [J. Mol. Spectrosc. 152 (1992) 29] and the present work, correction factors for predominant vibrational motions are determined and compared. Normal modes were decomposed into five non-redundant motions (P-O stretch, P=O stretch, P-O-P bend, P-O-P wag, and P=O wag). Standard deviations found for the HF, MP2, and DFT corrected frequencies compared to experiment are 9, 5, and 4 cm(-1), respectively. Electron distribution for selected molecular orbitals is considered.     

Scaling factors for carbon NMR chemical shifts obtained from DFT B3LYP calculations

A linear scaling of the calculated chemical shifts is used in order to improve the accuracy of the DFT predicted ¹³C NMR chemical shifts. The widely applied method of GIAO B3LYP/6-311+G(2d,p) using the B3LYP/6-31G(d) optimized geometries is chosen, which allows cost-effective calculations of the ¹³C chemical shifts in the molecular systems with 100 and more atoms. A set of 27 ¹³C NMR chemical shifts determined experimentally for 22 simple molecules with various functional groups is used in order to determine scaling factors for reproducing experimentally measured values of ¹³C chemical shifts. The results show that the use of a simple relationship (δ_scalc = 0.95 δ_calc + 0.30, where δ_calc and δ_scalc are the calculated and the linearly scaled values of the ¹³C chemical shifts, respectively) allows to achieve a three-fold improvement in mean absolute deviations for 27 chemical shifts considered. To test the universal applicability of the scaling factors derived, we have used complex organic molecules such as taxol and a steroid to demonstrate the significantly improved accuracy of the DFT predicted chemical shifts. This approach also outperforms the recently recommended usage of the Hartree-Fock optimized geometries for the GIAO B3LYP/6-311+G(2d,p) calculations of the ¹³C chemical shifts.     

Anharmonic vibrational frequency shifts upon interaction of phenol+ with the open shell ligand O2. The performance of DFT methods versus MP2

Anharmonic vibrational frequency shifts of the phenol(+) O-H stretching mode upon complex formation with the open-shell ligand O(2) were computed at several DFT and MP2 levels of theory, with various basis sets, up to 6-311++G(2df,2pd). It was found that all DFT levels of theory significantly outperform the MP2 method with this respect. The best agreement with the experimental frequency shift for the hydrogen-bonded minimum on the potential energy surfaces was obtained with the HCTH/407 functional (-93.7 cm(-1) theoretical vs -86 cm(-1) experimental), which is a significant improvement over other, more standard DFT functionals (such as, e.g., B3LYP, PBE1PBE), which predict too large downshifts (-139.9 and -147.7 cm(-1), respectively). Good agreement with the experiment was also obtained with the mPW1B95 functional proposed by Truhlar et al. (-109.2 cm(-1)). We have attributed this trend due to the corrected long-range behavior of the HCTH/407 and mPW1B95 functionals, despite the fact that they have been designed primarily for other purposes. MP2 method, even with the largest basis set used, manages to reproduce only less than 50% of the experimentally detected frequency downshift for the hydrogen-bonded dimer. This was attributed to the much more significant spin contamination of the reference HF wave function (compared to DFT Kohn-Sham wave functions), which was found to be strongly dependent on the O-H stretching vibrational coordinate. All DFT levels of theory outperform MP2 in the case of computed anharmonic OH stretching frequency shifts upon ionization of the neutral phenol molecule as well. Besides the hydrogen-bonded minimum, DFT levels of theory also predict existence of two other minima, corresponding to stacked arrangement of the phenol(+) and O(2) subunits. mPW1B95 and PBE1PBE functionals predict a very slight blue shift of the phenol(+) O-H stretching mode in the case of stacked dimer with the nearly perpendicular orientation of oxygen molecule with respect to the phenolic ring, which is entirely of electrostatic origin, in agreement with the experimental observations of an additional band in the IR photodissociation spectra of phenol(+)-O(2) dimer [Patzer, A.; Knorke, H.; Langer, J.; Dopfer, O. Chem. Phys. Lett. 2008, 457, 298]. The structural features of the minima on the studied PESs were discussed in details as well, on the basis of NBO and AIM analyses.