Experimental and Theoretical Studies of the Vibrational and Electronic Properties of (2e)-2-[3-(1h-imidazol-1-yl)-1-phenyl- Propylidene]-n-phenylhydrazinecarboxamide: an Anticonvulsant Agent

In the current investigation, the molecular structure of the anticonvulsant agent


Introduction
Epilepsy as a disease can be described as a group of serious brain disorders affecting people of all ages.It is characterized by the onset of spontaneous convulsive and non-convulsive seizures.Epilepsy affects approximately 50 million people with a high prevalence in the first decade of life and in individuals more than 60 years old [1,2].Conventional antiepileptic agents (i.e., phenytoin and carbamazepine) and recent antiepileptic drugs (i.e., gabapentin, vigabatrin, remacemide, and loreclezole) are clinically effective for different seizure types.However, only 50% of sufferers are adequately treated with the currently available antiepileptic drugs (AEDs), and a few AEDs have severe side effects [3][4][5][6][7].Therefore, the development of safer and more effective AEDs is necessary.The rational design and discovery of new antiepileptics are difficult due to incomplete data concerning the molecular modeling pathway of human epilepsy and the sophisticated mechanism of action for the majority of antiepileptic drugs [8].Therefore, a ligand-based approach is one of the most important strategies for the development of new AEDs using the available biological information from well known and novel drugs.(2E)-2-[3-(1H-imidazol-1-yl)-1-phenylpropylidene]-N-phenylhydrazinecarboxamide ((2E)-HIPC) is a hybrid structure incorporating both imidazole and arylsemicarbazone moieties that act as pharmacophore scaffolds in some candidate anticonvulsants [9,10].This compound has exhibited anticonvulsant activity with 16% seizure protection at a dose level of 750 μmol/kg in a subcutaneous pentylenetetrazole (scPTZ) screen with no neurotoxicity [11].
The molecular structure of (2E)-HIPC has not previously been explored using computational methods.Therefore, we used ab initio Hartree-Fock (HF) and density functional theory (DFT/B3LYP) approaches to model the molecular structure of (2E)-HIPC.DFT/B3LYP/6-311++G(d,p) is a hybrid DFT standard computational method that gave sufficiently accurate results for organic molecules without consuming significant amounts of computing resources.Fourier transform infrared (FT-IR) and Raman spectra were also collected, and the data were compared with the computed spectra.The molecular electrostatic potential (MEP) and frontier molecular orbitals (HOMO-LUMO) of (2E)-HIPC were computed to predict the structural changes and reactive sites.In addition, the NMR ( 13 C and 1 H) chemical shifts of the title molecule were computed by the gauge-including atomic orbital method using B3LYP/6-311++G(d,p) calculations, and the theoretical results were compared with the experimental values.

General
The melting point was determined on a Gallenkamp melting point apparatus and is uncorrected.The FT-IR spectrum was obtained using KBr pellets on a DR/Jasco FT-IR 6300 spectrometer with a spectral resolution of 2 cm −1 (JASCO Inc., Easton, MD, USA).The FT-Raman spectrum was recorded on a Bruker RFS 100/s, Ettlingen, Germany.Excitation wavelength was 1064 nm with a maximal power of 150 mW and the emission of Nd:YAG laser (Goettingen, Germany) was used.The measurements were performed on solid samples.NMR spectra were measured in DMSO-d6 on a Bruker NMR spectrometer (Ettlingen, Germany) operating at 500 MHz for 1 H and 125.76 MHz for 13 C at the Research Center of the College of Pharmacy of King Saud University, Arabia Riyadh, Saudi Arabia.Chemical shifts are expressed in δ-values (ppm) relative to TMS as an internal standard.Mass spectra were measured on an Agilent Triple Quadrupole 6410 QQQ LC/MS with an ESI (electrospray ionization) source (Santa Clara, CA, USA).The X-ray diffraction measurements of the title compound were performed using a Bruker SMART APEXII CCD diffractometer (Ettlingen, Germany).Crystallographic data have been deposited with the Cambridge Crystallographic Data Centre (supplementary publication number CCDC-933941) (Cambridge, UK).

