IJOC  Vol.7 No.2 , June 2017
Combined Experimental and Computational Investigation of 2-(2-Hydroxyphenylimino) Phenolic Derivatives: Synthesis, Molecular Structure and NLO Studies
A series of substituted 2-(2-hydroxyphenylimino) phenolic (salen) derivatives (1-4) have been synthesized and their structures of obtained compound were characterized by analytical, FT-IR, UV-Vis and 13C{1H}-NMR experimentally. The geometry structure optimization, frequencies (IR), NMR, electronic character, frontier molecular orbital (HOMO-LUMO) and first static hyperpolarizability (βtot) studies of reported compounds were calculated using DFT with B3LYP/6-311G(d,p) level of theory. The calculated HOMO and LUMO energies showed that charge transfer occurs within the molecule and from the MEP, the molecular stability and bond strength have been explained. In addition to that influence of energy gap (?Egap) between the HOMO-LUMO orbitals on the first static hyperpolarizability (βtot), we calculated the ?Egap for all the salen compounds. These results reveals that the smaller HOMO-LUMO ?Egap is, larger the βtot is.

1. Introduction

Schiff bases are important class of organic compounds have long attracted attention, owing to their remarkable biological and pharmacological properties, such as antibacterial, antiviral, antineoplastic and antimalarial activities. The functional applications of the azomethine group of Schiff base derivatives enable their use in numerous fields, they have

Ø nucleophilic imine group,

Ø an imine carbon that has both electrophilic and nucleophilic character,

Ø configurations isomerism from the presence of the C=N double bond.

These structural features of the Schiff base compounds give its physical and chemical properties [1] [2] [3] . Schiff bases constitute some of the most valuable groups of biomolecules. First reported in 1864 by Hugo Schiff , these compounds gained notoriety due to the ease way of preparation from commercially available inexpensive aldehydes/ketones and primary amines. The azomethine linkage (>C=N-) allows rapid access to vast libraries of structurally diverse molecular hybrids with interesting biological properties, including antifungal, antibacterial, antimalarial, anti-inflammatory, antiviral, antioxidant, pesticidal and in-vitro/ in-vivo inhibitory effects against experimental tumor cells. The electrophilic car- bon and nucleophilic nitrogen in (>C=N-) core confer to Schiff bases the possibility to interact with several nucleophilic and electrophilic biological species, which can lead to enzymes inhibition or DNA replication impairment. Then, Schiff bases are promising as lead compounds for the rational design of novel cytotoxic and cytostatic small molecules with a mechanism of action that may differ from that of clinically approved anticancer agents [3] [4] [5] [6] .

In addition, Schiff base compounds high potential applications in non-linear optical (NLO), optical communication, optical signal processing and transmission, optical data acquisition and storage, optical computing, and especially op- tical limiting effects utilized in the protection of optical sensors and human eyes from high-intensity laser beams [7] [8] [9] [10] . In the present work, we have synthesized certain salen based Schiff base compounds, 2-(2-hydroxyphenyli- mino) phenolic derivatives viz. salcylaldehyde (SA), o-hydroxyacetophenone (AA), o-vanillin (VA) and 2-hydroxy-1-naphthaldehyde (NA), with o-aminophenol [11] [12] [13] [14] and characterized by FT-IR, UV-Vis, 13C{1H}-NMR techniques, Cyclic Voltammetry studies. Moreover, theoretical studies were carried out on the molecular structure using density functional methods (B3LYP) invoking 6-311G(d,p) basis set. The energy of the highest occupied molecular (HOMO) orbital and lowest unoccupied (LUMO) molecular orbital have been predicted.

2. Experimental

Materials and Methods

All the chemicals and solvents used were purified and dried by standard methods. FT-IR spectra were recorded as KBr pellets with a PerkinElmer FT-IR spectrometer in 4000 - 400 cm−1 range. Microanalyses were carried out with a Vario El AMX-400 elemental analyzer at STIC, Cochin University of Science and Technology, Kerala, India. Electronic spectra were recorded in CH3CN as solvent with an Ocean optics spectrophotometer USB 4000. NMR spectra were also recorded on Bruker AVANCE III 500 MHz (AV 500) spectrometer; chemical shifts are expressed in ppm (δ units) relative to TMS signal as internal reference in DMSO.d6.

