Structure of Imidazolium-N-phthalolylglycinate Salt Hydrate: Combined Experimental and Quantum Chemical Calculations Studies

: An organic supramolecular salt hydrate (imidazolium:N-phthalolylglycinate:H 2 O; IM + NPG − -HYD ) has been examined for its charge-transfer (CT) characteristics. Accordingly, IM + – NPG − –HYD has been characterized thoroughly using various spectroscopic techniques. Combined experimental and quantum chemical studies, along with wave function analysis, were performed to study the non-covalent interactions and their role in CT in the supramolecular salt hydrate. Notably, IM + –NPG − –HYD crystalizes in two configurations ( A and B ), both of which are held together via non-covalent interactions to result in a three-dimensional CT supramolecular assembly. The through-space CT occurs from NPG – (donor) to IM + (acceptor), and this was mediated via noncovalent forces. We demonstrated the role of π – π stacking interactions (mixed-stacking donor-acceptor interactions) in the presence of charge-assisted hydrogen bonds in the regulation of CT properties in the self-assembly of the IM + –NPG − –HYD salt hydrate.


Introduction
The self-assembly characteristics of non-covalent interactions, such as hydrogen bonds, π-π stacking, and other weak van der Waals forces participating in (bio)chemical processes, have inspired supramolecular chemists to employ them effectively in crystal engineering [1,2]. Crystal engineering is defined as the study of intermolecular interactions in terms of crystal packing and exploiting such kinds of non-covalent interactions for designing multifunctional materials having the desired properties [3,4]. It has been proven that different intermolecular interactions and/or crystal packing could highly modulate the optical and electronic properties of molecular crystals [5]. For example, the formation of multi-component molecular materials (co-crystals, salts, and salt hydrates) has been extensively investigated with the development of crystal engineering or supramolecular chemistry [6][7][8].
Hence, a profound understanding of non-covalent interactions and the prediction of their strength are critical to designing the desired functional materials. Remarkably, the physical and chemical properties of co-crystals are not simply the sum of the molecular properties of their constituent molecular moieties [26,27]. In this regard, co-crystallization or crystallization of multicomponent molecular systems (co-crystals, salts, and salt hydrates) opens the door not only for the development of new multifunctional materials but also for the examination of novel phenomenon such as charge-transfer along the molecular crystallographic axis [28,29].
In 1952, Mulliken laid the foundation stone of the concept of charge transfer (CT) [30]. He observed new absorption bands for multicomponent molecular complexes; these bands were absent from the spectra of any of the constituent components individually. In addition, these molecular complexes contain an electron-donor component and an electron-acceptor component with low ionization energy and high electron affinity, respectively. Such type of multicomponent molecular complexes has an ionic excited state (ψ1 = D + ···A -) with comparatively low excitation energy, and mixing this state with the neutral ground state (ψ0 = D···A) results in ψ0 + λψ1 and ψ1 -λψ0 states. The transition dipole moment between these states is very large, resulting in strong electronic absorption; this transition dipole moment cannot occur in isolated components, be it an isolated donor or an isolated acceptor. Strikingly, the mixing of the molecular components results in color change because of the low excitation energy that may lie in the visible region of the spectrum. Mulliken observed that the mixing must lower the energy of the complex relative to the energy of the state ψ0: a chargetransfer contribution to the interaction energy. Hence, the CT interaction is the consequence of the charge density being transferred from the electron-donor component to the electron acceptor component, and the charge-transfer energy is the lowering in the energy associated with this transfer. Therefore, regarding the above-described Mulliken's CT concept and crystal-engineering concept, we combine both the concepts to achieve the desired charge-transfer complex.
In this study, we prepare a charge-transfer IM + -NPG − -HYD with mixed-stacking and chargeassisted hydrogen bonds between IM + and NPGmoieties. The proton transfer occurred from the carboxylic group of NPG to the amine group of IM, resulting in a charge-assisted hydrogen bond between both the moieties [31]. Generally, charge-assisted hydrogen bonding is stronger than normal hydrogen bond because of the presence of additional electrostatic interactions; consequently, the molecular components act as charge carriers. Nevertheless, charge-assisted hydrogen bonds have not been studied much as a designing factor in the self-assembly process of supramolecular entities [32][33][34][35][36][37]. There is still much scope for developing, understanding, and implementing charge-assisted hydrogen bonds in the context of crystal engineering.

