Skip to main content

Synthesis, crystal structure, Hirshfeld surface investigation and comparative DFT studies of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate


The ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1), a thiazole ester, was synthesized by refluxing 1-(2-nitrobenzylidene)thiosemicarbazide and ethyl bromopyruvate. The compound is characterized by spectrometric, spectroscopic and single crystal (SC-XRD) techniques. Non-covalent interactions that are responsible for crystal packing are explored by Hirshfeld surface analysis. All theoretical calculations were performed by DFT quantum chemical methods using 6-311G(d,p) and cc-pVTZ basis sets and compared. Theoretical harmonic frequencies of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1) were optimized. Confirmation of hydrogen bonding sites was analyzed by molecular electrostatic potential (MEP) and Mulliken population analysis. The vibrational frequencies of characteristic functional groups and chemical shifts were found in good agreement with experimental assignments. Frontier molecular orbital (FMO) revealed relatively small HOMO–LUMO (highest occupied molecular orbital-lowest unoccupied molecular orbital) gape, which speaks off the nearly planar geometry and extended conjugation, as compared to the substituents with no conjugation possible. It has also been observed that –NO2 substituent plays a vital role for this relatively small HOMO–LUMO gape and overall electronic properties when compared with similar thiazole carboxylates (2–6, Table 6). Ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1) was also evaluated for its anti-oxidant and anti-microbial activities.

Graphical Abstract


1,3-Thiazole is an important sulphur and nitrogen containing azole [1]. Binding capability of sulphur and nitrogen in 1,3-thiazole, to induce sigma and pi interactions with receptor moieties, makes it an interesting class of heterocycles in drug designing [2]. It is very significant motif in pharmaceutical chemistry [3,4,5]. Thiazole derivatives are known for their diverse biological significance like anticancer [6], antiallergic [7], antiviral [8], antibacterial [9], antifungal [10], antitubercular [11, 12] and antiinflammatory [13,14,15] activities. Besides these, antiparasitic agents viz, nitazoxanide, tenonitrozole and aminitrozole have 1,3-thiazole ring in their structure [16].

Computational chemists use both experimental data and theoretical models in order to determine the structural features like bond lengths, bond angles and torsion angles besides accounting for spectroscopic properties like harmonic frequencies, vibrations and chemical shifts at molecular level. An important tool in the hand of theoretical chemists to compare and analyze the properties like geometry, vibrational frequency and dipole moment [17,18,19,20,21,22,23,24] of the molecules, is Density Functional Theory (DFT). The thiazole ester, ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1) was synthesized for the purpose of biological studies. However, it was found interesting to investigate the structural features by comparing theoretical and experimental data and also to see the effect of the substituent i.e. nitro group (–NO2), on the conjugation leading to difference in the HOMO–LUMO gape with extended conjugation. The analysis of theoretical results revealed that –NO2 substituent has played a vital role in decreasing the HOMO–LUMO gape, furthermore; it has also exerted its effect on overall electronic properties when compared with compounds 2–6.

Materials and methods

Ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1) was synthesized by the reaction of 1-(2-nitrobenzylidene)thiosemicarbazide and ethyl bromopyruvate (Scheme 1). The pure reagent grade chemicals were purchased and used as such. Thiosemicarbazide, 4-nitrobenzaldehyde and ethyl bromopyruvate were purchased form Merck Germany. Pre-coated Silica 60 HF254 Aluminum sheets (Merck, Germany) were used to monitor the reaction with thin layer chromatography (TLC). The melting point determination, functional group identification, NMR chemicals shifts and single crystal analysis were carried out as reported elsewhere [25]. Dimethyl sulfoxide-deuterated (DMSO-d6) was used to record the proton and carbon NMR (nuclear magnetic resonance), using 300 and 75 MHz frequency, respectively. High resolution mass spectrometry (HRMS) was carried out on Bruker Micro TOF-ESI (time of flight-electrospray ionization) spectrometer, positive targeted mode.

Scheme 1
scheme 1

Synthesis of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate

Synthesis of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1)

The desired compound (1) was obtained by refluxing 1-(2-nitrobenzylidene)thiosemicarbazide (0.5 g, 0.00237 mol) and ethyl bromopyruvate (0.3 mL, 0.00238 mol) for 4 h in absolute ethanol (20 mL). The reaction progress assessed by TLC. When the reactants consumed, as shown by TLC, quenching of reaction by adding ice cold water resulted into precipitates, which on filtration and washing with excess of water yielded pure compound [25].

Colour: Yellow solid; Yield: 82%; Melting point: Above 300 °C; Rf: 0.52 (acetone/n-hexane, 1:2); FTIR (ATR, cm−1): ν̅ 1076, 1338, 1436, 1517, 1571, 1684, 2981; 1H–NMR: δ (ppm) 1.28 (3H, t, –CH3, J = 7.2 Hz), 4.24 (2H, q, –OCH2–CH3, J = 7.0 Hz), 7.60 (1H, dt, Ar–H, J = 1.5 Hz, J = 8.4 Hz), 7.78 (1H, t, Ar–H, J = 7.7 Hz), 8.02 (2H, m, Ar–H), 7.82 (1H, s, 1,3-thiazole ring C-5), 8.37 (1H, s, –CH=N–), 12.63 (1H, s, –N–NH–C–); 13C–NMR: δ (ppm) 14.6 (–CH3), 60.9 (–OCH2–), 120.2 (C-5, 1,3-thiazole ring), 125.2, 128.1, 128.7, 130.4, 137.4, 147.9 (Ar–C), 137.6 (–CH=N–), 143.3 (C-4, 1,3-thiazole ring), 161.7 (C=O), 168.2 (C-2, 1,3-thiazole ring); HRMS: 321.9680 [M + H]+, 343.0700 [M + Na]+.

Single Crystal X-Rays Diffraction (SC-XRD) studies

The crystallization of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1) from tetrahydrofuran (THF) yielded the crystals suitable for diffraction collection using single crystal X-ray diffractions. Table 1 presents the experimental conditions, different parameters and structure refinements. The data reduction and cell refinement carried out using SAINT and data collection by APEX 2 [26]. SHELXS97 [27] and SHELXL97 [28] were used for solving and refining the structure, respectively. PLATON [29], Mercury 3.6 software and ORTEP-3 for Windows [30] were used to draw molecular graphics. The CIF is deposited at Cambridge Crystallographic Data Centre (CCDC = 1,951,603).

