Abstract

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.

1. Introduction

A strip between the microwave and infrared region of electromagnetic spectrum is named terahertz band spanning from 0.1 to 10 THz [1]. Terahertz time domain transmission spectroscopy has been applied tremendously in the fields of materials characterization [25], defense and security [611], telecommunication [12], 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) [13]. Moreover, high intake of acrylamide increases the risk of endometrial cancer and possibly ovarian cancer in women [14]. 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) [15] 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) [16]. 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 [17] 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” [20] 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 [21]. 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 [22] with Ground State method B3LYP along with 6-311G basis set [23] together with (3df, 3pd) and ++ diffuse functions [24]. Raman activities are calculated at incidence frequency of 785 nm which is the frequency of excitation source of our Raman spectrometer [17]. 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.

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 [25], 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 [25], 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 [26]. 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.

7. Conclusion

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.

Acknowledgment

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.

Supplementary Materials

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.

  1. Supplementary Material