Spectroscopic Analysis
A shift in the vibrational frequency of the multi-component molecular crystal with respect to the individual molecular entities denotes the formation of new hydrogen bonds and other weak noncovalent forces [38][39][40]. As depicted in Figure S1, from the IR spectrum of NPG, carbonyl stretching frequencies of the carboxylic acid are observed at 1724.17 and 1415.18 cm −1 . The carbonyl stretching bands were the most distinctive IR features to distinguish the carboxylic acid group from the conjugate carboxylate moiety [41,42]. The asymmetric and symmetric stretching frequencies of the carboxylate group (COO -) were observed at 1708.60 and 1419.98 cm −1 , respectively. Significant shifts in these vibrational frequencies were observed in IM + -NPG − -HYD (see Figure S1 in supplementary materials), signifying the proton transfer from the carboxylic-acid group [43]. A further relevant observation is that the IR spectrum of the IM + -NPG − -HYD salt displayed a characteristic band at the 3502.27 cm −1 , revealing the υOH stretching frequency of hydrogen-bonded lattice water [44]. Furthermore, the vibrational frequencies have also been calculated to support and assign the mode of vibration, observed experimentally. The IR spectrum has been simulated using DFT/ M062X/ 6-31G** method given in Figure S2. The calculated frequencies and other spectral features were found to be within the range, shown in Table 1. Three factors could be responsible for the deviation in the computed spectrum: 1) the environmental factor as DFT calculations were performed with solvation effect (liquid phase, EtOH as solvent) while experimental data was obtained at solid-state; 2) the calculated frequencies are included only harmonic while experimental have both harmonic and anharmonic effect; and 3) basis set and DFT functional discrepancies. However, the pattern and trend of spectra were quite similar in both the cases which support the experimental observations. To examine IM and NPG, as well as the self-assembly and molecular recognition in the solution phase, tandem mass-spectroscopy (ESI-MS/MS) was performed. A m/z 294.1 molecular ion peak, which is precisely matched with the sum of the IM + -NPG − -HYD salt's hydrate, [C13H13N3O5 + 3H + ], was designated as the primary ion in the first stage ( Figure S3).
The UV-vis spectra of IM, NPG, and IM + -NPG − -HYD salt hydrate were recorded in the solution state (EtOH, 1 × 10 −4 M) to explore the possibility of CT interaction between the NPGand IM + molecular entities within the IM + -NPG − -HYD salt hydrate ( Figure 1). The absorption spectrum of the IM + -NPG − -HYD salt hydrate is characterized by the appearance of a new absorption band in the visible region (380-550 nm). Neither NPG nor IM displayed any absorption band in this region. Therefore, this new absorption band in the visible region (380-550 nm) is attributed to the formation of donor-acceptor charge transfer IM + -NPG − -HYD salt hydrate [45][46][47][48]. In comparison to UV-vis spectra of IM, NPG, and IM + -NPG − -HYD salt hydrate, a remarkable bathochromic shift in the absorption bands appeared in the region of 250-350 nm is also observed in IM + -NPG − -HYD salt hydrate. The observed significant bathochromic shift in the 250-350 nm absorption band is also indicated the charge transfer nature of the IM + -NPG − -HYD salt hydrate [45][46][47][48]. Furthermore, TDDFT calculation has been carried out at the cam-B3LYP/cc-pVTZ level of theory to calculate the vertical excitation energies and oscillator strengths of the electronic transition. The TDDFT simulated UV-vis spectrum (vertical excitation energies and oscillator strengths) of IM + -NPG − -HYD salt hydrate is in good agreement with the experimental spectrum. The TDDFT results suggested that a broad CT band in the range of 400-550 nm, comprises CT0 → CT1 , CT0 → CT2, CT0 → CT3, and CT0 → CT4 electronic transitions located at 430, 446,513, and 534 nm, respectively ( Figure 1). The vertical excitation energies for electronic transitions CT0 → CT1, CT0 → CT2, CT0 → CT3, and CT0 → CT4 were found to be 1.72, 1.79, 2.06, and 2.14 eV, respectively. Moreover, the molecular orbitals involved in the electronic transitions CT0 → CT1, CT0 → CT2, CT0 → CT3, and CT0 → CT4 were found to be HOMO → LUMO, HOMO-2 → LUMO, HOMO → LUMO+1, and HOMO-1 → LUMO+1 excitations, respectively [48]. The detail description of the CT transfer nature of the IM + -NPG − -HYD salt hydrate is given in Section 2.5. The 1 H NMR spectra of NPG, IM, and IM + -NPG − -HYD were recorded in DMSO-d6 solution and depicted in Figure 2. The 1 H NMR spectrum of NPG exhibited signals at 4.31, 7.85-7.93, and 13.25 ppm, attributed to the protons of the -CH2, aromatic ring, and -COOH group, respectively. In the 1 H NMR spectrum of IM, signals exhibited at 7.01-7.64 ppm are attributed to the aromatic protons of the ring, while an upfield signal appeared at 12.07 ppm is attributed to the -NH proton. Compared with the free NPG and IM 1 H NMR spectra, a significant peak shift was observed in the spectrum of IM + -NPG − -HYD, indicating the existence of hydrogen bonds and other noncovalent interactions also in the solution phase. A broad signal with low intensity at 11.90 ppm in IM + -NPG − -HYD is attributed to the -NH proton. Additionally, 1 H NMR spectrum for salt hydrated simulated theoretically using DFT/B3LYP/6-311 + G(2d,p) and PBE0/Def2-TZVPP methods to assign the experimentally observed 1 H NMR signals. The simulated 1 H NMR spectrum at PBE0/Def2-TZVPP level of theory is found to be in good match with the experimental one with a slight deviation. This deviation could be attributed to the absence of a solvation effect in the calculated shielding values (gas-phase) whereas the experimental values recorded in the solution phase.

