Academic Editor:Dr. Praveen Kumar Sharma, Lovely Professional University, Phagwara, Punjab, India-144411
Checked for plagiarism: Yes
Review by: Single-blind
Vibrational, Ultra Violet, Natural Bond Orbital Analysis of E-1 using Quantum Mechanical Computations and Experimental Spectra.
The FT-IR, FT-Raman and UV-Vis spectra of E-[1-(3'-methylthienyl)-5-Phenyl-2,4-Pentadiene-3-one (MPPO) were recorded. The optimized molecular bond parameters, harmonic frequencies were calculated using B3LYP method with 6-311++G (d,p) basis set.The various normal modes were precisely assigned with thehelp ofTED calculation. The theoretical spectrograms for FT-IR, FT-Raman and Ultra Violet visible. Spectra of the title molecule had been constructed. The ICT was calculated by means of Natural Bond Orbital analysis. The Non Linear Optical properties related to polarizability and hyperpolarizability based on the finite-field approach were calculated.The band gap energy was calculated using HOMO-LUMO analysis. Furthermore, the Molecular Electrostatic Potential, Mulliken atomic charges and thermodynamic properties of MPPO were also calculated.
Chalcones belong to the flavonoid families which are synthesized in factories to preserve the health of plants against infections and parasites. They have attractedincreasing attention due to numerous pharmacological applications 1, 2, 3, 4, 5, 6, 7, 8.It possess various multipronged activities such as anti-microbial 9, anti-depressants 10, anti-plasmodial 11, anti-aids 12, insect anti-feedant activities 13, 14, biological treatment due to its good anti-malaria 15, and in vitro anti-tumor activity 16. It exhibits radical quenching and hydroxyl adducts formation 17. Whilst chalcones are active against various protein targets, modification of the privileged core could lead to novel compounds with specifically targeted inhibitory activity 18.
Chidan Kumar et al., 19 synthesized the high quality single crystal of efficient novel NLO chalcone derivative (2E)-1-(5-chlorothiophen-2-yl)-3-(2,3,4-trimethoxyphenyl) prop-2-en-1-one and its structure was characterized by FTIR, FT-Raman and single crystal XRD techniques. The vibrational wavenumbers were computed using DFT and were assigned on the basis of PED analysis. The geometrical parameters obtained from XRD studywere compared with the calculated values by applying DFT/6-31G (d,p) basis set. Stability of the molecule, hyper conjugative interactions, charge delocalization and intra-molecular hydrogen bond had been studied using NBO analysis. Karunakaran et al., 20 reported the FTIR and FT-Raman spectra of trans-3-(o-hydroxyphenyl-1-phenyl)-2-propen-1-one (or simply 2-hydroxychalcone) were recorded in the regions 4000–400 cm-1 and 3500–100 cm-1, respectively in the solid phase. The vibrational wavenumbers were calculated by HF and DFT/B3LYP methods with 6-311++G (d, p) basis set, using Gaussian 09W program package. A detailed interpretation of the IR and Raman spectra of 2-hydroxychalcone was reported.
In light of the above literature survey, it was clear that inspite of many important applications of MPPO not many works were carried out on this molecule; particularly the complete vibrational analysis supported by quantum computations was not carried out. Hence in this study, FTIR, FT-Raman and UV-Vis. spectral analysis of the compound MPPO was carried out. The NBO analysis, with emphasis on HOMO-LUMO, NLO, MEP, Mulliken’s charges, and various thermodynamic parameters were also calculated using quantum computations by DFT/B3LYP method.
Mono benzal acetone (1.46g, 0.01mol) was dissolved in 15 mL of ethanol with slight warning. To this hot solution methyl thiophen-2-carboxaldehyde (1.4 mL, 0.01mol) and few drops of solution hydroxide solution 10%) were added. The solution gradually turned red on warming and yellow crystals separated. The product was filtered off and recrystalised twice from ethanol. The yield and melting point of the crystal is 76% and 148 respectively.
The FT-IR spectrum on an IFS 66v spectrophotometer of MPPO was recorded in the spectral region between 400–4000 cm-1 using the KBr pellet technique. The spectrum was recorded at room temperature with a scanning speed of 10 cm-1 per minute and at the spectral resolution of 2.0 cm-1 in the Department of Chemistry, Jamal Mohamed College, Trichy, Tamilnadu, India. The FT-Raman spectrum of title compound was recorded using the 1064 nm line of an Nd: YAG laser as excitation wavelength in the region 50–3500 cm-1 on Bruker model IFS 66V spectrophotometer equipped with an FRA 106 FT-Raman module accessory and at spectral resolution of 4 cm-1. The FT-Raman spectral measurement was carried out from SAIF Laboratory, IIT Madras, Tamilnadu, India. The UV-Vis absorption spectrum of MPPO was recorded in the range of 200–500 nm using a Shimadzu – 2600 spectrometer in the Department of Chemistry, Jamal Mohamed College, Trichy-20. The UV pattern was taken from a 10-5 molar solution of MPPO dissolved in benzene.
For meeting the requirements of both accuracy and computing economy, theoretical methods and basis sets should be considered. DFT had proved to be extremely useful in treating electronic structure of molecules. The density functional theory parameter hybrid model DFT/B3LYP/6 basis set was adopted to calculate the properties of the title molecule in this work. 50All the calculations were performed using the Gaussian 03W program package 21 with the default convergence criteria without any constraint on the geometry 22. It should be noted that Gaussian 03W package did not calculate the Raman intensities. The Raman activities were transformed into Raman intensities using Raint program 23.
The optimized geometrical parameters for MPPO are calculated using B3LYP/6-311++G(d,p) level of calculation and are presented in Table 1 in accordance with atom numbering given in Figure 1. The title molecule consists of thiophene and phenyl ring fused by carbonyl group. In MPPO the C=O group plays an important role. The carbonyl (C16=O17) bond length is observed at 1.190 Å24, whereas the calculated bond length is about 1.235Å.In thiophen ring, the two C-S bond lengths are found at 1.707 and 1.728Å.In this study, the C1-S5(1.762Å) and C4-S5(1.717Å) bond lengths are differ by 0.045Å and the bond angles C2-C1-S5(110.34˚) and C3-C4-S5(112.86˚) also differ by 2.52˚,which is due to the shortening of bond distance between atoms C1 and S5 atoms. The bond angle O17-C16-C18(124.37˚) is positively (~2.98˚) from C14-C16-O17(121.39˚). This may be due the intra- or inter-molecular interaction S---O---H. From the theoretical values, it is found that most of the optimized bond parameters are slightly greater than the experimental values, since the calculation has been done on a single molecule in gaseous state.Although the differences, calculated geometrical parameters represent a good approximation and they can be used as foundation to calculate the other parameters, such as vibrational frequencies and thermodynamics properties.The calculated bond parameters are compared with experimental values24, 25. The optimized structure is planar as it is evident from the dihedral angles C12C14C16O17=0.00° & C12C14C16C18=-180.00°; C14C16C18C20=180.00° & O17C16C18C20=-0.00°;C2C1C12C14=179.99˚& S5C1C12C14=-0.00˚ and C18C20C22C23=-179. 99˚ & C18C20C22C24=0.00˚ are listed in Table 1.Table 1. The optimized bond parameters of MPPO using B3LYP/6-311++G(d,p) basis set
|Parameters||B3LYP/ 6-311++G(d,p)||XRD a,b,c,d|
|Bond Lengths (Å)|
|C4-S5||1.717||1.717a, 1.707b, 1.714d|
|Bond Angles (°)|
|C1-C2-C3||112.64||112.20a, 114.72b, 112.5d|
|C1-S5-C4||91.27||92.00a, 91.81b, 92.2d|
|Dihedral Angles (°)|
The molecule MPPO belongs to C1 point group symmetry. This molecule has 32 atoms and hence90 normal modes of fundamental vibrations are possible, which span the irreduciable representation:61A'+29A''. All the 90 fundamental modes are active in both IR and Raman. The vibrational analysis of MPPO is performed on the basis of the characteristic vibrations of carbonyl, methyl, thiophene, pentadiene and phenyl ring modes. In this study, the harmonic vibrational frequencies are calculated at B3LYP level using 6-311++G(d,p)basis set have been collected in Table 2. Some discrepancies could be identified in between harmonic and experimental frequencies, which are scaled down by proper scale factor 26. To find the exact vibrational behavior of this molecule the TED analysis was performed on the gas phase. The combined experimental and theoretical spectra are shown in Figure 2,Figure 3.Table 2. The experimental and calculated frequencies of MPPO using B3LYP/6-311++G(d,p) level of basis set harmonic frequencies (cm−1) IR Raman intensities (Km/mol) reduced masses (amu) and force constants (mdynA°−1)
|Mode No||Calculated||Observed Frequencies (cm-1)||IR Intensity||Raman Intensity||Reduced Masses||Force Consts||Vibrational Assignments≥10% (TED)d|
|4||3188||3063||3064 vw||3059 vw||4.7||1.04||1.1||6.57||νC8H9(83)|
|13||3073||2953||2919 w||2921 vw||2.57||0.27||1.1||6.13||νC27H31(100)|
|15||1696||1629||1647 s||1647 vw||12.94||3.2||7.03||11.91||νC16O17(85)|
|16||1642||1578||1617 m||1616 w||17.67||1.55||5.53||8.78||νC24C27(18)+νC23C25(25)+νC25C29(14)+βH26C23C25(15)|
|17||1624||1560||1586 vs||1584 vs||100||21.83||5.42||8.43||νC4C3(32)+νC24C27(12)|
|18||1604||1541||1550 m||1551 vw||98.37||3.24||5.33||8.08||νC18C20(46)+βH32C29C27(20)|
|20||1534||1474||1492 m||1494 vw||5.18||0.67||3.93||5.44||νC16C18(45)|
|26||1467||1409||1404 m||1406 vw||2.2||0.77||1.82||2.31||βH19C18C20(58)|
|35||1216||1168||1183 m||1185 w||3.88||0.24||1.25||1.09||βH26C23C25(55)|
|39||1119||1075||1090 w||1090 m||36.76||1.15||1.92||1.42||νC18C20(12)+νC16C18(24)+βH19C18C20(12)+βH7C4S5(16)|
|62||718||689||697 w||691 vvw||2.72||0.02||1.34||0.41||τH19C18C20C22(44)+τC1C12C14C16(14)+τH30C25C29C27(26)|
The carbonyl stretching frequency has been most extensively studied by IR spectroscopy 27. The multiply bonded group is highly polar and therefore gives rise to an intense IR absorption band. The carbon–oxygen double bond is formed by Pπ–Pπ bonding between carbon and oxygen atoms. The loan pair of electrons on oxygen also determines the nature the carbonyl group. The νC=O mode of carboxylic acids is identical to the νC=O mode in ketones, which is expected in the region 1660-1740 cm-128.The stretching vibration of carbonyl group C=O can be observed as a very strong band in both FTIR and FT-Raman spectra at 1665cm-129, 30.According to above literature strong band observed in FTIR spectrum at 1647 cm-1 and weak band in FT-Raman spectrum at 1647 cm-1 is assigned to C=O stretching vibration.The value of this band is calculated 1629cm-1 (mode no:15) with a TED of 85%. This assignment is in moderate agreement with literature 31 in the case of thiophene-2-carohydrazide. The observed very weak band in FT-Raman at 724cm-1 is ascribed to βC16=O17 mode, while its correspondingharmonic frequency is 733cm-1 (mode no:60). These assignments are also supported by the literature31in addition to TED (36%) output.
The C-S stretching vibration is well known to mix with neighboring modes 32.Literature survey revealed that ν(C-S) vibration of thiophene ring was observed in the region 719-842cm-133. The harmonic wave numbers 811,687 (mode nos:56,63) and FTIR band observed at 812 cm-1 are assigned to νC-S modes. These modes are mixed with βCCC modes and also have considerable TED (>32%) values. Furthermore, these assignments also find support from the literature34.
Methyl groups are generally referred to as electron donating substituents in the aliphatic and aromatic ring system 35. For the assignments of CH3 group, basically nine fundamentals can be associated with each CH3 group namely: the symmetrical stretching in CH3 (CH3sym. stretching); anti-symmetrical stretching (CH3asym. stretching); the symmetrical (CH3 sym. deformation) and anti-symmetrical (CH3asy. deformation) deformation modes; the in-plane rocking (CH3ipr), out-of-plane rocking (CH3opr) and twisting (CH3) bending modes. The asymmetric/symmetric CH stretching modes of CH3 are expected in the region around 2980/2870 cm-1, respectively 36, 37, 38, 39. The CH3asy/sym frequencies are established at 3068, 3063/3052 cm-1 (mode nos: 3, 4/5), respectively. This assignment is further supported by observed bands (FTIR: 3064/FT-Raman: 3059 cm-1) in addition to TED values (>83%). The asymmetric/symmetric CH3 bending vibrations usually occur in the regions 1430-1470cm-1/1360-1400 cm-1, respectively 40. In accordance with the above literature, the harmonic bands 1445, 1432/1357cm-1 (mode nos: 22, 24/28) are respectively assigned to βasy CH/βsy CH of CH3 group. Further, these assignments are supported from observed FTIR band: 1447 cm-1 in addition to TED values (>60%).
The CH3 rocking modes usually appeared39in the region 1010-1070 cm-1. The bands established at 983 and 966 cm-1 in FT-Raman and FTIR were assigned to CH3 in-plane and out-of-plane rocking modes, respectively 41. The value of these bands are calculated at 1010 and 998 cm-1 (mode nos: 43 and 44) with considerable TED values (68% & 40%). The mode no: 44 coincideswith observed FT-Raman band (1000 cm-1). These assignments find support from the literatureSubramanian et al., 42 and were within the frequency intervals given by Varsanyi et al.,43. The CH3 torsional mode was assigned at 168 cm-142. It should be mentioned here, that this mode is pure with 95% TED value. Hence the mode no: 85 is undoubtedly assigned to τCH3 mode.
The ring C=C and C-C stretching vibrationsknown as semicircle stretching usually occur in the region 1450–1625 cm-144, 45. Varsanyi et al., 46have assigned the νC-C vibrations in the regions of 1590-1625,1575-1590,1470-1540,1430-1460 and 1280-1380cm-1 with variable intensity of bands. The ring νC-Cmodes occur in the region 1530-1625cm-142.In this region there are two (or) three bands due to skeletal vibrations, the strongest usually being at about 1500cm-1 for aromatic six member rings. Hence, the bands of strong to medium intensity in FTIR at1617 (m),1586 (vs) & 1012 (m) cm-1 and FT-Raman at 1616(w),1584(vs),1314(w) & 1258 (w) cm-1 are assigned to C-C stretching vibrations of benzene moiety. The frequencies assigned to νC-C modes agree well with the literature 44, 45 and also in line with the harmonic wavenumbers in the range 1578-976cm-1(mode nos:16,17,31,33,42, 47) in addition to TED values.
The βCCC and ΓCCCmodes are associated with smaller force constant than the stretching one and hence assigned to lower frequencies. The βCCC modes of benzene ring are expected to appear with considerable intensity under the reduced symmetry.The C-C-C in-plane bending vibrations can be observed in the FTIR spectrum at 1012cm-1 and in FT-Raman spectrum at 634 cm-1.Similarly the C-C-C out-of-plane bending modes can be observed as a medium band in FTIR spectrum at 455cm-1. These results are supported by computed wavenumbers 1011, 976, 637 and 670, 505, 470 cm-1 (mode nos:42, 47, 65 and 64, 71,72), respectively in addition to literature values 39.
Literature survey reveals that the C-C stretching vibrations in substituted thiophen rings were reported in ranges of 1329-1431,1420-1501 and 1419-1519cm-128, 30, 31. Inour case the very strong bands at 1586/1584 cm-1 in FTIR/FT-Raman and harmonic bands at 1560, 1370 and 1148 cm-1( mode nos:17,27 &37) are assigned to νC-C modes in MPPO. These assignments are having considerable TED values(>24%). As it is evident from the Table 2, the βCCC modes: 811,687 cm-1(mode nos:56,63) are mixed with νC-S modes. Mode no: 56 is in good agreement with observed FTIR band at 812cm-1 in addition to TED output (>15%). These assignments are also in accordance with literature28, 30, 31. The mode nos:67,79 (601, 228cm-1) are designated as ΓCCC modes of thiophen moiety with considerable TED values (>58%).
The harmonic frequencies established at 1541,1474,1330,1075 cm-1(mode nos: 18,20,29,39) and 1158cm-1, (mode no:36) are attributed to νCC modes of pentadiene moietyand νC20-C22,νC1-C12 modes, respectively. These assignments are having considerable TED values (>25%) and also find support from the observed bands:1550,1492, 1334,1090 cm-1(FTIR)/1551,1494,1090 cm-1 FT-Raman. The mode no:34 (1212/1193cm-1:FTIR) is attributed νC2-C8 mode. The harmonic wavenumbers:894,607,528, 445 cm-1(mode nos:52,66,69,73) and 689,737,336cm-1(mode nos:62,59,77) are respectively assigned to βCCC and τCCC modes of pentadiene, in which mode nos:69,73 and 62 are agreeable with observed bands: (FT-Raman/FTIR)531/424cm-1 and 691/697cm-1. These assignment are also having considerable TED values.
The νC-H modes of hetero aromatic structure are expected to occur in the range of 3000-3100cm-1 with some weak bands. The vibrational band in this region are not affected appreciably by the nature of the substituent’s 43, 44. In MPPO,the weak bands at 2919 (2921:FT-Raman), 2843cm-1 in FTIR spectrum have been assigned to νC-H vibrations. The present theoretical calculation by B3LYP/6-311++G(d,p) method places this modes at 3104,3041,2984,2953,2907cm-1 (mode nos:1,7,12,13,14). This shows that the bands have appeared at the expected position of the spectra, except the last band (2843cm-1) which is slightly less.The reduction in the frequency value is naturally due to the presence of CH3 group inMPPO.
In aromatic compounds the C-H in-plane /out-of-plane bending frequencies appear in the ranges of 1000-1300/750-1000cm-1, respectively45, 46. In this work the βC-H modes are assigned in the range 1307-1011cm-1 (modenos:31,33,35,38,42), which are in agreement with observed bands in FTIR at 1183,1012 cm-1 and in FT-Raman at 1314,1258,1185 cm-1. The calculated frequencies 991-777cm-1 (mode nos: 45, 49, 51, 55, 58) for C-H out-of-plane bending modes fall in the FTIR values of 985,760cm-1. These assignments are find support from the Varsanyi 46 and also find support from TED values. Hence, among the CH vibrations, only the CH stretching modes are found influenced by the substitutional group, whereas the βCH and ΓCHmodes are remain unaffected.
The C-H stretching modes are expected to appear with multiple weak bands in the frequency range 3000-3100cm-1. These bands are not affected appreciably by the nature of the substituent’s 43, 47.The νC-H modes were observed at 3011,3062 and 3072 cm-1 in FTIR spectrum by Balachandran et al., 31. Based on the above conclusion, the calculated frequencies 3027,2989cm-1 (mode nos: 9, 11) are assigned to νC-H vibrations of thiophen moiety. These assignments are aggreable with observed FTIR band at 3021 cm-1 and also find support from TED values (>75%).The C-H in-plane bending modes appeared by sharp but weak to medium bands in the range 1100-1500cm-1 and the bands were not sensitive to the nature of substitutent’s. The ΓC-H out-of-plane bending modes were expected to occur in the region 800-1000cm-147.With reference to this, the observed medium band at 1193cm-1 in FTIR is attributed to βC-H of thiophen ring. For the same mode the theoretical frequencies are:1370,1212cm-1 (mode nos:27,34). The harmonic wavenumbers:1066 and 862, 705cm-1 (mode nos:40 and 53,61) are designated as βHCS and τHCCS modes, respectively and theseassignments have considerable TED values (>50%).
In this work, it has been established well and the calculated wavenumbers in the range 3098-3016cm-1 (mode nos:2,6,8,10) are designated as νC-H modes of pentadiene moiety46, 47. For this mode, the corresponding vibrational bands are missing in the experimental spectra, which may be dueto isomer. Hence, the above assignments are justified with the help of TED values (>72%).The βC-H modes are attributed to the FTIR/FT-Raman bands present at 1404/1406cm-1 and for the same mode the harmonic frequencies are 1464,1437,1428,1409cm-1(mode nos: 21,23,25,26). This assignment is aggreable with literature 47in addition to TED results (>42%). The τCCCH vibrations are characterized by the FT-Raman/FTIR bands at 834,691/697cm-1,while the calculated wavenumbers are: 979, 958,831,825,689cm-1(mode nos:46,48,54,55,62) and also have considerable TED values (>42%).
The total dipole moment and first-order hyperpolarizability (β0) of MPPO is calculated by using B3LYP/6-311++G(d,p) basis set and are listed in Table 3.The calculated dipole moment and the mean first hyperpolarizability (β0) values are 0.7071 Debye and 10.2328x10-30 esu, respectively, and the values are presented in Table 3. Total hyperpolarizability of the title molecule is approximately twenty seven times greater than that of urea. The above results show that MPPO can be best material for NLO applications.Table 3. The NLO measurements of MPPO.
|Dipole moment ( μ ) Debye|
|Polarizability ( α0) x10-30esu|
|Hyperpolarizability ( β0 ) x10-30esu|
The NBO analysis is performed on MPPO using B3LYP/6-311++G(d,p) basis set. The E(2) energies and types of interactions are listed in Table 4.The EDs of the conjugated single as well as double bond of the aromatic/thiophene rings are 1.979& 1.606e/1.979&1.773e,respectively.This exhibits a strong delocalization inside the molecule MPPO. The non-bonding atoms transfer more energy to the acceptor, during the intra-molecular interactions. The LP(2) S5 and LP(2) O17 atoms transfer more energy (27.15 and 14.61 KJ/mol) to anti-bonding orbitals C3-C4 and C14-C16, respectively, when compare with σ-σ* transition for the same. These interactions causes a pronounced and moderate decrease of the lone pair orbitals occupancies; 1.581e and 1.886e, respectively, and there is a possibility for hyperconjugativeinteraction between them.The donor and acceptor interactions of πC22-C24→π*C23-C25(20.9 KJ/mol); πC27-C29→π*C23-C25(18.55KJ/mol) reveal the maximum hyperconjugative interaction in benzene ring of MPPO.Table 4. The second order perturbation theory analysis of Fock Matrix in NBO basis for MPPO
|Type||Donor NBO (i)||ED/e||Acceptor NBO (j)||ED/e||E(2) KJ/mol||E(j)-E(i) a.u.||F(i,j) a.u.|
|σ-σ*||BD ( 1) C 1 - C 2||1.97397||BD*( 1) C 1 - C 12||0.0274||3.73||1.22||0.06|
|BD*( 1) C 12 - C 14||0.01462||1.94||1.96||0.055|
|π-π*||BD ( 2) C 1 - C 2||1.73778||BD*( 2) C 3 - C 4||0.32831||16.62||0.27||0.06|
|BD*( 1) C 8 - H 10||0.00953||1.87||0.66||0.034|
|BD*( 1) C 8 - H 11||0.00953||1.87||0.66||0.034|
|BD*( 2) C 12 - C 14||0.17953||20.25||0.33||0.074|
|σ-σ*||BD ( 1) C 1 - S 5||1.97486||BD*( 1) C 4 - H 7||0.01311||2.76||1.22||0.052|
|BD*( 1) C 12 - H 13||0.01378||1.76||1.08||0.039|
|σ-σ*||BD ( 1) C 1 - C 12||1.97917||BD*( 1) C 12 - C 14||0.01462||1.53||1.96||0.049|
|BD*( 1) C 14 - H 15||0.01324||1.42||1.16||0.036|
|BD*( 1) C 16 - C 18||0.05908||34.4||0.24||0.082|
|BD*( 1) C 18 - H 19||0.01254||120.55||0.08||0.09|
|BD*( 1) C 22 - C 24||0.02627||27.23||0.32||0.083|
|BD*( 2) C 22 - C 24||0.36662||36.29||0.13||0.067|
|σ-σ*||BD ( 1) C 2 - C 3||1.97199||BD*( 1) C 1 - C 12||0.0274||4.38||1.19||0.065|
|BD*( 1) C 3 - C 4||0.01512||2.83||1.36||0.056|
|BD*( 1) C 4 - H 7||0.01311||3.02||1.25||0.055|
|σ-σ*||BD ( 1) C 2 - C 8||1.97883||BD*( 1) C 3 - C 4||0.01512||1.28||1.32||0.037|
|σ-σ*||BD ( 1) C 3 - C 4||1.98319||BD*( 1) C 4 - H 7||0.01311||1.34||1.3||0.037|
|π-π*||BD ( 2) C 3 - C 4||1.8072||BD*( 2) C 3 - C 4||0.32831||0.65||0.28||0.012|
|σ-σ*||BD ( 1) C 3 - H 6||1.97779||BD*( 1) C 3 - C 4||0.01512||1.33||1.21||0.036|
|σ-σ*||BD ( 1) C 4 - S 5||1.97845||BD*( 1) C 1 - C 12||0.0274||4.22||1.18||0.063|
|σ-σ*||BD ( 1) C 4 - H 7||1.98573||BD*( 1) C 3 - C 4||0.01512||1.3||1.23||0.036|
|BD*( 1) C 16 - C 18||0.05908||13.02||0.08||0.03|
|BD*( 1) C 22 - C 24||0.02627||7.54||0.16||0.031|
|σ-σ*||BD ( 1) C 8 - H 9||1.98812||BD*( 1) C 12 - C 14||0.01462||3.95||1.87||0.077|
|BD*( 1) C 16 - O 17||0.01724||4.38||1.27||0.067|
|BD*( 1) C 24 - H 28||0.02271||2.98||1.1||0.051|
|BD*( 2) C 27 - C 29||0.33144||2.97||3.85||0.104|
|BD*( 1) C 27 - H 31||0.01409||2.52||4.81||0.098|
|σ-σ*||BD ( 1) C 8 - H 10||1.9793||BD*( 1) C 22 - C 24||0.02627||74.35||0.12||0.086|
|σ-σ*||BD ( 1) C 12 - H 13||1.96197||BD*( 1) C 14 - H 15||0.01324||1.05||0.96||0.029|
|BD*( 1) C 14 - C 16||0.05431||5.72||1.06||0.07|
|σ-σ*||BD ( 1) C 12 - C 14||1.98129||BD*( 1) C 1 - C 12||0.0274||4.63||1.25||0.068|
|BD*( 1) C 14 - C 16||0.05431||2.42||1.28||0.05|
|BD*( 1) C 16 - C 18||0.05908||8.93||0.27||0.044|
|π-π*||BD ( 2) C 12 - C 14||1.81895||BD*( 2) C 16 - O 17||0.27328||21.41||0.29||0.072|
|σ-σ*||BD ( 1) C 14 - H 15||1.97581||BD*( 1) C 1 - C 12||0.0274||7.09||1.02||0.076|
|BD*( 1) C 16 - O 17||0.01724||3.76||1.16||0.059|
|σ-σ*||BD ( 1) C 14 - C 16||1.9799||BD*( 1) C 12 - H 13||0.01378||1.96||1.09||0.041|
|BD*( 1) C 12 - C 14||0.01462||2.04||1.91||0.056|
|BD*( 1) C 16 - C 18||0.05908||5.72||0.19||0.03|
|BD*( 1) C 18 - H 19||0.01254||4.58||0.03||0.011|
|σ-σ*||BD ( 1) C 16 - O 17||1.99368||BD*( 1) C 16 - C 18||0.05908||2.76||0.57||0.036|
|π-π*||BD ( 2) C 16 - O 17||1.96036||BD*( 2) C 12 - C 14||0.17953||4.25||0.43||0.039|
|BD*( 2) C 18 - C 20||0.13135||4.57||0.4||0.039|
|σ-σ*||BD ( 1) C 16 - C 18||1.98102||BD*( 1) C 12 - C 14||0.01462||1.59||1.91||0.049|
|BD*( 1) C 18 - H 19||0.01254||6.35||0.03||0.012|
|BD*( 1) C 20 - H 21||0.01242||1.68||1.09||0.038|
|σ-σ*||BD ( 1) C 18 - H 19||1.97184||BD*( 1) C 16 - O 17||0.01724||4.69||1.15||0.066|
|BD*( 1) C 20 - H 21||0.01242||1.26||0.94||0.031|
|BD*( 1) C 20 - C 22||0.02819||7.7||1||0.078|
|σ-σ*||BD ( 1) C 18 - C 20||1.98216||BD*( 1) C 16 - C 18||0.05908||14.3||0.27||0.056|
|BD*( 1) C 18 - H 19||0.01254||8.42||0.11||0.028|
|BD*( 1) C 20 - C 22||0.02819||3.69||1.23||0.06|
|π-π*||BD ( 2) C 18 - C 20||1.84033||BD*( 2) C 16 - O 17||0.27328||20.32||0.29||0.07|
|σ-σ*||BD ( 1) C 20 - H 21||1.97149||BD*( 1) C 16 - C 18||0.05908||73.62||0.05||0.056|
|BD*( 1) C 22 - C 24||0.02627||40.86||0.13||0.065|
|σ-σ*||BD ( 1) C 20 - C 22||1.97691||BD*( 1) C 16 - C 18||0.05908||9.02||0.2||0.038|
|BD*( 1) C 22 - C 23||0.02207||2.06||1.22||0.045|
|BD*( 1) C 22 - C 24||0.02627||32.81||0.27||0.085|
|BD*( 2) C 22 - C 24||0.36662||7.71||0.08||0.025|
|σ-σ*||BD ( 1) C 22 - C 23||1.97392||BD*( 1) C 20 - C 22||0.02819||2.17||1.17||0.045|
|BD*( 2) C 22 - C 24||0.36662||120.04||0.1||0.106|
|σ-σ*||BD ( 1) C 22 - C 24||1.9728||BD*( 1) C 23 - H 26||0.0143||3.64||1.12||0.057|
|BD*( 1) C 24 - C 27||0.01547||2.59||1.27||0.051|
|π-π*||BD ( 2) C 22 - C 24||1.5862||BD*( 2) C 18 - C 20||0.13135||13.05||0.29||0.059|
|BD*( 2) C 23 - C 25||0.31148||20.9||0.29||0.071|
|BD*( 2) C 27 - C 29||0.33144||2.84||3.47||0.09|
|σ-σ*||BD ( 1) C 23 - C 25||1.97917||BD*( 1) C 4 - H 7||0.01311||3.66||1.28||0.061|
|BD*( 1) C 20 - C 22||0.02819||3.26||1.21||0.056|
|BD*( 1) C 22 - C 23||0.02207||4.39||1.27||0.067|
|σ-σ*||BD ( 1) C 23 - H 26||1.98069||BD*( 1) C 22 - C 23||0.02207||0.69||1.07||0.024|
|σ-π*||BD ( 1) C 24 - C 27||1.97888||BD*( 2) C 3 - C 4||0.32831||3.77||0.9||0.056|
|BD*( 1) C 4 - H 7||0.01311||6.05||1.47||0.084|
|BD*( 1) C 16 - O 17||0.01724||4.9||1.55||0.078|
|BD*( 1) C 20 - C 22||0.02819||4.41||1.39||0.07|
|σ-σ*||BD ( 1) C 24 - H 28||1.97705||BD*( 1) C 16 - C 18||0.05908||27.9||0.02||0.023|
|BD*( 1) C 22 - C 23||0.02207||4.86||1.05||0.064|
|BD*( 1) C 22 - C 24||0.02627||27.94||0.1||0.048|
|σ-π*||BD ( 1) C 25 - C 29||1.97941||BD*( 2) C 3 - C 4||0.32831||4.09||0.69||0.052|
|BD*( 1) C 12 - C 14||0.01462||10.22||1.94||0.126|
|BD*( 1) C 16 - O 17||0.01724||9.45||1.34||0.101|
|BD*( 1) C 23 - H 26||0.0143||6.31||1.13||0.075|
|BD*( 1) C 24 - H 28||0.02271||8.21||1.17||0.087|
|σ-σ*||BD ( 1) C 25 - H 30||1.98016||BD*( 1) C 4 - H 7||0.01311||3.33||1.1||0.054|
|BD*( 1) C 22 - C 23||0.02207||5.6||1.08||0.069|
|BD*( 1) C 24 - H 28||0.02271||3.22||1||0.051|
|BD*( 1) C 29 - H 32||0.01371||2.72||2.61||0.075|
|σ-π*||BD ( 1) C 27 - C 29||1.98062||BD*( 2) C 3 - C 4||0.32831||4.13||0.69||0.052|
|BD*( 1) C 12 - C 14||0.01462||9.11||1.94||0.119|
|BD*( 1) C 16 - O 17||0.01724||9.4||1.34||0.1|
|BD*( 2) C 27 - C 29||0.33144||16.2||3.92||0.245|
|π-π*||BD ( 2) C 27 - C 29||1.64275||BD*(2) C 23 - C 25||0.31148||18.55||0.29||0.066|
|σ-σ*||BD ( 1) C 27 - H 31||1.98023||BD*(1) C 4 - H 7||0.01311||3.3||1.08||0.053|
|BD*(1) C 24 - H 28||0.02271||5.55||0.99||0.066|
|σ-π*||BD ( 1) C 29 - H 32||1.98091||BD*( 2) C 3 - C 4||0.32831||4.65||0.48||0.046|
|BD*( 1) C 12 - C 14||0.01462||26.51||1.73||0.191|
|BD*( 1) C 16 - O 17||0.01724||16.8||1.13||0.123|
|BD*( 1) C 24 - H 28||0.02271||9.57||0.96||0.086|
|BD*( 2) C 27 - C 29||0.33144||23.61||3.71||0.288|
|BD*( 1) C 27 - H 31||0.01409||34.65||4.66||0.359|
|n-σ*||LP ( 1) S 5||1.98351||BD*( 1) C 3 - C 4||0.01512||2.25||1.32||0.049|
|n-π*||LP ( 2) S 5||1.58098||BD*( 2) C 3 - C 4||0.32831||27.15||0.24||0.074|
|n-σ*||LP ( 1) O 17||1.969||BD*( 1) C 14 - C 16||0.05431||2.03||1.23||0.045|
|BD*( 1) C 16 - C 18||0.05908||5.78||0.22||0.032|
|BD*( 1) C 24 - H 28||0.02271||2.57||1.16||0.049|
|n-π*||LP ( 2) O 17||1.88598||BD*( 1) C 14 - C 16||0.05431||14.61||0.81||0.098|
|BD*( 1) C 24 - H 28||0.02271||2.64||0.74||0.04|
|π*-σ*||BD*( 2) C 1 - C 2||0.39784||BD*( 1) C 12 - C 14||0.01462||8.28||33.56||1.036|
|BD*( 1) C 16 - O 17||0.01724||8.11||32.96||1.013|
|BD*( 2) C 27 - C 29||0.33144||21.31||35.54||1.286|
|BD*( 1) C 27 - H 31||0.01409||8.56||36.49||1.099|
|BD*( 1) C 29 - H 32||0.01371||4.11||34.39||0.74|
|π*-π*||BD*( 2) C 3 - C 4||0.32831||BD*( 2) C 16 - O 17||0.27328||0.71||0.02||0.005|
|π*-σ*||BD*( 2) C 12 - C 14||0.17953||BD*( 1) C 29 - H 32||0.01371||2.9||2.02||0.22|
|π*-π*||BD*( 2) C 16 - O 17||0.27328||BD*( 2) C 12 - C 14||0.17953||35.93||0.05||0.077|
|BD*( 2) C 18 - C 20||0.13135||62.96||0.02||0.074|
|π*-π*||BD*( 2) C 18 - C 20||0.13135||BD*( 2) C 27 - C 29||0.33144||1.45||3.19||0.126|
|π*-π*||BD*( 2) C 22 - C 24||0.36662||BD*( 2) C 3 - C 4||0.32831||5.26||0.58||0.084|
|BD*( 1) C 16 - O 17||0.01724||3.86||1.23||0.141|
|BD*( 2) C 23 - C 25||0.31148||6.92||0.62||0.101|
|BD*( 1) C 24 - H 28||0.02271||3.37||1.06||0.121|
|π*-π*||BD*( 2) C 23 - C 25||0.31148||BD*( 2) C 27 - C 29||0.33144||6.83||3.19||0.232|
|π*-σ*||BD*( 2) C 27 - C 29||0.33144||BD*( 1) C 27 - H 31||0.01409||367.87||0.96||1.273|
Molecular orbital and their properties like energy are very useful to the physicists and chemists 44. The HOMO acts as an electron donor and the LUMO is the electron acceptor, and the gap between HOMO and LUMO characterizes the molecular chemical stability. The energy gap between the HOMO and the LUMO MOs is a critical parameter in determining molecular electrical transport properties because it is a measure of electron conductivity. Figure 4 shows, the HOMO-LUMO plots of MPPO. The HOMO is located over the thiophen ring and partly on CH3 group. The LUMO located all over the molecule except methyl group, since methyl group has large electronegativity. The energy gap of MPPO molecule is calculated as 3.4327 eV in gas phase using B3LYP/6-311++G(d,p) basis set.HOMO implies the energy of about -6.0551 eV and LUMO energy is -2.6224 eV. The various physico-chemical parameters are listed in Table 5 and the frontier molecular orbital energies are also listed in Table 6.Table 5. The Physico-chemical properties of MPPO
|Energy gap||3.432 eV|
|Ionization potential (IP)||6.055 eV|
|Electron affinity (EA)||2.622 eV|
|Electrophilicity Index (ω)||2.741|
|Chemical Potential (µ)||4.338|
|Occupancy||Orbital energies (a.u)||Orbital energies (eV)||Kinetic energies (a.u)|
The electronic spectra of the title moleculeis calculated using TD-DFT/B3LYP/6-311++G(d,p) basis set. The electronic transitions, positions of experimental absorption peaks, calculated absorption peaks (λmax), vertical excitation energies, oscillator strengths (f) are listed in Table 7. The observed and UV-Vis absorption maximum of MPPO are shown in Figure 5 In this study, the electronic transition is calculated at three excited states (ES) like as ES1, ES2 and ES3, respectively,in which ES1 appeared at 421.19nm (2.9436eV), ES2 at 398.42 nm (3.1119eV) and ES3 lies at 343.32nm (3.6113eV). In ES3 state there are three transitions from HOMO 64 to Lumo 68, HOMO66 to LUMO68 and HOMO 67 to LUMO 69. The molecule MPPO has sulphur and oxygen as lone pairs, hence π-π* transition is also possible.The electronic absorption calculation shows that the visible absorption maximafor MPPO corresponds to the electron transition between frontier orbitals such as transition from HOMO to LUMO. It should be mentioned here that the low energy absorptions is found at 303 nm and 352 nm belong to the dipole allowed σ-σ* and π-π* transition from HOMO to LUMO, respectively. The intense bands calculated at 398nm can be assigned to high delocalization of π-electrons in the thiophen moiety of MPPO. Similarly the intense band at 343nm can be distributed to the excitation of σ-electrons which are localized on the whole molecular system.Table 7. The electronic transitions of MPPO
|Calculated at B3LYP/6-311++G(d,p)||Oscillator Strength||Calculated Band gap(ev/nm)||Experimental Type Band gap (ev/nm)|
|Excited State-1||Singlet-A(f=0.0000)||2.9436 eV/421.19 nm|
|65 -> 68||0.68512||4.164218|
|Excited State-2||Singlet-A(f=0.6489)||3.6113 eV/343.32 nm||352 π–π*|
|67 -> 68||0.64129||3.432542|
|Excited State-3||Singlet-A(f=0.0507)||3.6113 eV/343.32 nm||303 π –π*|
|64 -> 68||0.28196||4.289929|
|66 -> 68||0.58917||3.960688|
|67 -> 69||0.13934||4.689976|
MEP indicates the net electrostatic effect produced at that point by the total charge distribution (nuclei + electron) of the molecule and correlates with dipole moments, electronegativity, partial charges and chemical reactivity of the molecules. It provides a visual method to understand the relative polarity of the molecule. The different values of the electrostatic potential represented by different colors; red represents the regions of the most negative electrostatic potential, blue represents the regions of the most positive electrostatic potential and green represents the regions of zero potential.The MEP is calculated with B3LYP/6-311++G(d,p) basis set and the MEP plot is shown in Figure 6.
The negative electrostatic potential corresponds to an attraction of the proton by the aggregate ED in the molecule (shades of red), while all the hydrogen atoms have positive electrostatic potential corresponds to the repulsion of the proton by the atomic nuclei (shades of blue). The MEP map shows that the negative region is mainly localized on oxygen and sulphur atoms while the positive region is localized on the surface of all the hydrogen atoms in MPPO. These regions indicate the active charge sites of the molecule. A region of zero potential envelopes the π-system of the aromatic rings 48.
The calculation of atomic charges plays an important role in the application of quantum mechanical calculation to molecular system 49.The Mulliken charges of MPPO was calculated at B3LYP/6-311++G(d,p) level of theory and are listed in Table 8. The mulliken atomic charge plot is shown in Figure 7. The carbon atom (C22) has more positive charge (0.9652 a.u), which isfused linkage with highelectronegative carbonyl group. Similarly,the methyl group carbon atom C8 has maximum negative charge (-0.9825 a.u), sinceit isattached with thiophene ring. In MPPO all the hydrogen atoms have positive charge.Table 8. The Mulliken atomic charges of MPPO
The various thermodynamic parameter are computed by B3LYP/6-311++G(d,p) basis set and are listed in Table 9. The thermodynamic functions like heat capacity (C0p,m), entropy (S0m) and enthalpy changes (ΔH0m) are obtained from the theoretical harmonic frequencies. From Table 10, it can be seen that these thermodynamic functions are increasing with temperature ranging from 100 to 1000 K due to the fact that the molecular vibrational intensities increase with temperature 51. The correlation equations between entropy, heat capacity, enthalpy changes and temperatures were fitted by quadratic formulas and the correspondingare fitting factors (R2) for those thermodynamic properties are 0.9999, 0.9992 and 0.9994, respectively. The corresponding fitting equations are as follows and the correlation graphics are shown in Figure 8.
S0m = 0.98075 + 0.00414T - 3.65824 x 10-5 T2 (R2 = 0.9999)
C0p,m = 6.04475 + 0.02552T +2.25471x 10-5 T2 (R2 = 0.9992)
ΔH0m = 4.00201 + 0.169T + 1.49276x 10-5 T2 (R2 = 0.9994)Table 9. The calculated total energy (a.u), zero point vibrational energies (Kcal/mol), rotational constants (GHz) and entropy (cal/mol K-1) for MPPO
|Zero-point Energy||158.198 (Kcal/Mol)|
|Rotational constants (GHZ)||1.013|
|T (K)||S (J/mol.K)||Cp (J/mol.K)||ddH (kJ/mol)|
The above data can be used to compute other thermodynamic energies according to the well-known relationships of thermodynamic functions and to predict the directions of chemical reactions 52. All the thermodynamic calculations were done in gas phase and they could not be used in solution.
A complete vibrational analysis had been carried out for the first time to the molecule MPPO. The bond parameters and vibrational wavenumbers agreed well with experimental results. The first order hyperpolarizability (β0 = 10.2328x10-30esu) of MPPO was calculated and found to be twenty seven times greater than that of urea and hence the molecule had moderate NLO activity. The considerabledecrease of the lone pair orbital occupancy (1.581e) of LP(2)S5was due to more E(2) energy transfer to anti-bonding orbital C3-C4. The Homo-Lumo energy gap was calculated about 3.4327 eV, which explain the eventual charge transfer occur within the molecule and also enhanced the biological activity. The recorded UV-Vis. spectral values agreed with calculated values. The λmax (352) was assigned to π-π* type. Furthermore, the MEP and Mulliken atomic charges had been calculated and also plotted. The good correlations between the statistical thermodynamics and temperature were also established.