Accuracy and Applicability of the New Exchange Correlation Functionals for Reproduction of the Infrared Spectra of Butyl Acrylate and Butyl Methacrylate Molecules

The butyl acrylate and butyl methacrylate were optimized by seven functionals. All the structures found are local minima and belong to the Cs symmetry. The calculated frequencies are scaled and ranked according to their square errors. The scaling factors of the B972 and B98 functionals fail to reproduce the infrared spectra. The calculated and scaled frequencies with G96LYP, OLYP, and HCTH functionals give acceptable correlations with the experimental spectra. The scaling factors for O3LYP/6-31G(f,p) and O3LYP/6-311+G(df,p) levels of theory reproduce very well the infrared spectrum of butyl acrylate, and the scaled frequencies at VSXC functional with Pople’s double zeta basis sets show the best accuracy in the case of butyl methacrylate.


Introduction
The density functional theory is one of the most efficient and promising methods of quantum physics and chemistry.It is a theory of electronic structure formulated in terms of electronic density as the basic unknown function instead of the electron wave function.According to the fundamental theorem of Hohenberg and Kohn, the electron density carries all the information that one might need to determine any propriety of the electron system.This methodology is the most used by scientists for predicting the molecular structures and reproducing the experimental data.It is well established from the large amount of the published paper on the calculated molecular properties that the DFT is more accurate and less time consuming than the ab initio methods.Unfortunately, the exchange correlation functional is known as an approximate form, and there is no systematic way to improve its accuracy as in the conventional ab initio methodology.In the last decade, many approximate exchange-correlation functionals have been introduced to the scientific community [1,2].To improve the accuracy of these functionals, one should apply them to medium and large systems.A considerable amount of reports on the vibrational frequency calculations are available in the literature using the earlier exchange correlation functions, namely, BLYP, BP86, B3LYP, and B3PW91 [3][4][5][6][7][8].Nonetheless, they fail to correctly predict the experimental vibrational spectra.This is due to the deficiency of the electron-electron interaction considered and the neglect of anharmonicity when calculating the vibrational frequencies, and the experimental data are almost adjusted with the simulated one using the scaling factors [9][10][11].
To the best of our knowledge, no one has tested the accuracy of this set of functionals on medium-sized molecules.
Organic Chemistry International However, there have been some works which made an assessment on the prediction of the infrared spectra.Jaramillo and Scuseria report that the performance of VSXC is comparable to the experimental and B3LYP data [19].In the vibrational analysis of thiosaccharin and thiosaccharinate anion, Jovanovski et al. [20] found that HCTH and OLYP gave comparable results as the B3LYP and they have been used all together for the assignment of bands.Raissi et al. [21] report the vibrational assignment of trifluoroacetyl acetone using B3LYP, B3PW91, BLYP, BPW91, and G96LYP functionals.They did not see the need for use of scaling factors.The G96LYP gave the smallest wavenumbers within this set of functionals.

Experimental
Butyl acrylate (BuA) and butyl methacrylate (BuMA) (purity greater than 99%) were purchased from Sigma-Aldrich, and they were used as they are without further purifications.The infrared spectra of liquid Butyl acrylate and butyl methacrylate were recorded within 1800-700 cm −1 range with a Perkin Elmer FTIR System-2000 model.

Computational
Molecular geometries were determined by gradient optimizations using the Kohn-Sham density functional theory [22].Seven functionals with different levels of complexity were investigated: B972, B98, G96LYP, HCTH, OLYP, O3LYP, and VSXC.The OLYP, HCTH, and G96LYP functionals are the generalized gradient approximation (GGA) functionals.The B972 and O3LYP are the hybrid functionals.The VSXC functional is the meta-GGA functionals, which depend on the kinetic energy density term.The B98 functional is hybrid meta-GGA functional incorporating 21.98% HF exchange.The O3LYP functionals are derived from Lee-yang-Parr correlation functional, the Handy's OPTX exchange functional which reparameterized the Becke88 exchange functional with a flexible formula, and the HF exact exchange [1].For each one of these functionals, the following four basis sets were applied: two of Pople's split valence basis sets 6-31G(d,p) and 6-311+(df,p) [23][24][25] and two flexible polarized (p) and correlation consistent (cc) valence basis set [26] (cc-p VDZ, AUG-cc-p VDZ).AUG denotes augmentation by diffuse functions.This selection represents a good deal of commonly used basis sets.The butyl acrylate and butyl methacrylate molecules were first optimized at a given methodology, and then harmonic vibrational normal modes were calculated.The obtained frequencies were scaled using the scaling factors proposed by Tantirungrotechai et al. [11].The scaled frequencies were compared with the experimental ones in the finger print region 1800-700 cm −1 .The scaling factors  for vibrational frequencies were evaluated by minimizing the square error between the calculated and the experimental values: where   and ] are the th calculated harmonic and the th experimental fundamental frequencies in cm −1 , respectively.The used scaling factors are represented in Table 1.
All the calculations were performed by using the Gaussian 03 rev.B01 program running on PCs running the windows operating system [18].
The Figure 1 shows clearly that the scaling factors give more accurate spectra than the unscaled ones.The difference between the scaled and unscaled frequencies rises up to 35 cm −1 for high frequencies.

