Synthesis and structural properties of 2-((10-alkyl-10H-phenothiazin-3-yl)methylene)malononitrile derivatives; a combined experimental and theoretical insight

Background Donor acceptor moieties connected through π-conjugated bridges i.e. D-π-A, in order to facilitate the electron/charge transfer phenomenon, have wide range of applications. Many classes of organic compounds, such as cyanine, coumarin carbazole, indoline, perylene, phenothiazine, triphenylamine, tetrahydroquinoline and pyrrole can act as charge transfer materials. Phenothiazines have been extensively studied as electron donor candidates due to their potential applications as electrochemical, photovoltaic, photo-physical and DSSC materials. Results Two phenothiazine derivatives, 2-((10-hexyl-10H-phenothiazin-3-yl)methylene)malononitrile (3a) and 2-((10-octyl-10H-phenothiazin-3-yl)methylene)malononitrile (3b) have been synthesized in good yields and characterized by various spectroscopic techniques like FT-IR, UV–vis, 1H-NMR, 13C-NMR, and finally confirmed by single crystal X-ray diffraction studies. Density functional theory (DFT) calculations have been performed to compare the theoretical results with the experimental and to probe structural properties. In order to investigate the excited state stabilities the absorption studies have been carried out experimentally as well as theoretically. Conclusions Compound 3a crystallises as monoclinic, P2 (1)/a and 3b as P-1. The X-ray crystal structures reveal that asymmetric unit contains one independent molecule in 3a, whereas 3b exhibits a very interesting behavior in having a higher Z value of 8 and four independent molecules in its asymmetric unit. The molecular electrostatic potential (MEP) mapped over the entire stabilized geometries of the molecules indicates the potential sites for chemical reactivities. Furthermore, high first hyperpolarizability values entitle these compounds as potential candidates in photonic applications.Graphical abstract Phenothiazines; a comparison of experimental and theoretical analysis Electronic supplementary material The online version of this article (doi:10.1186/s13065-016-0158-z) contains supplementary material, which is available to authorized users.


Background
In few years, a great interest has developed in molecules having electron donor-acceptor (D-A) properties and their modern applications as dye sensitized solar cells (DSSC) [1], photosensitizers [2] and redox sensitizers [3]. The metal based donor-acceptor (D-A) complexes are well known where a metal atom behaves as an electron acceptor and ligands as electron donor species [4][5][6]. Ruthenium metal is a key contributor in the synthesis of such complexes. To avoid the cost of metal and its environmental hazards there is a space for the synthesis of new organic donor-acceptor molecules. A salient feature of such organic based (D-A) molecules is that donor acceptor moieties are connected through π-conjugated bridges i.e. D-π-A, in order to facilitate the electron/ charge transfer phenomenon [7]. The classes of organic compounds that have been evaluated as (D-A) candidates include cyanine [8], coumarin [9], carbazole [10], indoline [11], perylene [12], phenothiazine [13], triphenylamine [14], tetrahydroquinoline [15] and pyrrole [16].
Molecules containing phenothiazine as electron donor part have been extensively studied due to their electrochemical [17], photovoltaic [18], photo-physical [19] and DSSC applications [1]. The synthesis of phenothiazine derivatives and their DSSC applications were claimed by many investigators, and the best results were produced in the solar cells where phenothiazine was used as electron donor and boradiazaindacene as electron acceptor candidates [19]. In addition to their physical applications, phenothiazine derivatives have been recognized as potent anti-psychotic [20], anti-infective [21], antioxidant, anti-cancer [22] and anti-Parkinson agents [23]. These were also qualified as valuable MALT1 protease [24], cholinesterase [25], and butyryl-cholinesterase enzyme inhibitors [26].
In addition to our recent work [27][28][29][30][31][32], here we report the synthesis and structural properties of two new phenothiazine derivatives (Fig. 1). Both compounds have been synthesized in high yields and characterized by spectroscopic as well single crystal diffraction studies. The DFT investigations have been performed to validate the spectroscopic results, and to investigate other structural properties like frontier molecular orbital (FMO) analysis, molecular electrostatic potential (MEP), natural bond orbital (NBO) analysis (intra and inter molecular bonding and interaction among bonds), and first hyperpolarizability analysis (nonlinear optical response).

