One-Dimensional Organic–Inorganic Material (C6H9N2)2BiCl5: From Synthesis to Structural, Spectroscopic, and Electronic Characterizations

The new organic–inorganic compound (C6H9N2)2BiCl5 (I) has been grown by the solvent evaporation method. The one-dimensional (1D) structure of the allylimidazolium chlorobismuthate (I) has been determined by single crystal X-ray diffraction. It crystallizes in the centrosymmetric space group C2/c and consists of 1-allylimidazolium cations and (1D) chains of the anion BiCl52−, built up of corner-sharing [BiCl63−] octahedra which are interconnected by means of hydrogen bonding contacts N/C–H⋯Cl. The intermolecular interactions were quantified using Hirshfeld surface analysis and the enrichment ratio established that the most important role in the stability of the crystal structure was provided by hydrogen bonding and H···H interactions. The highest value of E was calculated for the contact N⋯C (6.87) followed by C⋯C (2.85) and Bi⋯Cl (2.43). These contacts were favored and made the main contribution to the crystal packing. The vibrational modes were identified and assigned by infrared and Raman spectroscopy. The optical band gap (Eg = 3.26 eV) was calculated from the diffuse reflectance spectrum and showed that we can consider the material as a semiconductor. The density functional theory (DFT) has been used to determine the calculated gap, which was about 3.73 eV, and to explain the electronic structure of the title compound, its optical properties, and the stability of the organic part by the calculation of HOMO and LUMO energy and the Fukui indices.


Introduction
Halometallates (III) organic-inorganic materials have generated significant interest because of their special structural features and high potential applications, such as electronic and optical materials [1][2][3][4][5][6]. Among these classes of materials, halogeno bismuthate(III) with general formula A a B b X 3b+a , where A = organic amine and X = I, Br, Cl, have emerged as an up-and-coming class of luminescent and semiconducting materials and have stimulated considerable interest among researchers [7][8][9][10][11][12]. The anion substructure of these materials comprises MX 6 octahedra that may be isolated (zero-dimensional-0D) or linked by edges, corners or faces (1D, 2D, 3D architectures) [13][14][15]. Low dimensional (1D and 0D) substructures expand physical properties associated with organic and inorganic constituents. The anionic parts act as semiconductors and the organic cations act as potential barriers, contributing to original electronic and optical behavior [16][17][18]. Thus, in order to increase the dimensionality of the inorganic substructures and decrease the band gap, an additional successful scheme is the incorporation of multifunctional organic cations capable of influencing the bonding structures within the inorganic network and increasing the orbital interaction across adjacent chains. The organic cations control the structure through their size, in addition to hydrogen bonds and/or van der Waals interactions. Generally, the trend is that non-bulky cations with a delocalized charge, such as imidazolium groups, can form hydrogen bonds with the terminal nitrogen atoms (N) and lead to anions with high N/M ratios (M: post-transition metals), as in the case of 1D M III X 5 2− . In the present work, we describe the synthesis and experimental characterization supported with theoretical studies of a material of formula (C 6 H 9 N 2 ) 2 BiCl 5 . The results of the vibrational (Raman-infrared), optical (absorption-band gap and HOMO-LUMO transitions) properties, and Hirshfeld surface analysis are discussed. The density functional theory DFT calculation was used to explain the optical properties by determining the calculated gap, the dielectric function, and the refractive index. The title compound stability is discussed in terms of the total and partial density of states and the HOMO/LUMO quantum mechanical descriptors of the organic part. The calculated Fukui indices give an insight into the nucleophilic and electrophilic sites of the organic component.

Experimental and Computational Methods
2.1. Synthesis of (C 6 H 9 N 2 ) 2 BiCl 5 Bismuth(III) oxide (Bi 2 O 3 ), hydrochloric acid (HCl; 37% in water), and 1-Allylimidazole (C 6 H 8 N 2 ) were purchased from Sigma-Aldrich and used without further purification. The synthesis of (C 6 H 9 N 2 ) 2 BiCl 5 material was carried out by slow evaporation method at room temperature. Values of 1.0 mmol (0.65 g) of Bi 2 O 3 and 2.0 mmol (0.30 g) of 1-Allylimidazole were dissolved together in 20 mL of an aqueous hydrochloric acid solution (pH ≈ 3). The mixture was stirred with heating (50 • C) for one hour. Then, it was covered with paraffin film and allowed to evaporate at room temperature. Colorless block crystals with suitable dimensions for crystallographic study were formed after 10 days.

