Research Article | Open Access
Ramzan Ullah, Salah Ud-Din Khan, Muhammad Aamir, Rasheed 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, 7 pages, 2013. https://doi.org/10.1155/2013/148903
Terahertz Time Domain, Raman and Fourier Transform Infrared Spectroscopy of Acrylamide, and the Application of Density Functional Theory
We present terahertz time domain spectra of acrylamide in the frequency range from 0.2 to 2 THz with nearly constant refractive index having an average value of 1.33 and an absorption coefficient. Raman (95–3000 cm−1) and FTIR (450–4000 cm−1) spectra also show good agreement with density functional theory (DFT) B3LYP 6-311G++ (3df 3pd) calculations except C-H and N-H stretching frequencies even after scaling with scale factor of 0.9679. We use MOLVIB to rescale such frequencies to match experimental values.
A strip between the microwave and infrared region of electromagnetic spectrum is named terahertz band spanning from 0.1 to 10 THz . Terahertz time domain transmission spectroscopy has been applied tremendously in the fields of materials characterization [2–5], defense and security [6–11], telecommunication , and biomedical sciences.
Acrylamide (C3H5NO) is a heat generated toxicant found mostly in starchy foods cooked or processed at high temperature and may be a possible carcinogen according to the World Health Organization (WHO) . Moreover, high intake of acrylamide increases the risk of endometrial cancer and possibly ovarian cancer in women . This alarming situation demands an immediate and indepth exploration of acrylamide, and the most challenging task in this regard is to detect the acrylamide in foods in all forms accurately and efficiently. Acrylamide is the primary component of polymers and polyacrylamides out of which the latter is not only being used in genetic and biological laboratories as gel electrophoresis but also as flocculator and soil-conditioning agent in other fields. Molecular structure of acrylamide is shown in Figure 1.
2. Experimental Setup
Ultrafast Fiber Laser System (Near-Infrared System FemtoFiber Pro NIR from TOPTICA Photonics)  is used to generate subpicosecond THz pulses. Wavelength of the laser is set to 780 nm with an average power of greater than 140 mW. Pulse width is less than 100 fs and repetition rate is 80 MHz. Our THz setup is THz Time-Domain Spectroscopy System with fast scanning (EKSPLA) . A schematic diagram of our setup is shown in Figure 2. Beam splitter 1 (BS1) splits the laser beam in ratio 80 : 20 to reduce the optical power. Beam Splitter 2 (BS2) divides the laser beam into pump and probe beam. Split ratio is 50 : 50 by adjusting the half wave plate. Our THz emitter and detector is LTG-GaAs photoconductive switch. THz emitter is biased at 50 V and Si lens is mounted on the back of emitter and detector for enhanced collimation. Probe laser pulse hits the THz emitter after passing through the focusing lens L1 which focuses it to the gap of photoconductive antenna. THz signal is guided towards the sample and then to the detector by four off axis parabolic mirrors. Fast delay line is used to acquire the signal in real time while slow delay line is used to adjust the pump and probe laser beams temporally. It has an optional advantage to increase the resolution. THz signal is detected by digital signal processor built in a central controller which is connected with slow delay line, fast delay line, and THz emitter while detector is connected via operational preamplifier. Central controller displays the THz signal on personal computer in time domain. Average of 1024 signals was taken for each measurement to get a very high signal to noise ratio of more than 105 : 1. The whole setup is inside a box whose relative humidity was dropped below 5% by nonheating sorbent dryer, and temperature was maintained at 23°C. THz system is shown in Figure 2. RamanStation 400F  and Spectrum 100 FTIR [18, 19] from PerkinElmer are used for Raman and FTIR spectroscopy, respectively.
3. Sample Preparation
“Acrylamide” CAS no. powder is bought from “Sinopharm Chemical Reagent Co. Ltd”  and is ground to very fine powder. This powder is then pressed homogeneously between two plates of polystyrene and sealed from the sides. Refractive index and absorption coefficient of polystyrene are determined first with air as reference. Equation (14) was used for refractive index and (15) for absorption coefficient by setting . A constant refractive index of 1.58 and nearly zero absorption coefficient of polystyrene (0.2–2 THz) are calculated which is consistent with the literature . Thickness is determined by a digital micrometer capable of measuring the thickness up to 0.001 mm precision. Many samples are prepared with varying thicknesses with one reference. Reference cell is prepared by joining two polystyrene plates as shown in Figure 3.
4. Data Analysis
We derive the formulae for refractive index and absorption coefficient as follows. We know that absorption coefficient and phase can be represented as where is the imaginary part and is the complex refractive index of the material under consideration, is angular frequency, is the sample thickness, and is speed of light whose value is taken as 299 792 458 m/s.
We represent the electric fields as follows: Here, represents the complex transmission Fresnel’s coefficients with subscripts being the interface of two media and is the phase accumulated while passing through the media. is the phase accumulated while passing through the air through a distance equal to the sample width where the value of can be determined using (2). is the incident electric field: where and are representing two media and is the interface of two media. Dividing (4) by (3), Putting the values of and using (2) and adding appropriate subscripts,
As we assume that there is no absorption in air, so As is related to absorption coefficient “”, so by using (1) and adding appropriate subscripts,
Entering the values of and by using (5),
As polystyrene (PS) has no absorption and for the sample, we also assume that absorption is negligible. So,
Comparing the above result with Euler’s form of the equation, We can evaluate the values of and from (13):
Sample amplitude , reference amplitude , sample phase , and reference phase are put from experimental data after performing fast Fourier transform (FFT) on time domain traces in the previous formulae. is angular frequency, is speed of light, is sample thickness, and is the refractive index of polystyrene whose value was taken as 1.58. Error bars are showing 95% confidence interval calculated from the standard error of the mean. Error bars are almost constant for the refractive index until 1.8 THz and minimum for absorption coefficient in the range 0.2 to 1.2 THz. Error bars are large in the high frequency range because of low signal to noise ratio caused by weak input signal and large absorption in high frequency region. Time domain reference and sample pulses along with their FFTs are shown in Figure 4. Terahertz time domain signals were padded with zeros from 21 ps before performing FFT.
5. DFT Calculations
DFT calculations are carried out in Gaussian 09 package  with Ground State method B3LYP along with 6-311G basis set  together with (3df, 3pd) and ++ diffuse functions . Raman activities are calculated at incidence frequency of 785 nm which is the frequency of excitation source of our Raman spectrometer . Optimization and frequency calculations are done, and vibrational spectra are compared with experimental results. On the basis of these DFT calculations, experimental peaks are assigned to relevant molecular vibrations. List of all modes along with FTIR, Raman, and potential energy distribution (PED) analysis is given in Table 1.
|Ob: out of plane bending, t: torsion, Unscaled: raw DFT values without any scaling, Scaled: scaled with MOLVIB. |
Note: animation of all of these vibrational modes can be found in the Supplementary Material available online at http://dx.doi.org/10.1155/2013/148903.
6. Results and Discussion
As shown in Figure 5, refractive index of acrylamide is almost constant from 0.2 to 2 THz with an average value of 1.33, and absorption coefficient is also constant in 0.2–1.2 THz range after which it slowly starts rising up without any significant peak in this region. This agrees well with DFT calculation which predicts no intramolecular vibrational mode in this range. According to DFT calculations, first vibrational mode appears at 2.786 THz (92.94 cm−1 Unscaled in Table 1) which is obviously beyond the range of our THz spectrometer. Spectra of both FTIR and DFT are shown in Figure 6. As 3100 to 3000 cm−1 region is characterized by carbon-hydrogen (C-H) stretching , hence, the peaks at 2813 and 3034 and 3102 cm−1 are assigned to C-H stretching as shown in Figure 6. Two absorption lines at 3191 and 3347 cm−1 are due to hydrogen-nitrogen stretching as revealed by DFT calculations. After this region, vibrational motion is complex and multiple bonds contribute to the overall vibrational motion of the molecule as shown in the PED analysis in Table 1. The vibrational frequencies in the region 1400–1625 cm−1 are usually assigned to C=C and C-C stretching vibrations , and with the aid of PED analysis, it is clearly revealed that the vibrations in this range in FTIR and Raman experimental spectra are partially contributed by C=C and C-C interactions although there is contribution from other bonds too. Peaks are reasonably in good agreement with DFT values in the region (450–1800 cm−1) both for FTIR and Raman after scaling with the scale factor of 0.9679 for vibrational frequencies calculated by Andersson and Uvdal . However, DFT values are deviating over a large extent from experimental values in the C-H and N-H stretching regions of FTIR even after scaling because of the involvement of intramolecular hydrogen bonding. We use MOLVIB [27, 28] to further scale these values and are given in Table 1. PED analysis is also done by MOLVIB. Raman and DFT spectra are shown in Figure 7. Experimental Raman spectrum actually consists of 95–3000 cm−1 range, but because there was no peak after 1800 cm−1 in the experimental Raman spectrum, so graph in Figure 7 is cut off at 2000 cm−1.
We have calculated (0.2–2 THz) the relevantly constant refractive index of acrylamide with average value of 1.33, and absorption coefficient is also constant initially from 0.2 to 1.2 THz, and after that it starts rising up exponentially. No absorption peak was observed in the THz range which is also predicted by DFT calculations according to which first intramolecular vibrational mode appears at 2.78634 THz which is beyond the range of our spectrometer. FTIR (450–4000 cm−1) and Raman (95–3000 cm−1) spectra have been presented with different absorption lines and their association with corresponding vibrational modes on the basis of DFT calculations. Scaling was done by MOLVIB to match with the experimental values especially for C-H and N-H stretching regions.
The authors would like to extend their sincere appreciation to the Deanship of Scientific Research at King Saud University for its funding of this research through the Research Group Project no. RGP-VPP-255.
All the 24 vibrational modes of Acrylamide molecule as calculated by DFT B3LYP 6-311G++ (3df 3pd) are given in conjunction with the research paper. There is one animated .gif file for each vibrational mode with the name of the file being the frequency of the respective vibrational mode in cm-1 as given in Table 1 of the research article.
- Y.-S. Lee, Principles of Terahertz Science and Technology, Springer Science+Business Media, 2009.
- D. Grischkowsky, S. Keiding, M. van Exter, and Ch. Fattinger, “Far-infrared time-domain spectroscopy with terahertz beams of dielectrics and semiconductors,” Journal of the Optical Society of America B, vol. 7, no. 10, pp. 2006–2015, 1990.
- R. Piesiewicz, C. Jansen, S. Wietzke, D. Mittleman, M. Koch, and T. Kürner, “Properties of building and plastic materials in the THz range,” International Journal of Infrared and Millimeter Waves, vol. 28, no. 5, pp. 363–371, 2007.
- A. Katsounaros, M. Mann, M. Naftaly, K. Z. Rajab, Y. Hao, and W. I. Milne, “Terahertz time-domain spectroscopy characterization of vertically aligned carbon nanotube films,” Carbon, vol. 50, no. 3, pp. 939–942, 2012.
- B. Yu, F. Zeng, Y. Yang et al., “Torsional vibrational modes of tryptophan studied by terahertz time-domain spectroscopy,” Biophysical Journal, vol. 86, no. 3, pp. 1649–1654, 2004.
- J. S. Melinger, S. S. Harsha, N. Laman, and D. Grischkowsky, “Temperature dependent characterization of terahertz vibrations of explosives and related threat materials,” Optics Express, vol. 18, no. 26, pp. 27238–27250, 2010.
- M. C. Kemp, P. F. Taday, B. E. Cole, J. A. Cluff, A. J. Fitzgerald, and W. R. Tribe, “Security applications of terahertz technology,” in Terahertz for Military and Security Applications, vol. 5070 of Proceedings of SPIE, pp. 44–52, April 2003.
- J. F. Federici, B. Schulkin, F. Huang et al., “THz imaging and sensing for security applications—explosives, weapons and drugs,” Semiconductor Science and Technology, vol. 20, no. 7, pp. S266–S280, 2005.
- W. H. Fan, A. Burnett, P. C. Upadhya, J. Cunningham, E. H. Linfield, and A. G. Davies, “Far-infrared spectroscopic characterization of explosives for security applications using broadband terahertz time-domain spectroscopy,” Applied Spectroscopy, vol. 61, no. 6, pp. 638–643, 2007.
- J. Chen, Y. Chen, H. Zhao, G. J. Bastiaans, and X.-C. Zhang, “Absorption coefficients of selected explosives and related compounds in the range of 0.1-2.8 THz,” Optics Express, vol. 15, no. 19, pp. 12060–12067, 2007.
- A. Rahman, “Dendrimer based terahertz time-domain spectroscopy and applications in molecular characterization,” Journal of Molecular Structure, vol. 1006, no. 1–3, pp. 59–65, 2011.
- R. Piesiewicz, T. Kleine-Ostmann, N. Krumbholz et al., “Short-range ultra-broadband terahertz communications: concepts and perspectives,” IEEE Antennas and Propagation Magazine, vol. 49, no. 6, pp. 24–39, 2007.
- K. M. Wilson, L. A. Mucci, B. A. Rosner, and W. C. Willett, “A prospective study on dietary acrylamide intake and the risk for breast, endometrial, and ovarian cancers,” Cancer Epidemiology Biomarkers and Prevention, vol. 19, no. 10, pp. 2503–2515, 2010.
- P. D. Cunningham, N. N. Valdes, F. A. Vallejo et al., “Broadband terahertz characterization of the refractive index and absorption of some important polymeric and organic electro-optic materials,” Journal of Applied Physics, vol. 109, no. 4, Article ID 043505, 2011.
- M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 09, Revision A. 01, Gaussian, Wallingford, Conn, USA, 2009.
- V. 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.
- T. Clark, J. Chandrasekhar, G. W. Spitznagel, and P. v. R. 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.
- M. Arivazhagan and R. Meenakshi, “Vibrational spectroscopic studies and DFT calculations of 4-bromo-o-xylene,” Spectrochimica Acta Part A, vol. 91, pp. 419–430, 2012.
- M. P. Andersson and P. Uvdal, “New scale factors for harmonic vibrational frequencies using the B3LYP density functional method with the triple-zeta basis set 6-311+G(d,p),” Journal of Physical Chemistry A, vol. 109, no. 12, pp. 2937–2941, 2005.
- T. Sundius, “Molvib—a flexible program for force field calculations,” Journal of Molecular Structure, vol. 218, pp. 321–326, 1990.
- T. Sundius, “Scaling of ab initio force fields by MOLVIB,” Vibrational Spectroscopy, vol. 29, no. 1-2, pp. 89–95, 2002.
Copyright © 2013 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.