DFT Study of the Structure, Reactivity, Natural Bond Orbital and Hyperpolarizability of Thiazole Azo Dyes

The structure, reactivity, natural bond orbital (NBO), linear and nonlinear optical (NLO) properties of three thiazole azo dyes (A, B and C) were monitored by applying B3LYP, CAM-B3LYP and ωB97XD functionals with 6-311++G** and aug-cc-pvdz basis sets. The geometrical parameters, dipole moments, HOMO-LUMO (highest occupied molecular orbital, lowest unoccupied molecular orbital) energy gaps, absorption wavelengths and total hyperpolarizabilities were investigated in carbon tetrachloride (CCl4) chloroform (CHCl3), dichloromethane (CH2Cl2) and dimethlysulphoxide (DMSO). The donor methoxyphenyl group deviates from planarity with the thiazole azo moiety by ca. 38°; while the acceptor dicyanovinyl, indandione and dicyanovinylindanone groups diverge by ca. 6°. The HOMOs for the three dyes are identical. They spread over the methoxyphenyl donor moiety, the thiazole and benzene rings as π-bonding orbitals. The LUMOs are shaped up by the nature of the acceptor moieties. The LUMOs of the A, B and C dyes extend over the indandione, malononitrile and dicyanovinylindanone acceptor moieties, respectively, as π-antibonding orbitals. The HOMO-LUMO splittings showed that Dye C is much more reactive than dyes A and B. Compared to dyes A and B, Dye C yielded a longer maximum absorption wavelength because of the stabilization of its LUMOs relative to those of the other two. The three dyes show solvatochromism accompanied by significant increases in hyperpolarizability. The enhancement of the total hyperpolarizability of C compared to those of A and B is due to the cumulative action of the long π-conjugation of the indanone ring and the stronger electron-withdrawing ability of the dicyanovinyl moiety that form the dicyanovinylindanone acceptor group. These findings are facilitated by a natural bond orbital (NBO) technique. The very high total hyperpolarizabilities of the three dyes define their potent nonlinear optical (NLO) behaviour.


Introduction
In recent decades, researchers have become interested in the fabrication of metal surfaces that are functionalized with organic chromophores for tailoring their electrical, magnetic, optical and electroptical properties [1][2][3][4][5]. The adjustment of their key surface properties could be achieved by proper molecular design and/or by fine control of their film structure at the molecular level [6]. The proper chromophores suitable for these criteria could occur as organic second-order nonlinear optical (NLO) compounds having an electron donor (D) and electron acceptor (A) separated by a π-conjugated spacer (D-π-A) [7]. These structural rearrangements facilitate asymmetrical ground-state charge transfer emanating from the donor group (D) through the π-linker to the acceptor moiety (A) under the influence of an electric field [8]. The last few decades have witnessed the fabrication of thermally and photochemically stable NLO benzenoid chromophores with potent hyperpolarizabilities [9][10][11]. Toward the end of last century, a series of single-substituted thiazole ring and donor-acceptor thiazole-containing chromophores were synthesized, characterized and their superior hyperpolarizabilities compared to oxazoles, imadazoles and thiophenes were obtained [12,13]. In 2004, semi-empirical and ab initio calculations were performed on a series of push-pull π-conjugated styryl benzothiazoles dyes [14]. At the end of last decade, benzothiazolium salts having dimethylamino and diphenylamino electron-donating and nitro or cyano electron-withdrawing groups were synthesized and studied for their NLO properties both theoretically and experimentally by using Hyper-Rayleigh scattering [15]. In 2011, the linear and nonlinear properties of a series of triphenylamine-derived benzothiazoles were tuned by incorporating some electron-withdrawing groups [16]. Extremely polarizable NLO chromophores were manufactured also by incorporating the five-membered heteroaromatic thiazole, pyrrole or thiophene rings [17][18][19] and/or thiazole-annulated heteroaromatics, such as benzobisthiazole [20]. It has been established theoretically and experimentally that the location of the heteroatoms within heteroaromatic rings and the relative position of electron-donor and electron-acceptor groups could play vital roles in dictating the NLO activity of these chromophores [14][15][16][20][21][22][23].
El-Shishtawy et al. [31] have synthesized and investigated, experimentally, the structure and NLO properties of three thiazole azo dyes with a lateral methoxyphenyl donor group coupled with indandione, malononitrile and dicyanovinylindanone acceptor moieties. They showed that the three dyes are thermally stable and have strong absorption wavelengths whose maxima are dictated by the nature of the acceptor groups. In addition, they have also demonstrated that the dye with the dicyanovinylindanone acceptor group is an extremely promising NLO devise with a nonlinearity µβ coefficient amounting to three times those of the other two dyes. In this paper, we endeavor to further complement their findings theoretically and computationally. Density functional theory (DFT) and Time-dependent density functional theory (TD-DFT) with traditional hybrid and long-range corrected (LC) functionals and moderate basis sets will be used to monitor their photochromic and NLO properties. The intramolecular charge transfer (ICT) of the three dyes, from the methoxyphenyl donor toward the three-acceptor groups, will be studied by natural bond orbital (NBO) technique.

