International Journal of Spectroscopy

International Journal of Spectroscopy / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 841593 |

Diwaker, Abhishek Kumar Gupta, "Quantum Chemical and Spectroscopic Investigations of (Ethyl 4 hydroxy-3-((E)-(pyren-1-ylimino)methyl)benzoate) by DFT Method", International Journal of Spectroscopy, vol. 2014, Article ID 841593, 15 pages, 2014.

Quantum Chemical and Spectroscopic Investigations of (Ethyl 4 hydroxy-3-((E)-(pyren-1-ylimino)methyl)benzoate) by DFT Method

Academic Editor: Hakan Arslan
Received17 Feb 2014
Revised02 Apr 2014
Accepted05 Apr 2014
Published09 Jul 2014


In the present work we have reported the optimized ground state geometry, harmonic vibrational frequencies, NMR chemical shifts, NBO analysis, and molecular electrostatic potential surface map of the title compound using DFT/B3LYP/6-311++G(2d, 2p) level of theory. We have compared our calculated results with the experimentally obtained values and found that both are in close agreement with each other. We have used the gauge-invariant atomic orbital (GIAO) approach to calculate the NMR (13C and 1H) chemical shifts using Gaussian 09 package. TD-DFT (time-dependent DFT) approach has been used to simulate the electronic spectra of the title compound in order to account for excited states. Other molecular properties such as HOMO-LUMO energies, NBO analysis, and PED distribution analysis have been studied and reported using DFT/B3LYP/6-311++G(2d, 2p) level of theory.

1. Introduction

The title compound chosen for DFT studies to extract different molecular properties has been experimentally synthesized and prepared using 1 amino pyrene and (ethyl 3-formyl-4-hydroxybenzoate) at room temperature for six hours in the presence of dry MeOH [1]. The title compound shows sensing properties for selective detection of niobium ions in mixed aqueous media. In the literature survey we found that there are lot of research articles based upon fluorescent techniques for detection of various metal ions, however literature survey also reveals very few ab initio HF/MP2/DFT calculations of such type of compounds. Fluorescence is a very simple technique and acts as a convenient characterization tool for detection of very small amount (in ppm) of various metal ions in solutions [2]. From application point of view niobium metal is used in various kinds of applications such as superconducting magnets [3] and biological applications [4]. In interest of such applications the quantum mechanical calculations of the title compound are thoroughly investigated. The aim of this work is to predict the structural, electronic, vibrational, and spectral parameters and other molecular properties of the title compound using DFT approach [59].

2. Experimental Details

In this section we have reported short details about the methodology and characterization tools used for the title compound, however we advised the readers to consult [1] for more details. The chemical structure and fluorescent properties of the title compound are confirmed by single crystal X-ray diffraction, UV spectra, 1H and 13C NMR, and FTIR spectra. 1H and 13C NMR are recorded in chloroform using TMS as internal standard on a Varian Mercury 300 spectrometer operating at 300 MHz for 1H and 75 MHz for 13C. IR spectra are recorded on a Perkin-Elmer PE-983 infrared spectrometer as KBr pellets with absorption reported in cm−1. The ultraviolet absorption spectra were recorded on Shimadzu UV-2450 spectrophotometer. Fluorescent spectra measurements were performed on Agilent Technologies Cary Eclipse fluorescence spectrometer.

3. Computational Details

Using DFT/B3LYP/6-311++G(2d, 2p) level of theory [10] we have investigated the ground state optimized geometry of the title compound. The molecular geometry is fully optimized using tight convergence criteria along with redundant internal coordinates and Berny’s optimization algorithm. The optimized parameters obtained using DFT approach have been compared with the experimental values and are in close agreement with them. Further we have used the optimized ground state geometry of the title compound to study the different properties like NMR spectra, UV-Vis spectra, MEP surface mapping, PED analysis, and NBO analysis. Using DFT/B3LYP/6-311++G(2d, 2p) level of theory and GAIO (Gauge-Invariant Atomic Orbital) [11, 12] approach we have reported the NMR (13C and 1H) chemical shifts of the title compound and compared them with their experimental counterparts. To study the electronic transitions and excited states we have used the TD-DFT (time-dependent) method available in Gaussian 09 package. HOMO-LUMO energies are also calculated at DFT/B3LYP/6-311++G(2d, 2p) level of theory. Vibrational wavenumbers assignment is done by using VEDA 4 program. MEP surface mapping is investigated to comment upon the reactive nature of the title compound. In order to find out the various interactions between the filled and the vacant orbitals, NBO analysis [13] of the title compound has been done using NBO 3.1 program available in Gaussian 09 package at DFT/6-311++G(2d, 2p) level of theory. The unoptimized structure of our title compound is presented in Figure 1.

4. Results and Discussion

4.1. Molecular Geometry

We have used the DFT/B3LYP/6-311++G(2d, 2p) level of theory available in Gaussian 09 program to investigate the ground state geometry of the title compound. The geometry is fully optimized with tight convergence criteria and the structure is local minima on the PES. On comparison with the experimentally obtained parameters one can conclude that our calculation is successful, as the difference between calculated and experimental bond lengths, bond angles is of few Å. Figure 2 represents the stable conformation of the title compound using DFT calculations. The selected calculated bond lengths and angles (Å) for the title compound along with their corresponding experimental values are listed in Table 1.

Serial number (bond length/angstrom)CalculatedExperimental

1 (3, 12)1.361.35
2 (7, 8)1.241.23
3 (13, 14)1.301.29
4 (9, 7)1.371.34
5 (9, 10)1.481.45
6 (14, 15)1.411.41