Surface Investigation by SEM/EDX
SEM-EDX was used to investigate the morphology and the elemental composition of the crystals of the title compound. SEM images were acquired using a JEOL, JSM-6380-LA SEM-EDXS spectrometer, typically operated with an acceleration voltage of 15 kV coupled with an energy-dispersive X-ray spectrometry (EDX) detector system. In order to limit detrimental charge effects, a platinum-rich tape was used to partially cover the sample.

X-Ray Diffraction Analysis
Single crystal X-ray data were collected on a Bruker D8 Quest diffractometer (MoKα radiation λ = 0.71073 Å) at 298 K using Bruker APEX3 software package [19]. Data reduction was performed using SAINT [20]. Multi-scan absorption correction was performed using SADABS [21]. The structures were solved by direct methods with SHELXS using SHELXTL package and refined using full-matrix least-squares procedures on F2 via the program SHELXL-2014 [22]. All non-hydrogen atoms were refined with anisotropic displacement parameters. Hydrogen atoms were placed at calculated positions using a riding model. Molecular graphics were prepared using Diamond 3 [23]. Details of the data collection and crystallographic parameters are given in Table 1. Selected interatomic distances and bond angles are presented in Table 2. Additional information on the crystal structure study can be obtained in the form of a crystallographic information file (CIF), which was deposited in the Cambridge Crystallographic Data Center (CCDC) database (deposition number 1903636). For the other physical characterizations of the title compound, the single crystals have been selected according to their morphologies and color under binocular glass. Then, the single crystals were ground by means of an agate mortar to obtain a polycrystalline powder. The powder's purity has been verified by the powder X-ray diffraction (PXRD) followed by Rietveld refinement.
The PXRD pattern was recorded, at room temperature, using a Bruker D8 ADVANCE diffractometer equipped with a Cu anode (CuKα radiation λ = 1.54056 Å). The measurements were performed in the range 5-80 • under Bragg-Brentano geometry with a step of 0.02 • and a counting time of 2 s per step. The Rietveld refinement was performed by using the GSAS computer program [24,25]. The crystallographic data of the single crystal have been used as a starting set. The results of the Rietveld refinement are depicted in Figure 1. The final agreement factors were R p = 0.042, R wp = 0.126, and χ 2 = 0.051. The structural model was fully consistent with the one obtained from the single crystal X-ray diffraction data.  Figure 1. The Rietveld refinement of the room temperature powder diffraction pattern of (I).

Physical Measurements
The solid-state infrared spectrum of (C6H9N2)2BiCl5 was obtained using a Nicol NXR FTIR spectrometer, at room temperature, on KBr pellets ranging from 400 to 40 cm −1 . We recorded the Raman spectrum using a HORIBA Lab RAM HR Evolution Surfac Enhanced Raman Scattering Microscope in the 400-50 cm −1 region. Room-temperature d fuse absorption and diffuse reflectance spectra of polycrystalline powder (C6H9N2)2BiCl5 were measured using a JASCO V-770 spectrophotometer in the wav length range 300-700 nm, and BaSO4 plates were used as a reference.

Computational Methods
The structural model obtained from the X-ray single-crystal refinement has bee used as a starting model in the density functional theory (DFT) calculation. The electron structure and the optical property parameters were obtained using the electron exchang correlation functional proposed by Perdew-Burke-Ernzerhof (PBE) [26] and the project augmented-wave (PAW) pseudopotential plane-wave method, as implemented in th CASTEP code [27]. A plane-wave basis set was adopted with an energy cutoff of 600 e The atomic force, maximum displacement, and total energy convergence criteria were 0.

Physical Measurements
The solid-state infrared spectrum of (C 6 H 9 N 2 ) 2 BiCl 5 was obtained using a Nicolet NXR FTIR spectrometer, at room temperature, on KBr pellets ranging from 400 to 4000 cm −1 . We recorded the Raman spectrum using a HORIBA Lab RAM HR Evolution Surface-Enhanced Raman Scattering Microscope in the 400-50 cm −1 region. Roomtemperature diffuse absorption and diffuse reflectance spectra of polycrystalline powder of (C 6 H 9 N 2 ) 2 BiCl 5 were measured using a JASCO V-770 spectrophotometer in the wavelength range 300-700 nm, and BaSO 4 plates were used as a reference.

