A comparison between observed and DFT calculations on structure of 5-(4-chlorophenyl)-2-amino-1,3,4-thiadiazole

The crystal and molecular structure of 5-(4-chlorophenyl)-2-amino-1,3,4-thiadiazole 3 was reported, which was characterized by various spectroscopic techniques (FT-IR, NMR and HRMS) and single-crystal X-ray diffraction. The crystal structure 3 (C8H6ClN3S) crystallized in the orthorhombic space group Pna21 and the unit cell consisted of 8 asymmetric molecules. The unit cell parameters were a = 11.2027(2) Å, b = 7.6705(2) Å, c = 21.2166(6) Å, α = β = γ = 90°, V = 1823.15(8) Å3, Z = 8. In addition, the structural geometry (bond lengths, bond angles, and torsion angles), the electronic properties of mono and dimeric forms of compound 3 were calculated by using the density functional theory (DFT) method at B3LYP level 6-31+ G(d,p), 6-31++ G(d,p) and 6-311+ G(d,p) basis sets in ground state. A good correlation was found (R2 = 0.998) between the observed and theoretical vibrational frequencies. Frontier molecular orbitals (HOMO and LUMO) and Molecular Electrostatic Potential map of the compound was produced by using the optimized structures. The NBO analysis was suggested that the molecular system contains N-H…N hydrogen bonding, strong conjugative interactions and the molecule become more polarized owing to the movement of π-electron cloud from donor to acceptor. The calculated structural and geometrical results were in good rational agreement with the experimental X-ray crystal structure data of 1,3,4-thiadiazol-2-amine, 3. The compound 3 exhibited n→π* UV absorption peak of UV cutoff edge, and great magnitude of the first-order hyperpolarizability was observed. The obtained results suggest that compound 3 could have potential application as NLO material. Therefore, this study provides valuable insight experimentally and theoretically, for designing new chemical entities to meet the demands of specific applications.

antibacterial activities 21 . In addition, these thiadiazole compounds have also been a valuable source as reaction intermediates in the synthesis of commercially available drugs (megazol and furidiazine) 22 .
In view of the broad ranging applications of 2-amino-1,3,4-thiadiazole scaffold and in extension to our efforts to develop various biologically active conjugates 23 , we conducted the comparative studies between experimental and theoretical (DFT) data of 5-(4-chlorophenyl)-2-amino-1,3,4-thiadiazole, a molecule with anti-proliferative activity, which was earlier reported by Guan et al. 24 . The crystal structure and theoretical calculations of this compound have not been reported yet. Here, we have synthesized the 5-(4-chlorophenyl)-1,3,4-thiadiazol-2-amine (3) according to the reported literature, and characterized by using crystal X-ray, IR, NMR, HRMS. In addition, the optimized geometry parameters of the title compound was carried out at DFT and TD-DFT methods with the B3LYP functional and two basis sets 6-31+ G(d,p) and 6-311+ G(d,p), and conceivable correlations were investigated between observed and theoretical data.