Serial number (angle/degrees)CalculatedExperimental

1 (13-14-15)121.30120.98
2 (12-3-4)121.10120.37
3 (8-7-6)124.98124.94
4 (9-10-11)106.83106.84
5 (7-9-10)117.12117.78
6 (12-3-2)118.82118.86
7 (8-7-9)122.09123.01
8 (9-7-6)112.92112.95

: Bond lengths (in Angstrom); : Angles (in Degrees).

Correlation between [14] the calculated and the experimental parameters of bond lengths and bond parameters for the title compound are shown in Figure 3. Bond length and bond angles correlation values are 0.9802 and 0.9921, respectively.

4.2. Chemical Shifts

NMR spectroscopy is considered as a valuable tool for the structural and functional characterization of molecules. 1H and 13C NMR chemical shifts of the title compound are investigated using DFT/B3LYP/6-311++G(2d, 2p) level of theory with GIAO (gauge-invariant atomic orbital) approach in DMSO. The calculated 1H and 13C NMR chemical shifts of the title compound together with the corresponding experimental values are shown in Tables 2 and 3 as values relative to tetramethylsilane. The 1H NMR spectra of the title compound in DMSO show a triplet peak in the range from 1.36 to 1.40 ppm for C–CH3 and quartet peak in the range from 4.32 to 4.38 ppm for O–CH2. Corresponding calculated values are in the range from 1.31 to 1.46 ppm and 4.27 to 4.32 ppm, respectively. Aromatic protons in pyrene appeared in the range from 8.04 to 8.48 ppm with calculated values ranging from 8.32 to 9.03 ppm. Aromatic protons in substituted benzene showed a singlet peak at 8.66 ppm and doublet at 8.50 ppm and 7.16 ppm. Corresponding calculated values are 8.66 ppm, 8.50 ppm, and 7.44 ppm, respectively. Proton of Schiff base shows a singlet peak at 9.37 ppm while calculated value is at 9.34 ppm. The 13C NMR spectra of the title compound in DMSO showed a peak at 14.44 ppm for C–CH3 (calculated value as 13.48 ppm), at 60.99 ppm for O–CH2 (calculated value as 60.56 ppm), at 165.88 ppm for C=O (calculated value as 180.00 ppm), and at 162.81 ppm for C=N (calculated value as 169.94 ppm). Aromatic carbons of pyrene and substituted benzene showed peaks in the range from 115.67 to 165.12 ppm while calculated values are in the range from 121.53 to 175.80 ppm. We have also reported the hydrogen and carbon NMR chemical shifts by IGAIM and CGST methods using B3LYP/6-311++G(2d, 2p) basis sets [15]. From the comparison table we can conclude that theoretical values for carbon and hydrogen NMR chemical shifts calculated by GIAO method are in close agreement as compared to other CGCST and IGAIM methods. Linear correlation coefficients for linear regression analysis of theoretical and experimental values of 1H and 13C NMR chemical shifts using GIAO method are found to be 0.9959 and 0.9958 respectively.

AtomExperimental value6311++G6311++G6-311++G
(2d, 2p)/GIAO (2d, 2p)/CGCST

H3914.15(s, 1H, Ar OH)15.2215.6915.68
H409.37(s, 1H, N=CH)9.349.569.56
H41, H42, H43, H448.04–8.48(9H, Ar CH, Pyrene)8.32–9.038.46–9.468.45–9.46
H45, H46, H47, H48
H318.50(d, 1H, Ar CH)8.508.778.77
H327.16(d, 1H, Ar CH)7.347.707.70
H338.66(s, 1H, Ar CH)8.668.998.99
H34, H354.32–4.38 (s, 2H, O–CH2)4.27–4.324.75–4.824.74–4.81
H36, H37, H381.36–1.40(T, 3H, C–CH3)1.31–1.461.60–1.891.60–1.89

(2d, 2p)/GIAO (2d, 2p)/CGST

C15–C30117.52–141.50(16C, Ar C, pyrene)115.60–148.50118.77–167.66118.76–167.66
C1–C6115.67–165.12(6C, Ar C, substituted benzene)121.53–175.80120.60–176.80120.59–176.80
C7165.88(1C, Ar C=O)180.00179.0179.07
C13162.81(1C, Ar C=N)169.94167.66167.66
C1060.98(1C, OCH2)60.5668.6668.66
C1114.44(1C, CH3)13.4813.5113.50

4.3. Frontier Molecular Orbitals

Analysis of the HOMO-LUMO band gap helps us in understanding many molecular properties of a molecule like chemical reactivity, UV-Vis spectra, and stability of the molecule [16] along with optical and electrical properties. Chemical reactivity of a molecule can be determined from the HOMO-LUMO band gap. A small band gap implies low kinetic stability of the molecule. HOMO-LUMO separation is a result of significant degree of intermolecular charge transfer from the electron donor groups to the electron acceptor groups through conjugated paths. Energy gap between HOMO and LUMO has also been used to prove the bioactivity from intramolecular charge transfer (ICT). We have reported the HOMO-LUMO analysis of the title compound using DFT/B3LYP/6-311++G(2d, 2p) level of theory [17]. In our analysis we found that the title compound has a total of 986 orbitals out of which 103 are occupied and the remaining 883 are virtual orbitals. The orbitals numbered as 103 and 104 account for HOMO and LUMO orbitals. The HOMO-LUMO energies of the title compound have also been calculated using ab initio calculations and are found to be −5.63 electron volts and −2.50 electron volts, respectively. The HOMO-LUMO for the title compound has been shown in Figure 4 and the gap is found to be 3.13 electron volts. The HOMO LUMO distribution is mostly localized on the rings which show that they are type orbitals. transition implies an ED transfer between rings transition. From this value of band gap we can predict that the title compound can be used for organic solar cell applications, title compound has high kinetic susceptibility and low chemical reactivity. Using HOMO and LUMO energies along with equations as which is electronegativity, as chemical hardness with as chemical softness has been calculated for the title compound. The terms and are equivalent to and and are referred to as ionization potential and electron affinity, respectively. In addition to HOMO/LUMO energies the HOMO−1/LUMO+1 energies of the title compound have been calculated using B3LYP/6-311++G(2d, 2p) level of theory and are found to be −5.61 eV and −2.41 eV, respectively. Electron donating and electron withdrawing ability of the title compound are expresses in terms of , and and come out to be 4.065, 1.565, and 0.3194 [18].