Geometrical Analysis
The structure of the three dyes A, B and C were optimized to global minima using B3LYP, CAM-B3LYP and ωB97XD functionals with aug-cc-pvdz and 6-311++G** basis sets (See Figure 1). Some selected bond lengths and dihedral angles of these optimization procedures together with those extracted from the crystal structures of dyes B and C [31] were listed in Table 1. In excellent agreement with the available crystallographic data [31], all the geometry optimizations of the dyes indicated that the donor methoxyphenyl group deviates from planarity with the thiazole azo moiety by ca. 38 • ; while the acceptor dicyanovinyl and dicyanovinylindanone groups diverge by ca. 6 • . As Table 2 shows, an overall excellent agreement between the measured [31] and calculated values is met. This is indicated by the small absolute errors. These findings agree satisfactorily with similar theoretical computations [32]. The N6-N7 crystallographic [31] and calculated bond lengths of dyes B and C are in overall good agreement with each other, having average errors maxima of 0.031 and 0.038 Å, respectively. Dye B crystallographic S4-C1 bond length of 1.744 Å was exactly reproduced by ωB97XD/6-311++G** level of theory but with an error of 0.012 Å in dye C. Generally, the bond lengths computed by the different DFT functionals with the aug-cc-pvdz basis set are slightly longer and much closer to the crystallographic values [31] compared to those estimated by the 6-311++G** basis set.
The multiple bond character of C1-C33 bond with an average value of 1.426 Å is indicated by being shorter than that of 2-methylthiophene [33] by ca. 0.079 Å. That is, C1 and C33 are nearly sp 2 hybridized. The sp 2 hybridization environment around the bridge between the acceptor groups and the thiazole azo moiety is also facilitated by the value of the angle C33-C1-C2 being ca. 124 • . An unexpected agreement between the crystallographic [31] and the computed dihedral angles occurred satisfactorily, although the two measurements were carried out in different phases [34]. hybridized. The sp 2 hybridization environment around the bridge between the acceptor groups and the thiazole azo moiety is also facilitated by the value of the angle C33-C1-C2 being ca. 124°. An unexpected agreement between the crystallographic [31] and the computed dihedral angles occurred satisfactorily, although the two measurements were carried out in different phases [34].
Dye A Dye B

Frontier Molecular Orbitals (FMOs)
The frontier molecular orbitals (FMOs) are formed mainly by the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). The FMOs participate strongly in investigating the electrical and chemical properties of substrates [35]. They affect these properties through forming their polarities, together with their abilities for absorbing light. This means that they act as donor and acceptor orbitals [36], respectively. The HOMO and LUMO orbitals, the chemical hardness (η), the electronic chemical potential (μ) and the global electrophilicity index (ω) using the elected DFT functionals and basis sets are listed in Table 3. Figure 2 depicts the HOMO and LUMO orbitals for the studied dyes (A, B and C). On the one hand, the HOMOs for the three dyes are identical. They spread over the methoxyphenyl donor moiety, the thiazole and the benzene rings as π-bonding orbitals. On the other hand, the nature of the acceptor moieties dictated the energy and the shape of the LUMOs. The LUMOs of the A, B and C dyes extend over the indandione, malononitrile and dicyanovinylindanone acceptor moieties, respectively, as π-antibonding orbitals. The strengths of these acceptor groups are reflected in the stabilization of the LUMOs [31] in excellent agreement with our results shown in Table 3. The kinetic stability of the three dyes can be monitored by the HOMO-LUMO energy gaps [37]. This means that smaller HOMO-LUMO splittings lead to lower kinetic stability and higher chemical reactivity. The reverse is equally true. In addition, it is energetically favorable to move electrons easily between high-lying HOMOs and low-lying LUMOs [38]. All our results showed that Dye C is much more reactive than dyes A and B. However, dyes A and B are of comparable reactivity. With the exception of ωB97XD/6-311++G** model chemistry, all other elected levels of theory revealed that Dye A is slightly more reactive than Dye B.