Results and Discussions
X-ray crystal structure description and synthesis of compound 3. Figure 1 depicts the synthesis of 5-(4-chlorophenyl)-1,3,4-thiadiazolo-2-amine (3) compound. It was achieved by the reaction between the corresponding 4-chlorobenzoic acid (1) and thiosemicarbazide (2). In this process, using phosphorous oxychloride as an acid activating agent, the reaction mixture was refluxed for 4 h as per reported procedure 25 . The structural assignment of the purified compound (3) was established and confirmed by employing various spectroscopic techniques. The identified title molecule (3) was subjected to single-crystal X-ray diffraction analysis. The comprehensive X-ray crystallographic evidences were placed in Cambridge crystallographic data center (CCDC 1879792). A single colorless plank-shaped crystal of compound,3 was recrystallized through sluggish evaporation of hot ethanol over a period of five days. For the crystallographic analysis, Fig. 2 represented the 50% probability of ORTEP image and the structural established crystal data are illustrated in Table 1. The synthesized title molecule (3) develops in the orthorhombic space group Pna2 1 by Z = 8, and Table 2 represented the measured chosen torsion angles, bond distances and angles. The crystal compound embraces of two independent similar molecules in the asymetric unit and the heteroaromatic rings were planar. The dimers with 4-chlorophenyl and 1,3,4-thiadiazol-2-amine rings were bonded through C2-C3 and C10-C11 junction with a bond length of 1.464(5) and 1.469(5) Å. These were approximately equal bond lengths, which were within estimated range for C-C (1.467(4) Å) 26,27 (Fig. 2). The bond angles of these rings exhibited 121.0(3)° for C4-C3-C2 and 120.6(3)° for C12-C11-C10. Moreover, the C=N double bonds were conjugated with the unshared electron pairs on the N1 and S1 atoms of compound 3, and the corresponding endocyclic (N=C) bond lengths were almost similar with N2-C1 (1.316(4)), N3-C2 (1.296(4)), N5-C9 (1.324 (4)) and N6-C10 bonds (1.303(4)) Å. The carbon and carbon bonds (C-C) in the 4-chlorophenyl ring were of almost equal lengths in the range of 1.371 to 1.399 Å. In addition, similar conformations of four sets of asymmetric units were identified in the unit cell along the b axis of crystal packing diagram, which were linked over intermolecular symmetric hydrogen bond interactions (N-H…N) (Figs. 2 and 3). The strongest intermolecular hydrogen bonding existed between the protons of the exocyclic amino moiety and the nitrogen atoms of the 1,3,4-thiadiazole ring with the bond distance of 2.11 Å N(4)-H(4 A)…N(2) (167.3°)  DFT optimization investigation. The most appropriate density functional theory based function is B3LYP with 6-31+ G(d,p) and 6-311+ G(d,p) basis sets were used for creating the optimized geometries of monomer and dimeric form of molecule 3 in the gas phase (Fig. 5). The computed geometrical parameters (torsion angles bond lengths and angles) are represented in Table 2 and these parameters were compared with the found experimental www.nature.com/scientificreports www.nature.com/scientificreports/ parameters. The identified bond lengths were good reliable with the experimental parameters through an average divergence of 0.01 Å to 0.03 Å for both basis sets. The deviations could take place by the intramolecular hydrogen bonding interactions. It was concluded that the sequences of B3LYP functional with both 6-31+ G(d,p) and 6-311+ G(d,p) basis sets produce the best compatible results with the experimental values ( Table 3).
The calculation of the vibrational frequencies in IR-spectrum of monomer and dimeric forms of molecule 3 was carried out by using B3LYP function with 6-31+ G(d,p) basis set. The chosen characteristic vibrational band values and along with observed values are presented in Table 4. In IR spectrum, the N-H stretching frequency was appeared at 3255 cm −1 , this same band was identified at 3219 cm −1 and 3228 cm −1 for dimeric and monomeric of compound 3. The computed vibrational value was in excellent agreement with the observed value. In addition, the benzene ring of C-H stretching IR band was seemed at 3082 cm −1 (calculated as monomer 3189 cm −1 and dimer 3188 cm −1 ). The most prominent vibrational stretching absorbance was found at 1629 cm −1 due to the endocyclic group C=N group of the thiadiazole ring and the computed as 1641 cm −1 and 1652 cm −1 for monomer and dimer forms. The strong stretching band observed for C=N group, this may be due to the neighbouring atoms of electronic effects (unshared electron pairs on the S, N atoms) and intermolecular hydrogen bonding 28 Table 3. Intermolecular hydrogen bonding geometry of the compound 3 (Å, °). Symmetry codes: a + X,−1 + Y, +Z; b + X,1 + Y, +Z. www.nature.com/scientificreports www.nature.com/scientificreports/ at 735 cm −1 . The computed C-Cl stretching absorbance band (725 cm −1 for monomer and 724 cm −1 for dimer) was in good agreement with the observed value (706 cm −1 ). The correlation studies were summarized in Fig. 6 and these results indicates the good correlation founded (R 2 = 0.998) amongst the experimental and calculated vibrational absorbance values. The acquired vibrational frequency values were consistence with results obtained in the literature 29,30 .
The calculated 1 H and 13 C NMR spectrum of the structure 3 was analyzed by using B3LYP/6-31+ G(d,p) in DMSO and the chemical shift ( 1 H and 13 C NMR) values are presented in Table S1. A good agreement between the observed and calculated chemical shift values, except for nitrogen atom bonding with H1A and H1B protons. The calculated deviation was observed due to the hydrogen bonding interactions (N-H…N) and was given in Table 3.