4.4. MEP Surface Mapping

We have reported and plotted the MEP surface mapping, alpha density, and total density of the title compound using Gaussian 09 program. The molecular electrostatic potential surface along with Alpha density and total density for the title compound is represented in Figure 5.

MEP surface mapping is useful in understanding hydrogen bonding interactions as well as sites for electrophilic and nucleophilic attacks [19, 20]. The MEP surface provides us with net electrostatic effect caused due to total charge distribution. It can be considered as a fruitful quantity to understand the various molecular properties like hydrogen bonding and reactivity. It also provides a useful tool to know the relative polarity of the molecule [21]. Portion of the molecule which has –ve electrostatic potential will be susceptible to electrophilic attack. The surface is color coded as per the electrostatic potential (red is more electron rich area and blue is more electron poor area.). The total electron density plot of the title compound shows a uniform distribution. The order in the increase of the electrostatic potential as per color code will follow as red < orange < yellow < green < blue [22]. At last we conclude that the investigated molecule has several sites for electrophilic as well as nucleophilic attacks as shown in MEP surface mapping.

4.5. UV-Vis Studies and Electronic Properties

To find the electronic absorption spectrum including singlet and triplet states of the title compound the calculations were performed on fully optimized ground state geometry using DFT/B3LYP/6-31++G(2d, 2p) level of theory. THF is used as a solvent to simulate the electronic absorption. Figure 6 represents the computed electronic spectra of the title compound. The electronic spectra are recorded within a range of 200 nm–800 nm. Using TDDFT theory the oscillator strength along with excitation energy for the triplet and the singlet states has also been calculated. The different values for excitation energy along with oscillator strength as well as CI expansion coefficients are listed in Table 4. For the title compound the maximum absorption value obtained using TD-DFT/B3LYP/6-311++G(2d, 2p) basis set are 485 nm, 332 nm, and 285 nm, respectively, with THF as solvent in CPCM model. Corresponding experimental values as reported are 383 nm and 258 nm, respectively. The calculated band at 485 nm is intense and accounts for a type of transition. absorption band in the calculated spectrum indicates a transition and is close to experimentally calculated values.

Excited state (triplet/singlet) CI expansion coefficientEnergy (eV)Wavelength calc. (nm)Oscillator Strength (f)

5(T)−0.258293.2510381.37 (383)0.0000
20(S)−0.300754.3829282.88 (238)0.3427

4.6. Vibrational Spectra

In the present study we have reported the molecular vibrations of the title compound by means of FTIR spectroscopy. Our title compound is asymmetric top with C1-symmetry and is characterized by 141 normal modes of vibration. We have used DFT/B3LYP/6-311++G(2d, 2p) [23] method to investigate the normal modes of vibration of our title compound. The main reason for selecting this computational scheme is that it reproduces experimental frequencies with high accuracy and the same can be predicted from the comparison of the calculated values with the experimental ones. The calculated and experimental FTIR spectra of the title compound are shown in Figure 6. On comparison we found that the calculated values using the above method are found to be in close agreement with the experimental values. Calculated C–H stretching vibrations of aromatic rings appeared in the wavenumber range 2800–3200 cm−1. The same has been confirmed with the experimental IR where the wavenumber range for aromatic rings ranges from 3000 to 3200 cm−1. The bands observed in the wavenumber range from 3250 to 2850 cm−1 in the calculated IR spectra of the title compound are assigned to the alkyl C–H stretching vibrations and the same is confirmed with the experimental values. C=O (ester) stretching vibrations are predicted at 1649 cm−1 while for the same functional group experimental values are at 1711 cm−1. C=N (Schiff base) stretching vibrations are predicted at 1611 cm−1, while experimental values are at 1610 cm−1. We have also analyzed and reported our modes of vibrations in terms of PED. PED analysis is done by using VEDA 4 program [24]. This program generally uses the Gaussian output file in formatted checkpoint form as its input files for PED analysis. These input files contain information about orientation of coordinates, force constants (F-matrix), and frequencies with atom displacement matrix. The information on F-matrix must start form the line “Force Constants in Cartesian coordinates” (Figure 7).

We have repeated our PED analysis few hundred times to achieve maximum value of PED contributions. In PED interpretation each fundamental normal mode coordinate is expressed in terms of internal mode coordinates which is a combination of stretchings, bendings, or torsions. This transformation basically results in the nondiagonality of the force constant matrix, which means that PED contributions of different modes are mutually related to each other by nondiagonal terms. Further we explain how this procedure works as a normal mode coordinate is replaced by an internal set of coordinates and PEDs are calculated. A parameter EPm is used to express the maximum PEDs and is basically considered as optimization of the PED analysis. If our title compound consists of a large number of modes, then it will result in an increase of optimization time. Theoretically calculated and experimental wavenumbers (available) are summarized in Table 5. Detailed vibrational assignments, IR intensities, and computed wavenumbers along with the percentage of PED are given in Table 5. The spectra were analyzed in terms of the PED contributions by using the VEDA program.