Computational Methods
The structural model obtained from the X-ray single-crystal refinement has been used as a starting model in the density functional theory (DFT) calculation. The electronic structure and the optical property parameters were obtained using the electron exchangecorrelation functional proposed by Perdew-Burke-Ernzerhof (PBE) [26] and the projector augmented-wave (PAW) pseudopotential plane-wave method, as implemented in the CASTEP code [27]. A plane-wave basis set was adopted with an energy cutoff of 600 eV. The atomic force, maximum displacement, and total energy convergence criteria were 0.05 eV Å −1 , 0.002 Å, and 10 −6 eV, respectively. The limited memory BFGS method has been used for energy minimization. We carried out the Brillouin-zone integrations using the Monkhorst-Pack scheme [28], with a regularly spaced mesh of 2 × 2 × 2 points in the reciprocal unit cell.
The HOMO, LUMO, and Fukui indices of the organic part of the title compound have been obtained using DMol3 code [29]. The calculations were carried out through Mulliken population analysis [30] and using the PBE functional method with DNP basis set [31]. The convergence parameters were as follows: maximum displacement of 0.005 Å, SCF tolerance 1 × 10 −6 eV/atom, convergence energy tolerance 1 × 10 −6 Ha. After the geometry optimization convergence, the following parameters were calculated: the HOMO (highest occupied molecular orbital energy), LUMO (lowest unoccupied molecular orbital energy) [32], and Fukui indices (FI) [33].
We drew the Hirshfeld Surfaces [34] and their relative 2D fingerprint plots [35] using Crystal Explorer 17 software [36] with a final refined crystallographic information file as the input. The quantifying and decoding of the intercontact in the molecular packing were created using d norm (normalized contact distance) and 2D fingerprint plots, respectively. The dark-red spots on the d norm surface resulted from the short interatomic contacts, while the other intermolecular interactions appeared as light-red spots. The d i (inside) and d e (outside) represents the distances to the Hirshfeld surface from the nuclei, with respect to the relative van der Waals radii. The proportional contribution of intercontact over the surface was visualized by the color gradient (blue to red) in the fingerprint plot. The enrichment ratios Exy were obtained from the actual contacts between the different chemical species (x, y) and equiprobable proportions calculated from the surface chemical content [37,38]. An enrichment ratio larger than unity indicated that the contact type was favored and made the largest contribution to the crystal packing.

Energy-Dispersive X-ray Analysis (EDX)
The SEM image and the characteristic EDX spectrum are shown in Figure 2. The SEM micrograph ( Figure S1) shows a good crystal quality, with a width of 4.5 mm. Figure 2 shows the total area-averaged EDX result obtained for the synthesized (C 6 H 9 N 2 ) 2 BiCl 5 . Characteristic X-ray peaks associated with C, Bi, O, N and Cl are all identified and labeled in Figure 2. Each element's weight percentage was determined to be 32.51% carbon, 32.98% bismuth, 1.37% oxygen, 8.42% nitrogen, and 24.72% chlorine. The small weight percent, accounted for by the small oxygen peak visible in the EDX spectrum, is attributed to the impurity acquired during the synthesis or the coating processes.

Crystal Structure Description
The single crystal X-ray diffraction analysis showed that the title compound (C 6 H 9 N 2 ) 2 BiCl 5 , (I), crystallizes in the centrosymmetric space group C2/c and the asymmetric unit comprises one Bi 3+ cation (located on a crystallographic two-fold axis), three chloride anions (one of which is located in a crystallographic inversion center), and one 1-allylimidazolium cation ( Figure 3). The cations are related by glide planes and an inversion center. The Bi 3+ ion is coordinated by six chloride anions in a distorted octahedral geometry.
micrograph ( Figure S1) shows a good crystal quality, with a width of 4.5 mm. Figure 2 shows the total area-averaged EDX result obtained for the synthesized (C6H9N2)2BiCl5. Characteristic X-ray peaks associated with C, Bi, O, N and Cl are all identified and labeled in Figure 2. Each element's weight percentage was determined to be 32.51% carbon, 32.98% bismuth, 1.37% oxygen, 8.42% nitrogen, and 24.72% chlorine. The small weight percent, accounted for by the small oxygen peak visible in the EDX spectrum, is attributed to the impurity acquired during the synthesis or the coating processes.