Frontier Molecular Orbitals (FMOs)
The frontier molecular orbitals (FMOs) are formed mainly by the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). The FMOs participate strongly in investigating the electrical and chemical properties of substrates [35]. They affect these properties through forming their polarities, together with their abilities for absorbing light. This means that they act as donor and acceptor orbitals [36], respectively. The HOMO and LUMO orbitals, the chemical hardness (η), the electronic chemical potential (µ) and the global electrophilicity index (ω) using the elected DFT functionals and basis sets are listed in Table 3. Figure 2 depicts the HOMO and LUMO orbitals for the studied dyes (A, B and C). On the one hand, the HOMOs for the three dyes are identical. They spread over the methoxyphenyl donor moiety, the thiazole and the benzene rings as π-bonding orbitals. On the other hand, the nature of the acceptor moieties dictated the energy and the shape of the LUMOs. The LUMOs of the A, B and C dyes extend over the indandione, malononitrile and dicyanovinylindanone acceptor moieties, respectively, as π-antibonding orbitals. The strengths of these acceptor groups are reflected in the stabilization of the LUMOs [31] in excellent agreement with our results shown in Table 3. The kinetic stability of the three dyes can be monitored by the HOMO-LUMO energy gaps [37]. This means that smaller HOMO-LUMO splittings lead to lower kinetic stability and higher chemical reactivity. The reverse is equally true. In addition, it is energetically favorable to move electrons easily between high-lying HOMOs and low-lying LUMOs [38]. All our results showed that Dye C is much more reactive than dyes A and B. However, dyes A and B are of comparable reactivity. With the exception of ωB97XD/6-311++G** model chemistry, all other elected levels of theory revealed that Dye A is slightly more reactive than Dye B. Table 1. Some selected bond lengths (Å) and dihedral angles (degrees) of the optimized structures of A (top 6 lines), B (line 7 to 13) and C (bottom 7 lines) dyes which have been estimated by using different DFT functionals with 6-311++G** and aug-cc-pvdz basis sets. The crystal structures for B and C dyes are listed for comparison purposes.  Table 2. Absolute and average absolute errors of some selected bond lengths (in Å) and torsional angles (in deg.) for the thiazole azo dyes B (top 7 lines) and C (bottom 7 lines), which were estimated by using different DFT functionals with 6-311++G** and aug-cc-pvdz basis sets as compared to the experimental crystal data 1 .
The crystal structure of A is not available.  Table 3. The HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) orbitals energies (eV), the energy gaps (E.G./eV), the dipole moments (D.M./Debye) the chemical hardness (η/eV), electronic chemical potential (µ/eV), the global electrophilicity index (ω/eV) and the total hyperpolarizabilities (β tot /au) for gas-phase dyes A, B and C, which were calculated by applying B3LYP, CAM-B3LYP and ωB97XD functionals with 6-311++G** and aug-cc-pvdz basis sets. For comparison, the values for gas-phase p-nitroaniline (pNA) are given.

B3LYP
CAM-B3LYP ωB97XD 6-311++G** aug-cc-pvdz 6-311++G** aug-cc-pvdz 6-311++G** aug-cc-pvdz Likewise, the chemical hardness (η) is useful in studying the stability and reactivity of compounds. It is formulated in terms of the energies of the HOMOs and LUMOs [40]: This formula indicates that soft compounds have small chemical hardness, while hard ones have large chemical splittings. In other words, soft compounds have small excitation energies, that is, their electron densities are easily altered, while hard ones have large excitation energies or their electronic densities are difficult to modify [40]. Referring back to Table 3, it can be seen that Dye C is the softest, while Dye B is the hardest, among the three dyes.
Moreover, the electronic chemical potential (µ) shows the escaping tendency of electrons in compounds [41] and given by Equation (2) [42,43]: As Table 3 shows, the order of the electronic chemical potential for the three dyes is as follows: B > C > A. The global electrophilicity index (ω) estimates the stabilizing energy when a surrounding environment supplies a chemical entity with an additional electronic charge. The index (ω) relates to the electronic chemical potential (µ) and the chemical hardness (η) through Equation (3) [41]: The values of the global electrophilicity indexes shown in Table 3 for the three dyes indicate that Dye A is the strongest nucleophile, while Dye C is the strongest electrophile among the three substrates.
The values of the global electrophilicity indexes shown in Table 3 for the three dyes indicate that Dye A is the strongest nucleophile, while Dye C is the strongest electrophile among the three substrates.