Experimental wavenumbers (cm−1)Calculated wavenumbers (cm−1) (cm−1)Raman (cm−1)PEDInterpretation

3224.314.48183.79s11 93 (CH2), ring 1
3212.19.1987.51s13 94 (CH2), ring 4
32003205.532.772.41s3 94 (CH2), ring 1
3191.4560.42530.75s19 97 (CH3), ring 5
3189.2827.8114.52s14 92 (CH2), ring 4
3183.941.5736.38s4 99 (CH), ring 1
31813181.1551.58278.89s18 96 (CH2), ring 3
3171.6524.78229.97s12 −13 s15 −83 (CH2), ring 5 + (CH), ring 4
3167.972.7545.05s12 78 s15 −15 (CH2), ring 4 + (CH2), ring 5
3167.2411.8860.89s7 92 (CH2), ring 2
3161.690.2447.53s17 93 (CH3), ring 5
3158.711.7440.37s16 93 (CH2), ring 3
3119.248.5916.84s5 −13 s6 −10 s8 −17 s10 60Not defined
3099.9138.21156.66s6 −81 (CAH2)
30903090.570.9395.26S8 78 s10 −12 (CAH2)+ (CH2)
3059.2120.4988.61s2 99 (CAH)
3052.3815.95154.12s9 96 (CAH2)
3029.9926.36191.4s5 76 s10 21 (CAH3)
2894.86785.16225.28s1 96 (CAO)
1675.14219.1561.69s26 −60 (CAO)
17111667.7912.271522.05s39 42 (C–C), ring 3, 4, 5
1648.5371.86211.55s20 32 s23 −10 (CAO) + (C–C), ring 1, 4
1643.3351.26229.63s20 −16 s23 −28 (CAO) + (C–C), ring 1, 4
1638.2339.051129.96s23 −16 s28 −22 (C–C), ring 3 + (C–C), ring 4
1629.6122.253026.04s21 11 s28 −23 (C–C), ring 1 + (C–C), ring 5
1622.0716.12824.95s27 13 s81 11 (C–C), ring 4 + CCC, ring 5
16101610.778.328784.86s35 27 (CAN)
1583.638.54225.08s36 21 s39 −10 s60 11 (C–C), ring 2 + (C–C), ring 3 + HCC, ring 3
1554.366.423.52s62 70 s101 15 HCC + HCOC, (OCA)
1546.9913.71548.31s33 16 s68 −24 (C–C), ring 3 + HCC, ring 2
1539.548.0729.88s61 70 s62 −10 s103 10 HCH,0020(HCA) + HCH, (HCA) + HCOC, (OCA)
1526.219.123.7s59 73 s100 −12 HCH, (HCA) + HCOC, (OCA)
1524.2337.79260.32s25 12 s72 −10 (C=C), ring 5 + HCH, ring 1
1518.0871.682628.38s25 11 s55 −14 s56 −10 (C=C), ring 2 + HCC, ring 1
1495.7521.77369.81s67 −25 HCC, ring 3
1493.0789.532975.2s35 −13 s54 17 (CAN) + HOC
1476.414.864.47s60 22 HOC
1465.713.62246.14s71 −65 HCC, ring 5
1457.6612.615.74s63 84 HCC, (HCA)
1451.136.8516.4s54 −11 HOC, (OCA)
1444.016.01106.44s29 −15 s69 10 (C–C), ring 5 + HCC, ring 3
1422.523.561171.11s22 34 (C=C), ring 2
1418.968.01160.33s62 11 s101 67 HCH, (HCA) + HCOC, (OCA)
1407.1311.57189.62s31 −23 s64 −17 (C–C), ring 5 + HCN, (NCA)
1387.024.94692.54s31 −11 s64 31 (C–C), ring 5 + HCN, (NCA)
1371.353.01744.64s24 −51 s42 13 (C–C), ring 1 + (C–C), ring 5
1365.373.1720.73s85 10 CCC, ring 4
1347.538.3230.53s37 −41 s66 −10 (C=C), ring 2 + HCC, ring 4
1322.4132.6186.95s56 −45 HCC, ring 1
1301.850.3218.48s58 76 HCC, (CAH)
1293.461.1619.8s33 12 s60 10 s68 21 (C–C), ring 3 + HCC, ring 4 + HCC, ring 5
1290.84115.58764.31s30 58 (CO), (CAO)
1284.222.96455.48s41 20 s69 −17 (C–C), ring 5 + HCC, ring 5
1264.63151.2724.03s26 −16 s57 −14 (C–C), ring 5 + HCC, ring 1
1263.1562.223895.41s41 −11 s44 −11+ (C–C), ring 1 + (C–C), ring 5
1252.17107.36376.98Not defined
1231.74194.37485.07s32 −10 s72 21 (C–C), ring 2 + HCC, ring 5
1223.49154.6440.05s70 −52 HCC, ring 2
1211.31266.52575.9s70 −13 HCC, ring 2
1196.5852.3947.76s25 17 s66 40 (C=C), ring 4 + HCC, ring 5
1190.050.9829.65s39 14s65 49 (C=C), ring 3 + HCC, ring 3
1188.465.774.66s58 −14 s100 54 HCC (CAH) + HCOC, (OCA)
1163.1162.34.06s29 25s57 55 (C=C), ring 1 + HCC, ring 1
1147.163.0519.94s48 22 s61 −18 s94 −12 s103 −32 HCH (CAH) + CCO (CAO) + HCCO (OCA)
1143.320.5114.59s77 32 CCC, ring 4
1132.5518.83225.65s69 11 HCC, ring 5
1118.58268.594.61s46 −31 s55 15 (C=C), ring 1 + HCC, ring 1
1085.9316.5930.98s27 14 s34 22 (C=C), ring 3 + (C=C), ring 5
1039.7188.124.02s47 72 (C=C), (CAH)
1037.692.76.89s102 −80 HCCC, ring 4
1015.71119.59s112 −52 HCCC, ring 5
1014.4316.39160.22s99 −66 HCNC, (CAN)
1005.480.180.64s113 61 HCCC, ring 3
1004.515.0579.88s106 74 HCCC, ring 1
995.676.44.17s32 10 s40 −12 (C–C), ring 2 + (C=C), ring 3
982.591.2130.45s109 66 HCCC, ring 2
966.875.99189.31s97 38 HCCC, ring 1
964.3622.393.01s97 44 HCCC, ring 1
957.2780.571.93s96 75 HOCC (CAO)
938.591.654.09s110 72 HCCC, ring 5
925.2533.74166.09s44 −11 s79 −27 (CAN) + CCC, ring 4
881.43111.218.99s108 −60 HCCC, ring 2
876.5942.0310.19s48 −12 s107 25 (C–C), ring 1 + HCCC, ring 4
867.6138.532.37s98 60 HCCC, ring 1
866.725.462.94s48 −16 s98 −21 s107 −12 s108 −14Not defined
846.4110.782.34s50 −54 CCC, ring 4
833.533.46.97s107 14 s111 31 HCCC, ring 1 + HCCC, ring 3
830.770.382.5s104 70 HCCO, (CAO)
828.699.2166.05s111 21 HCCC, ring 2
815.143.697.4s128 −18 CCCC, ring 5
804.431.774.67s73 −27 OCO, (CAO)
793.291.0859.46s43 20 s82 18 (C–C), ring 1 + CCC, ring 4
778.978.824.11s114 67 HCCC, ring 5
764.695.889.65s82 17 CCC, ring 1
760.3419.230.69s136 −51 CCCC, ring 1
741.7817.230.25s107 −12 s111 11 s116 −11 HCCC, ring 1 + HCCC, ring 4 + CCCC, ring 5
724.926.876.4s141 55 OCOC, (CAO)
714.440.1918.14s28 −11 s49 16 (C–C), ring 3 + CCC, ring 4
702.648.9232.91s120 −15 s138 −10 s140 −24 CCCC, ring 1 + CCCC, ring 2 + CCCC, ring 3
676.9730.4711.99s41 11 s93 14 (C–C), ring 1 + CCO, (CAO)
651.4514.7211.48s51 −35 CCC, ring 1
640.9713.6140.41s126 12 s129 −11 CCNC, ring 1 + CCCC, ring 1
622.12.328.46s38 −13 s80 17 (C–C), ring 3 + CCN, (CAN)
578.440.119.98s41 10 s49 10 (C–C), ring 3 + CCC, ring 5
568.572.7512.66s121 −11 s122 −12 s124 −11s139 15Not defined
559.840.7113.11s122 23 s124 10 s134 19 CCCC, ring 1 + CCCC, ring 3 + CCCC, ring 5
546.1711.742.14Not defined
537.813.676.23s52 24 CCO, (CAO)
530.021.0612.8s117 24 s139 13 CCCN, (CAN) + CCCC, ring 3
521.453.0574.21s125 22 s127 14 CCCC, ring 3 + CCNC, (CAN)
512.460.831.44s88 23 CCC, ring 3
497.393.1911.62s76 −14 s83 11 s85 −15 CCC, ring 3 + CCN (CAN) + CCC, ring 5
483.982.272.37s123 −24 s135 24 CCCC, ring 1 + CCCC, ring 2
460.51.96.03s87 20 CCC, ring 4
444.531.193.6s124 10 CCCC, ring 1
432.840.6821.66s91 41 CCO (CAO)
415.62.4689.07s22 −11 s90 −14 (C–C), ring 4 + CCC, ring 3
386.890.51.35s86 −10 CCC, ring 2
372.3211.471.34s75 −16s94 56 COC, (CAO) + CCO, (CAO)
355.131.384.26s123 19 s135 10 CCCC, ring 3 + CCCC, ring 5
346.913.473.7s115 −11 NCCC, (CAN)
311.174.6420.43s45 −15 s95 31 (C–C), ring 4 + COC, (CAO)
289.375.168.61s115 24 NCCC, (CAN)
2730.462.42s105 −10 s121 17 s139 15 CCCO, (CAO) + CCCC, ring 3 + CCCC, ring 5
267.890.587.14s75 13 COC, (CAO)
252.251.2310.57s103 12s105 55 CCCO, (CAO)
241.042.819.87s75 12 COC, (CAO)
223.338.2610.27s138 −43 CCCC, ring 4
186.815.987.28s89 12s116 −17 CCC, (CA) + CCCN, (CAN)
171.60.3410.17Not defined 24 s133 37 CCCC, ring 1 + COCC, (CAO)
134.243.0426.07s75 −10 COC (CAO)
122.640.0811.46s119 15 s129 −11 CCCN, (CAN) + CCNC, (CAN)
100.651.286.44s119 10 s127 14 CCCN, (CAN) + CCNC, (CAN)
92.610.677.41s118 11 s131 −15 s133 −20 s134 13 COCC, (CAO)
79.90.643.96s92 23 s131 12 CCC, ring 4 + COCC, (CAO)
59.620.882.9s130 −16 s131 −37 COCC, (CAO)
39.590.122.58s130 48 s131 13 s132 −11 COCC, (CAO)
31.930.145.25s74 −42 CNC (CAN) 12 s132 47 CCNC, (CAN) + COCC, (CAO)
16.340.38.44s126 18 s132 −22 s135 12 CCNC, (CAN) + COCC, (CAO) + CCCC, ring 3