Results and discussion
The synthesis of two phenothiazine derivatives 3a and 3b has been accomplished in three steps beginning from 10-phenothiazine resulting in good yields (details are given in the experimental section). These compounds have been characterized by 1 H-NMR, 13 C-NMR, FT-IR and UV-vis. spectroscopic techniques, and finally their structures have been confirmed by X-ray diffraction analysis. Computational studies have been carried out to compare the theoretically calculated spectroscopic properties with the experimental results, and to investigate some structural properties as well.

X-ray diffraction analysis
Both compounds 3a and 3b have been recrystallized in methanol under slow evaporation method in order to grow suitable crystals to ensure the final structures, and to study their three dimensional interactions. The compound 3a, bearing a hexyl group at nitrogen, is crystallized in a monoclinic system having space group P2 1/a and 3b containing an octyl substituent at nitrogen has been crystallized in a triclinic system having space group P-1. Complete crystal data parameters for both compounds have been provided in Table 1. The ORTEP views of both 3a and 3b are shown in Fig. 2.
While analyzing the crystal structure it is observed that compound 3a exists as single independent molecule in an asymmetric unit. On the other hand, an interesting behavior has been observed for 3b which shows a high Z value of 8   The thiazine rings are not planar having the root mean square (rms) deviation values of 0.1721 (1) Å, 0.1841 (2) Å, 0.2184 (3) Å, 0.1392 (2) Å and 0.1593 (2) Å for compounds 3a and 3b (molecule A, molecule B, molecule C, molecule D) respectively. In compound 3a, the two aromatic rings are oriented at a dihedral angle of 24.80(1)°, while the thiazine ring is oriented at dihedral angles of 13.33 (1)° and 12.56 (1)° with reference to ring 1 (C1-C6) and ring 2 (C7-C12), respectively.

Geometry optimization
In the past decade, methods based on DFT have got the attention of researchers because of their accuracy and wide applications. The DFT investigations of both compounds 3a and 3b have been performed not only to validate X-ray results, but also to compare and investigate other spectroscopic and structural properties. The structures of both 3a and 3b have been optimized by using B3LYP/6-31G (d, p) level of theory, and the the optimized geometries are shown in Fig. 3. A comparison of bond angles and bond lengths for both compounds are listed in Additional file 1: Tables S2,  S3. Although the packing diagram of 3b shows four molecules in asymmetric unit, yet only molecule A has been considered for comparison. The experimental and simulated bond lengths/bond angles of all atoms for compounds 3a and 3b (A) are correlated nicely. A  deviation of 0.001-0.036 Å in bond lengths has been appeared for both compounds. Maximum deviations of 5.4° and 4.2° in dihedral angles from C14-C13-C5 bonds in 3a and from C23-C22-C21 bonds in 3b have been observed.

Vibrational analysis
The experimental vibrational spectra of phenothiazine derivatives 3a and 3b have been recorded as neat, and both the experimental as well as simulated spectra are shown in Fig. 4. The vibrational frequencies of both were computed at the same level as was used for energy minima structures and assignments were accomplished by using Gauss-View 05 program. A comparison of experimental and calculated vibrational frequencies is given in Table 2.
The simulated vibrations above 1700 cm −1 have been scaled by using a scaling factor of 0.958 and for less than 1700 cm −1 scaling factor is 0.9627 [37]. In the table only those simulated vibrations are given whose intensities are more than ten. For both compounds, the vibrations arise mainly from aromatic C-H, double bond C=C, C-N, C-S, nitrile, CH 2 , and CH 3 functional groups. From Table 2, it is clear that there exists an excellent agreement between the experimental and theoretical vibrations.