Crystal Structure of IM + -NPG ̶ -HYD (1:1:1)
The IM + -NPG − -HYD salt was block-shaped and crystallized in the monoclinic C2/c space group. The asymmetric unit consists of an NPG -, an IM + and a water molecule (see Figure 3), and the crystal packing is depicted in Figure S4. The details of bond distances, hydrogen bonds, structure refinement, and crystal data were presented in Tables 2, and 3. The C-O distances in the C10-O3 and C10-O4 of the carboxylic acid group in NPG are 1.243(3) and 1.264(3) Å, respectively, indicating the presence of a deprotonated form of NPG in the IM + -NPG − -HYD salt hydrate. There are two configurations, A and B, present in the supramolecular assembly of IM + -NPG − -HYD. Both these configurations are held together through the water molecule via hydrogen bonds. In configuration A, IM + and NPG − molecules are arranged alternatively and stacked on each other through π-π interactions (centroid-centroid separation: 3.760 Å) to form a one-dimensional donor-acceptor (···DADADA···) CT-chain ( Figure 4a). In configuration B, NPGis linked to IM + and H2O by hydrogen bonds, further selfassembling alternatively to each other to form a one-dimensional donor-acceptor CT-chain.   Furthermore, these one-dimensional CT-chains are perpendicular to each other and connected with a water molecule via N2-H2···O4, N3-H3A···O1W, O1W-H1WB···O3, and O1W-H1WA···O2 hydrogen bonds, resulting in 3D supramolecular CT assembly ( Figure 5). Furthermore, geometry optimizations of both the configurations have been carried out at the DFT level of theory using the M06-2X/6-31G** method and in good agreement with the crystal structure data. The calculated bond distances and the distances between IM + and NPG − are given in Table 1. Slight deviation observed in the calculated structure observed due to the fact that the DFT calculation performed in solvent phase EtOH as a solvent using the CPCM solvation method while experimental data collected from the solid-state material ( Figure S5). We have also calculated the interfragmentary interaction energy from the Hirshfeld charges and depicted in Figure 4.

