Abstract

The Fourier transform infrared (FTIR) spectrum of 2-hydroxy-2-methyl-1-phenylpropan-1-one has been measured in the region 4000–700 . The most stable conformation of title molecule was found after a careful potential energy surfaces study. The molecular geometry, vibrational frequencies, and infrared intensities have been calculated by using ab initio HF and density functional theory calculation B3LYP with 6-311+ basis set. Scaled frequencies and potential energy distribution were calculated for band assignment. We found an excellent agreement between the experimental and the simulated spectra. Energy gap between HOMO and LUMO explains the eventual charge transfer interactions taking place within the molecule.

1. Introduction

The radiation curing industry is one of the most rapidly developing fields in the entire coating industry. The low toxicity, cheapness, speed, control, and ease of formulation and operation are some of the main advantages of this growing technology. UV and/or visible light radiation is used to induce photochemical polymerization or cross-linking of a monomer, oligomer, or prepolymer formulation containing a certain type of unsaturation, such as an acrylic group, and an appropriate initiator [1, 2]. Among the large group of known photoinitiator structures, hydroxyalkylphenones, such as 2-hydroxy-2-methyl-1-phenylpropan-1-one (Figure 1), have gained much interest due to their high reactivity and universal applicability [3, 4].

Despite the high practical value of title molecule, its physicochemical properties have been relatively little looked into [812].

To our best knowledge no structural data and detailed interpretation of the vibrational spectra of the title molecule are presented in the literature. This prompted us to look into the vibrational spectroscopy of the title molecule more carefully.

2. Experimental

2-Hydroxy-2-methyl-1-phenylpropan-1-one (purity 97%) was purchased from Sigma-Aldrich and used as such without further purification. The infrared spectra of liquid films placed between KBr pellets were recorded within 4000–700 cm−1 range with a Perkin Elmer FTIR System-2000 model.

3. Computational

Potential energy curves were obtained by calculating the variation of total energy of title molecule with dihedral angles. For this purpose the calculations were performed with B3LYP functional and 6-311+G** basis set. The dihedral angles were varied from 0° to 360° with an interval of 10°.

All DFT calculations were performed using Becke’s 3 parameter (local, nonlocal, Hartree-Fock) hybrid exchange functionals with Lee-Yang-parr correlational functionals (B3LYP) [13, 14].

The harmonic vibrational frequencies of the stable conformation were calculated at the same level of theories used for the calculated optimized geometry. Two levels of theory were used for vibrational analysis B3LYP/6-311+G** and HF/6-311+G**; the calculated frequencies were scaled down by the wavenumber linear scaling (WLS) procedure of Yoshida et al. [15, 16] using the following equation, for DFT: and the scaling factor proposed by Baker et al. [17] for HF. The scaling procedures are used due to the fact that the anharmonic effects are neglected when calculating wavenumbers at HF and DFT theories, which make the calculated value higher than the experimental one. All the calculations were carried out with the Gaussian 03 program [18]. The vibrational assignments of normal modes were provided on the basis of PED calculation using the program GAR2PED [19].

4. Conformational Analysis

Due to the flexibility of 2-hydroxy-2-methyl-1-phenylpropan-1-one and possible rotational isomerism, we examined the variation of energy with respect to , , and dihedral angles. From Figure 2 we see that the energy is increasing when increasing dihedral angle to a maximum in energy of −3.25 105 kcal/mol at 180°. There is only one minimum at 0° with an energy of −3.38 105 kcal/mol.

It is well seen from Figures 3 and 4 that the molecule takes two minima at 0° and 180°. In Figure 3 the difference in energy between the two minima is 419.7 kcal/mol, and it is 8.785 kcal/mol in Figure 4.

The potential energy surfaces revealed the existence of tree conformers. The energy of the most stable conformation is −3.38 105 kcal/mol.

5. Geometry Optimizations