Molecular electrostatic potential (MEP) investigation. The MEP considered as an interaction energy
at a certain zone of a structure was involved in the electrical charge distribution from the proton, nuclei and electrons, which were positioned at r.
The V(r) was obtained from the following the Eq. (1) 31 .  Table 4. Most characteristic experimental and calculated vibrational frequencies (cm −1 ) and assignments of the compound 3 at B3LYP/6-31+G(d,p). www.nature.com/scientificreports www.nature.com/scientificreports/ Where Z A signifies the charge of the nucleus A, which is positioned at R A , (r′) denote the dummy integration variable and ρ(r′) represent the electron density function of the structure.
The probable nucleophilic (blue region) and electrophilic (red region) attack sites were established. The MEP was computed at the B3LYP functional and 6-31+ G(d,p) level of theory of monomeric and dimeric forms, and the images displayed in Fig. 7, which gives the valuable evidence about the intramolecular hydrogen bonding interactions and the possible nucleophilic and electrophilic sites 32 . The higher negative region was placed on the N-N atoms of the thiadiazole ring (red, Fig. 7), which was regarded as the feasible position for electrophilic attack with V(r)of −0.046 a.u, due to the N-H-N hydrogen bonding interfaces around this zone. The supreme positive zone was positioned at the amino group of nitrogen and sulfur atoms of thiadiazole ring of the molecule with V(r) of 0.046 (Fig. 7). This shows the probable site for nucleophilic attack on unshared electron pair atoms (N and S).
Frontier molecular orbital (FMO) analysis. The electric and optical parameters of the structure in the FMO theory was computed by using B3LYP/6-31+G(d,p) and B3LYP/6-311+G(d,p) basis sets. FM orbitals are the important orbitals in a structure, and it is considered as the lowest molecular orbitals (LUMO) and the highest molecular orbital (HOMO) 29,33 . The outermost orbital occupied by electrons called the HOMO, which can act as an electron donors. The energy of the LUMO is the foremost vacant innermost orbital unoccupied with electrons, it performs as an electron acceptor. The electrons donating capacity of molecule is allied to the E HOMO , and the over-all the greater the HOMO energy (smaller negative value), the superior the capacity to donate electrons 29,34 . The energy difference (ΔE) between the HOMO and LUMO, softness, electronegativity, chemical hardness and  www.nature.com/scientificreports www.nature.com/scientificreports/ electrophilicity index values of monomer and dimer of compound 3 are presented in Table 5 and Fig. 8. The HOMO orbital was displayed charge density localized around the 2-amino-1,3,4-thiadiazole ring, while a charge density localized on the 4-chlorophenyl ring in the LUMO orbitals of the both monomer and dimeric forms of both basis sets, signifying that these moieties can incentive the electron transition. The energy gap between the HOMO and LUMO orbitals indicates the chemical strength and reactivity of the molecule 35,36 . The energy gap difference between the HOMO and LUMO is ΔE of 3.03 and 3.01 eV for monomer of and ΔE of 2.59 and 2.93 eV for dimer of B3LYP/6-31+ G(d,p) and B3LYP/6-311+ G(d,p) basis sets. The energy of the smaller band space gap increases the reactive nature and stability of compound 3. In addition, the compound characterized based on reactivity, chemical hardness (η) electronegativity (μ), softness (ζ) and electrophilicity index (Ψ) were determined according to Eqs. (2)(3)(4)(5).
Where A and I is electron affinity and ionization potential.
As can be seen in the Table 5, the title compound can be referred as a soft molecule, since lesser energy gap systems are termed as soft molecules. The chemical behavior of the compound was identified by electronegativity factor. The greater value of the electronegativity of molecule 3 is a designated of its chemical reactivity.
Electronic absorption spectrum analysis. The electronic spectral analysis of the title structure was computed in both gas and various organic solvent phases by using time dependent-DFT (TD-DFT) IEFPCM model at the B3LYP/6-31++G(d,p) level theory on the ground state optimized geometry. Although, the observed absorption spectra were investigated at room temperature in the diverse solvents with 5×10 −5 M of a concentration, and the plots are presented in the Figs. 9a,b. The experimental maximum wavelength (λ max ) and absorbance was  Table 5. The calculated energies of the compound 3. www.nature.com/scientificreports www.nature.com/scientificreports/ measured at room temperature and the values are summarized in Table 6. Dipole moment (μ, D), major absorption energy E, oscillator strength (f), maximum wavelength (λ max ), electronic transition of excitation energy and atomic orbital contribution were calculated in solvent and gas phases, and the values are demonstrated in Table 7.
The observed absorption spectrum of the title compound was recorded in various solvents such as hexane, dichloromethane (DCM), acetonitrile, ethanol and DMSO, concluded that the very low solubility was observed in hexane solvent. As results from Fig. 9a, the molecule was showed high intensity band at 290, 297, 300, 307 and 314 nm in all the five solvents, which can be attributed to the n→π* electronic transition. This was occurred from the molecular orbitals of the auxochrome (NH 2 ) of exocyclic nitrogen lone pairs and the conjugated π-bond of azomethine (-C=N-N=C-) of 1,3,4-thiadiazole moiety. With the increasing the solvent polarity from non-polar to polar, n→π* transition get slightly shifted towards longer wavelength, which was probably due to the solvatochromic effect (Tables 6 and 7) 37 . Compound 3 displayed a positive absorption solvatochromism, which signifies that the absorption was redshifted with cumulative solvent polarity. The highest energy transition was