Quantum Chemical Calculations
The structure of (2E)-HIPC was first optimized using the AM1 method and then refined with optimization by DFT/B3LYP and ab initio HF methods at the 6-311++G(d,p) level as implemented in the Gaussian 09 software [15].The X-ray diffraction structure was used as the origin for the initial geometric optimization of the title compound.The vibrational frequencies at the HF and B3LYP levels of theory were evaluated for the optimized structures by assuming C1 point group symmetry.The B3LYP and HF methods provide good predictions of the structural parameters and binding energies of (2E)-HIPC.No negative frequencies were obtained during the optimization.To improve the agreement of the computed values with the experimental data, a spectral uniform scaling factor was used to balance the systematic errors caused by basis set incompleteness, neglect of electron correlation and vibrational anharmonicity.Therefore, the vibrational wavenumbers computed at the HF level were scaled by 0.905 [16], and the wavenumbers computed using B3LYP were scaled by 0.958 for wavenumbers greater than 1700 cm −1 and by 0.983 for wavenumbers less than 1700 cm −1 [17,18].With a few exceptions, the differences from the experimental values were less than 10 cm −1 after rescaling using these scaling factors.The Gauss viewing program [19] was used for visualization and verification of the assigned normal modes.
The variations in the thermodynamic functional parameters (i.e., heat capacity, entropy, and enthalpy) were examined from 100 to 700 K using the computed wavenumbers of the title molecule.The electronic absorption spectra for the optimized molecule were calculated with TD-DFT at the B3LYP/6-311++G(d,p) level.To analyze the reactive sites of the title compound, the molecular electrostatic potential were computed using the HF and B3LYP methods in the gas phase.DMSO and acetonitrile solvents were used in the DFT/B3LYP/IEFPCM calculations to improve the prediction of the NMR values and achieve good accuracy.The chemical shifts were calculated and analyzed.

Synthesis
The title compound 2 was prepared in a three-step reaction sequence as depicted in Scheme 1.
Scheme 1. Synthetic pathway for the target compound 2.

Potential Energy Surface Scan Grid
A potential energy surface scan was used to determine the optimized molecular structure of the (2E)-HIPC molecule in its equilibrium state and to determine the most stable structure at the global minimum of the potential energy surface.The energy surface scan grid for the title molecule was performed using the AM1 method, and the N16-C17-N18-C20 and C9-C8-C7-C5 dihedral angles were rotated for various geometries.The rotation angles of the two benzene rings with respect to the imidazole ring were computed by varying the dihedral angles.A potential energy scan of the (2E)-HIPC molecular geometry was performed by freezing the N16-C17-N18-C20 (SCAN1) and C9-C8-C7-C5 (SCAN2) dihedral angles and varying the dihedral angles values in 18 steps of 20° from −179° to 179° and from 0° to 180° for N16-C17-N18-C20 and C9-C8-C7-C5, respectively.A 3D view of the (2E)-HIPC rotational conformer is shown in Figure 1.This scan revealed that the optimized molecular geometries at 170 out of 361 with N16-C17-N18-C20 and C9-C8-C7-C5 dihedral angles of −19.42° and 440.713°, respectively, were the ground state conformations.The potential energy surface (PES) scan grid and the calculated optimized conformers were determined in the gas phase.The PES scan of the dihedral angle rotations resulted in 361 conformers.Two rotational global minimum energy conformers were identified for the title compound.Rotational conformers I and II, which are shown in Figure 2, were optimized with the AM1 and B3LYP methods in the gas phase, DMSO, and acetonitrile.The global minimum energies and dipole moments for the conformers are presented in Table S1.Conformer I exhibited the lowest global minimum energies of 0.1695246 and −1084.70550Hartree in DMSO solution using AM1 and B3LYP, respectively, compared with the obtained values in the gas phase and in acetonitrile.This result was confirmed by the computed dipole moment values in DMSO, which were 7.1931 and 7.4136 Debye using the AM1 and B3LYP methods, respectively.Therefore, conformer I was computed to have the minimum global energy, and this optimized molecular structure was used for further analysis.

