Abstract
A test on calculated vibrational modes of 2-(2,6-dichlorophenyl)-N-(1,3-thiazol-2yl) acetamide using ab initio density functional method has been performed. The calculated harmonic vibrational frequencies are scaled via two schemes, uniform, , and nonuniform, . Scaled vibrational modes are compared with experimental FTIR bands. A good correlation is shown between scaled frequencies with the correlation coefficient, . This clearly shows that both schemes efficiently reproduce observed spectrum. However, a close investigation of individual normal modes reveals that nonuniform scaling performs much better than uniform scaling especially in the high frequency region.
1. Introduction
Ab initio quantum chemical methods are reliable tools for extracting valuable information about any chemical system. Currently available ab initio methods are capable of reproducing experimental IR spectra of rigid and small- as well as medium-size molecules. In many chemical systems, the theoretically calculated frequencies and intensities of vibrational modes closely resemble the experimentally observed values. Hence, it may be expected that the normal modes analysis from theoretically predicted force field reflects well the modes of vibrations in a real molecule. Normal-mode analysis is commonly employed in the interpretation of the vibrational spectra of molecules [1–4]. However, there are some limitations of the model used to calculate the vibrational frequencies. These limitations include the effect of electron correlation [5] and basis set deficiencies as well as inter- and intramolecular interactions. It needs to be emphasized that the calculated frequencies represent vibrational signature of the molecule in its gas (isolated) phase, unlike experimental spectra which are recorded in liquid or solid phase of sample with impurity. This is experimentally verified that [6]. On comparing the vibrational spectra, it is observed that, at places where calculated frequencies come out to be lower than the experimental frequencies, weak intermolecular interactions may be taking place. From a practical point of view, the main disadvantage of vibrational spectroscopy is the lack of a direct spectra-structure relation. Ab initio calculations follow the approximation that the chemical bonds are elastic, but it is not the case especially at higher range of frequencies due to anharmonicity of vibrations [7]. Thus, the computed frequencies within this Born-Oppenheimer approximation [8] need to be properly scaled. These anharmonic vibrations mainly include torsions and inversions, such as the umbrella modes of halomethyl radicals [9]. Similarly, in lower range of frequencies (where nearly all modes are rotational), the interaction of modes of vibration as well as molecular interaction takes place very strongly which certainly affect the calculated frequencies.
The motivation for predicting or simulating vibrational spectra is to make vibrational spectroscopy a more practical tool [10–12]. In this respect, it is not only desirable but also essential to have calculated frequencies very close to their observed values. This can usually be expected by appropriate scaling of calculated frequencies. Appropriate scaling method is still a fundamental problem existing in the literature, even today. Merrick et al. had calculated some scale factors using harmonic approximation but these factors did not work well in the region where mixing of several bands occurred [13]. Later on, Alecu et al. have calculated new scale factors within harmonic approximation for a number of theoretical methods [14]. They tested a lot of functional systems, for example, electronic model chemistries, single-level wave, doubly hybrid density, neglect-of-diatomic-differential-overlap, and so forth. On the other hand, Mauricio Alcolea Palafox derived some scaling equations by statistical correlation and suggested that these scaling equations give more satisfactory result as compared to constant scale factors [15–17].
In present paper, we perform a critical survey on the scaling of vibrational frequencies of 2-(2,6-dichlorophenyl)-N-(1,3-thiazol-2yl) acetamide by using two different schemes. The first scheme is based on a constant scale factor as suggested by Alecu et al. which we call “uniform scaling.” The second scheme is based on a scaling equation as recommended by Palafox et al., which we refer to as “nonuniform scaling.”
2. Computational Methods
Our ab initio calculations employ a hybrid form of density functional, B3LYP in conjunction with 6–31+G(d,p) basis set for geometry optimization of the test molecule. The B3LYP combines Becke’s three-parameter exchange term [18] with the functional devised by Lee-Yang-Parr to treat electron correlation [19]. Vibrational frequency calculations are performed using optimized geometry at the same level of theory. The calculated frequencies are scaled by two different schemes as below.
For uniform scaling [13],
For nonuniform scaling [15], Gaussian 09 package [20] is used to perform all computations and assignments of calculated vibrational modes are made with the help of visual animation generated by GaussView 5.0 program [21].
3. Results and Discussion
First, we discuss some salient features of test molecule in condensed as well as isolated (gas) phase. The single crystal of 2-(2,6-dichlorophenyl)-N-(1,3-thiazol-2-yl) acetamide is grown by slow evaporation method (melting point = 489–491 K) [22]. A photograph of this crystal along with X-ray spectrum can be found in supplementary information (see Figures S1 and S2 in Supplementary Material available online at http://dx.doi.org/10.1155/2014/649268). In crystal phase, molecules are linked by pairs of N−HN hydrogen bonds as shown in Figure 1.
The measured X-ray crystallographic parameters are used to model an initial structure of isolated molecule for the process of geometry optimization. After optimization, we find a sofa shaped structure in which the mean plane of dichlorophenyl ring is almost perpendicular to that of thiazole ring as shown in Figure 2.
In Figure 3, we have shown a linear correlation between calculated and observed bond-lengths. Note that C–H bond-lengths are not included here. We find a good correlation with correlation coefficient () of 0.9991 which suggests that our computational scheme is fully capable of reproducing the experimental geometry.
FTIR spectrum of test molecule is recorded by using Shimadzu-Model Prestige 21 spectrometer in the region 400–4000 cm−1 with the sample (purity of 98%) in KBr pellet. The vibrational analysis of isolated molecule (Figure 2) is performed at the same level of theory as mentioned earlier. In Figure 4, we plot calculated IR and observed FTIR spectra for wavenumber region of 1800 cm−1–600 cm−1. The calculated wavenumbers are scaled by using (1) and (2). The assignments of all vibrational modes have been made along with their directions of polarization. Table 1 lists calculated and scaled wavenumbers of all vibrational modes as well as their assignments. For the sake of comparison, FTIR values of the most prominent modes are also included.
(a)
(b)
In Figures 5 and 6, we have shown the correlation between experimentally observed wavenumbers and theoretically scaled wavenumbers. We find the correlation coefficient of 0.99638-0.99639 which ensures that both scaling schemes are equivalently appropriate, at least in the present case. In order to compare these two scaling schemes for particular vibration mode, we divide our vibrational analysis into two categories.
3.1. Higher Frequency Region (above 1400 cm−1)
This region covers C–H stretching vibrations of phenyl ring and methylene group along with N–H as well as C=O stretching modes of amide group. N–H stretching mode is observed experimentally at 3203 cm−1 and corresponding scaled frequencies are calculated at 3475 cm−1 and 3459 cm−1 by uniform scaling and nonuniform scaling, respectively. The apparent discrepancy between experimental and theoretical frequencies is a consequence of intermolecular hydrogen bonding present in condensed phase (see Figure 1) but absent in case of isolated molecule. The most prominent C–H mode associated with phenyl ring vibrations is seen at 3045 cm−1 in FTIR which is scaled at 3090 cm−1 and 3078 cm−1 by our theoretical schemes.
The C=O stretching observed at 1687 cm−1 is scaled at 1684 and 1688 cm−1, respectively. Similarly, C–C stretching of phenyl ring at 1550 cm−1 corresponds to scaled values of 1546 cm−1 and 1551 cm−1. The C–H stretching associated with methylene group is found to be coupled with N–H stretching mode at 1436 cm−1 which is theoretically scaled at 1422 cm−1 and 1428 cm−1 by uniform scaling and nonuniform scaling schemes, respectively.
3.2. Lower Frequency Region (below 1400 cm−1)
In this region, the bending modes associated with ring systems and various groups are generally observed. The in-plane bending of methylene group observed at 1286 cm−1 and 1141 cm−1 is coupled with phenyl ring vibration. The uniform scaling gives these modes at 1287 cm−1 and 1130 cm−1 while nonuniform scaling calculates at 1295 cm−1 and 1140 cm−1, respectively. Similarly, out of plane bending of methylene coupled with phenyl ring at 931 cm−1 in FTIR is scaled at 940 cm−1 and 952 cm−1 by our theoretical schemes.
The in-plane bending of thiazole ring is observed at 856 cm−1 against scaled values of 850 cm−1 and 863 cm−1. On the other hand, the bending of phenyl ring observed at 785 cm−1 is scaled at 759 cm−1 and 773 cm−1. Furthermore, the coupling of vibrations of both rings is observed at 761 cm−1 corresponding to scaled values of 755 cm−1 and 769 cm−1. Moreover, even lower normal modes at 692 cm−1 and 592 cm−1 associated with bending of thiazole rings are scaled at 692 cm−1 and 602 cm−1 by uniform scaling method, whereas at 707 cm−1 and 617 cm−1 by nonuniform scaling scheme.
The above discussion suggests that scaled wavenumbers by both schemes successfully explain observed bands in case of title molecule, giving a good correlation with experimental frequencies (see Figures 5 and 6). The comparative response of two scaling schemes is given in terms of wavenumber difference between observed and scaled bands plotted in Figure 7. A close investigation on Figure 7 reveals that nonuniform scaling performs much better than uniform scaling especially in the high frequency region. Thus, the calculated vibrational modes should be scaled by appropriate scaling equation to account for anharmonicity in vibrations. For lower frequency region, where discrepancies between calculated and observed wavenumbers are observed mainly due to mixing of vibrational modes and the effect of electron correlation, both scaling schemes perform equally well.
4. Conclusions
We have performed a test on the scaling of calculated vibrational bands in case of 2-(2,6-dichlorophenyl)-N-(1,3-thiazol-2yl) acetamide, adopting two different scaling schemes, namely, scale factor and scaling equation. Scaled normal modes are compared with observed vibrational bands. The validity of both schemes is established by the correlation between theoretically scaled frequencies and experimentally observed frequencies with the coefficient of 0.9964. The analysis of individual modes has revealed that both schemes perform equally well for the lower range of frequencies below 1400 cm−1 (torsion, wagging, etc.). In case of high frequency region (above 1400 cm−1) where most of the stretching modes occur, scaling equation provides more accurate vibrational bands as compared to scale factor. To be specific, C−H stretching of phenyl ring observed at 3045 cm−1 (in present case) is better represented by 3078 cm−1 (by scaling equation) rather than by 3090 cm−1 (by scale factor) at B3LYP/6–31+G(d,p) theory. Thus, the present assessment is supposed not only to assist further studies on the prediction or simulation of vibrational spectra of new molecules but also to assign more accurate vibrational modes for previously reported molecules in the literature.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
Ambrish K. Srivastava wishes to thank Council of Scientific and Industrial Research (CSIR), New Delhi, India, for providing a research fellowship.
Supplementary Materials
Experimental Synthesis: 2,6-Dichlorophenylacetic acid (0.240 g, 1 mmol) and 2-aminothiazole (0.1 g, 1 mmol) and 1-ethyl-3-(3-dimethylaminopropyl)-carbodiimide hydrochloride (1.0 g, 0.01 mol) were dissolved in dichloromethane (20 mL). The mixture was stirred in the presence of triethylamine at 273 K for about 3 h. The contents were poured into 100 ml of ice-cold aqueous hydrochloric acid with stirring, which was extracted thrice with dichloromethane. The organic layer was washed with saturated NaHCO3 solution and brine solution, dried, and concentrated under reduced pressure to give the title compound. Single crystals were grown from methanol and acetone mixture (1:1) by the slow evaporation method (M.P.: 489–491K).