S.No Solvent
Wavelength λ max (nm) Absorbance    Consequently, the calculation was further engaged with TD-DFT-IEP-PCM model 38 . A strong absorption band was observed at 280 nm (f = 0.5110) and 285 nm (f = 0.6398) in the gas and solvent phase (DMSO), which corresponds to the excitation of electrons occurred from HOMO to LUMO transition (96.71% and 97.96% orbital contribution). It could be owing to the solvent effects, which can impact the geometry, molecular properties and electrostatic structure, thus causing shift in absorption bands 39 . The electronic transitions occurred from dipolar interaction with solvents of different polarities suggest that the intramolecular charge transfer (ICT) of the emission state in which the HOMO and LUMO were apparently localized on the 1,3,4-thidiazole moiety and the phenyl ring.
The distribution of energy levels and electronic transitions of molecular orbitals in the title compound in both solution and gas phases are overlapping. The frontier MO analysis in calculating the energies and the nature of frontier MOs indicate that these orbitals are the maximum feasible sites of interchange density, directing to the interaction between orbitals. There was a small variance in the energy gap in gas and solvent phases (Fig. S1). The estimated energy gap (ΔE = HOMO-LUMO) value in solvent phase for DMSO and DCM were 3.45 and 3.47 eV respectively, which corresponds to the intramolecular charge transfer (S 0 →S 1 ) transition (Fig. 10a). It may be assumed that for DCM and DMSO could stabilize the compound, due their different polarities and dielectric constants by elevating the energy gap between their HOMO and LUMOs. The energy gap decreased with increasing solvent polarity, which shifted the absorption peak to longer wavelength than 290-314 nm, consequently leading to decline in their absorption capacity in the desired UV range. The HOMO orbital shows that π-bonding and non-bonding (n-type) nature orbitals were identified on the thiadiazole moiety, phenyl ring and nitrogen atom of exocyclic amino group, the HOMO-1 orbital was localized on the 2-amino-1,3,4-thiadiazole moiety and chlorine atom, and the HOMO-2 orbital was localized on the 2-amino-1,3,4-thiadiazole moiety only. Similarly, the LUMO shown the π* anti-bonding type orbitals on carbon, sulfur and nitrogen atoms were localized in the phenyl and thiadiazole ring moieties. Therefore, the results concluded that the solvent polarity enhances the probability of transition (n→π*) in the present compound. A very good correlation was observed between the TD-DFT and experimental absorption spectra results.