Geometry of the Title Molecule
The optimized molecular structure of (2E)-HIPC with labeled atoms is shown in Figure 3. Selected optimized theoretical parameters in the gas phase and X-ray diffraction geometrical parameters of (2E)-HIPC are listed in Table S2.The optimized bond lengths of the C-C bond in the phenyl ring of (2E)-HIPC fall in the range of 1.390-1.409Å and 1.380-1.391Å in the B3LYP and HF methods, respectively.The observed XRD data for these bond lengths in (2E)-HIPC fall in the range from 1.368 to 1.393 Å, which are in good agreement with the values obtained using the HF method.The mean aromatic ring C-C bond distance is 1.385 Å by the HF method, which is consistent with the experimental value (1.382 Å).The optimized C-C bond length in the imidazole ring of (2E)-HIPC is 1.353 Å using the HF method, which is slightly longer than the observed XRD value (1.350 Å).The optimized bond length for the C=O bond in the title compound is 1.221 Å by the HF method, which is in a good agreement with the experimental value (1.227 Å).It is interesting that the calculated N15-N16 bond length in the side chain of the title molecule is almost the same as the observed XRD value (1.384 Å).Based on the above comparisons, although there are some differences between our values and the literature data, the optimized geometric parameters can reproduce the literature values well and are the basis for subsequent discussion.The theoretical optimized C-H bond lengths show large deviation from the experimental values, which may originate from the low scattering factors of hydrogen atoms in the XRD.The C-H bond lengths determined using the HF method lead to geometric parameters that are much closer to the experimental data.The ring C-N bond distances are in the range of 1.353 to 1.412 Å using both the HF and B3LYP methods, while the observed XRD values are in the range of 1.347 to 1.418 Å.The longest bond length was found in C20-N18 due to the close C=O group.The influence of different groups in the title compound on the bond length and dihedral angle parameters seems to be negligibly small.The mean angle deviation of the phenyl ring carbon atoms C-C-C was found to be greater by 0.2° in the HF method and by 0.1° in the B3LYP method than the corresponding XRD values.A large deviation was observed between theoretical and experimental XRD values in the bond angles of H38-N16-N15 and H35-C11-C12.The hydrogen bonding geometrical parameters [14] are presented in Table S3.