, , , and denote the stretching, bending, torsion, and out ( ABCD means the angle between the AD vector and the BCD plane) modes. Indices notation: s: symmetric; as: asymmetric; A: aliphatic; ring 1: C1-C2-C3-C4-C5-C6; ring 2: C15-C16-C17-C18-C19-C20; ring 3: C21-C22-C23-C24-C16-C17; ring 4: C17-C18-C24-C25-C26-C27; ring 5: C23-C24-C25-C28-C29-C30.
4.6.1. Ring, C=O and C=N Vibrations

The C–H stretching vibrations in the range 2800–3200 cm−1 are for aromatic compounds. From the PED analysis we found that C–H stretching vibrations for ring 1 are assigned at 3183 cm−1. This mode is very pure mode as its PED analysis is about 99%. The values observed in the range (3158–3224) cm−1 are assigned to the stretching vibrations of methyl hydrogen’s while their experimentally obtained counterparts are at 3200 cm−1 and 3118 cm−1, respectively. The percentage of PED calculated for these modes by VEDA 4 program varies from 92 to 99% indicating that they are pure modes. C–N modes of vibrations are assigned on the basis of PED calculations. In the PED analysis we found that C–N modes of vibrations are at 1407 and 1387 cm−1 respectively, however these modes are not pure modes and are mixed with C–C stretching modes, while experimentally observed value is at 1611 cm−1. On the basis of PED analysis the wavenumbers at 1648, 1643, 1675, 1554, and 1526 cm−1 are assigned to C=O stretching modes, however again these modes are not pure modes and are mixed with other modes of vibrations. Experimentally obtained values for C=O stretching modes is at 1649 cm−1. The PED analysis for various modes of the title compound along with their percentage values are summarized in Table 5.