Crystal Structure Description
The single crystal X-ray diffraction analysis showed that the title compound (C6H9N2)2BiCl5, (I), crystallizes in the centrosymmetric space group C2/c and the asymmetric unit comprises one Bi 3+ cation (located on a crystallographic two-fold axis), three chloride anions (one of which is located in a crystallographic inversion center), and one 1allylimidazolium cation ( Figure 3). The cations are related by glide planes and an inversion center. The Bi 3+ ion is coordinated by six chloride anions in a distorted octahedral geometry.
The polymeric [BiCl 5 ] 2− anion comprises distorted BiCl 6 octahedra that form onedimensional zigzag chains by sharing the cis corners Cl3 and Cl3 i (generated by the two-fold symmetry operation; i: −x + 1, y, −z + 1/2). The chains propagate along the [001] direction ( Figure 4). The (C 6 H 9 N 2 ) + cations are inserted in the voids between the chains ( Figure 4). The Bi-Cl distances involving the terminal chlorides (Bi-Cl1 = 2.559 (6) Å and Bi-Cl2 = 2.698 (7) Å, are predictably shorter than those involving the bridging ones (i.e., Bi1-Cl3 = 2.968 (2) Å). These values are in the range characteristic of pentachlorobismuthates and are close to the values found for (C 2 H 7 N 4 O) 2   Anisotropic displacement parameters are shown at the 50% probability level. The grey dashed lines present bonds between atoms generated by symmetry elements (i) −x + 1, y, −z + 1/2. The intra-molecular hydrogen bond within the cation is shown in a green dashed line.  Examination of the geometrical features of the organic moiety shows that the 1allylimidazolium cation exhibits a non-regular configuration. The plane of the allyl group makes a dihedral angle of 84.63 • (16) with the plane of the imidazole ring and the twist of the allyl group (torsion N1-C3-C2=C1) is −2.1 • (5)). This configuration is due to intramolecular hydrogen-bonding interactions C1-H1A· · · N1. Such a result is not similar to the known data for allylimidazolium cations in homologous compounds [41,42]. Indeed, in Bi 4 I 16 .4(C 6 H 9 N 2 ). 2(H 2 O) [42], the allyl group has regular configuration where the dihedral angle between the allyl group and imidazole ring is 77.9 • and the torsion N1-C1-C4=C6 is 129.2 • (2).
Generally, the crystal structure cohesion is achieved via various types of non-covalent interactions. The principal feature of interest in this structure is hydrogen bonding ( Figure 5), which contributes to the stabilization of the crystal packing. As shown in Figure 5, the N-H and C-H moieties in the (C 6 H 9 N 2 ) + cation act as hydrogen-bond donors with chlorine vertices of [BiCl 5 ] 2− anions ( Table 3). The presence of π-π stacking between imidazolium rings and unusual π(C=C) ... π(C=C) interactions shows another kind of interaction [43]. The measurement of their forces is determined by the centroid distance between the neighboring imidazolium ring (d = 3.779 (2) Å) and the interacting distance measured from the middle of the aromatic C-C bond (x, y, z; (2) Å, which demonstrates the presence of weak π-π stacking (Figure 6a,b).   Table 3. Hydrogen-bond geometry (Å , °) for (I).

Hirshfeld Surface Analysis, Two-Dimensional Fingerprint Plots and Enrichment Ratios (E XY )
The Hirshfeld surfaces of (I) were mapped over d norm , curvedness, and shape index (Figure 7a-c).
type of contact corresponds to the π-π interactions between cationic molecules (the distance between the centroids of two antiparallel organic cations is 3.799 (2) Å and the distance between centroids of two parallel allyl groups in the cations is 3.287 (2) Å ). The [BiCl5] 2− octahedra also contributes to the molecular surface and displays enrichment values (EBi … Cl = 2.43). The Cl … Cl contacts are the third most frequent interactions (11.19%) ( Table 4) due to the abundance of chlorine on the molecular surface (33.45%). However, these contacts are highly diminished with an enrichment ratio around 0.07 (Table 4). Finally, the five types of contacts contribute significantly to the stability of the crystal structure. In conclusion, this analysis for intermolecular interactions is consistent with those observed by X-ray diffraction analysis.  Intense red regions in the surface mapped over d norm indicate that close contact interactions are apparent around the chlorine, nitrogen and carbon atoms participating in N-H···Cl and C-H···Cl hydrogen bonds (Table 1) [34]. The curvedness and shape-index provide further chemical understanding into the molecular arrangement. A surface with low curvedness designates a flat region, and it is indicative of π-π stacking interactions in the crystal. The donor and the acceptors of π-π stacking can be known as blue and red triangles around the participating atoms on the surfaces mapped over shape-index properties corresponding to H···H contacts, as shown in Figure 7b. The overall fingerprint plots are calculated, including all intermolecular contacts, and the decomposed fingerprint plots, which focus on specific interactions ( Figure 8). Relative contributions to the Hirshfeld surface for all intermolecular contacts in (I) are shown in Figure S2. Globally, the H· · · Cl contacts are the most favored interactions, and their relative contribution reaches 60.9% (Figure 8a). There are indeed three N2-H7···Cl2 iii , C4-H4···Cl1 iv and C5-H5···Cl3 v hydrogen bonds in the crystal structure (Table 3). These contacts are the most frequent interactions due to the abundance of chlorine and hydrogen on the molecular surface (% S Cl = 33.45% and % S H = 57.9%) ( Table 4). The fingerprint plot decomposition shows that H· · · H contacts comprise 23.1% of the total Hirshfeld surface area (Figure 8b) and are the second most frequent interactions due to the abundance of hydrogen on the molecular surface (57.9%). The C· · · H contacts represent the third most important interaction on the surface, with a percentage around 6% (Figure 8c) and an enrichment ratio higher than the unit E C · · · H = 1.13 ( Table 4). The N· · · C and C· · · C contacts are favored and display enrichment values (E N · · · C = 6.87 and E C · · · C = 2.85) ( Table 4), and they represent a major attraction in the crystal despite the small surface areas of carbon and nitrogen (4.6% and 1.7% of the total surface). This type of contact corresponds to the π-π interactions between cationic molecules (the distance between the centroids of two antiparallel organic cations is 3.799 (2) Å and the distance between centroids of two parallel allyl groups in the cations is 3.287 (2) Å). The [BiCl 5 ] 2− octahedra also contributes to the molecular surface and displays enrichment values (E Bi · · · Cl = 2.43). The Cl· · · Cl contacts are the third most frequent interactions (11.19%) ( Table 4) due to the abundance of chlorine on the molecular surface (33.45%). However, these contacts are highly diminished with an enrichment ratio around 0.07 (Table 4). Finally, the five types of contacts contribute significantly to the stability of the crystal structure. In conclusion, this analysis for intermolecular interactions is consistent with those observed by X-ray diffraction analysis.  Table 4. Hirshfeld contact surfaces (grey), derived "random contact" (blue) and "enrichment ratios" (green) for compound (I).

Vibrational Properties
All assignments are based on the spectra of the previous study of 1-allylimidazolium cation [41,[44][45][46]. The experimental FTIR and FT-Raman vibrational spectra are shown in Figures 9 and 10, respectively.
The Raman bands observed below 400 cm −1 for (I) (Figure 10) correspond to the anions' internal vibrational modes. The bands observed between 300 and 125 cm −1 on the Raman spectrum were assigned to the Bi-Cl stretching modes. The intense band observed at 275 cm −1 was assigned to the terminal Bi-Cl bonds' stretching mode, whereas the weak band located at 238 cm −1 was assigned to the stretching mode of the bridging Bi-Cl bonds. The lattice vibrations were found at wavenumbers lower than 100 cm −1 .  Table 4. Hirshfeld contact surfaces (grey), derived "random contact" (blue) and "enrichment ratios" (green) for compound (I).

Vibrational Properties
All assignments are based on the spectra of the previous study of 1-allylimidazolium cation [41,[44][45][46]. The experimental FTIR and FT-Raman vibrational spectra are shown in Figures 9 and 10, respectively.

Optical Absorption
Two absorption bands centered at 337 and 362 nm are observed in the room-temperature UV-Visible absorption spectrum ( Figure 11). These bands are typically assigned to metal-centered (MC) transitions [47][48][49] and to the ligand-to-metal charge transfer transitions [50,51]. We measured the UV-Visible diffuse reflectance spectra at room temperature to determine the optical band gap of (C6H9N2)2BiCl5. The absorbance as a function of reflectance is given by the Kubelka-Munk theory (Equation (1)): . Figure 9. Experimental infrared spectrum of (I).

Figure 10.
Single crystal Raman spectrum of (I).

Optical Absorption
Two absorption bands centered at 337 and 362 nm are observed in the room-temperature UV-Visible absorption spectrum ( Figure 11). These bands are typically assigned to metal-centered (MC) transitions [47][48][49] and to the ligand-to-metal charge transfer transitions [50,51]. We measured the UV-Visible diffuse reflectance spectra at room temperature to determine the optical band gap of (C6H9N2)2BiCl5. The absorbance as a function of reflectance is given by the Kubelka-Munk theory (Equation (1)): (1) Figure 10. Single crystal Raman spectrum of (I).
The Raman bands observed below 400 cm −1 for (I) (Figure 10) correspond to the anions' internal vibrational modes. The bands observed between 300 and 125 cm −1 on the Raman spectrum were assigned to the Bi-Cl stretching modes. The intense band observed at 275 cm −1 was assigned to the terminal Bi-Cl bonds' stretching mode, whereas the weak band located at 238 cm −1 was assigned to the stretching mode of the bridging Bi-Cl bonds. The lattice vibrations were found at wavenumbers lower than 100 cm −1 .

Optical Absorption
Two absorption bands centered at 337 and 362 nm are observed in the room-temperature UV-Visible absorption spectrum ( Figure 11). These bands are typically assigned to metal-centered (MC) transitions [47][48][49] and to the ligand-to-metal charge transfer transitions [50,51]. We measured the UV-Visible diffuse reflectance spectra at room temperature to determine the optical band gap of (C 6 H 9 N 2 ) 2 BiCl 5 . The absorbance as a function of reflectance is given by the Kubelka-Munk theory (Equation (1)): The energy band gap obtained by extrapolation of the linear portion of the absorption edges was estimated to be 3.26 eV, indicating a semiconductor nature of the title compound.

Density Functional Theory Calculations
To understand the title compound's optical properties, we have determined the electronic structure, the density of state, and optical properties such as the dielectric function and the refractive index. Stability of the title compound depends primarily on the organic part; therefore, we calculated the HOMO/LUMO energy gap and the Fukui indices of the imidazolium cation.

Band Structure and Density of States
The calculated band structure of (C6H9N2)2BiCl5 is shown in Figure 12. The maximum of the highest occupied valence bands and the minimum of the lowest unoccupied conduction bands lay at the same point (between A-and G-points) ( Figure 12); therefore, the title compound is a direct semiconductor. The band-gap value was about 3.73 eV ( Figure  12). This calculated value is higher than the experimental optic gap determined from the UV-Vis absorption spectrum ( Figure 11). This difference may be due to the self-interaction error of the exchange and correlation (XC) functional [52][53][54]. The energy band gap obtained by extrapolation of the linear portion of the absorption edges was estimated to be 3.26 eV, indicating a semiconductor nature of the title compound.

Density Functional Theory Calculations
To understand the title compound's optical properties, we have determined the electronic structure, the density of state, and optical properties such as the dielectric function and the refractive index. Stability of the title compound depends primarily on the organic part; therefore, we calculated the HOMO/LUMO energy gap and the Fukui indices of the imidazolium cation.

Band Structure and Density of States
The calculated band structure of (C 6 H 9 N 2 ) 2 BiCl 5 is shown in Figure 12. The maximum of the highest occupied valence bands and the minimum of the lowest unoccupied conduction bands lay at the same point (between A-and G-points) ( Figure 12); therefore, the title compound is a direct semiconductor. The band-gap value was about 3.73 eV ( Figure 12). This calculated value is higher than the experimental optic gap determined from the UV-Vis absorption spectrum ( Figure 11). This difference may be due to the self-interaction error of the exchange and correlation (XC) functional [52][53][54]. We illustrated the total and partial densities of states of (C6H9N2)2BiCl5 in Figure 13. The Fermi level (showed by a dotted line) is set to zero (EF). Here, we have treated H: 1s 1 ; C: 2s 2 2p 2 ; N: 2s 2 2p 3 ; Cl: 3s 2 3p 5 and Bi: 6s 2 4f 14 5d 10 6p 3 as valence electrons. Low-lying bands located within the energy range between −23 and −20 eV below the Fermi level mainly consist of Bi-5d and C-2s, and N-2s orbitals. The highest occupied crystal orbital is a mixture of Cl-3p and C-3p. With the majority from Cl-3p, conversely, the lowest unoccupied crystal orbital located above 3.73 eV is a mixture of Bi-6p and C-2p states with a smaller contribution of N-2p and Cl-3p states.  We illustrated the total and partial densities of states of (C 6 H 9 N 2 ) 2 BiCl 5 in Figure 13. The Fermi level (showed by a dotted line) is set to zero (EF). Here, we have treated H: 1s 1 ; C: 2s 2 2p 2 ; N: 2s 2 2p 3 ; Cl: 3s 2 3p 5 and Bi: 6s 2 4f 14 5d 10 6p 3 as valence electrons. Low-lying bands located within the energy range between −23 and −20 eV below the Fermi level mainly consist of Bi-5d and C-2s, and N-2s orbitals. The highest occupied crystal orbital is a mixture of Cl-3p and C-3p. With the majority from Cl-3p, conversely, the lowest unoccupied crystal orbital located above 3.73 eV is a mixture of Bi-6p and C-2p states with a smaller contribution of N-2p and Cl-3p states. We illustrated the total and partial densities of states of (C6H9N2)2BiCl5 in Figure 13. The Fermi level (showed by a dotted line) is set to zero (EF). Here, we have treated H: 1s 1 ; C: 2s 2 2p 2 ; N: 2s 2 2p 3 ; Cl: 3s 2 3p 5 and Bi: 6s 2 4f 14 5d 10 6p 3 as valence electrons. Low-lying bands located within the energy range between −23 and −20 eV below the Fermi level mainly consist of Bi-5d and C-2s, and N-2s orbitals. The highest occupied crystal orbital is a mixture of Cl-3p and C-3p. With the majority from Cl-3p, conversely, the lowest unoccupied crystal orbital located above 3.73 eV is a mixture of Bi-6p and C-2p states with a smaller contribution of N-2p and Cl-3p states. Figure 13. Density of states of (I). Figure 13. Density of states of (I).

Optical Properties
The title compound's optical properties can be described by the complex dielectric function ε(ω), which represents the linear response of the system to an external electro-magnetic field with a small wave vector. The optical properties are also associated with electron mobility and recombination rate of the electron [55]. We can express it as: The imaginary part of the complex dielectric function can be expressed as: where M is the dipole matrix element; i and j are the initial and final states, respectively; and f i is the Fermi distribution function for the ith state. E i is the energy of an electron in the ith states. By utilizing the Kramers Kronig relationship [56][57][58], the real part, ε1(ω), can be determined from the imaginary part as: where p indicates the principal value of the integral. The real part ε1(ω) of the complex function as a function of the used energy from 0 to 29 eV in three directions is shown in Figure 14a. The static functions ε1(0) of the real part is 1.87 for E||x and E||y and 1.94 for E||z. The two plots of the ε1(ω) along (100) and (010] directions are stackable, and they are different from the (001 direction, which means that the title compound is anisotropic. The maximum ε1(ω) peak is progressively increased by increasing the photon energy at around 4, 4, and 3.5 eV for E||x, E||y, and E||z, respectively. These peaks occur due to the gradual transition of mobile electrons from the topmost of the valance to the conduction band's visible bottom [59]. We illustrated imaginary part plots of the dielectric function ε2(ω) in Figure 14b. The variations are the same along [100] and [010] directions, and they are different from the [001] direction, which indicates that the title compound displays anisotropy. The maximum values were around 5.15, 5.15, and 5.04 eV for E||x, E||y and E||z, respectively. These peaks depicted in Figure 14(b) belong to an energy transition between some orbitals corresponding to certain energy because ε2(ω) is related to the density of state.
The optical constants such as refractive index η(ω) and the extinction coefficient k(ω), are calculated in terms of the real and the imaginary parts of the complex dielectric function as follows [60]: The variations of the refractive index η(ω) and extinction coefficient k(ω) are calculated according to the three crystallographic directions. The results are displayed in Figure 15.  The optical constants such as refractive index η(ω) and the extinction coefficient k(ω), are calculated in terms of the real and the imaginary parts of the complex dielectric function as follows [60]: The variations of the refractive index η(ω) and extinction coefficient k(ω) are calculated according to the three crystallographic directions. The results are displayed in Figure  15. The variation of η(ω) and k(ω) along the [100] and [010] directions are the same. Thus, for the following part, we discuss only [100] and [001] polarization directions.
The refractive index η(ω) (Figure 15a) represents the ratio of denser to the rare medium. The static refractive index values with considerable potentials are 1.37, and 1.39 for The refractive index η(ω) (Figure 15a) represents the ratio of denser to the rare medium. The static refractive index values with considerable potentials are 1.37, and 1.39 for E||x and E||z, respectively. The maximum peaks of refractive index η(ω) were obtained at 4.2 and 3.7 eV for E||x, and E||z, respectively. The magnitude of refractive spectra increased due to less band-gap and by increasing photon energy.
The extinction coefficients (k(ω)) ( Figure 15b) increased drastically to 1.15 and 1.03 at about 5.64 eV and 5.34 eV along the [100] and [001] directions and then decreased rapidly to the minimum value at about 8.35 eV, which is the indication of small absorption in this spectral region. It is worth pointing out that the extinction coefficient k(ω) is larger than the refractive index η(ω) in the spectral region from 5.64 eV until 8.35 eV, which means that light cannot propagate in this region.
The extinction coefficients (k(ω)) ( Figure 15b) increased drastically to 1.15 and 1.03 at about 5.64 eV and 5.34 eV along the [100] and [001] directions and then decreased rapidly to the minimum value at about 8.35 eV, which is the indication of small absorption in this spectral region. It is worth pointing out that the extinction coefficient k(ω) is larger than the refractive index η(ω) in the spectral region from 5.64 eV until 8.35 eV, which means that light cannot propagate in this region.

Frontier Molecular Orbital
The Frontier molecular orbitals provide the nature of reactivity and some of the molecules' physical and structural properties. Both HOMO and LUMO are the main orbitals considered in the study of the chemical reactivity of the molecule. The HOMO energy characterizes the electron-donating ability; however, the electron-accepting ability is characterized by LUMO energy. Obviously, the energy difference between HOMO and LUMO orbitals, called the energy gap, characterizes the molecule's chemical stability [61]. The

Frontier Molecular Orbital
The Frontier molecular orbitals provide the nature of reactivity and some of the molecules' physical and structural properties. Both HOMO and LUMO are the main orbitals considered in the study of the chemical reactivity of the molecule. The HOMO energy characterizes the electron-donating ability; however, the electron-accepting ability is characterized by LUMO energy. Obviously, the energy difference between HOMO and LUMO orbitals, called the energy gap, characterizes the molecule's chemical stability [61]. The energy gap is largely responsible for the chemical and spectroscopic properties of the molecules [62].
The optimized geometry of the organic part of the title compound obtained from the DFT calculations is shown in Figure 16a. The HOMO and LUMO orbitals are shown in Figure 16b

Fukui Function
The Fukui function (FI) is an important tool to determine the regioselectivity of the organic molecule and its reactive regions in terms of nucleophilic (f + ) and electrophilic attack (f − ) [66]. The Fukui indices as function of the atomic charges are given by: The low HOMO/LUMO energy gap of the organic cation is consistent with a high reactivity. At this level of calculation, the HOMO is mainly localized on the imidazolium ring, while the LUMO has a contribution from the allyl moiety only. The ionization potential (IP), electron affinity (EA), Mulliken electronegativity (χ) and absolute hardness (η) can be deduced from the values of E HOMO and E LUMO as follows. The ionization potential (IP) and electron affinity (EA) are related directly to E HOMO and E LUMO using Equations (7) and (8)  However, the Mulliken electronegativity (χ) and absolute hardness (η) can be approximated using Equations (9) and (10) [64,65]:

Fukui Function
The Fukui function (FI) is an important tool to determine the regioselectivity of the organic molecule and its reactive regions in terms of nucleophilic (f + ) and electrophilic attack (f − ) [66]. The Fukui indices as function of the atomic charges are given by: f + k = (N + 1) − (N) (for nucleophilic attack) (12) f 0 k = [q k (N + 1) − q k (N − 1)]/2 (for radical attack) (13) where q k is the electronic charge of atom k and N is the number of electrons. The condensed Fukui function values (f − k , f + k , f 0 k ) were calculated for electrophilic, nucleophilic, and radical attacks have been performed using Dmol 3 code. Table 5 shows that C(5), identified in Figure 7 by the red circle, has the highest values f + k , f − k and f 0 k . Thus, this carbon atom is the most favorite site for electrophilic, nucleophilic, and radical attack. The results obtained from the Fukui function and the analysis of the LUMO and HOMO orbitals (Figure 7b,c) are in good agreement, and the two methods lead to the same predictions of most electron-deficient site.

Conclusions
The (C 6 H 9 N 2 ) 2 BiCl 5 hybrid is obtained via slow evaporation at room temperature. This organic-inorganic material was crystallized in the monoclinic space group C2/c. Hydrogen bonds connected the organic and inorganic components. Scanning electronic microscopy (SEM) and energy-dispersive X-ray (EDX) were carried out. The Hirshfeld contact analysis implies that the enriched Cl· · · H, H· · · H, and C· · · H hydrogen bonds are the driving forces in the molecular arrangement and in the formation of the crystal packing. The optical absorption revealed that the band-gap of this compound is 3.26 eV. The characteristic vibrational peaks clearly appeared in FTIR and Raman spectra. The DFT calculation supports the experimental data and shows that the title compound cannot be used for photovoltaic cell because of its high optical gap (3.26 eV) and the low stability of the organic part (∆E gap = 0.2102 eV).

Supplementary Materials:
The following are available online at https://www.mdpi.com/1422-006 7/22/4/2030/s1. Figure S1: Micrograph SEM image of (I); Figure S2: Relative contributions to the Hirshfeld surface for all intermolecular contacts in (I); Table S1: Unit cell parameters of the single crystal and the optimized structure; Table S2: The atomic coordinates in the structure S1 and S2; Table  S3: The atomic displacements from the starting set to the optimized structure; Table S4: Evaluation of the structure similarity.
Author Contributions: H.F. designed the experiments, validation, computational study, formal analysis and wrote the paper; H.C., validation; M.F., single crystal analysis, structure resolution, and validation. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
All authors declare they have no conflict of interest.