Frontier Molecular Orbitals (FMOs) and Ultraviolet Spectral Analysis
The excitation energies obtained using the TD-DFT/B3LYP method were computed for five excited electronic states.To enable a more relevant comparison between the theoretical and experimental data, the calculated transition energies were uniformly blue shifted such that the calculated and measured energies of the most intense transitions matched.The UV calculations were performed for a single molecule in the gas phase, DMSO, and acetonitrile without interactions with the environment, and the computed results were directly compared to the measured UV data for (2E)-HIPC in acetonitrile solution.The major energy contribution was approximately 0.98 in HOMO to LUMO for all three UV analyses.The 100-400 nm wavelength regions of the theoretical and experimental UV spectra of (2E)-HIPC are illustrated in Figure 4.The observed electronic absorption wavelengths at 250-320 nm were assigned to the π-to-π* transitions.The computed electronic absorption wavelengths were calculated to be 347.13,304.34, 283.91, 274.33, and 267.55 nm in the gas phase; 318.40, 286.17, 278.08, 262.57, and 252.13 nm in DMSO; and 319.90, 287.47, 279.77, 265.00, and 254.05 nm in acetonitrile.The corresponding dipole moment values were 6.6951, 8.4925, and 8.6364 Debye.Based on a comparison of the dipole moments of DMSO and acetonitrile, more electronic transitions occur in acetonitrile due to its value.It is worth mentioning that the experimental and calculated absorption wavelengths are in good agreement, as shown in Figure 4.
HOMO-LUMO or frontier molecular orbitals (FMOs) are an important measure of the nonlinear optical and electric properties, which can be obtained from quantum chemistry, UV-Vis, and NMR spectra [20].The HOMO has the ability to donate an electron, and the LUMO, which acts as an electron acceptor, has the ability to receive an electron.GaussSum 2.2 software [21] was used to calculate the minor and major energy contributions to the HOMOs and LUMOs.The energy gap between the FMOs determines the kinetic stability, chemical reactivity, optical polarizability, and chemical hardness/softness of a molecule.Chemical hardness and softness can be used as an additive tool to describe the thermodynamic aspects of chemical reactivity.To assess the trends in the transition energies of (2E)-HIPC, we carried out calculations for (2E)-HIPC in acetonitrile, DMSO, and the gas phase.The energies of the FMOs for (2E)-HIPC (i.e., HOMO, HOMO−1, LUMO and LUMO +1) were calculated at the B3LYP/6-311++G(d,p) level, and the results are listed in Table S4.The values of the HOMO-LUMO energy gap were 4.3772, 4.3994, and 4.0621 eV for (2E)-HIPC in acetonitrile, DMSO, and the gas phase, respectively.
The energy gap of the FMOs explains the ultimate charge transfer and interactions of atoms within the molecule.This activity influences the biological functions of the molecule.Three-dimensional plots of the HOMO and LUMO orbitals were computed at the B3LYP/6-311++G(d,p) level.Positive charge is shown in violet, and negative charge is shown in green.Based on the results shown in Figure 5, although the HOMO was localized on one benzene ring, the LUMO was localized on a different benzene ring as well as on the chain linking the two rings.The values of electronegativity, chemical hardness, softness, and electrophilicity index for (2E)-HIPC in acetonitrile were 3.9125, 2.1886, 1.0943, and 5.2418 eV, respectively.These values are between the corresponding values for (2E)-HIPC in DMSO and in the gas phase.The lowest value of the dipole moment was calculated for the gas phase.The stronger intermolecular interactions in acetonitrile are due to the dipole moment value being larger than those in DMSO and the gas phase.
The electronic transitions and UV spectra of (2E)-HIPC were computed using the TD-DFT/IEFPCM approach at the B3LYP/6-311++G(d,p) level.The calculations were performed for (2E)-HIPC in acetonitrile, DMSO, and the gas phase.The frontier orbital energies, absorption wavelengths (λ), oscillator strengths (f), and excitation energies (E) for (2E)-HIPC were computed in the gas phase, DMSO, and acetonitrile; the data are presented in Table S5.TD-DFT calculations for (2E)-HIPC predicted an intense electronic transition at 279.7 nm (f = 0.4339) in acetonitrile, 278.0 nm (f = 0.4700) in DMSO, and 284.9 nm (f = 0.5114) in the gas phase.The observed differences are due to shifting of the bands in both DMSO and acetonitrile due to their permittivities.The results are better represented in graphical form, as shown in Figure 4. Experimentally, absorption peaks of (2E)-HIPC in acetonitrile were observed at 283.1, 234.9, and 199.4 nm.The observed values are shifted downfield compared with the corresponding computed values.The major energy contributions were computed.