Mulliken atomic charges. Mulliken atomic charge results takes an imperative position in the applications
of quantum chemical calculation, as the molecular polarizability, atomic charges affect the dipole moment and electronic structural features of molecular systems. The development of donor and acceptor pairs of charge transfer was succeeded, because of charge distributions over the atoms in the present molecule. Mulliken population investigation of monomer and dimeric forms of compound 3 was estimated by using function B3LYP with 6-31+ G(d,p) basis set and results are demonstrated in Tables S2 and S3 (ESI). Figure 11 illustrates the charge distribution structure of compound 3. The charge distribution presented that the all hydrogen atoms were concentrated on positively charged. The H16 (+0.322236e) and H17 (+0.298055e) atoms in particularly have higher positive atomic chargers than the other hydrogen atoms in title molecule. This was due to the presence of electron negative nitrogen (N15) atom of the exocyclic amino group. These hydrogen atoms appeal to the positive charge from the electron-withdrawing nature of the nitrogen atom. Hydrogen attached carbon atoms (C5 and C9) and quaternary carbon atoms (C1, C2, C4 and C11) exhibited negative chargers, while, the enduring carbon atoms appeared as positively charged. Another positive charge was noted on the sulfur atom, which was in the geometric center of 1,3,4-thiadiazole moiety. Nitrogen atoms (N15 and N19) were considered as more basic sites and displayed the negative atomic chargers. The existence of net positive charge on the hydrogen atom and high negative charge on the nitrogen atom, could facilitate the charge transfer over the formation of intermolecular hydrogen bonds (N-H-N). Similarly, the dimeric form of compound 3 was also calculated and same pattern of charge distribution as related to monomeric form was found (Figs. S2-4). NBO analysis. Furthermore, the natural bond orbital (NBO) study of the title molecule 3 was computed at B3LYP/6-31+ G(d,p) level of model. To understanding into delocalization of electron density, intra-and inter-molecular bonding relations between the bonds, and investigate the charge transfer (CT) or conjugative relations in molecular structure 40 . The structure of the calculated compound was exhibited a type of Lewis www.nature.com/scientificreports www.nature.com/scientificreports/ structure (97.7%), 2.12% of valence non-Lewis, and 0.16% of Rydberg non-Lewis ( Fig. 12 and Table S4). The natural Rydberg basis (NRB) populations and natural minimum basis (NMB) was also calculated in the reported compound, which was exhibited the 99.6% of maximum percentage contribution of the NMB to the molecular charge allocation.
The second-order perturbation theory study of the Fock matrix was used to assess the donor-acceptor interactions in the NBO investigation. The strength of interaction stabilization energy E (2) allied with electron delocalization among the donor (i) and acceptor (j) (i→j) was assessed by the following Eq. (6).
where q i is the donor orbital occupancy, F(i, j) is the off diagonal NBO Fock matrix element, and ε j and ε i are implies diagonal elements (orbital energies). The greater stabilization energy E (2) value indicates more intensive interaction amongst the electron acceptors and donors, i.e. the high electrons donating ability and superior degree of conjugation of the entire system 39 . Numerous donor and acceptor interactions were identified for the phenyl and 1,3,4-thiadiazole rings, and lone pairs (n) of the sulfur, chlorine and nitrogen atoms. The potential intensive interactions and electron dentistry transfer from lone pair electrons to the anti-bonding orbitals along with the stabilization energy values of the compound 3 were tabulated in Table 8. These results reveal that the intramolecular interactions have occurred due to the overlap between bonding lone pair of nitrogen (LP(1)N15) and sulfur (LP(2)S3) and antibonding C2-N19 orbitals, which results in the intramolecular charge transfer (ICT) initiated stabilization of the molecular  www.nature.com/scientificreports www.nature.com/scientificreports/ structure. Within the 1,3,4-thiadiazole fragment the n(LP(1)N15)→π*(C2-N19) and n(LP(2)S3)→π*(C2-N19) interactions facilitated strong stabilization to the molecular system of the studied compound by 156.64 kJ/mol and 114.72 kJ/mol respectively, which result in increase in electron density distribution and thus, failing of the respective bonds. This higher energy provided the stabilization to the molecular structure and these interactions suggestively impact the crystal packing with this molecule. In parallel technique calculations, the electron transfer Donor NBO (i) Acceptor NBO (j) E (2)a (kJ/mol) E(j) − E(i) b (a.u.) F(i,j) c (a.u.) π (C1-S3) π* (C2-N15)  Table 8. Second order perturbation theory analysis of Fock matrix in NBO basis of compound 3. a E (2) means energy of hyper conjugative interactions (stabilization energy in kJ/mol). b Ej − Ei is energy difference between donor and acceptor i and j NBO orbitals. c F(i,j) is the Fock matrix element between i and j NBO orbitals.
www.nature.com/scientificreports www.nature.com/scientificreports/ from π(C4-C5)→π*(C7-C11) gave a stabilization energy of 88.74 kJ/mol. This interaction was owing to the π electron delocalization in a molecular system. NBO analysis further provided an evidence for intermolecular hydrogen bonding, the stabilization of energy of 3.51 kJ/mol (n→σ*) obtained from lone pair of nitrogen (LP(1) N19) to σ*(N15-H16), signifying the existence of intermolecular N-H…N hydrogen bonding. The charge delocalization occurred from the lone pair of nitrogen atom to C-S bond through n→π* interaction. The involved N18 atom lone pair 1 and N19 atom lone pair 1 to π*(C1-S3) and π*(C2-S3) interactions were provided corresponding stabilization energies of 69.07 kJ/mol and 65.35 kJ/mol, respectively, which displayed hyper conjugative interaction. Thus, the NBO study suggested that presented molecule contains strong conjugative interactions and N-H…N hydrogen bonding, and the molecule become a highly polarized because of the exchange of π-electron cloud from donor to acceptor.

