Journal of Spectroscopy

Journal of Spectroscopy / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 2073613 |

Ramzan Ullah, Ishaq Ahmad, Yuxiang Zheng, "Fourier Transform Infrared Spectroscopy of “Bisphenol A”", Journal of Spectroscopy, vol. 2016, Article ID 2073613, 5 pages, 2016.

Fourier Transform Infrared Spectroscopy of “Bisphenol A”

Academic Editor: Eugen Culea
Received02 May 2016
Revised02 Jul 2016
Accepted19 Jul 2016
Published23 Aug 2016


FTIR (400–4000 cm−1) spectra of “Bisphenol A” are presented. Absorption peaks (400–4000 cm−1) are assigned on the basis of Density Functional Theory (DFT) with configuration as B3LYP 6-311G++ (3df 3pd). Calculated absorption peaks are in reasonable reconciliation with experimental absorption peaks after scaling with scale factor of 0.9679 except C-H and O-H stretching vibrations.

1. Introduction

“Bisphenol A” is a carbon based synthetic compound used primarily in the plastic industry and epoxy resins. BPA mimics hormone-like behavior which questions upon its use in food containers and baby bottles and so forth. BPA has been controversially associated with a number of adverse health effects including neurotoxicity, genotoxicity, and carcinogenicity [1]. BPA (C15H16O2) molecule is shown in Figure 1. BPA has been used tremendously in many areas like optical media, food containers, medical devices, electronics, protective coatings, and automotive industry are very few to mention. Besides that, the debate over the use of BPA for public health safety continues where US Food and Drug Administration (FDA) considers BPA totally safe for food containers and packaging [2] but some other researches oppose this idea [3]. In this state of perplexity, banning the use of BPA in some applications mostly related to foods and human contact is on the rise [4]. Situation further gets worse when some researches show that “Bisphenol S” is also toxic like BPA which is a common replacement for BPA to bypass the law [5]. Keeping in view its widespread use and significance, its optical properties are very important. Some environmental hormones have already been studied by THz, Raman, and FTIR spectroscopy [6]. BPA has also been studied by THz spectroscopy where one of its vibrational modes was observed but in FTIR spectroscopy, main focus was on CH and OH stretching [7]. Various molecular vibrations of BPA have been identified through Raman spectroscopy by our research group [8]. However, Raman spectroscopy and FTIR spectroscopy are complementary to each other. To get a complete picture of the molecular dynamics, both Raman and FTIR spectroscopic studies are essential. So, here we present a complementary FTIR spectroscopic study of BPA which will help a lot to understand the behavior of this molecule to a wide extent. The goal of this study is to find molecular vibrations of BPA in the Infrared regime which will pave the way to a better understanding of the interaction of this molecule and in finding the origin of its toxicity. This study may also help researchers and health agencies to devise better ways for the characterization and detection of BPA.

2. Experimental Set-Up

Spectrum 100 FTIR [9] from PerkinElmer is used to record the spectra from 400 to 4000 cm−1. BPA, CAS number (80-5-7), is purchased from “Sinopharm Chemical Reagent Co. Ltd.” [10]. Sample as well as Potassium Bromide (KBr) is ground to very fine powder and then a tablet is formed under high pressure after mixing them. This tablet is then put on the optical path of Infrared radiation to take the measurements. All the measurements were taken at room temperature and atmospheric pressure.

3. DFT Calculations