Electrostatic Potentials, Total Electron Densities, and Molecular Electrostatic Potentials
The molecular electrostatic potentials (MEPs) of the title compound were calculated using the B3LYP and HF methods in the gas phase, as shown in Figure 6.The MEP maps of (2E)-HIPC indicate that the regions of negative electrostatic potential lie outside of the compound and near the oxygen atom.For the imidazole ring of (2E)-HIPC, the negative electrostatic potential lies outside the ring.
Therefore, no favorable location exists for the formation of a normal cation above the ring.The molecular electrostatic potential (red color) was located on the N and O atoms belonging to the imidazole ring and C=O group, respectively.Less negative values were observed on the N-H groups (blue color).Both regions should be different because they are the principal acceptor and donor sites of the H bonds.In the vicinity of the nitrogen and oxygen atoms, the charge magnitude varies from −0.008 to −0.02.In the benzene ring, the charges vary from 0.008 to 0.8 from the inner to the outer contour.The electrostatic potential of different atoms on the surface is represented by different colors.
The value of the electrostatic potential increased from red (more negative) to blue (more positive).The coloring of these maps represents a range of −0.06023 a.u.(deepest red) to 0.06023 a.u.(deepest blue) at the B3LYP/6-311++G(d,p) level and −0.06321 a.u. to 0.06321 a.u. at the HF/6-311++G(d,p) level.The atomic sites colored in the MEP map indicate the strongest attraction in blue and the strongest repulsion in red.The electrostatic potential of C17 is more positive (0.718286) due to the high electronegativity of the neighboring O19 atom.N18 is more negative (−0.638291) due to the nearby oxygen and intramolecular bonding of H39 with N15.All of the hydrogen atoms have positive potentials.C7 and C8 have positive potentials, and C9 has a negative potential in the side chain carbon atoms.All carbons in the rings have negative electrostatic potentials except for C14 due to its bonding with two nitrogen atoms.

Infrared and Raman Spectral Analysis
In addition to three translational and three rotational degrees of freedom, the maximum number of potentially active normal modes of vibrations of a nonlinear molecule that contains N atoms is equal to (3N-6).All of the normal modes of vibration are both Raman and infrared active.The vibrational wavenumbers were calculated using the DFT/B3LYP method and the results, which were compared with the experimental values, are presented in Table S6.This comparison indicated that the calculated vibrational frequencies were overestimated due to the anharmonicity in the real system.The assigned vibrational modes of the experimental bands are based on normal mode analysis and comparison with the theoretically scaled wavenumbers obtained using the B3LYP method.The experimental and simulated FT-IR and FT-Raman spectra of (2E)-HIPC are shown in Figures 7 and 8, respectively.

C-H Vibrations
The aromatic rings exhibit C-H stretching as well as in-plane and out-of-plane bending vibrations.Aromatic compounds commonly exhibit multiple weak bands in the 3100-3000 cm −1 region [22] due to ring C-H stretching vibrations.These vibrations are not strongly affected by the nature of the substituents [23,24].All of the C-H stretching vibration modes exhibit very weak intensities.In the current study, (2E)-HIPC exhibits 10 C-H stretching vibrations.Among the four expected C-H stretching vibrations, three were observed at 3048, 2930, and 2908 cm −1 in the FT-IR spectra, and the same types of vibrations were observed at 3023, 2980, 2930, and 2880 cm −1 in the FT-Raman spectra.These vibrations are shifted downfield due to the imidazole moiety and the side-chain substitutions.All of the bands were observed as weak vibrations.
The vibrational modes due to C-H in-plane bending vibrations were observed in the 1000-1300 cm −1 region [25].For the title compound, the C-H in-plane bending vibrations were observed at 1198, 1180, 1165, and 1152 cm −1 .The wavenumbers obtained using the B3LYP/6-311++G(d,p) level of theory and subsequently rescaled were in good agreement with the experimentally observed data.The C-H out-of-plane bending vibrations appear in the 1000-700 cm −1 region [26].The vibrations identified at 740, 730, 713, and 700 cm −1 were attributed to C-H out-of-plane bending for (2E)-HIPC.