The Schiff base compounds (Scheme 1) were prepared by the reported literature procedure [11] [12] [13] [14] with modification of the substitutions and the purity of the Schiff bases were checked by TLC. Melting points were recorded with a Inlab, India micro heating table and were uncorrected.

To an ethanolic solution of salcylaldehyde (SA), o-hydroxyacetophenone (AA), o-vanillin (VA) and 2-hydroxy-1-naphthaldehyde (NA) (2 cm3, 20 mmol), with o-aminophenol (2 cm3, 20 mmol) were added with stirring. The mixture was then refluxed for 4 - 6 hours with the controlled conditions about 140˚C - 180˚C. On cooling the solution, a solid compound which separated out was filtered, dried and recrystallized from ethanol/DMSO (75:25). The purity of the ligand was checked by TLC. Yield, 85%, m.p. 143.6˚C - 245˚C.

All the salen compunds were stable at room temperature, non-hygroscopic and insoluble in water partially soluble in methanol, ethanol and soluble in CH2Cl2, CHCl3, DMF, DMSO, CH3CN, etc.

3. Computational Methods

All the computational studies have been carried out with the GAUSSIAN 03W program package [15] . Density functional theory (DFT) method has been applied because of its excellent compromise between computational time and description of electronic correlation and qualitative structure-activity relationships (QSAR) studies [16] . B3LYP, a hybrid functional of the DFT method, which consists of the Becke’s three parameters exact exchange functional B3 [17] combined with the non-local gradient corrected correlation functional of Lee-Yang- Parr (LYP) [18] has been used. The standard triple split valence basis set 6-311G [19] with a set of d, p polarization functions on heavy atoms and hydrogen atoms are used throughout the computational process.

A fully relaxed potential energy scan was carried out against the dihedral angle C2-C1-C7-N1 (Scheme 2) at B3LYP/6-311G(d,p) level, the minimum energy conformations from the energy scan, a further geometry optimization was performed at the same level of theory. Vibrational frequencies of the optimized structures were computed using the same level of theory and thermodynamic

Scheme 1. Preparation of Schiff base compounds.

Scheme 2. Potential energy scan (PES) for salen compounds.

corrections [20] [21] were obtained at 298 K and 1 atm, and added to electronic energies.

Calculated electronic properties such as dipole moment, first static hyper-po- larizability HOMO and LUMO energies, MEP, energy gap, electronic affinity (EA), electronegativity (χ), hardness (η), softness (S), electrophilic index (w) and ionization potential (IP) have been studied for all the four compounds under consideration. The DFT-based reactivity descriptors were obtained from the Equations (1)-(4) [22] [23] [24] [25] which play an important role in many areas of research.

Electronegativity (χ)


Hardness (η)


Softness (S)


The experimental polarizability is obtained as an average polarizability, given

by and the computed average polarizability is also


4. Results and Discussion

4.1. Experimental Results

4.1.1. Elemental Analysis and FT-IR Spectral Data

The elemental analysis data (Table 1) for the all Schiff base compounds are well agreed with the proposed molecular formulae.

The FT-IR spectroscopy is a powerful tool for the assignments of fundamental functional group determinations of organic compounds. In the title compounds frequencies of the functional groups viz., ν>C=N-, νph-OH, νph-C-O, νph-N=C- and ν-N=CH--N=C(CH3)- are of great importance in the infrared spectra. The infrared spectra of all the four Schiff bases (Figure 1) exhibited a strong band in the 1692 - 1642 cm−1 region due to the characteristic of azomethine n(C=N) group. All the title compounds were displayed a band in the region 3399 - 3308 cm−1, which could be due to ν(ph-OH). A strong band observed at 1373 - 1362 cm−1 in


Figure 1. (a) FT-IR spectrum of salen compounds (Experimental); (b) Vibrational frequency calculation by DFT method.

Table 1. Lemental analysis of Salen compounds.

the four compounds has been assigned to phenolic C-O stretching and a band due to ph-N=C- of aminophenol group exhibited in the region 1274 - 1296 cm−1. Further the spectra of all the Schiff bases showed a weak band in the 3020-3069 cm−1 region due to -N=CH-/-N=C(CH3)- groups of NA, SA, VA and AA [14] .