4.7. NBO Analysis

In order to understand the hyper conjugation as well as delocalization of the title compound we have investigated the natural bond orbital analysis of the title compound using NBO 3.1 program implemented in Gaussian 09 package [25]. We have used DFT/B3LYP/6-311++G(2d, 2p) level of theory in order to understand different kind of interactions between the filled and the vacant orbitals. We can investigate both intra- and intermolecular interactions using NBO analysis. In addition to this NBO analysis is also useful for understanding charge transfer conjugative interactions in different compounds. Using DFT/B3LYP/6-311++G(2d, 2p) level of theory the second-order perturbation theory analysis of Fock matrix in NBO basis [26] for title compound is listed in Table 6. For each donor and acceptor the stabilization energy associated with the delocalization is determined as Large value shows the intensive interaction between electron-donors and electron-acceptors groups and greater extent of conjugation of the whole system. The possible intensive interactions are also listed in Table 6. The second-order perturbation theory analysis of Fock matrix in NBO basis shows strong intramolecular hyper conjugative interactions of electrons. From Table 6 we can see that the intramolecular hyper conjugative interactions are formed by the orbital overlap between oxygen, nitrogen, and carbon-carbon bond orbitals. This orbital overlapping is responsible for ICT causing stabilization of the system under study. From the analysis of Table 6 we found that the strong intramolecular hyper conjugative interaction is of C7–O9 from n2 (C7–O9)which increases ED (0.10070 e) that weakens the respective bonds leading to stabilization of 32.78 kcal mol−1. Similarly another strong intramolecular hyper conjugative interaction of C7–O8 from n2(C7–O8) increases ED (0.28565 e) that weakens the respective bonds leading to stabilization of 46.67 kcal mol−1. We have also found another strong intramolecular hyper conjugative interaction of C21–C22 from n1 (C21–C22) which increases ED (0.18831 e) that also weakens the respective bonds leading to stabilization of 46.88 kcal mol−1. We predicted one more strong intramolecular hyper conjugative interaction of C21–C22 from n1(C21–C22) which increases ED (0.18831 e) that weakens the respective bonds leading to stabilization of 46.88 kcal mol−1, as well as strong intramolecular hyper conjugative interaction of C24–C25 from n1(C24–C25) which increases ED (0.46625 e) that weakens the respective bonds leading to stabilization of 71.49 kcal mol−1. These interactions are observed as an increase in electron density (ED) in C–C antibonding orbitals that weakens the respective bonds.

Donor ( )TypeED/eAcceptor ( )TypeED/e (2)a (kcalmol−1) (a.u) (a.u)