Table 1 SC-XRD experimental details of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1)

Results and discussion

The synthesis of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1) was achieved in two steps form thiosemicarbazide via thiosemicarbazone and cyclization with ethyl bromopyruvate. The synthesis was confirmed by spectroanalytical techniques. The spectroscopic data verified the proposed synthesis and the structure was later confirmed by the SC-XRD.

Single crystal X-ray

The SC-XRD inspection showed that ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1) (Fig. 1) exists in E-configuration, the benzylidenehydrazinyl moiety A (C1–C7/N1/N2) and the thiazole ring B (C8/N4/C9/C10/S1) are found to be planar with respective root mean square (r. m. s.) deviations of 0.0287 and 0.0032 Å. The dihedral angle between A/B is 8.41 (9)°. The planarity of the thiazole ring in two most closely related reported crystal structures, one having methyl substituted phenyl ring and second having bromo substituted phenyl ring was obvious by r.m.s. deviation of 0.004 and 0.0072 Å [25]. The nitro group C (N1/O1/O2) and carboxylate group D (C11/O3/O4) are oriented at a dihedral angle of 26.84 (10)° and 4.6 (5)° with the parent groups A and B, respectively. The C-atoms (C12/C13) are at a distance of -0.003 (11) and 1.329 (13) Å, respectively, from the carboxylate group as compared to their distance of 0.034 (16) and 1.440 (21) Å from the carboxylate group in closely related crystal structure having bromo substituted phenyl ring instead of nitro substituted phenyl ring [25]. Intramolecular H-bonding of type C–H…N and C–H…O are responsible for stabilization of molecular configuration. CH of imine functional group is found to engage with O-atom of nitro group (O2) to form S(6) loop through intra C–H…O bonding. Similarly, N-atom of imine functional group is found to engage with CH of phenyl ring to form S(5) loop through C–H…N bonding. Another S(5) is formed by C–H…O bonding, where O-atom is from carboxylate group and CH is from alkyl group (C12/C13). These intramolecular H-bonded loops are displayed in Fig. 1 and specified in Table 2. The molecules are primarily connected to each other in the form of dimers through strong N–H…O and comparatively weak C–H…N bonding to form \({\mathrm{R}}_{2}^{2}(9)\) loop as displayed in Fig. 2 and specified in Table 2 [32]. Infinite C7 chain is formed by N3-H3A…O3 bonding whereas C4 chain is formed by C–H…N bonding. These chains extend along b crystallographic axis in a zigzag fashion. The molecules are also interlinked by C–H…O bonding, where para CH of phenyl ring (C1–C6) acts as donor and O-atom is from nitro group. Due to this C–H…O bonding, another zigzag C6 chain is formed along c crystallographic axis. Besides the role of H-bonding in packing of molecules, weak non-covalent interaction named as off-set π–π stacking interaction also play significant role in the crystal packing. Five membered ring of one molecule located at asymmetric location is engaged in off-set π–π stacking interaction with the five membered ring of neighboring symmetry related molecule (-x,1-y,-z). Centroid to centroid distance and ring off-set for this particular interaction are found to be 3.92 Å and 1.798 Å, respectively as displayed in Fig. 3.

Fig. 1
figure 1

The ORTEP diagram drawn at probability level of 50% for (1). H-atoms are displayed by tiny circles of arbitrary radii

Table 2 Hydrogen-bond geometry (Å, º) for (1)
Fig. 2
figure 2

Packing diagram showing 2D sheet formed by N–H…O, C–H…N and C–H…O bonding for (1)

Fig. 3
figure 3

Graphical view of π–π stacking interaction for (1). Distance shown is measured in Å. H-atoms are omitted for clarity

Hirshfeld surface analysis (HSA)

For further exploration of non-covalent interactions that are accountable for the packing of molecules, Hirshfeld surface inspection is conducted on Crystal Explorer 17.5 [33]. Hirshfeld surface (HS) plotted over dnorm can be utilized to represent intermolecular interactions by colour coding [34,35,36]. The atomic contacts having distance shorter, equal and longer than summation of Van der Waals radii are represented by red, white and blue regions on HS, respectively. Bright red spots on the HS near NH group, carbonyl O-atom, para CH of nitro phenyl ring, one of the O-atom of nitro group indicate that these atoms are involved in H-bonding as displayed in Fig. 5. N-atom and CH of five membered ring are also engaged in H-bonding but no red spot is obvious in Fig. 4 due to graphical limitations. No red spot is found on the HS near imine N-atom indicating that this atom is not engaged in any kind of H-bonding. HS plotted over shape index is utilized to check that whether π–π stacking interaction is present in crystal packing or not. Consecutive triangular regions of red and blue on the HS around the five membered ring indicated that π–π stacking interaction is present in crystal packing as displayed in Fig. 5. As these triangular regions on the HS are not exactly in the center of five membered ring suggesting off-set π–π stacking interaction. In order to find out the percentage contribution of each interatomic contact involved in crystal packing, 2D fingerprint plots are utilized. In these plots, di and de stands for distance from the HS to nearest atom inside and outside the HS, respectively. 2D fingerprint plot for overall interactions in the crystal packing is displayed in Fig. 6a. The sky blue region in the center of this plot indicates the existence of π–π stacking interaction in crystal packing. The contribution of each individual interatomic contact is calculated including reciprocal contact. It is found that O…H, H…H, C…H, N…H and S…H interatomic contacts are the major contributors in crystal packing with percentage contribution of 29.9%, 28%, 12.4%, 7.4% and 5.5%, respectively as displayed in Fig. 6b-f. C…C, N…C, S…C, O…C, O…N, S…N, N…N, O…O, S…S, and S…O interatomic contacts are found to be minor contributor in crystal packing with percentage contribution of 3.5%, 3.4%, 2.4%, 2.3%, 1.6%, 1.1%, 1%, 0.3%, 0.2% and 0.1%, respectively as displayed in Fig. 7a–j. Furthermore, the interaction of an atom located inside the HS is computed with the atoms of molecules present in the surrounding of HS to have an idea which atom of a molecule is interacting strongly with the atoms of neighboring molecules. H-atoms located inside the HS interact strongly with the atoms of molecules in the neighbouring of HS, percentage contribution of this interaction is found to be 52.4%. The percent (%) contribution of all such interactions are displayed in Fig. 8.