Hirshfeld Surface Analysis
Hirshfeld surface analysis of charge transfer IM + -NPG − -HYD gives useful information about the role of different intermolecular interactions in determining the optical and electronic properties. The Hirshfeld surface of IM + and NPGseparately presented in Figure S5 and plotted over a dnorm range from −1.138 to 0.880 Å (IM + ) and −0.731 to 1.024 Å (NPD − ), shape index from −1.0 to 1.0 Å, and curvedness from −4.0 to 0.4 Å. All these three surfaces are transparently displayed to visualize IM + -NPG − -HYD in an analogous orientation, around which they were mapped. We have taken the asymmetric unit among the entire unit cell from the single-crystal X-ray data. The information present in the hydrogen-bonding tables is corroborated adequately with the red spots, where the large rounded depressions (deep red) visible on the dnorm surface are indicative of the hydrogenbonding interaction, and other visible spots are indicative of the H···H, C···C, C···N, C···H, O···O, O···N, and C···O interactions. The dominant O···H and C···H interactions can be viewed in the Hirshfeld surface map by the light red circular region in Figure S6a. These light red spots are the result of the N···H, C···O, C···N, and C···C interactions. Other visible spots on the surface indicate H···H and further weaker and longer interactions apart from the hydrogen bonds. The presence of π-π stacking interactions between the fused aromatic ring of NPGmoiety and the five-membered heteroaromatic ring of IM + moiety was reflected in the shape-index surface, as depicted in Figure S5b. As the inspection on the Hirshfeld surface show adjacent red and blue triangles on the shape-index surface, a characteristic feature is the presence of π-π stacking interactions. Moreover, a clear flat region in the direction of the bottom of both sides of the NPGand IM + moieties is evident on the curvedness surface, also depicting the presence of π-π stacking interactions between these two moieties, as depicted in Figure S5c.
Furthermore, we have also presented two-dimensional (2D) fingerprint plots for each observed interaction between the atom pairs in the salt hydrate. The fingerprint plots can be decomposed to highlight particular atoms of pairwise interactions ( Figure 6). This decomposition allows the parting of contributions from different types of interaction, which overlay the full fingerprint plot. Complementary regions are visible in the fingerprint plots, where one moiety acts as a donor (de > di) and the other as an acceptor (de < di). The H···O/O···H intermolecular interactions appear as distinct spikes in the 2D fingerprint plot (see Figure 6). The N···H intermolecular interaction appears as two small broad spikes on in the 2D fingerprint plot. The proportion of H···O/O···H interactions consists of 18.7% and 23.3% of the total Hirshfeld surface area. Such kind of dissimilar proportion for O···H interactions characterizes the oxygen of water interacting with the carbonyl oxygen (C=O) of the NPGmoiety and is accountable for the building of the supramolecular self-assembly. A similar type of unequal proportions was also observed for C···H, C···N, N···H, and C···O interactions, clearly indicating that these interactions participate in the growth of 3D supramolecular assembly of IM + -NPG − -HYD.

Quantum-Theory-of-Atoms-in-Molecules Analysis
The topology of calculated electron density at the DFT/ M062X/cc-pVTZ level of theory has been analyzed by employing the QTAIM approach. The quantum theory of atoms-in-molecules (QTAIM) has proved to be a great tool for investigating a variety of chemical interactions. According to the QTAIM, which was developed by Bader, a chemical bond is distinctively described by a line of maximum electron density linking two bonded atoms (bond path, BP) on which a bond critical point (BCP) (3, -1) is present [49,50]. Koch and Popelier [51] quantified that a hydrogen bond would be present in a molecular structure if the electron density is within the range of 0.002-0.040 atomic unit (a.u.), and the resultant Laplacian within the range of 0.024-0.139 a.u. The presence of hydrogen bonds was further supported by Rozas et al. According to Rozas et al. [52,53], these hydrogen-bond interactions may be categorized as follows: i) strong and covalent hydrogen bond, when (∇ 2 ρ) < 0 and H < 0; ii) medium and partially covalent hydrogen bond, when (∇ 2 ρ) > 0 and H < 0; and (iii) weak and typically electrostatic hydrogen bond, when (∇ 2 ρ) > 0 and H > 0.