C1–C2 1.68952C3–C4 0.4559525.080.260.075
C1–C2 1.68952C5–C6 0.3553115.880.280.060
C1–C6 1.63511C9–C10 0.341944.651.250.068
C3–C4 1.55498C1–C2 0.2746413.450.290.058
C3–C4 1.55498C5–C6 0.3553125.350.280.077
C3–C4 1.55498C13–N14 0.2194023.300.250.072
C5–C6 1.65905C1–C2 0.2746422.390.280.072
C5–C6 1.65905C1–C2 0.4559515.420.260.059
C5–C6 1.65905C7–O8 0.2856527.550.240.074
C13–C14 1.91336C3–C4 0.455956.890.940.023
C13–C14 1.91336C15–C16 0.443689.430.370.058
C15–C16 1.96604C13–N14 0.2194011.210.240.049
C15–C16 1.96604C17–C18 0.4691219.030.280.066
C15–C16 1.96604C19–C20 0.3286119.830.270.067
C15–C16 1.96604C21–C22 0.1883115.740.290.064
C17–C18 1.52506C15–C16 0.4436820.290.260.066
C17–C18 1.52506C19–C20 0.3286119.870.260.066
C17–C18 1.52506C24–C25 0.4662519.540.270.065
C17–C18 1.52506C26–C27 0.1933614.490.280.061
C19–C20 1.68814C15–C16 0.4436819.180.280.068
C19–C20 1.68814C17–C18 0.4691217.640.290.066
C21–C22 1.79067C15–C16 0.4436817.230.280.066
C24–C25 1.53004C17–C18 0.4691220.010.270.065
C24–C25 1.53004C26–C27 0.1933614.620.270.061
C26–C27 1.79745C17–C18 0.4691214.320.290.062
C26–C27 1.79745C24–C25 0.4662514.020.300.062
LP(2)O8 1.86096C6–C7 0.0589415.130.690.093
LP(2)O8 1.86096C7–O9 0.1007032.780.570.123
LP(1)O9 1.96939C7–O8 0.017696.521.080.075
LP(2)O9 1.80240C7–O8 0.2856546.670.300.108
LP(2)O9 1.80240C7–O8 0.2856546.670.300.108
LP(1)O12 1.97419C3–C4 0.040457.371.050.079
LP(2)O12 1.79696C3–C4 0.4559537.120.310.102
LP(1)N14 1.84279C3–C4 0.0327610.670.740.082
LP(1)C23 0.99827C21–C22 0.1883148.680.140.097
LP(1)C23 0.99827C24–C25 0.4662571.490.140.103
LP(1)C23 0.99827C21–C22 0.3307168.230.130.104

(2) means stabilization energy.
bEnergy difference between the donor and acceptor NBO orbitals.
is the Fock matrix element between and NBO orbitals.

5. Conclusions