Skeletal Vibrations C-C Vibrations
Aromatic C=C stretching vibrations appear in the 1430-1625 cm −1 region, and the ring stretching vibrations for C-C are expected to appear in the 1300-1000 cm −1 region [27][28][29].The ring vibrations can vary due to different substitutions.In our investigation, the band at 1425 cm −1 in the FT-IR spectrum and the bands at 1579, 1540, 1516, 1495, and 1406 cm −1 in the Raman spectrum were assigned to C=C stretching vibrations.The band at 1301 cm −1 was assigned to a C-C vibration.The C-C-C dihedral in-plane and out-of-plane bending vibrations were observed at 996, 980, 965, 936, 910, and 761 cm −1 and 405, 396, 345, 322, 285, and 160 cm −1 in the expected region, respectively.The theoretical values exhibited excellent agreement with the experimental ones.As anticipated, the in-plane deformations were observed at higher frequencies than the out-of-plane vibrations.

Carbon-Oxygen Vibrations
In the C=O bond, the lone pair of electrons in the oxygen determine the nature of the carbonyl group.The C=O bond functionality and the carbonyl stretching vibration have been previously studied by infrared spectroscopy [30].This multiple bonded group is highly polar and therefore gives rise to an intense infrared absorption band.Carboxylic acids are identified by the C=O vibration, which gives rise to a strong absorption between 1700 and 1800 cm −1 .In general, the C=O stretching modes of acids are more intense than those for ketonic stretching modes.The characteristic ketonic C=O band in the current study appear at 1620 cm −1 in the FT-Raman spectrum, and this peak is very intense [31].The π-π bonding between the carbon and oxygen results in the formation of a double bond, and intermolecular hydrogen bonding reduces the wavenumbers of the C=O stretching absorption to a greater degree than the presence of the molecular H.This C=O vibrational mode of the title molecule appears in the expected range, which indicated that it was not strongly affected by other vibrational modes.The computed C=O frequency, which was calculated at the B3LYP/6-311G(d,p) level of theory, was 1622 cm −1 , which strongly supports the assignment of the experimental data.

N-H Vibrations
The N-H stretching modes typically occur in the 3500-3300 cm −1 region.The asymmetric -NH stretching mode appears in the range from 3500 to 3420 cm −1 , and the symmetric -NH stretching mode occurs in the 3420-3340 cm −1 region.In the current investigation, the title compound exhibited one N16-H38 asymmetric stretch at 3590 cm −1 and one N18-H39 symmetric stretch at 3433 cm −1 , which are in agreement with the B3LYP/6-311G++(d,p) calculated values.The NH2 in-plane deformation vibrations occur in the 1650-1580 cm −1 wavelength region of the spectrum.Therefore, a very strong band was observed in the FT-IR spectra of (2E)-HIPC at 1315 cm −1 , which was assigned to the in-plane deformation mode of the -NH group.The -NH out-plane bending was observed at 655 cm −1 in the FT-IR spectrum.