Aromatic (CH), (C=C) and aliphatic (C=C) vibrations
The aromatic (CH) vibrations generally appear in the region 2800-3100 cm −1 [38]. The bands appeared in this region are normally of very low intensity, and not much affected by substituents. In the simulated spectra, the aromatic CH stretching vibrations of both compounds 3a and 3b have been predicted at 3086, 3077 cm −1 and 3085, 3077 cm −1 respectively. The calculated aromatic CH stretching vibrations coincide well with the experimental value appearing at 2916 cm −1 for both compounds. The symmetric and asymmetric stretching vibrational regions of aromatic ring (C=C) usually lie in between 1600-1200 cm −1 [39]. The experimental scans of 3a and 3b show aromatic C=C stretching vibrations at 1574, 1402 cm −1 and 1570, 1405 cm −1 respectively. The simulated aromatic stretching C=C peaks are found in strong correlation and appear at 1603, 1568, 1526, 1395 cm −1 for compound 3a, and 1594, 1526, 1395 cm −1 for compound 3b. An aliphatic C=C group in conjugation with aromatic ring is also present in both compounds and appears at 1559 cm −1 experimentally whereas this stretching vibration appears at at 1553 cm −1 for both 3a and 3b.
Aromatic in-plane and out of plane CH bending vibrational regions are usually weak and are observed in the range 1000-1300 cm −1 and 650-900 cm −1 respectively [40]. In the simulated spectra, in plane CH (aromatic) bending vibrations are observed in the range of 1428-1286 cm −1 for compound 3a, and in the region of 1352-1139 cm −1 for compound 3b.
These simulated values appear in nice agreement with the experimental values having appeared at 2848 cm −1 for compound 3a, and 2847 cm −1 for compound 3b. Along with the stretching vibrations, several scissoring, in-plane and out of plane bending, methylene (CH 2 ) and methyl vibrations are observed in the simulated and experimental spectra and a nice agreement is found between them. Both compounds 3a and 3b show the CH 2 scissoring vibrations in the range 1456-1448 cm −1 and 1453-1448 cm −1 respectively and these are correlated well with the experimental 1458 and 1462 cm −1 values respectively. The in-plane bending CH 2 vibrations are observed in the range 1337-1275 cm −1 and 1337-1287 cm −1 for 3a and 3b respectively. These bending vibrations are in agreement with the experimental counterparts having appeared at 1317 cm −1 , 1218 and 1323, 1228 cm −1 for 3a and 3b respectively.

Nitrile and C-N Group vibrations
The nitrile symmetric stretching vibrations of very high intensity appear at 2245 cm −1 in the simulated spectra for 3a and 3b. The nitrile asymmetric stretching vibrations of low intensity also appear at 2230 and 2231 cm −1 for both compounds. In the experimental scans, the nitrile vibrations appear at 2214 and 2215 cm −1 for 3a and 3b respectively, and are found in excellent correlation with the simulated values. The simulated C-N-C stretching frequency appear at 1483 cm −1 for both 3a and 3b and is in full agreement with its experimental counterpart observed at 1472 and 1474 cm −1 respectively.
The assignments of N-Ph stretching modes are difficult, as there are problems to discriminate them from other aromatic ring vibrations. For substituted aromatic rings, Silverstein et al. [41] defined the N-Ph stretching bands in the range 1200-1400 cm −1 . In the present study of compound 3a, the observed N-Ph symmetric stretching bands appear at 1338 and 1279 cm −1 in the simulated spectrum and are in very good agreement with the experimental 1363 cm −1 value. Similarly, the calculated N-Ph stretching frequencies of 3b appearing at 1337 and 1279 cm −1 also show good agreement with the experimental band at 1363 cm −1 .

Nuclear magnetic resonance (NMR) studies
For the last two to three decades, nuclear magnetic resonance spectroscopy has been unavoidable tool for structural investigations of organic and biological molecules. The 1 H and 13 C chemical shifts contain very important information about the structural environment of unknown compounds. Nowadays, a powerful method to predict and compare the structure of molecules is to combine the theoretical and experimental NMR methods. The DFT simulations using Gaussian software are playing very active role in this regard. A full and true geometry optimization of both compounds 3a and 3b has been performed by using B3LYP/6-311 + G (2d, p) basis set. An accurate optimization of molecular geometries is vital for reliable calculations of magnetic properties and their comparison with experimental results. The chemical shift calculations of both compounds have been performed by using the fully optimized geometries, adopting the GIAO method at the same level of theory and referred by using the internal reference standard i.e. trimethylsilane. Both the experimental as well as simulated NMR spectra have been recorded in CDCl 3 (for experimental 1 H and 13 C NMR see Additional file 1: Figs. S3-S6). The detailed simulated and experimental 1 H-NMR values are given in Table 3.
Both phenothiazine derivatives (3a and 3b) mainly have aromatic and aliphatic protons. In the experimental 1 H-NMR spectra, aromatic and double bonded protons appear in the range 7.74-6.83 ppm (compound 3a) and Scaling factor used 0.958 for vibrations between 3200 and 1700 cm −1 and 0.9627 used below 1700 cm −1 . Only those simulated values are given, those have shown intensity above 10 υ s symmetric streching, υ as asymmetric streching, β ın plane bending, γ out of plane bending, τ twisting, ρ scissoring, ω wagging