CT in IM + -NPG − -HYD
The DFT calculations have been performed to investigate the electronic properties of the supramolecular IM + -NPG − -HYD salt hydrate. In configuration A, the total Hirshfeld charges (atomic dipole moment corrected, ADC) for NPG -, IM + , and HYD fragments were calculated to be -0.921934 (−0.841933), 0.939092 (0.939091), and -0.017215 (−0.017215), respectively. While in configuration B, the Hirshfeld charges for NPG -, IM + , and HYD fragments were calculated to be -0.675391, 0.804130, and -0.128862, respectively. The atom wise Hireshfeld and atomic dipole moment corrected (ADC) charges are given in Table S3. Additionally, the static dipole moments of 13.284517 and 7.289192 a.u., with the vector pointing from NPGto IM + fragment, in configurations A and B, respectively, confirm the CT from NPGto IM + ( Figure S8) [48]. The interfragmentary CTs were further confirmed by the analysis of the frontier molecular orbitals (FMOs), i.e., HOMO and LUMO, of IM + -NPG − -HYD (see Figure 9). Furthermore, as a result of the intermolecular hydrogen bonds and π-π stacking in the ground state of this donor-acceptor salt hydrate, the CT process occurs from HOMO to LUMO, because of which the electron cloud reorganizes over the donor-acceptor moieties, thereby forming new molecular orbitals (MOs) [56]. Consequently, the electronic transitions occur between these newly formed MOs from the ground state to the excited CT states. Both HOMO and LUMO are symmetrically and spatially localized on NPG -(donor) and IM + (acceptor), respectively. The effective spatial separation between HOMO and LUMO, coupled with a small ∆EHOMO-LUMO as compared to individual NPG and IM, facilitates the induction of anion-to-cation CT in IM + -NPG − -HYD, through space.
The time-dependent density functional theory (TD-DFT) based calculations were performed to determine the oscillator strengths (f) and the corresponding energies of the vertical transitions for IM + -NPG − -HYD, and the major MOs contributing toward the electronic transitions are presented in Table S3. Therefore, the 20 excited states were calculated for IM + -NPG − -HYD using a closed-shell procedure at CAM-B3LYP/cc-pVTZ level and, as discussed hereafter. The calculated absorption spectrum for IM + -NPG − -HYD is depicted in Figure 1. The calculated spectrum reproduces the major spectral features and is in good agreement with the experimental absorption spectrum. The strongest absorption band, which is centered at 300 nm in the experiment, is also reproduced in the TD-DFT calculation which is mainly the results of HOMO-2 to LUMO+4 transitions from the singlet ground state to the 10th singlet excited state (S0-S10). Other transitions having low oscillator strengths are also observed in the calculated spectrum in this absorption region, which is presented in Table S4. By analyzing the FMOs responsible for these transitions, it is indicated that the character of this band corresponds primarily to π-π* transition. The other broadband with weak absorption appeared in the higher wavelength region (350-550 nm); this broadband was absent from the free IM and NPG absorption spectra and was reproduced in the TD-DFT calculation, which mainly involves multiple superimposed absorptions from the singlet ground state (S0) to the singlet excited states (Sn, where n = 1-4) with similar oscillator strengths, comprising the donor (NPG -)-acceptor (IM + ) CT transitions. Interestingly, the FMOs involved in these electronic transitions were localized mostly on either the donor or acceptor moieties, resulting in an effective CT [57]. When an electronic excitation occurs, we supposed that the charge density transferred from the donor to acceptor. The excitation can be attributed mainly to the transitions of the MOs from HOMOs to LUMOs. Nevertheless, this model of single orbital pair is not appropriate in most practical cases, and this analysis considers the IM + -NPG − -HYD transition processes to comprise the contributions of multiple pairs of MOs, and chargedensity transfers are considered as electrons and holes in the resulting text. The distribution isosurfaces of electron and hole for IM + -NPG − -HYD is depicted in Figures S9 and S10. Spatially distributed electrons and a hole in a particular excitation correspond to CT excitation, while alternatively or merged distribution of the electron and hole represents local excitation [58]. Hence, in IM + -NPG ̶ -HYD, mostly electrons and holes are spatially distributed at IM + and NPGfragments, respectively, indicating that both the absorption bands at 300 and 350-550 nm possess a CT character. Additionally, we have also calculated the extent of interfragmentary CT and charge-density difference (CDD) from the donor (NPG -) to acceptor (IM + ) fragment, as depicted in Figure 10. Hence, these theoretical results are in good agreement with the experimental results, and it can be concluded that IM + -NPG − -HYD behaves as a CT salt hydrate.