NMR Analysis
NMR parameters are used to determine the chemical environment of the molecule.The isotopic shifts provide essential information for understanding the microscopic environment of the molecules.DFT calculations were used to assign and analyze the experimental NMR spectra.DFT calculations were performed for (2E)-HIPC in different solvents to improve the accuracy of the predicted NMR values.The observed and simulated 1 H and 13 C NMR spectra of (2E)-HIPC are shown in Figure 9.The observed and theoretical values from the 1 H and 13 C NMR spectra of (2E)-HIPC are listed in Table S7.The chemical shifts were calculated for (2E)-HIPC in DMSO and acetonitrile, and the results were compared with the experimental data for (2E)-HIPC in DMSO.In the 1 H NMR spectrum, the signal at 145.1 ppm was due to the carbon atom (C17) attached to the oxygen atom.Similarly, the peaks located at 153.6 and 138.9 ppm were due to carbon atoms in the side chain, C7 and C14, respectively, located between the two nitrogen atoms in the imidazole ring.The chemical shifts for carbon atoms C20 and C5 in the benzene rings increased due to the attachment of the side chains.Of these two carbon atoms, C20 was upshifted due to its attachment to a nitrogen atom.In addition, the computed values were higher than the observed values.The chemical shifts of ring carbon atoms C2, C4, C6, C24, C22, C1, and C3 were observed at 128 ppm.The chemical shifts of the first two carbons (i.e., C2 and C4) were estimated to be smaller than the calculated values.The more upfield-shifted carbon atoms in the benzene ring were C21, C23, and C25.The observed chemical shifts of the imidazole ring carbon atoms (i.e., C11 and C12) were 119.4 and 126.4 ppm, respectively.
The 1 H NMR spectrum of (2E)-HIPC contains five sets of peaks between 3.33 and 10.27 ppm.The 1 H NMR signal splitting was caused by spin-spin coupling between nearest nuclei.The two triplets obtained at 3.33 and 4.13 ppm were due to the splitting of the signals of the side-chain hydrogen atoms and were assigned to the equivalent protons H31 and H32 and to the equivalent protons H33 and H34, respectively.The two triplets observed at 3.33 and 4.13 ppm were due to the splitting of the signals of the side-chain hydrogen atoms.The 1 H NMR splitting from 6.87 to 7.43 ppm was due to the signals from the hydrogen atoms (i.e., H35, H36, and H37) attached to the imidazole ring.The observed signal due to NH38 at 10.27 ppm exhibited a large chemical shift compared with the calculated values for (2E)-HIPC in DMSO and acetonitrile.The signal of the other NH39 exhibited an upfield chemical shift (8.88 ppm) due to intramolecular hydrogen bonding with N15.This intramolecular hydrogen bonding was confirmed by comparison to the calculated values for (2E)-HIPC in both DMSO and acetonitrile.The root mean square deviations (RMSDs) for 13 C NMR are 4.499 and 3.794 in DMSO and acetonitrile, respectively.However, a small deviation was observed for 1 H NMR in both solvents.The calculated chemical shift values for both the 1 H and 13 C NMR spectra of (2E)-HIPC in DMSO exhibited good linear fits with the experimental values.

Conclusions
A complete vibrational analysis of (2E)-HIPC was performed using the HF and DFT-B3LYP methods with the 6-311G++(d,p) basis sets.The stable rotational conformer of our molecule was identified via a potential energy surface grid scan using the AM1 and B3LYP methods.The global minimum energy and dipole moment were analyzed, and conformer I was used for further calculations.The standard statistical thermodynamic functions were obtained from the theoretical harmonic frequencies at the B3LYP and HF levels in the gas phase.The correlation equations between the statistical thermodynamic functions and temperature were also obtained.The thermodynamic functions fit well with regression values of 0.9999.The electronic transitions of the title compound were determined using TD-DFT calculations in acetonitrile, DMSO and the gas phase to simulate the electronic absorption spectrum of (2E)-HIPC.The electronic transitions are more intense in acetonitrile than in DMSO.In the acetonitrile solvent, the observed values were shifted downfield compared with the computed values.The observed and simulated spectra had good agreement when using the DFT-B3LYP/6-311G++(d,p) method.The RMSD value was lower in acetonitrile than in DMSO for 13 C NMR, while for 1 H NMR, there was not much deviation in either solvent.
Figure 1 illustrates the various levels of energies using colors where dark blue indicates the global minimum, light blue indicates a local minimum, red indicates a local maximum, and dark red indicates a global maximum.
The N-H•••N and C-H•••O bonds are identified as the NH and CO groups of (2E)-HIPC.These bonds are overestimated compared with the XRD values, whereas the H•••N and H•••O bonds are best fit with the HF and B3LYP values.Additionally, most of the intramolecular bond angles are in good agreement with the experimental values.

Figure 5 .
Figure 5. 3D plots of the HOMO and LUMO of (2E)-HIPC at the B3LYP level.