Results and Discussion
All the optimizations lead to the Cs symmetry of the butyl acrylate and butyl methacrylate molecules.For molecular structure and atom numbering, see Figure 2. It is too knotty to compare all the observed and calculated frequencies for each spectrum.We based this confrontation on the nine wellknown and prominent bands in the experimental spectrum of butyl acrylate and eight bands for butyl methacrylate.The observed bands values and assignments are shown in Table 2.
The first point that should be discussed is the influence of the basis sets used during the simulation for each level of theory.We notice a bit changes in each spectrum but this would not affect the general quality of the functionals used.The differences in the frequency for a given band do not exceed 30 cm −1 .
When exploring the different Figures 3,4,5,6,7,8,9,10,11,12,13,14,15, and 16, we notice that the accuracy of the simulated spectra for the mentioned functionals would be classified into three classes.So we calculate the percentage of the sum of square errors of each level of theory.The hundred percent is the highest value of the sum of square errors found within the 28 spectra calculated for each molecule.This is explicitly shown in where %Δ vib  is the percentage of the sum of square errors of each functional/basis set level of theory, Δ vib  is the sum of square errors of a given functional/basis set system, and Δ vib max is the maximum value of Δ vib  found through the 28 systems (functionals/basis set).Then the ranking is made on the basis of %Δ vib  values.When %Δ vib  < 1%, we have good correlated spectra, 1% ≤ %Δ vib  ≤ 4% acceptable spectra, and for %Δ vib  > 4%, we have bad correlate spectra.

Bad Correlated Spectra. Figures 3 to 6 represent a com-
parison between the experimental infrared spectrum and the simulated spectra for butyl acrylate and butyl methacrylate, respectively.In these figures, the B972 and B98 functionals are used.The spectra at these levels of theory result into a very bad reproduction of the infrared experimental spectrum of the title molecules.The calculated and scaled frequencies for both functionals less estimate the experimental bands.This deficiency overcomes the barrier of 200 cm −1 .The calculated %Δ vib  is superior to 80% for butyl acrylate and 75% for butyl    methacrylate.The B98 functional gives the maximum errors for both molecules.

Good Correlated Spectra.
The calculated and scaled frequencies at the O3LYP functional show a good correlations with the experimental infrared spectrum of the BuA molecule with Pople's basis sets see Figure 15.The majority of the calculated frequencies at O3LYP functional give an overestimation of the experimental bands.The frequencies go down about 15 cm −1 when the 6-311+G(fd,p) and augcc-p DVZ basis sets are used; however, they reproduce well the ] 2 and ] 3 bands.Regardless of inexactitudes related to the gas phase simulation environment, the O3LYP/6-31G(f,p) and O3LYP/6-311+G(df,p) levels of theory represent the best results for reproducing the infrared spectrum.The calculated %Δ vib  is of about 0.5%.
The percentage of square errors calculated for the scaled frequencies of BuMA revealed that the use of VSXC with Pople's basis sets gives the best correlations with the experimental spectrum (Figure 16) within all functionals used.The %Δ vib  calculated are VSXC 6-31g(d,p) (0.71%) < VSXC 6-311+g(df,p) (0.97%).

Conclusion
We have tested through this work the applicability and the capability of the scaling factors proposed in the literature to reproduce the experimental infrared spectra for two medium-sized molecules composed of about twenty-four atoms, ten of them are heavy.The B972 and B98 are not recommended for frequency calculation as they show a big deviation from the observed bands.The HCTH and OLYP functionals give acceptable results especially when we consider the environmental approximation between the gas phase simulation and the infrared spectra recorded in the liquid state.The O3LYP/6-31G(f,p) and O3LYP/6-311+G(df,p) levels of theory represent the best choice to study the infrared spectrum of butyl acrylate molecule.In the    Organic Chemistry International case of butyl methacrylate, VSXC/6-31g(d,p) shows the best results.Nevertheless, the hybrid DFT functional O3LYP is a good competitor for reproducing the infrared spectrum of the latter molecule.The G96LYP contains the Gill96 exchange functional, which has a rather simple form with only one parameter; thus, frequencies' scaling factors closest to 1 would be probably more concise for large molecules regarding time consuming.
The discrepancies between the obtained spectra within all the functionals used are probably due to the quality of the exchange correlation equations used.The method used to obtain the parameters of the functionals.The set of molecules used for parameterization and the fraction of HF exchange incorporated within the hybrid functionals play an important role to improve their quality.

Figure 3 :
Figure 3: Comparison between the experimental infrared spectrum of butyl acrylate and the simulated spectra using different basis sets at the B972 functional.

Figure 4 :Figure 5 :
Figure 4: Comparison between the experimental infrared spectrum of butyl acrylate and the simulated spectra using different basis sets at the B98 functional.

Figure 7 :
Figure 7: Comparison between the experimental infrared spectrum of butyl acrylate and the simulated spectra using different basis sets at the G96LYP functional.

Figure 8 :− 1 )
Figure 8: Comparison between the experimental infrared spectrum of butyl methacrylate and the simulated spectra using different basis sets at the G96LYP functional.

Figure 10 :
Figure 10: Comparison between the experimental infrared spectrum of butyl acrylate and the simulated spectra using different basis sets at the OLYP functional.

Figure 11 :
Figure 11: Comparison between the experimental infrared spectrum of butyl acrylate and the simulated spectra using different basis sets at the VSXC functional.

Figure 12 :
Figure 12: Comparison between the experimental infrared spectrum of butyl methacrylate and the simulated spectra using different basis sets at the HCTH functional.

Figure 13 :Figure 14 :
Figure 13: Comparison between the experimental infrared spectrum of butyl methacrylate and the simulated spectra using different basis sets at the OLYP functional.

Figure 15 :
Figure 15: Comparison between the experimental infrared spectrum of butyl acrylate and the simulated spectra using different basis sets at the O3LYP functional.

Figure 16 :
Figure 16: Comparison between the experimental infrared spectrum of butyl methacrylate and the simulated spectra using different basis sets at the VSXC functional.

Table 1 :
The vibrational frequency scaling factors evaluated for different exchange correlation functionals and basis sets.

Table 2 :
Observed bands for butyl acrylate and butyl methacrylate with their assignments.