Theoretically all the fundamental vibrations were active in IR. The results showed that the DFT (B3LYP) method applied in this work leads to vibrational wavenumbers which are in good agreements with the experimental data. It is noteworthy that the very important role of vibrational frequencies salen compounds, based on experimental data as well as theoretical calculations, the computed vibrational frequencies in Table 2 and Figure 1 for SA, AA, VA and NA at the DFT level of theory are in acceptable agreement with the experimental data.

4.1.2. UV-Vis Spectra

The UV-Visible spectra of the salen compounds (Figure 2) exhibit mainly two bands. The first band on the higher energy side, the range lmax= 350 - 385 nm, due to the excitation of the p-electrons (p®p* transitions) of the aromatic azo-

Figure 2. UV-Vis Spectra.

Table 2. Vibrational frequencies of salen compounds (Experimental & Theoretical).

methine group. The second band, lmax = 450 - 480 nm range is assigned to an intramolecular charge transfer (ICT transition) involving the salen compounds. This band observed in salicylaldimine compounds of salen derivatives is facilitated by the presence of intramolecular hydrogen bonding between the -OH group and the azomethine nitrogen [26] .

4.1.3. 13C{1H}-NMR Spectra

The formations of Schiff base were conveniently monitored by peak ratios in the 1H-NMR spectra. 1H-NMR spectra of all the four compounds (Table 3 and Figure 3) were taken in DMSO.d6 solvent. The aromatic region was a set of multiplets in the range 6.6 - 8.4 ppm for all the Schiff bases, while the azomethine proton of SA, VA and NA compounds were observed in the range 9.0 - 9.5 ppm. The phenolic-OH protons of all the four Schiff bases were observed as a singlet in the region 13.8-15.6 ppm. 1H-NMR spectra exhibited a strong peak at 3.8 and 3.5 ppm for -OCH3 protons of VA and -CH3 protons of AA salen compounds.

13C-NMR spectral data (Table 3 and Figure 3) were consistent with 1H-NMR spectral data. The methyl carbon of aliphatic substituents, azomethine (>C=N) for all the Schiff bases and for Ph-C-CH3 of AA and Ph-OCH3 of VA compounds were 45 and 56 ppm respectively. The resonance observed at 108 - 178 ppm was assigned to the phenyl group carbon of salen compounds.

4.2. DFT Study

4.2.1. Conformational Analysis and Optimized Geometries

At the minimum energy conformations (Figure 4) obtained from the energy scan, further geometry optimization was performed with the B3LYP/6-311G (d,p) basis set. No geometrical parameter constraint was imposed during the optimization, except those favoring the stabilizing effects due to hydrogen bonding between two adjacent -OH groups. The optimized structures of the most stable conformers of neutral form of salen compounds are shown in Figure 5.

The optimized geometrical parameters are shown in Table 4 (Supplement materials). From the data of bond distances and bond angles, it can be seen that no significant geometrical change has been observed for all the compounds. From the data of dihedral angle in Table 4, it can also be seen that compounds with o-hydroxyl group in A-ring are completely planar, while others have some


Figure 3. (a) 1H-NMR Spectra of salen compounds; (b) 13C-NMR Spectra of salen compounds.

Figure 4. PES scan optimization of salen compound AA.

Table 3. NMR Chemical Shifting of Salen compounds.

degree of deviation from the planarity due to the torsion between A-ring and the plane of aminophenol system B-ring (Scheme 2).

In Table 5 the dipole moment of four compounds at B3LYP/6-311G(d,p) level of computation have been reported. The Dipole moment values lie in the range of 1.7675D to 4.8823D for salen compounds. These values are rather high, reflecting the numerous polarized hydroxyl or carbonyl functions distributed over the structures. The three structures of salen compounds differ only by the presence of methyl at C7, methoxy at C3 and naphthyl groups in the A-ring and there is a little difference in the dipole moment value. The molecular dipole moment is used to measure bond polarities and charge densities in a molecule.


Figure 5. Optimized structure of Salen compounds with the B3LYP/6-311G(d,p) basis set.

Table 4. Geometry parameters of Schiff bases and its radicals at B3LYP/6-311G(d,p) (Bond distances and Bond Angle in Å).