Using DFT/B3LYP/6-311++G(2d, 2p) level of theory a detailed study of molecular structure, NMR chemical shifts, electronic properties, MEP surface mapping, NBO analysis, and vibrational and PED analysis of the title compound has been investigated and reported. On comparison with experimentally obtained parameters by one of coauthors of this paper we found that both of them are in agreement with each other. HOMO-LUMO analysis of the title compound shows that the electron charge distribution is mainly concentrated over the rings and there may be a charge transfer through system which accounts for bioactivity of the molecule. The title compound has also large band gap as reported in HOMO-LUMO analysis which accounts for its future applications as a useful material in solar cell devices. Molecular electrostatic surface maps give an idea about the chemical reactivity of the title compound. Our overall simulated results for different molecular properties of the title compound are obtained for the first time and we hope that they are helpful in the synthesis and design of new applications.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors thank Indian Institute of technology Mandi for providing the infrastructure required for computational studies as well MHRD scholarships. The authors also want to thank Dr. C. P. Parameswaran for allowing them to use some of the experimental data for comparison with simulated results along with useful and fruitful discussions for the completion of the paper.


  1. A. K. Gupta, A. Dhir, and C. P. Pradeep, “A fluorescence 'turn-on' chemodosimeter for selective detection of Nb5+ ions in mixed aqueous media,” Dalton Transactions, vol. 42, no. 36, pp. 12819–12823, 2013. View at: Publisher Site | Google Scholar
  2. J. R. Lakowicz, Topics in Fluorescence Spectroscopy. Volume 4: Probe Design and Chemical Sensing, Plenum Press, New York, NY, USA, 1994.
  3. J. L. H. Lindenhovius, E. M. Hornsveld, A. den Ouden, W. A. J. Wessel, and H. H. J. ten Kate, “Powder-in-tube (PIT) Nb3Sn conductors for high-field magnets,” IEEE Transactions on Applied Superconductivity, vol. 10, no. 1, pp. 975–978, 2000. View at: Publisher Site | Google Scholar
  4. D. Touati, “Iron and oxidative stress in bacteria,” Archives of Biochemistry and Biophysics, vol. 373, no. 1, pp. 1–6, 2000. View at: Publisher Site | Google Scholar
  5. B. Çatikkaş, E. Aktan, and Z. Seferoǧlu, “DFT, FT-Raman, FTIR, NMR, and UV-Vis studies of a hetarylazo indole dye,” International Journal of Quantum Chemistry, vol. 113, no. 5, pp. 683–689, 2013. View at: Publisher Site | Google Scholar
  6. A. Raj, Y. Sheena Mary, C. Yohannan Panicker, H. T. Varghese, and K. Raju, “IR, Raman, SERS and computational study of 2-(benzylsulfanyl)-3,5-dinitrobenzoic acid,” Spectrochimica Acta A: Molecular and Biomolecular Spectroscopy, vol. 113, pp. 28–36, 2013. View at: Publisher Site | Google Scholar
  7. V. Chiş, S. Filip, V. Miclǎuş et al., “Vibrational spectroscopy and theoretical studies on 2,4- dinitrophenylhydrazine,” Journal of Molecular Structure, vol. 744–747, pp. 363–368, 2005. View at: Publisher Site | Google Scholar
  8. M. Kumru, V. Küçük, and P. Akyürek, “Vibrational spectra of quinoline-4-carbaldehyde: combined experimental and theoretical studies,” Spectrochimica Acta A: Molecular and Bimolecular Spectroscopy, vol. 113, pp. 72–79, 2013. View at: Publisher Site | Google Scholar
  9. R. J. Xavier and P. Dinesh, “Conformational stability, vibrational spectra, HOMO-LUMO and NBO analysis of 1,3,4-thiadiazolidine-2,5-dithione with experimental (FT-IR and FT-Raman) techniques and scaled quantum mechanical calculations,” Spectrochimica Acta A: Molecular and Biomolecular Spectroscopy, vol. 113, pp. 171–181, 2013. View at: Publisher Site | Google Scholar
  10. C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Review B, vol. 37, no. 2, pp. 785–789, 1988. View at: Publisher Site | Google Scholar
  11. R. Dltchfield, “Molecular orbital theory of magnetic shielding and magnetic susceptibility,” The Journal of Chemical Physics, vol. 56, no. 11, pp. 5688–5691, 1972. View at: Google Scholar
  12. K. Wolinski, J. F. Hinton, and P. Pulay, “Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations,” Journal of the American Chemical Society, vol. 112, no. 23, pp. 8251–8260, 1990. View at: Google Scholar
  13. E. D. Glendening and A. E. Reed, NBO Version 3.1, TCL, University of Wisconsin, Madison, Wis, USA, 1998.
  14. M. Izadyar and M. Khavani, “Quantum chemistry aspects of the solvent effects on the ene reaction of 1-Phenyl-1,3,4-triazolin-2,5-dione and 2-methyl-2-butene,” International Journal of Quantum Chemistry, vol. 114, no. 10, pp. 666–674, 2014. View at: Publisher Site | Google Scholar
  15. H. Pir, N. Günay, Ö. Tamer, D. Avci, and Y. Atalay, “Theoretical investigation of 5-(2-Acetoxyethyl)-6-methylpyrimidin-2,4-dione: Conformational study, NBO and NLO analysis, molecular structure and NMR spectra,” Spectrochimica Acta A: Molecular and Biomolecular Spectroscopy, vol. 112, pp. 331–342, 2013. View at: Publisher Site | Google Scholar
  16. Diwaker, “Quantum mechanical and spectroscopic (FT-IR, 13C, 1H NMR and UV) investigations of 2-(5-(4-Chlorophenyl)-3-(pyridin-2-yl)-4,5-dihydropyrazol-1-yl)benzo[d]thiazole by DFT method,” Spectrochimica Acta A: Molecular and Biomolecular Spectroscopy, vol. 128, pp. 819–829, 2014. View at: Publisher Site | Google Scholar
  17. J. Fleming, Frontier Orbitals and Organic Chemical Reactions, John Wiley and Sons, New York, NY, USA, 1981.
  18. R. G. Pearson, “Absolute electronegativity and hardness correlated with molecular orbital theory,” Proceedings of the National Academy of Sciences, vol. 83, no. 22, pp. 8440–8441, 1986. View at: Publisher Site | Google Scholar
  19. P. Politzer, P. R. Laurence, and K. Jayasuriya, “Molecular electrostatic potentials: an effective tool for the elucidation of biochemical phenomena,” Environmental Health Perspectives, vol. 61, pp. 191–202, 1985. View at: Google Scholar
  20. A. Pullman, B. Pullman, and R. Lavery, “Molecular electrostatic potential versus field. significance for DNA and its constituents,” Journal of Molecular Structure: THEOCHEM, vol. 93, pp. 85–91, 1983. View at: Publisher Site | Google Scholar
  21. V. D. Vitnik, J. Ž. Vitnik, N. R. Banjac, N. V. Valentić, G. S. Ušćumlić, and I. O. Juranić, “Quantum mechanical and spectroscopic (FT-IR, 13C, 1H NMR and UV) investigations of potent antiepileptic drug 1-(4-chloro-phenyl)-3-phenyl-succinimide,” Spectrochimica Acta A: Molecular and Biomolecular Spectroscopy, vol. 117, pp. 42–53, 2014. View at: Publisher Site | Google Scholar
  22. B. Çatıkkaş, E. Aktan, and Z. Seferoǧlu, “DFT, FT-Raman, FTIR, NMR, and UV-Vis studies of a hetarylazo indole dye,” International Journal of Quantum Chemistry, vol. 113, no. 5, pp. 683–689, 2013. View at: Publisher Site | Google Scholar
  23. M. H. Jamróz, J. C. Dobrowolski, and R. Brzozowski, “Vibrational modes of 2,6-, 2,7-, and 2,3-diisopropylnaphthalene. A DFT study,” Journal of Molecular Structure, vol. 787, no. 1–3, pp. 172–183, 2006. View at: Publisher Site | Google Scholar
  24. Gaussian 09, “Gaussian,” Wallingford CT, 2004. View at: Google Scholar
  25. N. Günay, H. Pir, D. Avci, and Y. Atalay, “NLO and NBO analysis of Sarcosine-Maleic acid by using HF and B3LYP calculations,” Journal of Chemistry, vol. 2013, Article ID 712130, 16 pages, 2013. View at: Publisher Site | Google Scholar
  26. W. B. Person and J. H. Newton, “Dipole moment derivatives and infrared intensities. I. Polar tensors,” The Journal of Chemical Physics, vol. 61, no. 3, pp. 1040–1049, 1974. View at: Google Scholar

Copyright © 2014 Diwaker and Abhishek Kumar Gupta. 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