Frontier molecular orbital analysis and UV-vis absorption studies
Frontier molecular orbital analysis has proved very helpful in understanding the electronic transitions within molecules and analyzing the electronic properties, UV-vis absorptions and chemical reactivity as well [42]. The FMO analysis also plays an important role in determining electronic properties such as ionization potential (I. P.) and electron affinity (E. A.). The HOMO (highest occupied molecular orbital) represents the ability to donate electrons and its energy corresponds to ionization potential (I. P.), whereas the LUMO (lowest unoccupied molecular orbital) acts as electron acceptor and its energy corresponds to electron affinity (E. A.) [43]. Frontier molecular orbital (FMO) analysis is carried out at the same level of theory as used for the geometry optimization, applying pop = full as an additional keyword. The HOMO and LUMO surfaces along with the corresponding energies and energy gaps are shown in Additional file 1: Fig. S6. Compound 3a contains 93 filled orbitals, whereas 3b contains 103 filled orbitals. The HOMO-LUMO energy difference in both 3a and 3b has been found to be 2.96 eV. The kinetic stabilities of compounds can be assigned on the basis of HOMO-LUMO energy gap [44]. A low HOMO-LUMO energy gap means less kinetic stability and high chemical reactivity. It is clear that the HOMO-LUMO energy gaps in compounds 3a and 3b are very less, indicating that electrons can easily be shifted from HOMO to LUMO after absorbing energy. The experimental UV-vis absorption spectra of both compounds 3a and 3b in various solvents like dichloromethane, chloroform, methanol and dimethyl sulphoxide (DMSO) have been recorded within 250-700 nm range, and the combined spectra are shown in (Fig. 5). The theoretical absorption studies are also carried out by using TD-DFT method at B3LYP/6-31G (d, p) level of theory in gas phase, and polarizable continuum model (PCM) is applied to account for solvent effect (For simulated UV-vis spectra see Additional file 1: Fig. S7). A comparison of characteristic experimental and simulated UV-vis. absorption wavelengths (λ max ) of the both compounds in gas phase and different solvents (DCM, chloroform, methanol and DMSO) has been given in Table 4. As both the compounds have same chromophores; thus there is no significant difference in their absorption maxima.
Different solvents covering a wide range of polarity and dielectric constant have been selected in order to explore the solvent effect on the absorption maxima, but no significant difference has been observed. The experimental UV-vis. spectra of both compounds show mainly two absorption bands. In dichloromethane, λ max1 and λ max2 values for compound 3a appear at 320 and 474 nm corresponding to the π-π* and n-π* transitions respectively [45], and for 3b the values appear at 321 nm and 474 nm. In chloroform the absorption maxima of 3a are found at 321 nm (λ max1 ), 478 nm (λ max2 ) and for 3b they have been appeared at 321 nm (λ max1 ), 478 (λ max2 ). Similarly, the absorption maxima values appear at 317 nm (λ max1 ), 478 nm for compound 3a, and 317 nm (λ max1 ), 463 nm (λ max2 ), for compound 3b in methanol (polar protic) and DMSO (polar aprotic) respectively. The gas phase simulated spectrum of compound 3a show absorption maxima  Table 4.

Molecular electrostatic potential (MEP)
Molecular electrostatic potential (MEP) is associated with the electronic cloud. The electrophilic/nucleophilic reacting sites as well as hydrogen bonding interactions can be described in any compound on the basis of MEP [46,47]. Recognition process of one molecule by another, as in drug-receptor and enzyme substrate interactions, is related to electrostatic potential V(r), because the two species show interaction to each other through their potentials. The MEP analysis can be performed by using the following mathematical relation, described previously [48].
Here summation (Σ) runs over all nuclei A in a molecule, polarization and reorganization effects are ignored. Z A is charge of nucleus A, located at R A and ρ (r′) is the electron density function of a molecule. Usually, the preferred nucleophilic site is represented by red color and the preferred electrophilic site is represented by blue color. The electrostatic potential values at the surface are represented by different colors. The potential decreases in the order: red < orange < yellow < green < blue. The color code of the map is in the range between 0.0550 a.u. (deepest red) and 0.0550 a.u. (deepest blue), where blue corresponds to the strongest attraction and red corresponds to the strongest repulsion. Regions of negative V (r) are associated with lone pairs of electronegative atoms. According to the MEP analysis of compounds 3a and 3b, there are two negative regions at each molecule (red  coded region) shown in Fig. 6. These red coded regions are nitrile functional groups of the both compounds. As these two compounds differ only at the alkyl chain lengths located at the nitrogen in a heterocyclic ring, therefore the reactive sites are same. Apart from the nitrile groups the rest is lying between yellow and green regions. This shows that no strong electrophilic sites exist in both the compounds.