Table 5. Electronic energies for Schiff base neutral molecule and radicals in gas phase at B3LYP/6-311G(d,p) level of theory.


Bond polarity is one of the factors that determine the physicochemical property of molecules. The calculated values of the total dipole moments, which signifies the relatively polarized nature of the systems and they are soluble in polar solvents like CH3CN, DMSO, CHCl3, etc.

The computed NMR chemical shifts for SA, AA, VA and NA at the DFT level of theory are in acceptable agreement with the experimental data. Differences between the calculated and measured values may be a result of solvent interactions.

The 13C-NMR chemical shifts of selected carbons were calculated on the optimized structures of SA, AA, VA and NA using GIAO/DFT method with B3LYP/6-311(d,p) basis set for all atoms. Calculated and measured 13C chemical shifts of selected atoms are numbered in Figure 6.

4.2.2. HOMO-LUMO Orbital Distribution

4.2.3. Spin Density Distribution

The spin density is often considered to be a more realistic parameter which provides a better representation of the reactivity [6] [9] of salen moiety. It should be pointed out that the more delocalized the spin density in the salen compound is, the easier the radical is formed [5] . In the Figure 8, the spin densities of all the radicals mainly distribute on the phenolic oxygen atom and the phenyl ring A


Figure 6. NMR shielding of salen compounds.

HOMO-Neutral SALUMO-Neutral SAHOMO-Neutral AALUMO-Neutral AAHOMO-Neutral VALUMO-Neutral SAHOMO-NALUMO-Neutral NA

Figure 7. HOMO-LUMO diagrame of salen compounds.

Spin Density-SASpin Density-AASpin Density-VASpin Density-NA

Figure 8. Spin density distribution analaysis of salen compounds.

Table 6. Total energies, Frontier orbital energies, Softness, Egap and EI in gas phase calculated at B3LYP/6-311G(d,p) level of theory.

Table 7. All β(a.u.) components and β(tot) × 10−31 (esu) value calculated using DFT level of theory for salencompounds (1 a.u. = 8.6393 × 10−33esu).

and B. Comparing all the radicals, we find that the radicals have almost the same spin density distribution with each other, indicating the presence of additional hydroxyl or o,o-dihydroxyl group on A and B-ring has almost no influence on the spin density distribution [27] [28] .

4.2.4. Molecular Electrostatic Potential Surface (MEP) Study

Electrostatic potential surfaces are mainly used to study the reactive species of electrophilic or nucleophilic attacks/substitution in the chemical reactions, biological process, catalysis and also molecular modeling. Electrostatic potential mapped surface displays the molecular size, shape and potential values. In this study, 3 dimensional surfaces of mapped onto the constant electron density surface is as shown in Figure 9. Different values of electrostatic potential at the surface are represented by different colors [29] .

4.2.5. Computed Non-Linear Optical (NLO) Properties

The NLO response calculation was performed on the optimized geometry using the same level of theory. The first static hyperpolarizability is a third rank tensor that can be described by a 3 × 3 × 3 matrix. The 27 components of the 3D matrix

Figure 9. Molecular Electrostatic Potential surface diagram of AA.

can be reduced to 10 components due to the Kleinman symmetry [27] [30] (βxyy = βyxy = βyyx, βyyz = βyzy = βzyy, ..., likewise other permutations also take same value). It can be given in the lower tetrahedral format. The output from Gaussian 03W provides 10 components of this matrix as βxxx, βxxy, βxyy, βyyy, βxxz, βxyz, βyyz, βxzz, βyzz and βzzz respectively. Many types of hyperpolarizabilities have been discussed in the literature [31] . When reporting a single value of β, one of the common formats is to simply treat the three independent values for β as a quasi-Pythagorean problem and solve for the average β by Equation (5):


The complete equation for calculating the magnitude of the total first statichyperpolarizability from Gaussian 03W output is given as Equation (6):