Fig. 4
figure 4

HS plotted over dnorm in the range − 0.5540 to 1.4207 a.u. for (1)

Fig. 5
figure 5

HS plotted over shape index in the range − 1 to 1a.u. for (1)

Fig. 6
figure 6

2D plots for (1), a for overall interaction, bf for major contributor in crystal packing

Fig. 7
figure 7

a–j 2D plots for minor contributor in crystal packing for (1)

Fig. 8
figure 8

Percentage contribution of interaction of an atom located inside HS to atoms of molecules located in the surrounding of HS

Quantum chemical calculations

Guassian09 package [37,38,39] and Gauss-View molecular visualization software [40] was used in order to carry out calculation studies on personal computer. DFT (DFT/B3LYP) method [41, 42] with 6-311G(d,p) and cc-pVTZ basis sets were used for geometry optimization and further calculations (Fig. 9). The optimized structures obtained from these two methods were compared (Fig. 10) which showed almost same results. Optimized geometry was used in order to ascertain different aspects like vibrational spectra, chemical shifts, quantum chemical and thermodynamic properties of the compound, theoretically [43].

Fig. 9
figure 9

Optimized structure of compound (a) 6-311G(d,p) (b) ccPVTZ

Fig. 10
figure 10

Overlay structure of compound optimized by 6-311G (d,p) and cc-pVTZ method

Vibrational spectra

The vibrational spectra of the compound were computed theoretically (DFT/B3LYP method with 6-311G(d,p) and cc-pVTZ basis sets). The results obtained by two methods were subjected to scaling factors of 0.966 and 0.966 [44, 45] for 6-311G(d,p) and cc-pVTZ basis sets, respectively. The vibrational modes of the compound are assigned by animation of Guass View software, for both used methods. Total energy distribution (TED) calculations for compound (1) were executed with the scaled quantum mechanical (SQM) program [46, 47]. Several factors were taken into account for matching of the experimental and theoretical vibrational data of the compound (1), such as: animation option of the Gauss View package program, relative intensities of the vibrational bands and previously reported data in the literature [37]. The compound (1) contains 34 atoms and it shows 96 normal vibrational modes having a member of C point group with identity (E) symmetry operation. The conformer having the lowest optimization energy has been selected for calculating vibrational spectra of the molecule. At the end all the experimental results were compared to the theoretical results. Theoretical findings for infrared wavenumbers together with experimental results and TED values of the compound (1) are given in supporting information (Additional file 1: Table S1).

In the high wavenumber region bands appeared at around 3493 cm−1 are due to NH vibrations while those appeared at 3273–3036 cm−1 are due to symmetric C–H vibrations. The broadened shape of the band at this region could be attributed to the indication of residue of moisture in the sample and intra and intermolecular interactions. The band appeared at 1685 and 1571 cm−1 are due to C = O and –NO2 stretchings respectively [48]. The bands at 1298, 1251, 1244, 1160 and 1103 cm−1 are due to C–H stretching. The bands at 1336 and 742 cm−1 are due to C–C stretching [49, 50]. For the sake of better comparison experimental and calculated results were plotted as shown in Fig. 11.

Fig. 11
figure 11

Comparison of FT-IR spectra with experiment and calculated values (DFT/B3LYP method with the 6-311G(d,p) and cc-pVTZ basis sets)

NMR spectra

The 1H- and 13C–NMR spectra of the compound (1) were obtained in deuterated dimethyl sulfoxide and signals calibrated with respect to TMS and assigned the chemical shifts to respective protons and carbons on the basis of the chemical environment. The proton NMR spectrum exhibited a triplet and quartet at 1.28 and 4.24 ppm corresponding to methyl and methylene of ethyl ester moiety in the molecules, respectively. The ortho substitution of the benzene ring was confirmed by the presence of a triplet and a triplet of doublets at 7.60 and 7.78 ppm, respectively. The triplet of doublets indicates one small coupling constant (J = 1.5 Hz) corresponding to meta splitting and one large coupling constant (J = 8.4 Hz) for the ortho splitting. Two remaining aromatic protons were multiplet at 8.04–7.99 ppm. The singlet at 7.82 ppm, integrating for one proton, was assigned to thiazole ring proton. The benzylidene –CH– proton was observed as singlet at 8.37 ppm. The –NH proton was observed as a downfield singlet at 12.63 ppm.

The carbon NMR spectrum showed thirteen signals corresponding to total number of carbons in the compound. Two signals in aliphatic region at 14.6 and 60.9 ppm were assigned to the ethoxy carbons i.e. –CH3 and –OCH2, respectively. The six aromatic carbons of benzene ring were seen between 125.2 to 147.9 ppm. The benzylidene –CH– and thiazole ring C-4 was observed at 137.6 and 143.3, respectively. The carbonyl of ester and thiazole C-2 were observed at 161.7 and 168.2 ppm, respectively. The proton and carbon spectra confirmed the proposed structure for the synthesized compound.

GIAO approach was used in order to calculate chemical shifts using B3LYP/6-311G(d,p) and B3LYP/cc-pVTZ, which are then converted to the TMS scale as reported earlier [25] using the formula i.e. δ = Σo–Σ. IEF-PCM (Integral-Equation-Formalism Polarizable Continuum Model) model was employed in order to incorporate the solvent effects [51]. Correlations graphs were produced for comparison of theoretically calculated NMR results to experimental values (Figs. 12 and 13). The correlation co-efficients for B3LYP/6-311G(d,p) (Fig. 12) are 0.9964 for –1H and 0.9977 for –13C NMR, while correlation co-efficients for B3LYP/cc-pVTZ (Fig. 7) 0.9976 for –1H and 0.9982 for –13C NMR were observed which reflects better agreement in experimental and theoretical results. The experimental and theoretically calculated NMR data by both methods has been presented in supporting information (Additional file 1: Tables S2, S3).