Nonlinear optical properties (NLO).
We further extended the study for the computing of polarizability (α) and hyperpolarizability (β), which were significant to envisage the nonlinear optical properties of the molecules. It is an emerging tool in telecommunications and signal processing field. In recent years, organic non-linear materials have received a considerable attention in various optoelectronic technologies. Materials of nonlinear optical properties were instigated by intramolecular charge transfer and delocalization of the electrons from the donor to acceptor electron groups in conjugated π-electron systems, which have promising materials for NLO applications 41 .
Therefore, the nonlinear optical properties of the organic π-conjugated molecule 3 were studied at the B3LYP function with the 6-31 G(d) basis set in the gas phase. The first-order hyperpolarizability (β tot ), average linear polarizability (α) and total static dipole moment (μ) of this molecular system were predicted by using Eqs. (7, 8 and 9) and results are illustrated in Table 9. For molecule 3, the computed polarizability (α) and first hyperpolarizability (β tot ) were 18.079 × 10 −24 esu and 2.208 × 10 −30 esu respectively. The β tot of investigated compound was approximately 6 times more than for urea (0.3728 × 10 −30 esu). Urea is a prototype NLO material and most often used reference material for comparative purpose 42 . The small energy gap (ΔE) of compound 3 indicates the occurrence of molecular charge transfer, and it can be penetrating to the stimulus of the polar environment. Thus, it led to easy electron distribution and high polarization. Small energy offers maximal NLO response 43 , hence great magnitude of the first hyperpolarizability of the title molecule might potentially serve as NLO material.