Natural Bond Orbital Analysis
The NBO analysis was carried out to investigate the role of charge transfer in the charge-assisted hydrogen bond and π···π stacking interaction using the second-order perturbation theory of the Fock matrix within the framework of natural bond orbital (NBO) analysis. From the NBO analysis, the strength of the CT donor-acceptor interactions between NPGand IM + or H2O in IM + -NPG − -HYD can be obtained by probing the interactions between the filled (i, donor) Lewis-type NBOs and empty (j, acceptor) non-Lewis-type NBOs. In configuration A, both the NPGand IM + fragments behave as a donor and also as an acceptor. In addition, the charge transferred from the donor: NPGto the acceptor: IM + and vice versa are attributed to the weak non-classical O···C and N···C contacts, respectively. This CT mainly derives from the lone electron pairs of both oxygen and nitrogen in NPGand IM + , respectively, along with the empty antibonding orbitals of the carbon of the acceptor fragment. The energies of these interactions E(2) and donor-acceptor NBO orbital pairs are presented in Table 5 and Figures 11 and 12. Such behavior indicates that the charge facilitates the whole supramolecular assembly along the crystallographic molecular axis. While in configuration B, the CT resulting from the charge-assisted hydrogen bond is mainly attributed to the electron-donor (NPG -) and electron-acceptor (IM + ) interactions of the lone pair electrons of oxygen (from the carboxylate group) of NPGwith an empty antibonding orbitals of the hydrogen-bonding donor hydrogen atom of IM + LP(O) to BD* (D-H). Table 5. Results of NBO analysis of donor-acceptor interaction in IM + -NPG − -HYD, determined using the second-order perturbation theory

Conclusions
An organic supramolecular salt hydrate (imidazolium:N-phthalolylglycinate:H2O; IM + -NPG − -HYD) has been synthesized and characterized by single X-ray crystallography, 1 H NMR, IR, tandem mass-spectroscopy (ESI-MS/MS), elemental analysis, and UV-vis spectroscopy. The IM + -NPG − -HYD crystallizes in the monoclinic space group C2/c (Z = 8). The IM + -NPG − -HYD salt hydrate was formed via supramolecular self-assembly. The mixed π-π stacking donor-acceptor interactions and charge-assisted hydrogen-bonding interactions between IM + and NPGare the primary driving intermolecular interactions for the self-assembly, and they contribute to the CT characteristics of IM + -NPG − -HYD. The nature of the anion-to-cation CT was confirmed by comparing UV-vis spectra of IM and NPG, new absorption broadband ranging from 350 to 550 nm has appeared in IM + -NPG − -HYD, corresponds to the CT band , confirms the formation of donor-acceptor CT IM + -NPG − -HYD.
A series of DFT calculations have been performed to reproducing the experimental results. The DFT calculated structures and 1 H NMR spectrum are in good agreement with experimental results. The TDDFT calculated UV-vis spectrum reproduced major spectral features observed in the experimentally recorded spectrum. Also, electron-hole and Hirshfeld charge analysis validated the charge transfer (CT) nature of salt hydrate. Furthermore, the Hirshfeld surface, QTAIM, and NBO analysis explored the nature and strength of non-covalent interaction and their role in the charge transfer phenomenon.