UV-Visible Spectral Analysis
It has become general knowledge that the π→π* and n→π* electronic transitions in πconjugated organic compounds lead to UV-Vis. spectra [44]. They are due to electron motions between the FMOs; like the promotion of an electron from the HOMO to the LUMO. Dyes A, B and C have many π-bonds in the thiazole azo, the two benzene and indane rings, together with the nitrogen, oxygen and sulphur atoms lone pairs. The experimental and theoretical maximum absorption wavelength (λmax) for dyes A, B and C in chloroform are depicted in Table 4. The estimated values are computed by the time-dependent density functional theory (TD-DFT) [45] procedure using the polarizable continuum model (PCM) method [46] with TD-B3LYP/6-311++G** model chemistry. The experimental maximum wavelength bands for A, B and C dyes in chloroform showed up at 623, 619 and 687 nm, respectively [31]. The longest wavelength due to Dye C is ca. 10% longer than the average wavelength of the other two dyes. This difference in trend was nearly reproduced theoretically using the elected levels of theory. That is, our calculated maximum wavelengths, in chloroform, of Dye C are longer than those of dyes A and B by 5%-12%. In excellent agreement with El-Shishtawy et al. [31], the longer maximum wavelength of Dye C resulted from its more stabilized LUMO relative to the LUMOs of dyes A and B. It can be concluded that the extra stability of the LUMO of Dye C (see Table 3) originates from the stronger electron-withdrawing ability of the dicyanovinylindanone group compared to those of the indandione and dicyanovinyl moieties [31]. The UV-Vis. spectra of the three dyes (A, B and C) are simulated in CCl4, CHCl3, CH2Cl2 and DMSO solvents applying TD-CAM-B3LYP/6-311++G** level of theory. The results are shown in Table 5. The

UV-Visible Spectral Analysis
It has become general knowledge that the π→π* and n→π* electronic transitions in π-conjugated organic compounds lead to UV-Vis. spectra [44]. They are due to electron motions between the FMOs; like the promotion of an electron from the HOMO to the LUMO. Dyes A, B and C have many π-bonds in the thiazole azo, the two benzene and indane rings, together with the nitrogen, oxygen and sulphur atoms lone pairs. The experimental and theoretical maximum absorption wavelength (λ max ) for dyes A, B and C in chloroform are depicted in Table 4. The estimated values are computed by the time-dependent density functional theory (TD-DFT) [45] procedure using the polarizable continuum model (PCM) method [46] with TD-B3LYP/6-311++G** model chemistry. The experimental maximum wavelength bands for A, B and C dyes in chloroform showed up at 623, 619 and 687 nm, respectively [31]. The longest wavelength due to Dye C is ca. 10% longer than the average wavelength of the other two dyes. This difference in trend was nearly reproduced theoretically using the elected levels of theory. That is, our calculated maximum wavelengths, in chloroform, of Dye C are longer than those of dyes A and B by 5%-12%. In excellent agreement with El-Shishtawy et al. [31], the longer maximum wavelength of Dye C resulted from its more stabilized LUMO relative to the LUMOs of dyes A and B. It can be concluded that the extra stability of the LUMO of Dye C (see Table 3) originates from the stronger electron-withdrawing ability of the dicyanovinylindanone group compared to those of the indandione and dicyanovinyl moieties [31]. The UV-Vis. spectra of the three dyes (A, B and C) are simulated in CCl 4 , CHCl 3 , CH 2 Cl 2 and DMSO solvents applying TD-CAM-B3LYP/6-311++G** level of theory. The results are shown in Table 5. The interactions between these solvents and the dyes render them fairly solvatochromic, with red shifts of 0.063, 0.065 and 0.039 eV in the π→π* bands of A, B and C, respectively, on moving from the less polar CCl 4 to the highly polar DMSO. As Table 5 shows, the HOMOs and LUMOs of the dyes are further stabilized relative to the polarity of the solvents, with the latter being greatly affected. It seems, then, that the solvatochromic behaviour could have resulted mainly from the relatively strong interactions between these solvents and the indandione, malononitrile and dicyanovinylindanone acceptor moieties [47]. Table 4. The Dipole moments (D.M./Debye), the energy gaps (E.G./eV), the total hyperpolarizabilities (β tot × 10 −28 /esu) and the maximum absorption wavelength (λ max /nm) for the chloroform-solvated dyes (A, B and C) which were estimated by utilizing B3LYP, CAM-B3LYP and ωB97XD functionals with 6-311++G** and aug-cc-pvdz basis sets. Some experimental related values of the three dyes in chloroform are also listed for comparison purposes.