Density of states (DOS).
Density of states and partial density of states (PDOS) predictions were carried out by using GaussSum 43 . The TDOS plot, a population study per orbital reveals the composition of the molecular orbitals in a specific energy range very clearly (Fig. 13a) Table 9. The molecular electric dipole moments (μ), Polarizability (α) and hyperpolarizability (β tot ) values of compound 3.
compound 3 are illustrated in Fig. 13b and Fig. S5 (ESI). Five portions in the energy range from 0 to −20 eV were below the Fermi level of valence band (VB). The band (1) between −4 eV to 13.5 eV was formed by the hybridization of C 2p, N 2p, S 3p, Cl 2p and H 1 s orbitals. Likewise, the commencing energy band lying among −4 eV and the Fermi level was dominated by C 2p orbitals with partially contribution of N 2p and Cl 3p orbitals. One narrow band with a sharp peak (2) was about −14.5 eV due to the N 2p orbitals with partly aid of C 2p and S 3p orbitals. The lower parts (3)(4)(5) of the VB between −16, 18.5 and 20 eV were made up majorly of the C 2p and H 1 s orbitals as of strong hybridization of the C and H orbitals in the C 6 H 4 group. We detailed the experimental and theoretical characterizations of the organic compound, 1,3,4-thiadiazol-2-amine and analyzed its UV absorption as function of solvent effect. The interplay of several aspects, such as crystal structure, chemical bonding and thermodynamic properties of the compound was paramount. A good rational correlation between the suitable structural and geometrical response of experimental and theoretical parameters was observed. The compound 3 exhibited n→π* UV absorption of UV cutoff edge, and great magnitude of the first-order hyperpolarizability was observed. The obtained results suggest that compound 3 could have potential application as NLO material. Therefore, this study provides valuable insight both experimentally and theoretically, for designing new chemical entities to meet the demands of specific applications.

conclusions
Compound 5-(4-chlorophenyl)-2-amino-1,3,4-thiadiazole (3) was synthesized from 4-chlorobenzoic acid and thiosemicarbazide in POCl 3 . The structure was characterized by various spectroscopic techniques (FT-IR, NMR and HRMS) and single-crystal X-ray diffraction. XRD analysis revealed that 3 crystallizes in Pna2 1 orthorhombic space group. The electronic and geometrical parameters of monomeric and dimeric forms of title molecule 3 was computed by using DFT/B3LYP/6-31+ G(d,p), B3LYP/6-31++ G(d,p) and B3LYP/6-311+G(d,p) basis sets in ground state. The observed and computed values were compared. For mono and dimeric forms energy gap (HOMO-LUMO) of compound were 3.04 eV and 2.60 eV. The low energy gap of HOMO and LUMO showed the stability of compound 3. Moreover, the MEP shows that the positive potential sites were around the unshared electron pairs of sulfur and nitrogen atoms. Based on UV absorption the highest energy transition was found in polar solvents such as ethanol and DMSO, which resulted in increase in dipole moment (for ethanol 3.60 D and DMSO 3.65 D). In addition, the NBO analysis suggested that molecular system contains N-H…N hydrogen bonding and strong conjugative interactions and molecule become more polarized owing to the movement of π-electron cloud from donor to acceptor. The theoretical results for both basis sets were good linear agreement with the experimental X-ray crystal structure data of the title compound.