Natural bond orbital (NBO) analysis
Natural bond orbital analysis is an efficient method for studying intra-and intermolecular bonding and interaction among bonds, and provides a convenient basis to probe charge transfer or conjugative interaction [49]. The NBO approach describes the bonding anti-bonding interaction quantitatively and is expressed by means of second-order perturbation interaction energy E (2) [50][51][52][53]. This energy estimates the off-diagonal NBO Fock matrix element. The stabilization energy E (2) associated with i (donor) to j (acceptor) delocalization is approximated from the second-order perturbation approach as given below: where q i is the donor orbital occupancy, ε i and ε j are the diagonal elements (orbital energies) and F (i, j) is the offdiagonal Fock matrix element. The larger the E (2) value is, the greater is the interaction between electron donors and electron acceptors and the extent of conjugation of whole system. The various second-order interactions between the occupied Lewis type (bond or line pair) NBO orbitals and unoccupied (anti-bonding and Rydberg) non-Lewis NBO orbitals are investigated by applying DFT at the B3LYP/6-31G (d, p) level. As a result of our study, the compounds 3a and 3b are types of Lewis structures with 97.93 and 98.03 % character, valance-non Lewis character of 1.90 and 1.79 % respectively. Both the compounds share the same Rydberg non-Lewis character of 0.16 %.
The intramolecular hyperconjugative interactions result in the transfer of charge from donor (π) to acceptor (π*) orbitals. This charge transfer increases the electron density (occupancy) in antibonding orbitals and weakens the respective bonds [54]. From the significant entries in Table 5, it is clear that the occupancy of π bonds (C-C) for benzene rings of the title compounds (3a and 3b) lie in the range of ~1.59-1.71. On the other hand, the occupancy of π* bonds (C-C) for benzene rings range from ~0.33-0.42. This delocalization leads to the stabilized energy in the range of ~17. 15-25.19 kcal/mol.
The pi-bond of ethylenic moiety (C13-C14) also shows an average of ~20 kcal/mol stabilization energy when it is delocalized to either acetonitrile group. The strongest stabilization energy to the system by 31.28 kcal/mol is due to the lone pair donation of nitrogen atom N (1) to the antibonding π* (C2-C3) orbital. On the other hand, the same lone pair gives a stabilization energy of 24.09 kcal/ mol when it is conjugated with the antibonding π* (C11-C12) orbital of the aromatic ring. This clearly shows that the delocalization of lone pair of nitrogen N (1) is more towards that aromatic ring which has extended conjugation due to presence of electron withdrawing acetonitrile groups. The lone pair donation from sulfur atom (S1) to the antibonding π* (C1-C6) and (C7-C8) orbitals of both phenyl rings results in the stabilization energies of 12.09 and 11.23 kcal/mol respectively. The occupancy of lone pair electrons in sulfur atom (S1) is 1.84 as compared to 1.69 of lone pair on nitrogen atom (N1). As a consequence, the stabilization energies arising from the lone pair donation of sulfur atom to the antibonding π* (C-C) bonds of phenyl rings are comparatively smaller than those arising from lone pair donation of N1 atom. A plausible reason could be due to the deviation of sulfur atom from planarity because of its larger size. All σ to σ* transitions involving C-C bonds correspond to the weak stabilization energies in the range of ~2.53-4.58 kcal/mol.

Hyperpolarizability and non-linear optical properties
Recently, compounds having non-linear optical (NLO) properties have got appreciable attention of researchers because of their wide applications in optoelectronic devices of telecommunications, information storage, optical switching and signal processing [55]. Molecules containing donor acceptor groups along with pi-electron conjugated system are considered as strong candidates for possessing NLO properties [56].
In each 3a and 3b, the phenothiazine moiety is connected to a nitrile group through a conjugated double bond, and these molecules are anticipated to show nonlinear optical (NLO) properties. For the estimation of NLO properties, the first hyperpolarizability (βo) analysis for compounds 3a and 3b has been performed by employing same level of theory as for geometry optimization i.e. 6-31G (d, p) along with POLAR as an additional keyword. The first hyperpolarizability, a third rank tensor, is always described by a 3 × 3 × 3 matrix. The total 27 components of the 3D matrix can be reduced to 10 components as a result of Kleinman symmetry [57]. From the Gaussian output file ten components of 3D matrix have been identified as β xxx , β xxy , β xyy , β yyy , β xxz , β xyz , β yyz , β xzz , β yzz and β zzz respectively, and the values are given in Table 6.
Among all types of hyperpolarizabilities reported in literature, the more attractive is β tot . (First hyperpolarizability) [49] and it can be measured by using the following mathematical relation; respectively. These values are in excellent agreement with the reported values in literature [58,59], and this agreement proves that both compounds are strong candidates for NLO applications.