Nonlinear Optical (NLO) Properties
Experimentalists and theoreticians have adopted many different methods and conventions for the determination of hyperpolarizability. This situation has brought about some form of ambiguity when theoretical data are compared with measured ones [48]. In this contribution, we compute the total hyperpolarizability, β tot , by the relation: where The total hyperpolarizabilities in atomic units (a.u.) are related to the electrostatic units (esu) by the relation: 1 a.u. = 8.6393 × 10 −33 esu.
The computed dipole moments, energy gaps and the total hyperpolarizabilities for gas-phase dyes A, B and C, together with those of gas-phase p-nitroaniline (pNA) using B3LYP/6-311++G** model chemistry, are listed in Table 3. pNA is selected for comparison purposes because it is a typical example of a donor-acceptor charge-transfer species with high hyperpolarizability that was exposed to both experimental [39] and theoretical [49] investigations. In Table 3, it can be easily seen that dyes A and B have nearly equal total hyperpolarizabilities, which amount to 14-and 12-fold that of pNA, respectively, using the same level of theory. In contrast, the total hyperpolarizabilities of Dye C amount to twice those of dyes A and B and up to ca. 25-fold that of pNA. In Table 4 are listed the estimated dipole moments, energy gaps and total hyperpolarizabilities of the chloroform-solvated dyes (A, B and C) using the elected levels of theory together with their experimental counterparts in the same solvent. Apparently, the B3LYP functional overestimated the total hyperpolarizabilities of both the gas-phase and solution substrates compared to those obtained from CAM-B3LYP and ωB97XD counterparts, while the latter functionals yielded comparable values [49]. In addition, the total hyperpolarizabilities of the solvated dyes are ca. three-fold greater than their gas-phase peers [50]. This phenomenon could be related to the increase of dipole moment change between the ground and excited states [15] (cf . Tables 3 and 4). Our calculated total hyperpolarizabilities of the three dyes in chloroform are in good match with the measured ones in the same solvent [31]. In addition, the estimated total hyperpolarizabilities of the three dyes in CCl 4 , CHCl 3 , CH 2 Cl 2 and DMSO solvents using CAM-B3LYP/6-311++G** level of theory, listed in Table 5, are dictated by the polarity of these solvents, that is, the solvatochromic behaviours of the dyes have brought about pronounced NLO characteristics [47]. As Tables 3-5 show, the values of total hyperpolarizabilities of the three dyes (A, B and C) are inversely proportional to the magnitude of the energy gaps [51][52][53]. This is because narrow energy gaps enhance intramolecular charge transfer and hence produce larger total hyperpolarizabilities. In addition, the total hyperpolarizability of Dye A is slightly more than that of Dye B, in spite of the fact that the dipole moment of the latter is ca. twice that of the former. The enhanced total hyperpolarizability of dyes A and C could be due to the longer extension of the π-conjugation of the indandione and dicyanovinylindanone moieties compared to that of the dicyanovinyl group, while Dye B competes with them through the stronger electron-withdrawing ability of the dicyanovinyl moiety [31].