The most stable structure is recognized as true minima due to the lack of imaginary harmonic frequencies. The geometry parameters are represented in Table 1. To the best of our knowledge, exact experimental data of the geometrical parameters of the 2-hydroxy-2-methyl-1-phenylpropan-1-one are not available in the literature. Therefore our optimized structural parameters are compared with the XRD data of closely related molecules benzoic acid and butyric acid [5, 6]. The theoretical results are almost comparable with the experimental XRD values. The carbon-oxygen = 1.421 and carbon-carbon = 1.554, = 1.540, and = 1.540 bond lengths are in agreement with the XRD values [7] 1.38, 1.55, 1.54, and 1.54, respectively. In the phenyl ring, the bond lengths and the bond angles are in good agreement with the experimental one [7]. Thus the ring is slightly distorted and smaller than 120° at the point of substitution.

In Table 1, a small difference between the two levels of theory HF and DFT should be noticed for the calculated coordinates. The bond lengths are less estimated by HF with a maximum difference of 0.01A°. However, this difference is about 0.02A° for C14O23 bond. The B3LYP calculations lead to a planar structure of the phenyl ring, ketone and hydroxyl groups. Otherwise, these groups lie in a near planar form according to the HF optimization. This could be seen directly from dihedral angles (, , , , , and ) deviating about 10° from those obtained by B3LYP level of theory.

The zero point corrected energy for the most stable structure of the title molecule is −337.911 103 kcal/mol at B3LYP/6-311+G** and −335.786 103 kcal/mol at HF/6-311+G**. The energy calculated by DFT is lower showing more stability than that calculated at HF method.

6. Vibrational Analysis

The title compound contains 24 atoms leading to 66 normal vibrational modes. We compared the experimental IR spectrum with the computed IR spectrum, and we checked whether we can recognize all the different bands. Table 2 represents the calculated and scaled fundamental wavenumbers, intensities of vibrational peaks, and PED calculations along the internal coordinates obtained by B3LYP/6-311+G** level of theory. The corresponding experimental wavenumbers together with assignments and potential energy distribution are also reported in this table. Assignments have been made on the basis of PED, and they are given as per the internal coordinate system recommended by Pulay et al. using DFT [20]. The experimental and theoretical spectra are shown in Figure 5.

6.1. Phenyl Ring Vibrations

The substituted-like molecule gives rise to CH stretching, CH in-plane and out-of-plane bending vibrations. The heteroaromatic structure shows the presence of CH stretching vibration in the region 3100–3000 cm−1 which is the characteristic region for the ready identification of CH stretching vibration [15, 21]. In this region the bands are not affected appreciably by the nature of the substituent. The aromatic CH stretching frequencies rise from the modes observed at 3062, 3047, 3060, and 3080 cm−1 of benzene and its derivatives [22].

In this work, the CH stretching vibrations are observed at 3025, 3060, 3070, 3090, and 3108 cm−1 in the FTIR. The calculated values of these modes for the title molecule have been found to be 3031, 3041, 3051, 3066, and 3085 cm−1 at the B3LYP/6-311+G** level of theory. The calculated PEDs for these normal modes are greater than 85% (see Table 2).

The CH-in plane bending vibrations usually occur in the region 1430–990 cm−1 and are very useful for characterization purposes. The band due to CH-in plane ring vibrations interacting somewhat with CC stretching vibration is observed as a number of medium-to-weak sharp bands in the region of 1300–1000 cm−1 [23, 24]. In this study, the FTIR peaks at 1077, 1170, 1192, and 1315 cm−1 are assigned to CH in-plane bending vibrations as reported in Table 2. The strong peak below 900 cm−1 indicates clearly the aromatic aspect of the molecule. Substitution patterns on the ring can be judged from the CH out-of-plane bending in the region 900–650 cm−1, and these bands are highly informative [25]. The CH out-of-plane bending vibrations are strongly coupled vibrations in the region 900–667 cm−1 [26]. In the present study, the calculated frequencies at 994, 950, 857, 721, and 697 cm−1 may be assigned to the experimental bands at 934, 836, and 716 cm−1. The calculated PEDs show that these bands are not pure CH out-of-plane bending vibrations.