Fig. 12
figure 12

Correlation graphics of calculated and experimental chemical shifts of compound (a) 1H–NMR (b) 13C–NMR by B3LYP/6-311G(d,p)

Fig. 13
figure 13

Correlation graphics of calculated and experimental chemical shifts of compound (a) 1H–NMR (b) 13C–NMR by B3LYP/cc-pVTZ

HOMO and LUMO analysis

Frontier molecular orbitals (FMO) have unique importance for calculation of electronic states and stability of a molecule because these orbitals are at frontiers of electron occupation i.e. highest occupied molecular orbital (HOMO) and lowest lying unoccupied molecular orbital (LUMO) [44]. Theoretically calculated HOMO and LUMO are shown in Fig. 8. The compound possessed 83 occupied molecular orbitals. The respective energies for HOMO and LUMO are − 6.1093 and − 2.8798 eV for B3LYP/6-311G(d, p) levels and − 6.1148 and − 2.8763 for B3LYP/cc-pVTZ levels, respectively. Ionization potential (I), electron affinity (A) electronegativity (χ), chemical hardness (η) and chemical softness (S) have been obtained from FMO data (Table 3) [52]. Both methods produced almost the same results.

Table 3 The calculated FMO parameters using B3LYP/6-311G(d, p) level

FMO theory assists in the study of photo-electronic properties of the interesting compounds. The energy of HOMO’s and LUMO’s were computed by DFT method using B3LYP level along with 6-311G(d,p) and cc-pVTZ subset levels. The results are presented in Fig. 8 which shows the electronic contribution in HOMO-2, HOMO-1 and HOMO with LUMO, LUMO + 1 and LUMO + 2 of the compound. There is the delocalization of HOMO on all atoms of the molecule except ethyl group while concentration of LUMO at all atoms except ester group (Fig. 14). Theoretical data of HOMO (H) and LUMO (L) energies has been described in Table 4 which shows almost same values of energies for both methods i.e. 6-311G(d,p) and cc-pVTZ [53].

Fig. 14
figure 14

FMOs at B3LYP/cc-pVTZ level

Table 4 FMO energies calculated using B3LYP/6–31 + G(d, p) level and B3LYP/cc-pVTZ level

UV–Vis analysis

TDDFT studies (B3LYP method with 6-311G(d,p) and cc-pVTZ basis sets) were carried out in gas phase to get absorption properties. Calculated parameters i.e. wavelength of maximum absorption (λmax), oscillator strength (f), excitation energy (EDFT) and orbitals mainly involved for excitation are listed in Table 5 below [54].

Table 5 Wavelength, excitation energy and oscillator strength of compound

The value of computed absorption maxima (λDFT) were found to be 436.53, 340.35 and 437.92, 342.88 nm for B3LYP/6-311G(d, p) and B3LYP/ccPVTZ levels respectively with percentage of orbital contribution of [H → L (70%) and H → L + 2 (66%)]. Calculated wavelength (λDFT) 340.35 and 342.88 employing two different basis sets matched nicely with experimental (λExp) value of 350 (DMSO), 338 (DMF) and 360 (THF). It may be concluded that solvent plays its part for the difference between λDFT and λExp value (Additional file 1: Figure S1).

Molecular electrostatic potential (MEP)

The information about electron acceptor and electron-donor regions in a molecule can be easily visualized by MEP. Here the electrostatic potential decreases in the following order red < orange < yellow < green < blue. Molecular electrostatic map of compound is shown in Fig. 15 (B3LYP/cc-pVTZ basis set) where blue and red indicates the areas of strongest attraction and strongest repulsion, respectively. It is concluded that area of strong repulsions i.e. red resides on O1, O2, O3 and N4 atoms while that of strong attractions i.e. blue resides at N2 and N3 atoms, respectively. Therefore, it is assessed that N2 and N3 atom are electron-donor while O1, O2, O3 and N4 atoms are electron-acceptor [55].

Fig. 15
figure 15

MEP map of compound in gas phase usingB3LYP/cc-pVTZ

Mulliken population analysis

The charge distribution on the molecule was calculated theoretically with the help of B3LYP methods using 6-311G(d, p) and cc-pVTZ basis sets (Fig. 16). It is the charge distribution within the molecule, there by, rendering important information about hydrogen bonding within the molecule [54, 55]. It can also be used to characterize the electronic charge distribution in a molecule and the bonding, antibonding, or nonbonding nature of the molecular orbitals for particular pairs of atoms. Mullikan charges of hydrogen atoms are in positive number while negative for carbons except carbon atoms that are attached to oxygen and nitrogen atoms (due to strong electronegativity of O and N atoms). Oxygen and Nitrogen atoms play the vital role in inter- and intra-molecular hydrogen bonding.

Fig. 16
figure 16

The Mullikan charges diagram of compound (a) B3LYP/6-311G(d,p) (b) B3LYP/ccPVTZ

Thermodynamic properties

The vibrational analysis (B3LYP method with 6-311G(d,p) and B3LYP/cc-pVTZ basis sets) was used to calculate the thermodynamical parameters and properties of the compound i.e. entropy (\({S}_{m}^{o}\)), heat capacity (\({C}_{p,m}^{o}\)) and enthalpy (\({H}_{m}^{o}\)). Thermodynamical parameters i.e. entropies, standard heat capacities and enthalpies exhibit a direct relation with temperature. Their correlation equations for both B3LYP/6-311G(d,p) (Eqs. 13) and B3LYP/cc-pVTZ (Eq. 46) levels have been calculated (Fig. 11) and may be used for further studies of the compound and may serve as guideline for similar compounds. From these correlations it is evident that entropy, enthalpy and heat capacity of the compound (1) increases with increase in temperature.

Heat capacity is due to the rotational and vibrational energy that a molecule possesses. At low temperature there is not enough energy in order to excite many rotations and vibrations, so the slope is small. At very high temperature many energy levels are populated, so the rate of increase becomes constant. The same phenomenon has been observed in case of molecule (1) (Fig. 17). As enthalpy is sum of total energies of a molecule, so, it follows the same pattern of slow increase at low temperature but at high temperature it increases steadily due to the dominance of entropy. Entropy of a molecule increases with increase in temperature. The same phenomenon of enthalpy and entropy has also been observed in case of compound (1) also.