General Methods and Instrumentation
All the starting materials were purchased from well-known chemical suppliers (Sigma-Aldrich St. Louis, Missouri MO, United States) and were used without further purification. The solvents used were of analytical or chromatographic grade.
The infrared spectra were obtained (KBr disk, 400-4000 cm −1 ) on a Perkin-Elmer Model 1320 spectrometer (Waltham, Missouri MO, United States). The electronic spectrum was recorded on PerkinElmer UV-vis spectrophotometer (Waltham, Missouri MO, United States).. The NMR spectra were obtained on a Bruker DRX−400 spectrometer (Billerica, Missouri MO, United States) with Me2SO−d6 as the solvent.

Synthesis of IM + -NPG − -HYD (1:1:1)
A total of 333 mg (2 mmol) of NPG and 104 mg (2.0 mmol) of IM were milled with 5-10 drops of aqueous methanolic solution for 30 min at room temperature, and salt formation was confirmed using the IR, NMR, and ESI techniques. Then, the compound (50 mg) was dissolved in an aqueous solution of methanol and left to evaporate slowly at room temperature. Good-quality crystals suitable for single-crystal X-ray crystallography analysis were obtained after 15 days.

Single-Crystal X-Ray Crystallography
A Bruker SMART APEX CCD diffractometer (Billerica, Missouri MO, United States) was used to collect the single-crystal X-ray data at 100 K by employing graphite monochromatic MoK radiation (0.71073). The linear absorption coefficient, scattering factors for the atoms, and anomalous dispersion correction were referred from the International Table for X-ray crystallography [59]. The data integration and reduction were performed using SAINT software [60]. Empirical absorption correction was applied to collected reflections by using SADABS [61], and the space group was determined using XPREP [62].
The structure was solved via direct methods using SHELXTL-2016 and refined on F via fullmatrix least-squares using the SIR -97 program package [63]. Only a few H atoms could be located in the difference Fourier maps of the structure. The positions of the remaining were calculated using idealized geometries (riding model) and assigned fixed isotropic displacement parameters. All the non-H atoms were refined anisotropically. The structure-refinement and crystal data are presented in Table 5. The selective bond distance and angles are presented in Tables S1 and S2, respectively.

Computational Details
The single-point calculation, full geometry optimization, and NBO analysis [64] were performed at DFT level of theory by using the M06-2X functional [65] (this functional was specifically developed to describe weak dispersion forces and noncovalent interactions) with the help of the Gaussian-09 program package [66]. The calculations were performed using cc-pVTZ basis sets [67] for all the atoms. All the geometry optimizations and TDDFT calculations were performed by employing the polarizable continuum model, CPCM (EtOH as solvent). No symmetry restrictions were applied during geometry optimization. Furthermore, the Hessian matrix was calculated analytically for the optimized structures to prove the location of the correct minima (no imaginary frequencies). The 1 H NMR spectrum of IM + -NPG − -HYD salt hydrated simulated using DFT/B3LYP/6-311+G(2d,p) and PBE0/Def2-TZVPP methods [68,69]. We have simulated the UV-vis spectrum (vertical excitation energies and oscillator strengths) of IM + -NPG − -HYD using the CAM-B3LYP/cc-pVTZ method. The CAM-B3LYP functional is working well to calculate CT excitations by adopting the Coulombattenuating method as long-range corrected modification [70]. The electron-hole distribution, charge density difference (CDD), and inter-fragment CT analysis between the singlet ground state (S0) and singlet excited states (Sn, where n = 1-20) have been carried out with Multiwfn program [58] on the time-dependent wavefunction solution of IM + -NPG − -HYD generated at the CAM-B3LYP/cc-pVTZ level of theory. In addition to the topological analysis of the electron-density distribution by using the QTAIM and ELF method, it has been performed employing the Multiwfn program. The Hirshfeld surface was mapped using Crystal Explorer [71] software using crystal structure coordinates of CIF files. The Cartesian atomic coordinates of the calculated equilibrium structures in ethanol are provided in the ESI-material.