The ring carbon-carbon stretching vibration occurs in the region 1625–1430 cm−1. In general, the bands are of variable intensity and are observed at 1625–1590, 1590–1575, 1540–1470, 1460–1430, and 1380–1280 cm−1 from the frequency ranges given by Varsanyi [27] for the five bands in the region. In this work, the frequencies observed in the FTIR spectrum at 1598, 1578, 1489, and 1447 cm−1 are well reproduced by the theoretical frequencies at 1607, 1587, 1497, and 1452 cm−1. The calculated PEDs show that these modes are not pure CC stretching vibrations. The normal mode no. 31 is attributed to the ring breathing vibration; the calculated PED for this mode is 61%.

6.2. The CC Vibrations

The CC stretching vibrations of the title molecule occur in the region from 1000 to 1600 cm−1. The majority of these modes of vibration are combined with CH and OH bending vibrations. The calculated PEDs for CC stretching are of about 13% for normal modes #37, #38, and #39. However the normal mode #40 is attributed to the CC stretching vibrations with a potential of 37% and less than 10% for normal modes #51 and #53. The CCC in-plane bending frequencies are found at 372, 443, 529, 618, 701, 860, and 923 cm−1.

The CC torsions are likely to occur in the region below 500 cm−1. For the title molecule, we found five normal modes representing the torsion about the bond . The calculated wavenumbers are 22, 79, 163, 265, and 325 cm−1. The calculated PEDs are about 50%. These modes are not pure thereby they are generally mixed with other CO and CH3 torsions.

6.3. The C=O Vibrations

The C=O stretch of carboxylic acids is identical to the C=O stretch in ketones, which is expected in the region 1740–1660 cm−1 [28]. Krishnakumar et al. [29] observed very strong band at 1661 cm−1 in the IR spectrum for 1-naphthyl acetic acid, and it is assigned to the C=O stretching vibrations. In the present work a strong band observed in the FTIR spectrum at 1683 cm−1 is assigned to the C=O stretching vibrations. The calculated peak is at 1677 cm−1 with PED value of 82% as reported in Table 2.

6.4. The COH Vibrations

The CO stretching vibration is assigned to the band at 959 cm−1 in the FTIR spectrum. The predicted frequency shows a good correlation. The PED value is 83% for mode #27 and 22% for mode #28; most of the CO vibrations are mixed vibrations as shown in Table 2.

The OH group gives rise to three vibrations, namely, stretching, in-plane bending, and out-of-plane bending. The OH group vibrations are likely to be the most sensitive to the environment, so they show pronounced shifts in the spectra of the hydrogen-bonded species. In the case of unsubstituted phenol it has been shown that the frequency of OH stretching vibration in the gas phase is 3657 cm−1 [27]. Subramanian et al. [30] report the infrared spectrum of 2,5-di-tert-butylhydroquinone and assign the strong band at 3402 cm−1 to the OH stretching vibration. Koczon et al. [22] studied the experimental and theoretical IR and Raman spectra of picolinic, nicotinic, and isonicotinic acids and assigned the bands at 3437, 3447, and 3436 cm−1 for these molecules to the OH stretching vibration. In the FTIR spectrum of the title compound, the band observed at 3458 cm−1 is assigned to the OH stretching mode. The computed wavenumber is at 3482 cm−1 with a PED value of 85%. The carboxylic acids show in-plane bending band of OH at 1134 cm−1 (FTIR). This band is well reproduced by the calculated and scaled one at 1151 cm−1. The PED calculation predicts the existence of another OH in-plane bending at 1365 cm−1. This band is well correlated with the one at 1363 cm−1 in the experimental infrared spectrum. The calculated frequencies at 416, 418, 445, 506, and 807 cm−1 are assigned to the OH out-of-plane bending. The last band is well correlated with the experimental one at 796 cm−1. These normal modes are mixed with the skeleton bending vibrations and torsions. The predicted PEDs are about 50%.

6.5. The Methyl Group Vibrations