Fig. 17
figure 17

Thermodynamic properties of compound (a) B3LYP/6-311G(d,p) (b) B3LYP/cc-pVTZ

$${S}_{m}^{o}= -8{E}^{-5}{T}^{2}+0.3079T+71.92 \left({R}^{2}=0.9911\right)$$
$${H}_{m}^{o}= 7{E}^{-5}{T}^{2}+0.0363T+151.81 \left({R}^{2}=0.9901\right)$$
$${C}_{p,m}^{o}=-1{E}^{-4}{T}^{2}+0.2593T+7.5953 \left({R}^{2}=0.9917\right)$$
$${S}_{m}^{o}= -8{E}^{-5}{T}^{2}+0.3081T+71.98 \left({R}^{2}=0.9916\right)$$
$${H}_{m}^{o}=7{E}^{-5}{T}^{2}+0.0362T+151.92 \left({R}^{2}=0.9924\right)$$
$${C}_{p,m}^{o}=-1{E}^{-4}{T}^{2}+0.2588T+7.5854 \left({R}^{2}=0.9913\right)$$

Comparison of HOMO and LUMO analysis of compound with other derivatives

Careful analysis of molecular orbitals of molecules (16) reveals that except –CO– and OEt groups, HOMO and LUMO are completely delocalized on all atoms. The energy gap (∆E) between HOMO and LUMO is summarized in Table 6 below shows lowest energy gap (3.2385 eV) in case of compound (1) having electron withdrawing –NO2 substituent while highest energy gap (3.982 eV) in case of molecule (4) having methyl substituent. The energy gaps of other molecules are in between (3.819 eV-3.9149 eV) although there is a little difference. This data clearly indicates different substituents (electron donor/electron acceptor) has great effect on electronic properties of the molecule (Table 6). Therefore compound (1) possess better capability for electron transfer (lowest ∆E) as compared to other compounds (26) [25].

Table 6 Comparison of frontier orbital energies, electronegativity, hardness and softness of compounds (1–6)

The calculated data of EHOMO and ELUMO can be used in order to determine other molecular parameters which are also summarized in Table 6. There is an inverse relationship between energy gap, softness and reactivity of species while a direct relation towards hardness and stability of the molecule. Therefore, compound (4) which possess highest energy gap will resist any change indicating low reactivity while compound (1) will have high reactivity. The calculated energy gap of compound (1) is lowest but ionization potential (6.1148 eV) and electron affinity (2.8763 eV) values of compound (1) are greater than others. The hardness of molecule (1) is lowest (1.6192 eV) while softness is greatest (0.3088 eV) thus confirming the inverse relationship of energy to softness and direct relationship to hardness, respectively. The calculated values of softness (0.2511 eV–0.2618 eV) and hardness (1.9095 eV–1.9574 eV) of other molecules are comparable (Table 6) [58, 59].

The chemical potential (μ) is related to the stability of a molecule and there is a direct relationship. The value of μ for compound (1) is least (− 4.4955 eV) while higher (− 3.5402 eV) for compound (3). The order of chemical potential is as under:

$$\begin{aligned} & \left[ {{\mathbf{1}}\left( {\mu = - {4}.{\text{4955 eV}}} \right)} \right] < \, \left[ {{\mathbf{6}}\left( {\mu = - {4}.0{\text{795 eV}}} \right)} \right] \, < \, \left[ {{\mathbf{5}}\left( {\mu = - {4}.0{1}0{\text{9 eV}}} \right)} \right] \, \\ & \quad < \, \left[ {{\mathbf{2}}\left( {\mu = - {3}.{\text{9274 eV}}} \right)} \right] \, < \, \left[ {{\mathbf{4}}\left( {\mu = - {3}.{\text{873 eV}}} \right)} \right] \, < \, \left[ {{\mathbf{3}}\left( {\mu = - {3}.{54}0{\text{2 eV}}} \right)} \right] \\ \end{aligned}$$

The data also indicates that –CH3, –Br and –OH substituent has almost comparable effect whereas –NO2 substituent has an immense effect on overall electronic properties of the compound. Based on calculated data of hardness it can be inferred that the compound (2) (1.9574 eV) is the most stable and less reactive compared to compound (6) (1.6192 eV). Although hardness values of compounds (26) are almost comparable (Table 6).

Antioxidant assays [57]

Percent free radical scavenging assay (%FRSA), total antioxidant capacity (TAC) and total reducing power (TRP) studies of compound (1) exhibit 77.52 ± 0.34 percent scavenging, 144.07 ± 0.16 and 62.13 ± 0.84 µgAAE/mg values, respectively. The %FRSA, TAC and TRP results support the use of compound (1) as an antioxidant agent. However, the antifungal and antibacterial studies did not exhibit any significant activity.


The synthesis of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1) and its characterization by spectroscopic techniques (UV–Vis, FT-IR, 1H- and13C–NMR), HRMS and SC-XRD is accomplished. The structure was optimized by computational methods. X-rays analysis revealed that molecular configuration is stabilized by C–H…N and C–H…O bonding whereas crystal packing is stabilized by intermolecular N–H…O, C–H…N, C–H…O bonding and weak off-set π–π stacking interaction. Non-covalent interactions that are responsible for crystal packing are further explored by Hirshfeld surface analysis. Moreover, the experimental and theoretical spectroscopic results were in excellent agreement with each other. The MEP map shows that N2 and N3 atom are electron-donor while O1, O2, O3 and N4 atoms are electron-acceptor. The atoms which are electron rich or electron deficient are responsible for intra- and intermolecular attractive or repulsive interaction. The nearly planar geometry is also reflected by the low HOMO–LUMO gape, due to extended conjugation in the molecule. Theoretical data of compounds (2–6) has also been compared with that of (1) that indicates that –CH3, -Br and –OH substituents has almost comparable effect whereas –NO2 (Compound 1) substituent has great effect on overall electronic properties of the compound. The larger value of hardness than softness reflets greater stability of the compound that is important for technological applications. The compound (1) was also tested for antioxidant (% FRSA, TAC and TRP studies) and antimicrobial (antibacterial and antifungal) potential. It is concluded based on antioxidant results that compound (1) could be used as UV absorber in sunscreen lotions and may protect the skin damage from UV radiations.