Gaussian 09 package [11] is utilized to carry out DFT related calculations which are optimized to minimum and Ground State method B3LYP with 6-311G basis set [12] is used together with (3df, 3pd) and ++ diffuse functions [13]. MOLVIB [14, 15] is used to calculate Potential Energy Distributions (PEDs) given in Table 1. Another table given in the Supplementary Material contains more information about simulation and values of important parameters obtained like energy, dipole moment and so forth (see Supplementary Material available online at

Sr. numberExperimental FTIR (cm−1)Simulated (B3LYP) cm−1PED (%)

1531543ν CCC4 (36), γ COob (12), ν CCCR2 (9), ν CCCR1 (9), γ CCob (5)
2552563γ COob (28), γ CCob (15), ν CCC4 (14), γ CHob (12), τ R1t (7), τ R2t (6)
3564572γ COob (34), γ CHob (15), γ CCob (13), ν CCC4 (12), τ R2t (9), τ R1t (8)
4721734τ R1t (23), τ R2t (23), γ COob (20), γ CCob (19), γ CHob (7)
5734768ν CC4 (26), ν CCC4 (20), ν CO (8), δ CCHb (6), ν CCHa (6), ν CCCR1 (5)
6758825γ CHob (95)
7758826γ CHob (95)
8827844ν CO (23), ν CCR2 (14), ν CCR1 (14), ν CCCR2 (12), ν CCCR1 (12), ν CC4 (6)
9827851γ CHob (43), ν CO (9), ν CCR2 (9), ν CCR1 (8), γ COob (6)
10827853γ CHob (70), γ COob (11)
11827856γ CHob (52), γ COob (8), ν CCR1 (5), ν CO (5), ν CCR2 (5)
1210131026ν CCHa (28), ν CCHb (28), ν CCCR1 (8), ν CCCR2 (7), ν CCR1 (7), ν CCR2 (6)
1310131029ν CCCR1 (21), ν CCCR2 (21), ν CCR1 (15), ν CCR2 (15), ν CCHR1 (12), ν CCHR2 (11)
1410131033ν CCHa (16), ν CCHb (16), ν CCCR2 (13), ν CCCR1 (11), ν CCR2 (11), ν CCR1 (9)
1510831122ν CC4 (26), ν CCR1 (12), ν CCR2 (12), ν CCHR1 (9), ν CCHR2 (9), ν CCC4 (7)
1611021132ν CCR2 (15), ν CCR1 (15), ν CCHR2 (15), ν CCHR1 (15), ν CC4 (10), ν CCHb (7)
1711131139ν CCHR2 (16), ν CCHR1 (16), ν CCR2 (14), ν CCR1 (14), ν CCHa (9), ν CCHb (9)
1811491164ν CC4 (29), ν CCC4 (17), ν CCHR2 (12), ν CCHR1 (12), ν CCR2 (7), ν CCR1 (7)
1911771190ν COH (55), ν CCR1 (12), ν CCR2 (8), ν CO (7), ν CCHR1 (7)
2011771191ν COH (56), ν CCR2 (13), ν CCR1 (10), ν CCHR2 (7), ν CO (6), ν CCHR1 (5)
211218 1253 (Tentative)ν CC4 (45), ν CCC4 (13), ν CCHa (8), ν CCHb (8), ν CCR1 (6), ν CCR2 (6)
2212181281.70 (Tentative)ν CO (55), ν CCR2 (10), ν CCR1 (9), ν CCHR2 (8), ν CCHR1 (7)
2312181281.91 (Tentative)ν CO (54), ν CCR1 (9), ν CCHR1 (9), ν CCR2 (8), ν CCHR2 (8)
2412961326ν CCR1 (31), ν CCR2 (29), ν CCC4 (12), ν CCHR1 (9), ν CCHR2 (8)
2513621370.15ν CCHR2 (32), ν CCHR1 (20), ν COH (16), ν CCR2 (16), ν CCR1 (9)
2613621370.41ν CCHR1 (33), ν CCHR2 (20), ν COH (17), ν CCR1 (16), ν CCR2 (10)
2713841401δ HCHb (23), δ CCHb (23), ν HCHa (22), ν CCHa (22), ν CC4 (5)
2814351461ν CCR2 (26), ν CCHR2 (24), ν CCR1 (11), ν CCHR1 (11), ν CCC4 (10), ν COH (9)
2914461462ν CCR1 (26), ν CCHR1 (25), ν CCR2 (11), ν CCHR2 (11), ν COH (9), ν CCC4 (8)
3014631510.8ν HCHa (44), ν HCHb (44)
3114631511.55ν HCHb (46), ν HCHa (45)
3215101545ν CCHR2 (26), ν CCHR1 (24), ν CCR2 (16), ν CCR1 (14), ν CO (8)
3315101548ν CCHR1 (26), ν CCHR2 (24), ν CCR1 (16), ν CCR2 (15), ν CO (8)
3415981628ν CCR1 (33), ν CCR2 (33), ν CCHR1 (7), ν CCHR2 (7), ν COH (6)
3516121653ν CCR2 (39), ν CCR1 (22), ν CCHR2 (14), ν CCHR1 (8), ν CCCR2 (6)
3616121654ν CCR1 (39), ν CCR2 (22), ν CCHR1 (14), ν CCHR2 (8), ν CCCR1 (6)
3728703029ν CH3a (59), ν CH3b (41)
3828703034ν CH3b (59), ν CH3a (41)
3929333093ν CH3a (57), ν CH3b (42)
4029643098ν CH3b (56), ν CH3a (44)
4129753104ν CH3a (77), ν CH3b (22)
4229753105ν CH3b (79), ν CH3a (20)
4333373833.87ν OH (100)
4433373834.41ν OH (100)

Multiple assignments.
Ob: out of plane bending; t and τ: torsion. R1: Ring 1, R2: Ring 2, a: Label a, b: Label b, ν: stretching, δ: in-plane deformation, γ: out-of-plane deformation, and R: ring.
Scale factor for DFT B3LYP 6-311++G (3df, 3pd) is 0.9679.

4. Results and Discussion

Figure 2 shows experimental spectra of “Bisphenol A” in the range 400–1700 cm−1 along with calculated DFT spectra in absorbance mode for easy comparison. Absorption peaks are reasonably in good agreement in the region 400 to 1700 cm−1 after scaling with scale factor of 0.9679 for vibrational frequencies calculated by Andersson and Uvdal [16]. However, the large deviation of experimental C-H and O-H stretching vibrations with the DFT peaks as shown in Figure 3 has been addressed by Ullah et al. [7].

Experimental peak of 531 cm−1 is assigned to 543 cm−1 B3LYP which is a combination of different motions including CO, CC out of plane bending, and CCC interactions in both the rings and otherwise. 531 cm−1 experimental absorption peak is not clearly visible in Figure 2 which lies behind 552 cm−1. It is clear from Table 1 that 552 cm−1 and 564 cm−1 experimental peaks are mainly due to CO, CC, and CH out of plane bending along with little amount of torsions in both the rings and hence are assigned to 563 cm−1 and 572 cm−1 (B3LYP), respectively, as shown in Figure 2. FTIR absorption peak of 721 cm−1, which, according to Table 1, is assigned to simulated value of 734 cm−1 and corresponding PED analysis shows that the major contribution in this vibrational mode comes from the torsions in both the rings and CC out of plane bending. Absorption peak of 734 cm−1 is assigned to 768 cm−1 which is caused by CC stretching. Next absorption peak, 758 cm−1 which is clearly marked with its value in Figure 2, is really interesting to see that it is due to the CH out of plane bending motion and we can clearly see that there is only one experimental peak with this value but B3LYP shows two distinct peaks of 825 cm−1 and 826 cm−1 in Table 1. Hence, we assign 758 cm−1 to both of these values. Moreover, this CH out of plane bending does not come from single carbon and hydrogen atoms. It is actually the combination of all the CH out of plane bending motions in both the rings. A relatively high absorption peak of 827 cm−1 (Figure 2) is assigned to 4 simulated vibrational modes of 844, 851, 853, and 856 cm−1 due to their closeness. The major contribution in 844 cm−1 mode is from CO stretching and CC stretching in both the rings. The other three come mainly from CH out of plane bending movements. So the experimental peak 827 cm−1 might belong to any of them. In a similar fashion, 1013 cm−1 is assigned to 3 simulated values of 1026, 1029, and 1033 cm−1. Interestingly, the first and last peaks (1026 and 1033 cm−1) have similar contributions as given by PED analysis in Table 1 but the middle one 1029 cm−1 has different motions. However, they are too close to each other that DFT also suggests a single experimental peak in such circumstance. So 1013 cm−1 has once again multiple assignments. 1083 cm−1 and 1102 cm−1 are related to B3LYP values of 1122 cm−1 and 1132 cm−1, respectively, which are contributed mostly by CC stretching in both the rings. 1113 cm−1 and 1149 cm−1 are hereby assigned to 1139 cm−1 and 1164 cm−1, respectively, according to Table 1 and relative PED analysis shows that these modes are mainly because of CCH and CC interactions in the molecule. Multiple assignment of the experimental peak of 1177 cm−1 to both 1190 cm−1 and 1191 cm−1 is due to their nearness. The respective PED analysis shows that the interactions are similar in such a way that for one DFT peak 1190 cm−1, it is the motion in one ring and for the other 1191 cm−1, it is the same motion but in the second ring. Since their frequencies are bit different, DFT gives two peaks but experimentally, it shows only one as shown in Figure 2. Immediately next to the 1177 cm−1 peak, there is rather a broadband peak with small fluctuations in Figure 2. Its value is given in Table 1 as 1218 cm−1 which is the highest point along the broad peak. If we compare it with DFT spectra, a small peak of 1253 cm−1 might be assigned to it. However, due to broadband nature of this curve, it is hard to assign it with certain degree of confidence. Moreover, DFT spectra also show a rather high peak adjacent to 1253 cm−1, as given in Figure 2, which is actually two peaks (1281.70 cm−1 and 1281.91 cm−1) according to Table 1. These two peaks can also be assigned to 1218 cm−1. However, the shape of the experimental curve does not show a good behavior which might be due to some defect and hence tentative assignment is given here. PED analysis shows the former as a result of CC stretching and the last two from CO stretching. 1296 cm−1 is assigned to 1326 cm−1 which is from CC stretching in two rings. 1362 cm−1 is assigned to two absorption peaks of 1370.15 cm−1 and 1370.41 cm−1 resulting from CCH interactions in the rings but with different configurations. In a similar fashion, 1384 cm−1 can be assigned to 1401 cm−1 originating from HCH and CCH interactions. 1435 cm−1 and 1446 cm−1 are attributed to 1461 cm−1 and 1462 cm−1, respectively, resulting from CC and CCH in the rings. 1463 cm−1 which is not very clear in Figure 2 is assigned to two DFT peaks, 1510.8 cm−1 and 1511.55 cm−1, both of which are resulting from HCH interactions (not occurring in the rings) in different ways as shown in PED analysis. A sharp absorption peak of 1510 cm−1 in Figure 2 goes to 1545 cm−1 and 1548 cm−1 where CCH and CC modes are prevalent. 1598 cm−1 is assigned to calculated frequency of 1628 cm−1 and similarly experimental peak of 1612 cm−1 to 1653 cm−1 and 1654 cm−1 (multiple). Spectra end here and after a long pause, peaks start appearing again near 3000 cm−1 where C-H and then O-H stretching interactions exist which are clearly marked in Figure 3. The detailed analysis and assignment of these modes can be found in [6]. Spectra in the range 1700 to 2800 cm−1 are not shown due to nonexistence of any peak.

Simulated values of B3LYP given in Table 1 are raw values. These values can be corrected by using a scale factor of 0.9679 as discussed earlier in this section. DFT values are rounded off to nearest whole number. Only major peaks are labeled in Figures 2 and 3. The detail of the absorption peaks both experimental as well as calculated is given in Table 1. Assignment is done on the basis of comparison of the relative intensities of the experimental and calculated spectra.

According to DFT calculations, there are 14 C-H stretching vibrations in different combinations between 3000 and 3300 cm−1. Only few of them are visible in Figure 2. As “Bisphenol A” is relatively large molecule and many peaks are very near to each other, so it is hard to assign the individual peaks distinctly without multiple assignments as indicated at the end of Table 1.

5. Conclusion

We have presented FTIR (400–4000 cm−1) spectra of “Bisphenol A.” Peaks (400–4000 cm−1) have been assigned on the basis of Density Functional Theory (DFT) with configuration as B3LYP 6-311G++ (3df 3pd). Calculated absorption peaks are in reasonable agreement with experimental peaks after scaling with scale factor of 0.9679 except C-H and O-H stretching vibrations.

Competing Interests

The authors declare that they have no competing interests.


This work has been partially supported by the National Natural Science Foundation of China (no. 61275160). The authors would like to thank Professor S. Y. Wang for providing computer resources.

Supplementary Materials

Supplementary material contains information about Density Functional Theory (DFT) calculations done for Bisphenol A (BPA). Values of important parameters like Dipole Moment, Energy and others, are also given.

  1. Supplementary Material


  3. A. M. Hormann, F. S. Vom Saal, S. C. Nagel et al., “Holding thermal receipt paper and eating food after using hand sanitizer results in high serum bioactive and urine total levels of bisphenol A (BPA),” PLoS ONE, vol. 9, no. 10, Article ID e110509, 2014. View at: Publisher Site | Google Scholar
  5. R. Viñas and C. S. Watson, “Bisphenol S disrupts estradiol-induced nongenomic signaling in a rat pituitary cell line: effects on cell functions,” Environmental Health Perspectives, vol. 121, no. 3, pp. 352–358, 2013. View at: Publisher Site | Google Scholar
  6. R. Ullah, S. U.-D. Khan, M. Aamir, and R. Ullah, “Terahertz time domain, Raman and fourier transform infrared spectroscopy of acrylamide, and the application of density functional theory,” Journal of Spectroscopy, vol. 2013, Article ID 148903, 17 pages, 2013. View at: Publisher Site | Google Scholar
  7. R. Ullah, H. Li, and Y. Zhu, “Terahertz and FTIR spectroscopy of ‘Bisphenol A’,” Journal of Molecular Structure, vol. 1059, no. 1, pp. 255–259, 2014. View at: Publisher Site | Google Scholar
  8. R. Ullah and Y. Zheng, “Raman spectroscopy of ‘Bisphenol A’,” Journal of Molecular Structure, vol. 1108, pp. 649–653, 2016. View at: Publisher Site | Google Scholar
  11. M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 09, Revision A.01, Gaussian, Inc, Wallingford, Conn, USA, 2009.
  12. V. A. Rassolov, J. A. Pople, M. Ratner, P. C. Redfern, and L. A. Curtiss, “6-31G basis set for third-row atoms,” Journal of Computational Chemistry, vol. 22, no. 9, pp. 976–984, 2001. View at: Publisher Site | Google Scholar
  13. T. Clark, J. Chandrasekhar, G. W. Spitznagel, and P. V. Schleyer, “Efficient diffuse function-augmented basis sets for anion calculations. III. The 3-21+G basis set for first-row elements, Li-F,” Journal of Computational Chemistry, vol. 4, no. 3, pp. 294–301, 1983. View at: Publisher Site | Google Scholar
  14. T. Sundius, “Molvib—a flexible program for force field calculations,” Journal of Molecular Structure, vol. 218, pp. 321–326, 1990. View at: Publisher Site | Google Scholar
  15. T. Sundius, “A new damped least-squares method for the calculation of molecular force fields,” Journal of Molecular Spectroscopy, vol. 82, no. 1, pp. 138–151, 1980. View at: Publisher Site | Google Scholar
  16. M. P. Andersson and P. Uvdal, “New scale factors for harmonic vibrational frequencies using the B3LYP density functional method with the triple-ζ basis Set 6-311+G(d,p),” Journal of Physical Chemistry A, vol. 109, no. 12, pp. 2937–2941, 2005. View at: Publisher Site | Google Scholar

Copyright © 2016 Ramzan Ullah et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.