Since these β values of the first static hyperpolarizability (β) tensors of the output file of Gaussian 03W are reported in atomic units (a.u.), the calculated values were converted into electrostatic units (1 a.u. = 8.6393 × 10−33 esu). The first static hyperpolarizability (β) values of these Schiff bases were calculated under static electronic field by the DFT (B3LYP/(6-311G(d,p)) method. From Table 7, compounds NA, VA and AA have larger the first static hyperpolarizability (βtot) values than those of compound SA. As mentioned above, the DFT results show that the LUMO orbitals of all compounds were obtained from the linear combination of the orbitals of phenyl moiety, the compounds AA, VA and NA with an electron-donating and withdrawing group (CH3, OCH3 and naphthyl) have a larger βtot value compare with compound SA. Hence the compound VA was predicted to have larger NLO property. The energy difference (∆Egap) between the HOMO and the LUMO orbitals 3.5506, 4.3538, 3.8368, 3.4014 eV for the compounds SA, AA, VA, NA respectively has a larger influence on the βtot value [32] [33] . To understand the relationship between the static hyperpolarizabilities and the HOMO-LUMO energy gap, viz. the compound VA (3.4014 eV) lower the HOMO-LUMO energy gap show that the larger the βtot (13.455 a.u) value.

4.3. Conclusions

Schiff bases (1-4) have been successfully synthesized and characterized by elemental analysis, FT-IR, UV-Vis, NMR spectroscopy and cyclic voltammetry. The evaluation of the quantum mechanical studies reveal significant activity in structure optimization, PES, vibrational study, NMR Chemical shifts, spin density, hardness, electronegativity, dipole moment, EHOMO/ELUMO, ∆Egap and Softness which provide the evidence for a very strong positive correlation between experimental and theoretical predictions. These compounds also have towards considerable antioxidant activity along with potential to prevent DNA oxidative damage by free radicals.

Comparison between the four considered molecules indicates compound VA that requires the lowest energy for both H atom and electron transfer mechanisms. This theoretical approach confirms the important role of A and B ring in exhibiting antioxidant properties. Inspection of deprotonation processes of dihydroxyl groups have shown that, π electron delocalization of phenyl ring A and B (evidenced from UV-Vis study) plays a major role in the stabilization of products and thus in the lowering of the associated energies. The variables related to the chemical potential allow classifying salen type of Schiff base compounds that has the tendency to give electrons more than to attract them, which demonstrates biological importance.

In addition to that in order to understand the relationship between the βtot values and the substitution groups of salen compounds, the frontier orbital compositions have been analyzed and the energy gaps between the HOMO and LUMO orbitals were also calculated. The compound VA with electron-donating group (-OCH3) will produce the larger βtot value (13.455 a.u) than lower βtot (10.003, 9.244, 4.359 a.u) for the compounds NA, AA, SA respectively. The energy gaps between the HOMO and LUMO orbitals show that the lower the HOMO-LUMO energy gap, larger the first static hyperpolarizability (βtot).

The findings of this work sustenance the view that some of these synthesized compounds are promising sources of potential drugs that may be efficient as preventive agent(s) in some diseases.

Cite this paper
Nadhiya, V. and Kumaresan, R. (2017) Combined Experimental and Computational Investigation of 2-(2-Hydroxyphenylimino) Phenolic Derivatives: Synthesis, Molecular Structure and NLO Studies. International Journal of Organic Chemistry, 7, 185-217. doi: 10.4236/ijoc.2017.72015.
[1]   Ayman, A., Aziz, A., Salem, M.A.N., Sayed, M.A. and Aboaly, M.M. (2012) Synthesis, Structural Characterization, Thermal Studies, Catalytic Efficiency and Antimicrobial Activity of Some M(II) Complexes with ONO TRIDENTATE SCHIFF BASE N-Salicylidene-o-Aminophenol (saphH2). Journal of Molecular Structure, 1010, 130-138.

[2]   Opstal, T. and Verpoort, F. (2003) Synthesis of Highly Active Ruthenium Indenylidene Complexes for Atom-Transfer Radical Polymerization and Ring-Opening- Metathesis Polymerization. Angewandte Chemie International Edition, 42, 2876-2886.

[3]   da Silva, C.M., da Silva, D.L., Modolo, L.V., Alves, R.B., Martins, C.V.B., de Resende, M.A. and de Fátima, A. (2011) Schiff Bases: A Short Review of Their Antimicrobial Activities. Journal of Advanced Research, 2, 1-8.

[4]   Sztanke, K., Maziarka, A., Osinka, A. and Sztanke, M. (2013) An Insight into Synthetic Schiff Bases Revealing Antiproliferative Activities in Vitro. Bioorganic & Medicinal Chemistry, 21, 3648-3666.

[5]   Hameed, A., Al-Rashida, M., Uroos, M., Ali, S.A. and Khan, K.M. (2017) Schiff Bases in Medicinal Chemistry: A Patent Review (2010-2015). Expert Opinion on Therapeutic Patents, 27, 63-79.

[6]   Sirajuddin, M., Uddin, N., Ali, S. and Tahir, M.N. (2013) Potential Bioactive Schiff Base Compounds: Synthesis, Characterization, X-Ray Structures, Biological Screenings and Interaction with Salmon Sperm DNA. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 116, 111-121.

[7]   Govindarasu, K. and Kavitha, E. (2014) Vibrational Spectra, Molecular Structure, NBO, NMR, UV, First Order Hyperpolarizability, Analysis of (S)-(-)-N-(5-Nitro-2- Pyridyl) Alaninol by Density Functional Theory. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 127, 498-510.

[8]   Chandra, S. and Saneetika, X. (2004) EPR, Magnetic and Spectral Studies of Copper(II) and Nickel(II) Complexes of Schiff Base Macrocyclic Ligand Derived from Thiosemicarbazide and Glyoxal. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 60, 147-153.

[9]   Finkel, T. and Holbrook, N.J. (2000) Review Article Oxidants, Oxidative Stress and the Biology of Ageing. Nature, 408, 239-247. https://doi.org/10.1038/35041687

[10]   Zhou, L., Cai, P., Feng, Y., Cheng, J., Xiang, H., Liu, J., Wu, D. and Zhou, X. (2012) Synthesis and Photophysical Properties of Water-Soluble Sulfonato-Salen-Type Schiff Bases and Their Applications of Fluorescence Sensors for Cu2+ in Water and Living Cells. Analytica Chimica Acta, 735, 96-106.

[11]   Priyarega, S., Prabhakaran, R., Aranganayagam, K.R., Karvembu, R. and Natarajan, K. (2007) Synthetic and Catalytic Investigations of Ruthenium(III) Complexes with Triphenylphosphine/Triphenylarsine and Tridentate Schiff Base. Applied Organometallic Chemistry, 21, 788-793.

[12]   Mazzi, U., Refosco, F., Tisato, F., Bandoli, G. and Nicolini, M. (1986) Mixed Tri- and Bi-Dentate Schiff-Base Complexes of Technetium(V) and Rhenium(V). The Crystal Structure of N-(2-Oxidophenyl)Salicylideneiminato-NOO’)Oxo(8-Quino-linolato-NO)Technetium(V). Journal of the Chemical Society, Dalton Transactions, No. 8, 1623-1628.