Availability of data and materials

CCDC (Cambridge Crystallographic Data Center) number 1951603 contains the supplementary crystallographic data for the ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate (1). This data can be obtained free of charge via, or from the Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, UK; fax: (+ 44) 1223-336-033; or e-mail: While the computational data is presented in Additional files 1, 2.



Cambridge Crystallographic Data Centre


Density Functional Theory

DMSO-d 6 :

Dimethyl Sulfoxide-deuterated


Frontier Molecular Orbital


Free Radical Scavenging Assay


Fourier Transform Infrared


Highest Occupied Molecular Orbital


Hight Resolution Mass Spectrometry


Hirshfeld Surface


Integral-Equation-Formalism Polarizable Continuum Model


Lowest Unoccupied Molecular Orbital


Molecular Electrostatic Potential


Nuclear Magnetic Resonance


Oak Ridge Thermal Ellipsoid Plot


Single Crystal X-Rays Diffraction


Total Antioxidant Capacity


Time Dependent Density Functional Theory


Total Energy Distribution


Thin Layer Chromatography


Time of Flight-Electrospray Ionization


Total Reducing Power


Ultra violet-visible


  1. Siddiqui N, Arshad MF, Ahsan W, Alam MS. Thiazoles: A valuable insight into the recent advances and biological activities. Int J Pharm Sci Drug Res. 2009;1(3):136–43.

    CAS  Google Scholar 

  2. Hudson BM, Nguyen E, Tantillo DJ. The influence of intramolecular sulfur–lone pair interactions on small-molecule drug design and receptor binding. Org Biomol Chem. 2016;14:3975–80.

    CAS  PubMed  Google Scholar 

  3. Cunico W, Gomes CRB, Vellasco WTJ. Chemistry and Biological Activities of 1,3-Thiazolidin-4-ones. Mini Rev Org Chem. 2008;5:336–44.

    CAS  Google Scholar 

  4. Verma A, Saraf SK. 4-thiazolidinone–a biologically active scaffold. Eur J Med Chem. 2008;43:897–905.

    CAS  PubMed  Google Scholar 

  5. Hamama WS, Ismail MA, Shaaban S, Zoorob HHJ. Progress in the chemistry of 4- thiazolidinones. Heterocycl Chem. 2008;45(4):939–56.

    CAS  Google Scholar 

  6. Wells G, Bradshaw D, Diana P, Seaton A, Shi DF, Westwell AD. Antitumour benzothiazoles. Part 10: The synthesis and antitumour activity of benzothiazole substituted quinol derivatives. Bioorg Med Chem Lett. 2000;10(5):513–5.

    CAS  PubMed  Google Scholar 

  7. Hargrave KD, Hess FK, Oliver JT. N-(4-Substituted-thiazolyl)oxamic acid derivatives, new series of potent, orally active antiallergy agents. J Med Chem. 1983;26:1158–63.

    CAS  PubMed  Google Scholar 

  8. El-Sabbagh OI, Baraka MM, Ibrahim SM, Pannecouque C, Andrei G, Snoeck R, Balzarini J, Rashad AA. Synthesis and antiviral activity of new pyrazole and thiazole derivatives. Eur J Med Chem. 2009;44:3746–53.

    CAS  PubMed  Google Scholar 

  9. Chavan AA, Pai NR. Synthesis and antimicrobial screening of 5-arylidene-2-imino-4-thiazolidinones. ARKIVOC. 2007;14:148–55.

    Google Scholar 

  10. Bharti SK, Nath G, Tilak R, Singh SK. Synthesis, anti-bacterial and anti-fungal activities of some novel Schiff bases containing 2,4-disubstituted thiazole ring. Eur J Med Chem. 2010;42:651–60.

    Google Scholar 

  11. Shiradkar MR, Murahari KK, Gangadasu HR, Suresh T, Kalyan CA, Panchal D, Kaur R, Burange P, Ghogare J, Mokalec V, Raut M. Synthesis of new S-derivatives of clubbed triazolyl thiazole as anti-Mycobacterium tuberculosis agents. Bioorg Med Chem. 2007;15:3997–4008.

    CAS  PubMed  Google Scholar 

  12. Mjambili F, Njoroge M, Naran K, Kock CD, Smith PJ, Mizrahi V, Warner D, Chibale K. Synthesis and biological evaluation of 2-aminothiazole derivatives as antimycobacterial and antiplasmodial agents. Bioorg Med Chem Lett. 2014;24:560–4.

    CAS  PubMed  Google Scholar 

  13. Helal MHM, Salem MA, El-Gaby MSA, Aljahdali M. Synthesis and biological evaluation of some novel thiazole compounds as potential anti-inflammatory agents. Eur J Med Chem. 2013;65:517–26.

    CAS  PubMed  Google Scholar 

  14. Holla BS, Malini KV, Rao BS, Sarojini BK, Kumari NS. Synthesis of some new 2,4-disubstituted thiazoles as possible antibacterial and anti-inflammatory agents. Eur J Med Chem. 2003;38:313–8.

    CAS  PubMed  Google Scholar 

  15. Kalkhambkar RG, Kulkarni GM, Shivkumar H, Rao NR. Synthesis of novel triheterocyclic thiazoles as anti-inflammatory and analgesic agents. Eur J Med Chem. 2007;42:1272–6.

    CAS  PubMed  Google Scholar 

  16. Cohen A, Verhaeghe P, Crozet MD, Hutter S, Rathelot P, Vanelle P, Azas N. Tandem synthesis and in vitro antiplasmodial evaluation of new naphtho[2,1-d]thiazole derivatives. Eur J Med Chem. 2012;55:315–24.

    CAS  PubMed  Google Scholar 

  17. Fitzgerald G, Andzelm J. Chemical applications of density functional theory: comparison to experiment, Hartree-Fock, and perturbation theory. J Phys Chem. 1991;95:10531–4.

    CAS  Google Scholar 

  18. Ziegler T. Density functional theory as a practical tool for the study of elementary reaction steps in organometallic chemistry. Pure Appl Chem. 1991;63:873–8.

    CAS  Google Scholar 

  19. Andzelm J, Wimmer E. Density functional Gaussian-type-orbital approach to molecular geometries, vibrations, and reaction energies. J Chem Phys. 1992;96:1280.

    CAS  Google Scholar 

  20. Scuseria GE. Comparison of coupled-cluster results with a hybrid of Hartree-Fock and density functional theory. J Chem Phys. 1992;97:7528.

    CAS  Google Scholar 

  21. Dickson RM, Becke AD. Basis-set-free local density-functional calculations of geometries of polyatomic molecules. J Chem Phys. 1993;99:3898.

    CAS  Google Scholar 

  22. Johnson BG, Gill PMW, Pople JA. The performance of a family of density functional methods. J Chem Phys. 1993;98:5612.

    CAS  Google Scholar 

  23. Oliphant N, Bartlett RJ. A systematic comparison of molecular properties obtained using Hartree-Fock, a hybrid Hartree-Fock density-functional-theory, and coupled-cluster methods. J Chem Phys. 1994;100:6550.

    CAS  Google Scholar 

  24. Zhang Y, Guo ZJ, You XZ. Hydrolysis theory for cisplatin and its analogues based on density functional studies. J Am Chem Soc. 2001;123:9378–783.

    CAS  PubMed  Google Scholar 

  25. Haroon M, Akhtar T, Yousuf M, Baig MW, Tahir MN, Rasheed L. Synthesis, spectroscopic characterization and crystallographic behavior of ethyl 2-(4-methyl-(2-benzylidenehydrazinyl))thiazole-4-carboxylate: Experimental and theoretical (DFT) studies. J Mol Struct. 2018;1167:154–60.

    CAS  Google Scholar 

  26. Bruker,. APEX2 and SAINT. Madison, Wisconsin, USA: Bruker AXS Inc.; 2007.

    Google Scholar 

  27. Sheldrick GM. SHELXS-97. Program for the Solution of Crystal Structures: University of Gottingen, Germany; 1997.

    Google Scholar 

  28. Sheldrick GM. SHELXL-97. Program for Crystal Structures Refinement: University of Gottingen, Germany; 1997.

    Google Scholar 

  29. Spek AL. Structure validation in chemical crystallography. Acta Crystallogr. 2009;D65:148–55.

    Google Scholar 

  30. Farrugia LJ. WinGX suite for small-molecule single-crystal crystallography. J Appl Crystallogr. 1999;30:837–8.

    Google Scholar 

  31. Bruker BA, Inc., Madison, Wisconsin, USA, 2007 CrossRef CAS PubMed; G.M Sheldrick, A short history of SHELX. Acta Crystallogr Sect A Fundam Crystallogr. 2008; 64: 112–122.

  32. Bernstein J, Davis RE, Shimoni L, Chang NL. Patterns in hydrogen bonding: functionality and graph set analysis in crystals. Angew Chem Int Ed Engl. 1995;34:1555–73.

    CAS  Google Scholar 

  33. Wolff MA, Grimwood SK, McKinnon DJ, Turner JJ, Jayatilaka MJ, Spackman D. Crystal Explorer 17.5. University of Western Austr. 2012.

  34. McKinnon JJ, Jayatilaka D, Spackman MA. Towards quantitative analysis of intermolecular interactions with Hirshfeld surfaces. Chem Commun. 2007;2:3814–6.

    Google Scholar 

  35. Ashfaq M, Tahir MN, Kuznetsov A, Mirza SH, Khalid M, Al A. DFT and single crystal analysis of the pyrimethamine-based novel co-crystal salt: 2,4-diamino-5-(4-chloro-phenyl)-6-ethylpyrimidin-1-ium:4-hydroxybenzoate:methanol:hydrate (1:1:1:1) (DEHMH). J Mol Struct. 2020;1199:127041.

    CAS  Google Scholar 

  36. Ashfaq M, Bogdanov G, Glebov V, Ali A, Tahir MN, Abdullah S. Single crystal investigation, Hirshfeld surface analysis and DFT exploration of the pyrimethamine-based novel organic salt: 2, 4-diamino-5-(4-chlorophenyl)-6-ethylpyrimidin-1-ium 3-carboxybenzoate hydrate (1:1:1). J Mol Struct. 2021;1224:129309.

    CAS  Google Scholar 

  37. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Montgomery JA, Vreven T, Kudin KN, Burant JC, Millam JM, Iyengar SS, Tomasi J, Barone V, Mennucci B, Cossi M, Scalmani G, Rega N, Petersson GA, Nakatsuji H, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Klene M, Li X, Knox JE, Hratchian HP, Cross JB, Bakken V, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Ayala PY, Morokuma K, Voth GA, Salvador P, Dannenberg JJ, Zakrzewski VG, Dapprich S, Daniels AD, Strain MC, Farkas O, Malick DK, Rabuck AD, Raghavachari K, Foresman JB, Ortiz JV, Cui Q, Baboul AG, Clifford S, Cioslowski J, Stefanov BB, Liu G, Liashenko A, Piskorz P, Komaromi I, Martin RL, Fox DJ, Keith T, Al-Laham MA, Peng CY, Nanayakkara A, Challacombe M, Gill PMW, Johnson B, Chen W, Wong MW, Gonzalez C, Pople JA. Gaussian 03, Revision E.01. Wallingford, CT: Gaussian Inc.; 2009.

    Google Scholar 

  38. Asghar A, Yousuf M, Mubeen H, Nazir R, Haruna K, Onawole AT, Rasheed L. Synthesis, spectroscopic characterization, molecular docking and theoretical studies (DFT) of N-(4-aminophenylsulfonyl)-2-(4-isobutylphenyl) propanamide having potential enzyme inhibition applications. Bioorg Med Chem. 2019;27(12):2397–404.

    CAS  PubMed  Google Scholar 

  39. Dennington R, Keith T, Millam J. Gauss View, Version 412. Shawnee Mission: Semichem Inc.; 2007.

    Google Scholar 

  40. Becke AD. Density-functional thermochemistry. III. The role of exact exchange. J Chem Phys. 1993;98:5648.

    CAS  Google Scholar 

  41. Ditchfield R, Hehre WJ, Pople JA. Self-Consistent Molecular-Orbital Methods. IX. An extended gaussian-type basis for molecular-orbital studies of organic molecules. J Chem Phys. 1971;54:724.

    CAS  Google Scholar 

  42. Dennington I, Keith RT, Millam J, Eppinnett K, Hovell W. GaussView. Shawnee Mission: Semichem Inc; 2003.

    Google Scholar 

  43. García-Melchor M, Braga AA, Lledós A, Ujaque G, Maseras F. Acc, Computational Perspective on Pd-Catalyzed C-C cross-coupling reaction mechanisms. Chem Res. 2013;46:2626–34.

    Google Scholar 

  44. Fleming I. Frontier Orbitals and Organic Chemical Reactions. London: Wiley; 1976.

    Google Scholar 


  46. Rauhut G, Pulay P. Transferable scaling factors for density functional derived vibrational force fields. J Phys Chem. 1995;99:3093–100.

    CAS  Google Scholar 

  47. Baker J, Jarzecki AA, Pulay P. Direct scaling of primitive valence force constants: an alternative approach to scaled quantum mechanical force fields. J Phys Chem. 1998;102:1412–24.

    CAS  Google Scholar 

  48. George S. Infrared and Raman Characteristics Group Frequencies. New York: Wiley; 2001.

    Google Scholar 

  49. Sen F, Ekici O, Dincer M, Cukurovali A. A comparative study on 4-(4-(3-mesityl-3-methylcyclobutyl)thiazole-2-yl)-1-thia-4-azaspiro[4.5]decan-3-one: Experimental and density functional methods. J Mol Struct. 2015;1086:109–17.

    CAS  Google Scholar 

  50. Sen F, Dincer M, Cukurovalı A, Yılmaz I. N-[4-(3-methyl-3-mesityl-cyclobutyl)-thiazol-2-yl]-succinamic acid: X-ray structure, spectroscopic characterization and quantum chemical computational studies. J Mol Struct. 2013;1046:1–8.

    CAS  Google Scholar 

  51. Cances E, Mennucci B, Tomasi J. A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics. J Chem Phys. 1997;107:3032.

    CAS  Google Scholar 

  52. Pearson RG. Absolute electronegativity and hardness correlated with molecular orbital theory. Proc Natl Acad Sci. 1986;83:8440.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Chattaraj PK, Sarkar U, Roy DR. Electrophilicity Index. Chem Rev. 2006;106:2065.

    CAS  PubMed  Google Scholar 

  54. Muthu S, Porchelvi EE, Karabacak M, Asiri AM, Swathi SS. Synthesis, structure, spectroscopic studies (FT-IR, FT-Raman and UV), normal coordinate, NBO and NLO analysis of salicylaldehyde p-chlorophenylthiosemicarbazone. J Mol Struct. 2015;1081:400–12.

    CAS  Google Scholar 

  55. Sen F, Dincer M, Cukurovali A. Structural and spectroscopic characterization of 4-(3-methyl-3-phenylcyclobutyl)-2-(2-propylidenehydrazinyl)thiazole: A combined experimental and DFT analysis. Spectrochemica Acta Part A: Mol Biomol Spec. 2015;150:257.

    CAS  Google Scholar 

  56. Yousuf M, Youn IS, Yun J, Rasheed L, Valero R, Shi G, Kim KS. Violation of DNA neighbor exclusion principle in RNA recognition. Chem Sci. 2016;7:3581–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Haroon M, Khalid M, Akhtar T, Tahir MN, Khan MU, Muhammad S, Al-Sehemi AG, Hameed S. Synthesis, crystal structure, spectroscopic, electronic and nonlinear optical properties of potent thiazole based derivatives: Joint experimental and computational insight. J Mol Struct. 2020;1202:127354.

    CAS  Google Scholar 

  58. Haroon M, Khalid M, Akhtar T, Tahir MN, Khan MU, Saleem M, Jawaria R. Synthesis, spectroscopic, SC-XRD characterizations and DFT based studies of ethyl2-(substituted-(2-benzylidenehydrazinyl))thiazole-4-carboxylate derivatives. J Mol Struct. 2019;1187:164.

    CAS  Google Scholar 

  59. Zahra SSAM, Qasim M, Gul B, Zia M, Mirza B, Haq IU. Polarity based characterization of biologically active extracts of Ajuga bracteosa Wall ex Benth and RP-HPLC analysis. BMC Compl Altern Med. 2017;17:443–58.

    Google Scholar 