Method
All analytical grade chemicals and solvents were purchased from BDH, and used without further purification. Stuart Scientific (SMP3, version 5.0, UK) melting point apparatus was used to record the melting point, and the reported m. p. were uncorrected. 1 H-NMR spectra were recorded on a Bruker-AVANCE-III 600 MHz at 300 K, and chemical shifts were reported in ppm with reference to the residual solvent signal. FT-IR spectra were Table 5 Significant donor-acceptor interactions of 3a/3b and their second order perturbation energies calculated at B3LYP level using 6-31G (d, p) basis set a E (2) means energy of hyperconjucative interactions (stabilization energy) b Energy difference between donor (i) and acceptor (j) NBO orbitals

Computational details
Theoretical studies were performed by using Gaussian 09 software at density functional theory (DFT) level, as instituted in program [66]. The visualization of the results/optimized geometries was achieved by using Gauss view 05 [67]. The energy minima optimization of both compounds was carried out at B3LYP/6-31G (d, p) and B3LYP/6-311 + G (2d, p) levels of theory (the later was used further for nuclear magnetic studies). Frequency simulations were performed at the same level, to confirm the optimized geometries as a true minimum (no imaginary frequency). In addition, frequency simulations at B3LYP/6-311G (d, p) level were used for vibrational analysis. Nuclear magnetic resonance studies were performed at B3LYP/6-311 + G (2d, p) level, by adopting GIAO method in chloroform solvent and applying polarizable continuum model (PCM) for the solvent consideration. Chemical shift values were referred by using the internal reference standard i.e., tetramethylsilane. UV-vis absorption studies were simulated by using TD-DFT method and at B3LYP/6-31G (d, p) level of theory. MEP, NBO, FMO and first hyperpolarizability analyses were simulated at B3LYP/6-31G (d, p) level of DFT.

Experimental
The synthesis of both phenothiazine derivatives was carried out in three steps starting from simple phenothiazine. First step was alkylation of nitrogen, followed by subsequent aldehyde formation and then conversion to final product (Fig. 1).

General procedure for the synthesis of N-alkylated phenothiazine (1a, 1b)
In a round bottom flask a mixture of potassium hydroxide (2.003 g, 0.0357 mol), 10-phenothiazine (2.91 g, 0.0119 mol), 1-bromohexane (for 1a) or 1-bromooctane (0.0179 mol) (for 1b) and potassium iodide (in catalytic amount) in 50 ml dimethyl sulfoxide (DMSO) were taken. The reaction mixture was stirred for 5 h at room temperature and water (200 ml) was added. The crude product was extracted with CHCl 3 (3 × 50 ml) and the organic layer was washed with saturated ammonium chloride solution and then with water. The organic layer was dried over anhydrous sodium sulfate and filtered, after removing the solvent under reduced pressure, crude product was purified by flash column chromatography (eluent: n-hexane) to obtain colorless oil 1a in 88.68 % yield, and 1b in 86.15 % yield.

Conclusions
In this study, two novel phenothiazine derivatives 2-((10-hexyl-10H-phenothiazin-3-yl)methylene)malononitrile (3a) and 2-((10-octyl-10H-phenothiazin-3-yl) methylene)malononitrile (3b) have been synthesized and characterized by using FT-IR, UV-vis, 1 H, 13 C-NMR spectroscopic techniques and finally their structures are confirmed by single crystal X-ray diffraction studies. The DFT studies have shown a strong agreement between the simulated and experimental results. The optimized geometries of the both compounds at 6-31G (d, p) level have been used further for investigating structural properties. Frontier molecular orbital analysis shows that both the molecules have very low HOMO-LUMO energy gap, and therefore are kinetically less stable. The molecular electrostatic potential investigations reveal that electronegative region in both the compounds is spread over the nitrile groups. The high first hyperpolarizability values signify that these compounds can have very good nonlinear optical responses. The phenothiazine derivatives have very wide applications not only in dye sensitized solar cells but also in clinical field, and hopefully the results of this study will increase the interest of researchers working in this field.