[13]   Tisato, F., Refosco, F., Mazzi, U., Bandolin, U. and Nicolini, M. (1987) Technetium(V) and Rhenium(V) Complexes with N-(2-Mercaptophenyl)Salicylideneiminate. Crystal Structure of Chlo-ro(oxo)[N-(2-Sulphidophenyl)Salicylideneiminato-NOS]Technetium(V). Journal of the Chemical Society, Dalton Transactions, No. 7, 1693-1699. https://doi.org/10.1039/dt9870001693

[14]   Priyarega, S., Tamizh, M.M., Karvembu, R., Prab-hakaran, R. and Natarajan, K. (2011) Synthesis, Spectroscopic Characterization and Catalytic Oxidation Properties of ONO/ONS Donor Schiff Base Ruthenium(III) Complexes Containing PPh3/ AsPh3. Journal of Chemical Sciences, 123, 319-325.

[15]   Frisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Zakrzewski, V.G., Montgomery Jr., J.A., Stratmann, R.E., Burant, J.C., Dapprich, S., Millam, J.M., Daniels, A.D., Kudin, K.N., Strain, M.C., Farkas, O., Tomasi, J., Barone, V., Cossi, M., Cammi, R., Mennucci, B., Pomelli, C., Adamo, C., Clifford, S., Ochterski, J., Petersson, G.A, Ayala, P.Y., Cui, Q., Morokuma, K., Malick, D.K., Rabuck, A.D., Raghavachari, K., Foresman, J.B., Cioslowski, J., Ortiz, J.V., Baboul, A.G., Stefanov, B.B., Liu, G., Liashenko, A., Piskorz, P., Komaromi, I., Gomperts, R., Martin, R.L., Fox, D.J., Keith, T., Al-Laham, M., Peng, C.Y., Nanayakkara, A., Gonzalez, C., Challacombe, M., Gill, P.M.W., Johnson, B., Chen, W., Wong, M.W., Andres, J.L., Gonzalez, C., Head-Gordon, M., Reploge, E.S. and Pople, J.A. (2005) Gaussian 03, Revision D.01. Gaussian, Inc., Pittsburgh, PA.