For the assignment of CH3 group frequencies, nine fundamentals can be associated to each CH3 group [3133]. The CH stretching in CH3 occurs at lower frequencies than those of aromatic ring (3100, 3000 cm−1). Moreover, the asymmetric stretch is usually at higher wavenumber than the symmetric stretch. In the present work, in the experimental spectrum the CH3 asymmetric stretching frequency is assigned to the bands at 2991 and 2979 cm−1, whereas the symmetric stretching frequency is assigned to the band at 2935 cm−1. These bands are well reproduced by the calculated ones at 2995, 2992, 2979, and 2977 cm−1 for the asymmetric stretching and 2920, 2917 cm−1 for the symmetric stretching vibrations. The PED values for these normal modes are greater than 70%.

For methyl substituted benzene derivatives, the antisymmetric and symmetric deformation vibrations of methyl group normally appear in the region 1465–1440 cm−1 and 1390–1370 cm−1, respectively [21, 31, 34]. We attributed the observed doublet in the FTIR spectrum at 1465 and 1457 cm−1 to the asymmetric deformation of CH3 vibrations. These bands are well correlated with those calculated at 1480, 1473, and 1462 cm−1 with PED values greater than 80%. The bands at 1374 and 1387 cm−1 (FTIR) are assigned to the symmetric CH3 deformation. The simulated frequencies at 1397 and 1377 cm−1 represent good correlations for the symmetric CH3 deformations. The normal modes no. 28, #30, #33, #38, and #41 are assigned to the CH3 rocking vibrations. These vibrations find location with the CH ring bending vibrations in the range from 984 to 1280 cm−1. This is in concordance with Subramanian et al. [30] assignments of the infrared bands of 2, 5-di-tert-butyl-hydroquinone. They reported that the rocking vibrations of methyl group for this molecule take location in the range from 837 to 1247 cm−1, and they are mixed with the CH ring vibrations.

7. HOMO-LUMO Analysis

The total energy, energy gap, and dipole moment have an effect on the stability of a molecule. In this work we calculate these parameters for title compound in gas phase and solvent. For this purpose we used B3LYP/6-311+G** level of theory. Obtained results are reported in Table 3.

The energy gap between the HOMO and LUMO is an important parameter in determining molecular electrical transport proprieties and molecular stability. The energy gap between the HOMO and the LUMO is a critical parameter in determining the molecular electrical transport proprieties because it is a measure of electron conductivity. The energy values of HOMO are computed −7.456, −7.374, and −7.401 eV and LUMO are −2.231, −2.204, and −2.258 eV, and the energy gap values are 5.225, 5.17, and 5.143 eV in gas phase, DMSO and chlorophorm for 2-hydroxy-2-methyl-1-phenylpropan-1-one molecule, respectively. Lower value of the energy gap explains the eventual charge transfer interactions taking place within the molecule. In order to understand the bonding scheme of the title molecule we draw the surfaces of the orbitals frontier. We examine the four important molecular orbitals: the highest and the second highest occupied MOs and the lowest and the second lowest unoccupied MOs. These molecular orbitals are represented in Figure 6.

According to Figure 6 the charge distribution for the HOMO is localized on the C=C ring, OH, and CH of methyl group. The LUMO is characterized by a charge distribution at CH bond on the ring and bond.

8. Conclusions

A meticulous conformational search of the title molecule leads to a unique stable conformation. For this purpose the calculations were conducted with reliable DFT and HF levels of theory using 6-311+G** basis set. The fully optimized geometries at HF and DFT (B3LYP) were compared with the experimental data of related compounds presented in the literature. The comparison shows a good agreement. Infrared spectra were recorded, and the vibrational bands were assigned on the basis of the PED calculations obtained from B3LYP/6-311+G** calculations. A very good agreement between the observed and the calculated wavenumbers was observed. The scaled frequencies at HF/6-311+G** level of theory are in general greater than those calculated at B3LYP/6-311+G**. The latter ones are more accurate. Through the molecular orbital analysis we found that the energy gap between HOMO and LUMO is lower; this would explain the eventual charge transfer in the title molecule.