Download references


The corresponding author is grateful to Prof. Dr. Shahid Hameed, Department Chemistry, QAU, Islamabad-45320, Pakistan for providing NMR facility.


No funding available for this research work.

Author information




MH has conducted the research, compiled the experimental data and write up of synthetic part. TA: Concept idea and supervisor of this research and paper write up. MY and LR performed DFT Quantum chemical calculation and write up of that part. MNT and MA collected XRD single crystal data and Hirshfeld surface analysis and write up. SZ and IH carried out biological studies and write up. All authors read and approved final manuscript.

Corresponding authors

Correspondence to Tashfeen Akhtar or Muhammad Yousuf.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

No competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Synthesis, crystal structure, Hirshfeld surface investigation and comparative DFT studies of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate. Table S1. Comparison of the experimental and calculated vibrational wavenumbers (cm-1) of compound 1 in gas phase. Table S2. Experimental and theoretically calculated 1H NMR. Table S3. Experimental and theoretically calculated 13C NMR.

Additional file 2.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Haroon, M., Akhtar, T., Yousuf, M. et al. Synthesis, crystal structure, Hirshfeld surface investigation and comparative DFT studies of ethyl 2-[2-(2-nitrobenzylidene)hydrazinyl]thiazole-4-carboxylate. BMC Chemistry 16, 18 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Thiazole
  • DFT calculations
  • Anti-oxidant
  • XRD
  • Spectroscopy