[16]   Kohn, W. (1999) Nobel Lecture: Elec-tronic Structure of Matter—Wave Functions and Density Functionals. Reviews of Modern Physics, 71, 1253.

[17]   Becke, A.D. (1998) Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Physical Review A, 38, 3098.

[18]   Lee, C., Yang, W. and Parr, R.G. (1988) Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Physical Review B, 37, 785.

[19]   Frisch, M.J., Pople, J.A. and Binkley, J.S. (1984) Self-Consistent Molecular Orbital Methods 25. Supplementary Functions for Gaussian Basis Sets. The Journal of Chemical Physics, 80, 3265.

[20]   Hehre, W.J., Radom, L., Schleyer, P.V.R. and Pople, J.A. (1986) Ab Initio Molecular Orbital Theory. Wiley, New York.

[21]   Kohn, W., Becke, A.D. and Parr, R.G. (1996) Density Functional Theory of Electronic Structure. The Journal of Physical Chemistry, 100, 12974-12980.

[22]   Parr, R.G. and Pearson, R.G. (1983) Absolute Hardness: Companion Parameter to Absolute Electronegativity. Journal of the American Chemical Society, 105, 7512-7516.

[23]   Parr, R.G., Szentpaly, L.V. and Liu, S. (1999) Electrophilicity Index. Journal of the American Chemical Society, 121, 1922-1924. https://doi.org/10.1021/ja983494x

[24]   Chattaraj, P.K., Lee, H. and Parr, R.G. (1991) HSAB Principle. Journal of the American Chemical Society, 113, 1855-1856. https://doi.org/10.1021/ja00005a073

[25]   Sadasivam, K. and Kumaresan, R. (2011) Antioxidant Behavior of Mearnsetin and Myricetin Flavonoid Compounds—A DFT Study. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 79, 282-293.

[26]   Gabr, A.A. (1990) Spectrophotometric Studies on Some Schiff Bases deriveD from Benzidine. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 46, 1751-1757.

[27]   Wright, J.S., Johnson, E.R. and Dilabio, G.A (2001) Predicting the Activity of Phenolic Antioxidants: Theoretical Method, Analysis of Substituent Effects, and Application to Major Families of Antioxidants. Journal of the American Chemical Society, 123, 1173-1183.

[28]   Sadasivam, K. and Kumaresan, R. (2011) Theoretical Investigation on the Antioxidant Behavior of Chrysoeriol and Hispidulin Flavonoid Compounds—A DFT Study. Computational and Theoretical Chemistry, 963, 227-235.

[29]   Politzer, P. and Truhlar, D.G., Eds. (1981) Plenum, New York.

[30]   Kleinman, D.A. (1962) Nonlinear Dielectric Polarization in Optical Media. Physical Review, 126, 1977.

[31]   Pipek, J. and Mezey, P.Z. (1989) A Fast Intrinsic Localization Procedure Applicable for ab Initio and Semiempirical Linear Combination of Atomic Orbital Wave Functions. The Journal of Chemical Physics, 90, 4916. https://doi.org/10.1063/1.456588

[32]   Wang, P., Zhu, P., Wang, C. and Ye, C. (1999) Theoretical Investigation and Molecular Design of Pyrazine Derivatives with Large Hyperpolarizabilities (β). Journal of Molecular Structure: THEOCHEM, 459, 155-162.

[33]   Tang, G.D., Jiang, Z.J., Li, R.Q., Zhang, J.F., Zhang, Y. and Zhang, C. (2009) DFT Studies on Nonlinear Optical Properties of Neutral Nest-Shaped Heterothiometallic [MOS3Py5Cu3X] (M = Mo, W; X = F, Cl, Br, I) Clusters. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 74, 228-232.