Natural Bond Orbital (NBO) Analysis
The second-order perturbation energies (E (2) ) are commonly used as a quantitative tool for investigating bonding and antibonding interactions through natural bond orbital (NBO) technique [54][55][56][57]. The second-order perturbation energies (E (2) ) are given by Equation (6): where the term F 2 (i, j) represents the off-diagonal matrix elements, q i gives the donor orbital occupancy, ε i and ε j estimate the donor and acceptor orbital energies, respectively. The quantities of E (2) evaluate the magnitude of interaction between the donor and acceptor orbitals. They, therefore, show the extent of delocalization throughout the chemical species. Table 6 lists the extremely influential interactions between the bonding or lone pair Lewis Type NBO occupied orbitals with antibonding Non-Lewis NBO unoccupied orbitals of dyes A, B and C. They were estimated by using HF/6-31+G*//B3LYP/6311++G** level of theory. The n→π* and π→π* interactions define the charge transfer from the donor methoxyphenyl group towards the acceptor moieties and contribute effectively to the stabilization of the three dyes.  Table 6. Some selected most efficacious second-order perturbation (E (2) ) assessment of the hyperconjugative energies (kcal/mol) that trace the charge transfer from the methoxyphenyl donor group to the indandione, dicyanovinyl and dicyanovinylindanone acceptor moieties of the thiazole azo dyes A, B and C, respectively. They were computed applying HF/6-31+G*//B3LYP/6311++G** level of theory. The huge charge transfer that occurred in Dye C explains its larger total hyperpolarizability compared to those of A and B. On the one hand, this is because Dye C has both the long π-conjugation extension of the indanone moiety and the strong electron-withdrawing ability of the dicyanovinyl group. On the other hand, dyes A and B have either the long π-conjugation or the strong electron-withdrawing potency. The large dipole moment and the strong electron-withdrawing ability of Dye B yielded total hyperpolarizabilities comparable to those of Dye A that has a smaller dipole moment. The competitiveness of Dye A originates from the long π-conjugation extension associated with the indandione ring [31]. The intramolecular charge transfer associated with the π C1-C2 →π* C35-C50 interaction is effectively shown by the C1-C33 bonds being multiply bonded (1.419-1.431 Å). All these findings are in excellent agreement with those reported by El-Shishtawy et al. [31].

Computational Details
The Gaussian09 Suite of programs [58] were used to perform the quantum mechanical molecular orbital calculations of the three dyes, A, B and C. The GaussView [59] and Chemcraft [60] softwares were applied to monitor the structure and properties of the studied dyes. A number of density functional theory (DFT) [61] functionals with the triple zeta and polarization functions at the hydrogen and carbon atoms (6-311++G**) [62][63][64] and augmented correlation-consistent polarization with double zeta functions (aug-cc-pvdz) [65] basis sets were tested. They were used to optimize the geometrical structures of those dyes for investigating their reactivity, hyperpolarizability, linear and nonlinear optical (NLO) properties. The DFT functionals comprise the Becke, three-parameter, Lee-Yang-Parr exchange-correlation hybrid functional (B3LYP) [66,67], the Coulomb-attenuating method that includes the hybrid features of B3LYP and the long-range correction (CAM-B3LYP) [68], and the long-range corrected (LC) hybrid density functional with empirical atom-atom dispersion corrections (ωB97XD) [69]. The Gaussian NBO software [56,70] was applied to perform a natural bond orbital (NBO) analysis for the three dyes.
The time-dependent density functional theory (TD-DFT) [71] and the polarizable continuum model (PCM) [72] technique were applied to compute the UV-Vis. spectra of the CCl 4 , CHCl 3 , CH 2 Cl 2 and DMSO solvated dyes using the elected levels of theory. The frontier molecular orbitals (FMO) and the global chemical reactivity descriptors of the three dyes were examined by using all tested model chemistries

Conclusions
The geometry, reactivity, linear and nonlinear optical behaviour of three donor-acceptor thiazole azo dyes were monitored by DFT calculations. The donor methoxyphenyl group and the acceptor indandione, malononitrile and dicyanovinylindanone moieties, incorporated into dyes A, B and C, respectively, were not coplanar with the thiazole azo spacer group. The HOMO-LUMO analysis showed that Dye C is more reactive than both A and B. This property ties up nicely with the longer absorption wavelength of the former. This linear behaviour is assigned to both the long π-conjugation extension of the indanone ring and the strong electron-withdrawing ability of the dicyanovinyl moiety. The three dyes showed solvatochromism on moving from CCl 4 to DMSO solvents. The solvatochromic behaviours were reflected in pronounced NLO properties. The calculated total hyperpolarizabilities of Dye C were more than two-fold those obtained from A and B and ca. 25-fold of that from pNA. An NBO investigation supported these results. The enhancement of the linear and nonlinear behaviour of Dye C originates, in part, from the π C33-C50 →π* C51-C63 transition, which stabilized this dye by 37.99 kcal/mol. The delocalization energies of A and B are nearly comparable, with the former having a slightly higher amount. This is because the latter is devoid of the π C33-C50 →π* C51-O64 and π C33-C50 →π* C52-O63 transitions, which stabilized the former by 67.35 kcal/mol, while the former is destitute of the π C33-C50 →π* C51-N53 and π C33-C50 →π* C52-N54 transitions, which availed 59.17 kcal/mol for the latter. All our theoretical findings are in excellent agreement with experiment [31] and affirm the use of the three dyes as potential NLO devices.