Methods
All the chemicals (laboratory grade) and solvents were purchased from Sigma Aldrich and Merck and used without any further purification. NMR analysis was recorded on a Bruker AVANCE III 400 MHz spectrometer (399.995 MHz for 1 H and 100.4296 MHz for 13 C), and chemical shifts (δ) values were presented in parts per million (ppm). The δ values of DMSO-d 6 representing at 2.50 ppm for 1 H and 39.51 ppm for 13 C NMR, were used for the structure elucidation. Functional groups of molecules were identified by infrared spectrum (IR) spectrum, which was recorded on a Perkin Elmer Spectrum 100 FT-IR Spectrometer. Further, high-resolution mass data was done by using a Bruker microTQF-Q II ESI instrument operated at ambient temperature.
A colorless plank-shaped crystal (3) with a dimension 0.27 × 0.17 × 0.09 mm 3 was selected and mounted on a MITIGEN holder in paratone oil. The appropriate crystal was recorded on a 'Bruker SMART APEX-II CCD' single crystal diffractometer. The crystal was kept at T = 100 (2) K during data assortment. Data was measured using ω and φ scans using MoKα radiation. The total number of cycles and images were based on the approach calculation from the program COSMO (BRUKER) and the maximum resolution that was achieved was Θ = 28.415° (0.75 Å). The structure was solved with the ShelXS 44 structure solution program using direct methods and by using Olex2 45 as the graphical interface. The model was refined with the version 2016/6 ShelXL 46 using Least www.nature.com/scientificreports www.nature.com/scientificreports/ Squares minimization. All non-hydrogen atoms were refined anisotropically. Hydrogen atom positions were intended geometrically and refined using the riding model.
The molecular structure in ground state of the title compound was fully optimized using density functional theory (DFT). B3LYP functional 47 in combination with 6-31+ G(d,p), 6-31++ G(d,p), and 6-311+ G(d,p) basis set was employed. The optimization was performed both in gas and solvent phases. The molecular electrostatic potential map (MEP) was also predicted by using the Merz−Kollman procedure 48 at the 6-31+ G(d,p) level and with iso-surface of the electron density (0.004 electrons per Å 3 ). The 1 H and 13 C chemical shifts, derived on a δ-scale in relation to the TMS, (reference) were calculated via a method developed by Gregory et. al. 49 The Gauge-Independent Atomic Orbital (GIAO) at B3LYP/6-31+ G(d,p) was applied for estimating the 1 H and 13 C NMR of the DMSO solvent. SCRF-TD-DFT/B3LYP/6-31++ G(d,p) was used for calculating the UV-electronic absorption spectrum. Frontier molecular orbital (FMO), associated to the highest molecular orbital (HOMO) and the lowest molecular orbital (LUMO) calculations were done to determine the reactivity nature 50,51 of the compound 3. B3LYP/sto-3g basis set and GaussSum 43 was used for calculating the density of states spectra. The Gaussian03 program package and Gauss-View molecular visualization program was used for all the calculations 52 .
Procedure for the synthesis of 5-(4-chlorophenyl)-2-amino-1,3,4-thiadiazole (3). A mixture of 4-chlorobenzoic acid 1 (0.10 mmol) and thiosemicarbazide 1 (0.10 mmol) in 30 mL of phosphorous oxychloride was added dropwise with continuously stirring and refluxed gently for 1 h. The reaction mixture was permitted to cool and followed by carefully addition of water (50 mL). Again, the reaction mixture was refluxed for another 3 h. The movement of the reaction was monitored by TLC. After the completion of the reaction, the reaction mixture was basified with aqueous 30% NaOH (pH = 8). The separated solid was filtered and washed with cold water and the obtained white solid was recrystallized with hot ethanol.