N-Arylation of Protected and Unprotected 5-Bromo-2-aminobenzimidazole as Organic Material: Non-Linear Optical (NLO) Properties and Structural Feature Determination through Computational Approach

The interest in the NLO response of organic compounds is growing rapidly, due to the ease of synthesis, availability, and low loss. Here, in this study, Cu(II)-catalyzed selective N-arylation of 2-aminobenzimidazoles derivatives were achieved in the presence of different bases Et3N/TMEDA, solvents DCM/MeOH/H2O, and various aryl boronic acids under open atmospheric conditions. Two different copper-catalyzed pathways were selected for N-arylation in the presence of active nucleophilic sites, providing a unique tool for the preparation of NLO materials, C-NH (aryl) derivatives of 2-aminobenzimidazoles with protection and without protection of NH2 group. In addition to NMR analysis, all synthesized derivatives (1a–1f and 2a–2f) of 5-bromo-2-aminobenzimidazole (1) were computed for their non-linear optical (NLO) properties and reactivity descriptor parameters. Frontier molecular orbital (FMO) analysis was performed to get information about the electronic properties and reactivity of synthesized compounds.


Introduction
There are 30 derivatives of 2-aminobenzimidazole registered in the world as drugs that exhibit diverse pharmacological activities, e.g., antiparasitic, antifungal, antiviral, and antiallergic [1]. The 2-aminobenzimidazoles are interesting compounds with many biological effects, such as immunotropic, diuretic, antihistamine, and highly selective characteristics of p38aMaP inhibition [2][3][4]. Moreover, 2-aminobenzimidazole derivatives with antiviral activity against herpes simplex virus (HSV), human cyclomegalo virus (HCMV), and HIV have also been developed and patented [5][6][7]. The benzimidazole moiety is also a component of chemosensor receptors, which are utilized to recognize anions selectively and play an essential part in several biological activities [8]. The polyfunctionality of the cyclic guanidine moiety in the 2-aminobenzimidazole has made it a building block for the synthesis of a wide range of pharmacologically important benzimidazole derivatives [9,10]. Therefore, the selective N-arylation of amino-substituted N-heterocycles without protection and with protection of nucleophilic sites is important synthetically, as it gives an easy approach to a variety of different bio-actively potent N-arylated compounds.
Optimizing a catalytic pathway for substrates, where a couple of hetero atom sites are available for cross-coupling reactions, is an actual challenge for the organic chemist. The selective N-arylation of heteroaromatics, which has more than one nucleophilic site, could be precious because it quickly develops molecular complexity in target molecules with the fewest synthetic influences. However, a large number of transition metal-catalyzed protocols have developed as incredible tools for numerous types of regio/chemoselective N-(hetero) arylation reactions in organic synthesis [11]. Different Pd-and Cu-catalyzed arylations, such as Goldberg [12], Ullman [13], and Buchwald-Hartwig [14] procedures, were used for N-arylations of azoles. More recently, the copper-catalyzed Chan-Lam coupling method for arylation of N-nucleophile by using copper acetate and aryl boronic acids was introduced as the standard because of its mild reaction conditions [15,16]. The Chan-Lam coupling has been used for N-arylation of various amines, anilines, esters, imidazoles, and nitrogen-containing heterocycles [17].
These days computational chemistry [18] has emerged as a well-recognized partner of experimental chemistry. Computational chemistry is an extremely vast topic but herein we limited ourselves to density functional theory (DFT) [19]. Over the last 40 years, density functional theory (DFT) has become the most dominant and powerful tool in computational quantum chemistry for the modeling and simulation of chemical systems. Materials demonstrating great non-linear optical response play a key role in telecommunication, optical information processing, optical computing, etc. [20,21]. Nadeem and co-workers studied benzimidazole derivatives to tune the second-order nonlinear optical molecular switching by proton abstraction. Their results illustrate that substituted compounds have robustly large off-on NLO switching with a difference in β o values of 7, 63, 85, and 75 times larger than their neutral counterparts, respectively [22]. Tayade and Sekar examined Benzimidazole-Thiazole based NLOphoric Styryl Dyes both experimentally and theoretically. The nonlinear response based on α, and β was increased after the substitution of acceptor groups [23]. Thakare et al. experimentally and theoretically analyzed the NLO response of BODIPY-benzimidazole conjugate. The results obtained from the DFT method are in good accordance with those produced from solvatochromic method. The values obtained for nonlinear absorption coefficient (β) and third-order susceptibility χ(3) are 7.45 × 10 12 and 3.85 × 10 13 , respectively [24]. These reports shed light on the importance of substituted benzimidazole synthesis and their use in the field of optical and nonlinear optical materials using the DFT method.
In the present study, 2-aminobenzimidazole and its derivatives were synthesized via Chan-Lam cross-coupling reaction and analyzed theoretically for structural, spectroscopic [25], and NLO properties [26]. Density functional theory (DFT) calculations were used to investigate quantum chemical parameters, such as electron affinity (E A ), ionization potential (I), electronic chemical potential (µ), electrophilicity index (ω), and chemical hardness (η) [27] of synthesized derivatives of 5-bromo-2-aminobenzimidazole. We hope that electron-donating and withdrawing groups substitution can enhance the NLO response of the newly designed bezimidazoles, and we will obtain large first-order hyperpolarizability for these compounds.

Chemistry
For the acquisition of optimized protocol for selective C-NH arylation of 5-bromo-2aminobenzimidazole, we initiated our examination by adopting two different protocols in which we used 5-bromo-2-aminobenzimidazole (1) and its derivative N-(5-bromo-1Hbenzo[d]imidazol-2-yl)acetamide (2) with aryl boronic acids in presence of Cu(OAc) 2 as a catalyst, Triethylamine (Et 3 N) or Tetramethylethylenediamine (TMEDA) as a base, and the solvent DCM and methanol. The reaction was carried out at normal temperature under an open-air environment (Scheme 1).

N-(5-bromo-1H-benzo[d]imidazol-2-yl)acetamide
Generally, N-arylation of the heterocyclic compound with aryl boronic acid takes place under anhydrous conditions. In one protocol of direct N-arylation of 5-bromo-2aminobenzimidazole (1) to produce its derivatives 1a-1f (Figure 1), we used a TMEDA base and Cu(OAc)2·H2O catalyst with mixed parotic solvent (CH3OH/H2O 8:1). The desired products were obtained in good yield within only 2 h [29]. In contrast, in the second protocol for selective N-arylation of N- (5-bromo-1H-benzo[d]imidazol-2-yl)acetamide (2) to produce its derivatives 2a-2f (Figure 2), we used anhydrous Cu(OAc)2 instead of Cu(OAc)2·H2O with aryl boronic acid in presence of dry DCM instead of CH3OH/H2O, Et3N, and 4 Å molecular sieves [15,30]. This method was not found efficient, as we failed to isolate the desired product in good yield, even after 24 h, and various tedious steps, such as protection of 5-bromo-2-aminobenzimidazole (1) involved in this method. It was noticed that, by the addition of TMEDA base and a small amount of water, the yield of desired products is enhanced significantly. Generally, N-arylation of the heterocyclic compound with aryl boronic acid takes place under anhydrous conditions. In one protocol of direct N-arylation of 5-bromo-2aminobenzimidazole (1) to produce its derivatives 1a-1f (Figure 1), we used a TMEDA base and Cu(OAc) 2 ·H 2 O catalyst with mixed parotic solvent (CH 3 OH/H 2 O 8:1). The desired products were obtained in good yield within only 2 h [29]. In contrast, in the second protocol for selective N-arylation of N-(5-bromo-1H-benzo[d]imidazol-2-yl)acetamide (2) to produce its derivatives 2a-2f (Figure 2), we used anhydrous Cu(OAc) 2 instead of Cu(OAc) 2 ·H 2 O with aryl boronic acid in presence of dry DCM instead of CH 3 OH/H 2 O, Et 3 N, and 4 Å molecular sieves [15,30]. This method was not found efficient, as we failed to isolate the desired product in good yield, even after 24 h, and various tedious steps, such as protection of 5-bromo-2-aminobenzimidazole (1) involved in this method. It was noticed that, by the addition of TMEDA base and a small amount of water, the yield of desired products is enhanced significantly.    In the present study, it was observed that the base can significantly affect the Chan-Lam coupling. Et 3 N and pyridine are commonly used bases in Chan-Lam cross-coupling reactions. These bases contain amine additives which donate the electron pair to Cu(II) and allow it to oxidize Cu(III) intermediate and also capture the acetic acid which produces during the transmetallation step. For selective N-arylation of 5-bromo-2-aminobenzimidazole (1) without protection, we used TMEDA and Et 3 N as bases and obtained products in good yield with TMEDA base. It is observed that TMEDA act as a bidentate ligand and forms a more stable complex with copper which enhances its reactivity toward the NH group of imidazole ring as compared to Et 3 N for its selective N-arylation [31,32].

Density Functional Theory (DFT) Studies
All synthesized compounds (1a-1f and 2a-2f) were computationally studied to investigate the structural-properties relationships. First of all, compounds were optimized at B3LYP functional of DFT by using GAUSSIAN 09 software [33]. For molecular orbitals description, Pople'6-31+G(d,p) [34] was used ( Figure 3A,B). Frequency analysis were completed at the same level of theory for further confirmation of these structures as true minima energy structures on potential energy surfaces. A frontier molecular orbitals (FMOs) analysis and reactivity descriptor parameters were calculated at B3LYP/6-31+G(d,p) method [35]. Nonlinear optical (NLO) properties along with polarizabilities and first hyperpolarizability parameters were evaluated using CAM-B3LYP [36], LC-BLYP [37], and ωB97XD [38] density functionals with 6-31+G(d,p) basis set.
reactions. These bases contain amine additives which donate the electron pair to Cu(II) and allow it to oxidize Cu(III) intermediate and also capture the acetic acid which produces during the transmetallation step. For selective N-arylation of 5-bromo-2-amino benzimidazole (1) without protection, we used TMEDA and Et3N as bases and obtained products in good yield with TMEDA base. It is observed that TMEDA act as a bidentate ligand and forms a more stable complex with copper which enhances its reactivity toward the NH group of imidazole ring as compared to Et3N for its selective N-arylation [31,32].

Density Functional Theory (DFT) Studies
All synthesized compounds (1a-1f and 2a-2f) were computationally studied to investigate the structural-properties relationships. First of all, compounds were optimized at B3LYP functional of DFT by using GAUSSIAN 09 software [33]. For molecular orbitals description, Pople'6-31+G(d,p) [34] was used ( Figure 3A,B). Frequency analysis were completed at the same level of theory for further confirmation of these structures as true minima energy structures on potential energy surfaces. A frontier molecular orbitals (FMOs) analysis and reactivity descriptor parameters were calculated at B3LYP/6-31+G(d,p) method [35]. Nonlinear optical (NLO) properties along with polarizabilities and first hyperpolarizability parameters were evaluated using CAM-B3LYP [36], LC-BLYP [37], and ωB97XD [38] density functionals with 6-31+G(d,p) basis set.

Frontier Molecular Orbital (FMO) Analysis
FMOs analysis is performed to get information about the electronic properties and reactivity of compounds. The energies of HOMOs (EHOMOs) and LUMOs (ELUMOs) explain electronic properties and HOMO-LUMO gaps (GH-L) validate kinetic stability and reactivity of compounds. Large GH-L declared the less reactivity and more stability of compound and vice versa [39]. The energies of HOMOs (EHOMOs), energies of LUMOs, (ELUMOs), and HOMO-LUMO energy gaps (GH-L) of all compounds are summarized in Table 3.

Frontier Molecular Orbital (FMO) Analysis
FMOs analysis is performed to get information about the electronic properties and reactivity of compounds. The energies of HOMOs (E HOMOs ) and LUMOs (E LUMOs ) explain electronic properties and HOMO-LUMO gaps (G H-L ) validate kinetic stability and reactivity of compounds. Large G H-L declared the less reactivity and more stability of compound and vice versa [39]. The energies of HOMOs (E HOMOs ), energies of LUMOs, (E LUMOs ), and HOMO-LUMO energy gaps (G H-L ) of all compounds are summarized in Table 3. Among the compounds 1a-1f, 1d has the smallest G H-L (4.25 eV), whose E HOMO and E LUMO are −6.02 and −1.77 eV, respectively. Thus, 1d is considered as kinetically less stable and moderately reactive among 1a-1f series. Previously reported gaps for substituted benzimidazole-, thiazole-, and benzothiazole-based compounds range from 3 to 0.63 eV. They are reported as highly stable compounds [23,40,41]. On the other side, 4.69 eV of G H-L is calculated for compounds 1a, 1b and 1e, which represent more stability and less reactivity of these compounds. The other compounds have G H-L values of 4.25 to 4.52 eV, which represent their adequate stability with moderate reactivity. In the other series (2a-2f), 2d has the lower G H-L (4.79 eV), similar to the results of compound 1d. The rest of the compounds have G H-L values from 4.91 to 5.08 eV and represent high stability and less reactivity of these compounds. All compounds of this series have a larger G H-L and are ultimately less reactive compared to compounds 1a-1f (Table 3).
The FMOs isodensity distributions of compounds 1a-1f are shown in Figure 4A. The isodensity in HOMOs is distributed on benzimidazole rings and bromine atoms attached to these rings in compounds 1a and 1e. Bromine is a strong electron-withdrawing group so electronic densities move toward it. However, some of the isodensity also reside on oxygen, chlorine, nitrogen, and sulfur atoms of compounds 1b-d and 1f, respectively, because these atoms also have more electronegativity like bromine atoms. The result is the shift of electronic density from benzimidazole rings towards heteroatoms (C, O, N, S, and Cl) attached as substituents. This specific distribution explains the lower HOMO-LUMO gaps of these compounds (1b-d and 1f) compared to other compounds (1a and 1e). The isodensity shift is prominent for 1d compared to other compounds. The isodensity in LUMOs is mainly distributed on phenyl rings in all compounds (1a-1f, except 1d). These results clarify that the electron densities are migrating mainly from the donor (electron-donating groups terminal) to the acceptor unit (electron-withdrawing groups terminal). Compounds (2a-2f) also have similar isodensity distributions and their HOMO and LUMO densities are given in Figure 4B. Isodenisties are localized on the benzimidazole rings in compounds 2a-d and 2e. LUMO densities are present on the aromatic rings. In compounds 2c and 2f, HOMO densities shifted towards the chloro-substituted aromatic and thiophene rings and LUMO densities are on the benzimidazole rings. The gaps are mainly due to the smaller values of the LUMOs of 1d and 2d, and since the LUMOs are mainly centered on the N-arylated part, the gap is smaller due to the accepting nature of the pyridine. As shown by its highest electrophilicity value, the ω values of 1d and 2d are 3.56 and 3.80 eV, respectively. rings in compounds 2a-d and 2e. LUMO densities are present on the aromatic rings. In compounds 2c and 2f, HOMO densities shifted towards the chloro-substituted aromatic and thiophene rings and LUMO densities are on the benzimidazole rings. The gaps are mainly due to the smaller values of the LUMOs of 1d and 2d, and since the LUMOs are mainly centered on the N-arylated part, the gap is smaller due to the accepting nature of the pyridine. As shown by its highest electrophilicity value, the ω values of 1d and 2d are 3.56 and 3.80 eV, respectively.

Reactivity Descriptor Parameters
For further interpretation of the reactivity of compounds (1a-1f and 2a-2f), some other reactivity descriptor parameters were also analyzed and the results of their analysis are shown in Table 4. These parameters involve electron affinity (E A ), ionization potential (I), electronic chemical potential (µ), chemical hardness (η), and electrophilicity index (ω). According to Koopman's theorem, the negative values of HOMOs and LUMOs correspond to ionization potential (I) and electron affinity (E A ), respectively [42][43][44].
The value of chemical hardness (η) is mathematically calculated as follows: Chemical hardness (η) = (E HOMO − E LUMO )/2 The chemical hardness of 1a-2f compounds is from 2.13 to 2.35 eV. The highly unstable and most reactive compound is 1d based on its small chemical hardness value (η = 2.13 eV). Compounds 1a, 1b, and 1e are stable compounds which are also confirmed from their η values of 2.34, 2.33, and 2.35 eV, respectively. The chemical hardness of other compounds is between 2.24 and 2.26 eV representing their moderate stability and reactivity. Ionization energy is between 5.78 and 6.02 eV and electron affinity are in the range of 1.10-1.77 eV. Looking towards compounds 2a-2f, we see that 2d has lower chemical hardness (2.39 eV), and others are noticed as hard compounds; their η values range from 2.45 to 2.54 eV. The ionization potential and electron affinity values of 2a-2f are more than those of 1a-2f (Table 4). Table 4. Ionization potential (I), electron affinity (E A ), chemical hardness (η), electronic chemical potential (µ), and electrophilicity index (ω) of compounds 1a-1f and 2a-2f.
Compounds Electronic chemical potential (µ) tells about the charge transfer at their ground state inside the compounds, and it is mathematically represented as follows: Electronic chemical potential (µ) = (E HOMO + E LUMO )/2 where the µ of compounds 1a-1f is from −3.44 to −3.89 eV. The highest µ value (−3.89 eV) and lowest µ value (−3.44 eV) were calculated for 1d and 1e, respectively. In other series (1a-1f), compound 2d also has a high value of µ (−4.26 eV), which signifies that more charge transfer is from an electron donor to acceptor group.
Based on energy, the electrophilicity index (ω) represents the stability of compounds when an extra charge is transferred from the surrounding [44], and it is mathematically calculated as follows: Electrophilicity index (ω) = µ 2 /2η For compounds 1a-1f, the ω values range from 3.61 to 4.02 eV. The compound 1d has the highest ω value of 3.56 eV which indicates that it is less stable for the incoming charge. The presence of electron donor and acceptor groups linked through extended conjugation is the reason for less stability of this compound and enhance its reactivity as incoming charge causes more delocalization of electronic density. In the series of 2a-2f, the ω value of 2d is 3.80 eV and represents the more charge accepting capability of these compounds. Other compounds have a low ω value, i.e., up to 2.63 eV.

Non-Linear Optical (NLO) Properties
From the start of this century, scientists are following the methods to produce nonlinear optical (NLO) materials because of their potential consumptions in optical data storage, optical communication, optical computing, optical limiting, medical imaging, laser devices, etc. [45,46]. To increase the NLO response of the compounds, different strategies are adopted, such as push-pull mechanism [47], metal-organic framework [48], excess electron system [45,49], etc. The organic molecules develop a strong NLO response as the transformation of electrons from donor to acceptor group is responsible for the improved value of β o [50]. Polarizability (α o ) and first hyperpolarizability (β o ) parameters are used to measure the NLO response of respective compounds.
Polarizability (α o ) and hyperpolarizability (β o ) parameters of optical and nonlinear optical properties are studied at different density functionals, i.e., CAM-B3LYP, LC-BLYP, and ωB97XD (vide infra). CAM-B3LYP has a 0.65 fraction of nonlocal exchange at an asymptotic distance and is a highly reliable method for calculating hyperpolarizabilities [36]. Full range-separated functionals LC-BLYP [37] and ωB97XD [38] have the correct 1.00 fraction of nonlocal exchange they are also well-known funtionals in the field of NLO. Almost comparable results are obtained at all selected density functionals where polarizability ranges from 182 to 239 au and hyperpolarizability is between 1.73 × 10 −30 and 5.63 × 10 −30 esu for 1a-1f and 2a-2f at LC-BLYP. Polarizability values for all them (1a-1f and 2a-2f) is from 188 to 246 au and hyperpolarizability is between 5.31 × 10 −30 and 1.90 × 10 −30 esu at ωB97XD. Moreover, 188 to 247 au and 2.00 × 10 −30 to 5.66 × 10 −30 esu are the polarizability and the hyperpolarizability ranges for 1a-1f and 2a-2f at CAM-B3LYP (Table 3 and  Supplementary Materials Table S1). The highest polarizability (247 au) and hyperpolarizability (5.66 × 10 −30 esu) are calculated at CAM-B3LYP functional. Therefore, the results of CAM-B3LYP are given in the main manuscript other results are given in supplementary material (Table S1). The short-range intermolecular interaction plays important role in describing optical and nonlinear optical properties of all these compounds due to which the results of CAM-B3LYP are better than the other two functionals.
Polarizability is the distribution of electron density Benzimidazole ring act as electron-rich species when phenyl ring having different substitution are attached. The shifting of electronic density occurs due to an extended conjugated system between electron-withdrawing and donating groups, the corresponding intersystem charge transfer (ICT) increases the β o . The compound 1b displayed the largest β o value (5.66 × 10 −30 esu), because of the strong electron-donating group (MeO-) at the para position of the aromatic ring, and on the other side of this molecule, the electron-withdrawing bromine group is attached. According to the push-pull mechanism, shifting of electronic density (ICT) occur under extended conjugation. A β o value of (4.61 × 10 −30 esu) was observed for compound 1a where a lower electron-donating methyl group is attached to the phenyl ring compared to the methoxy group. The β o values of 1d, 1e, and 1f compounds are 4.40 × 10 −30 , 4.00 × 10 −30 , and 4.30 × 10 −30 esu, respectively. All these compounds have electron-donating property-based substations, i.e., pyridine, di-methylbenzene, and thiazole rings, which enhance β o. In compound 1c, chlorine and bromine groups are present on opposite terminals of the compound and result in a lower β o value (3.81 × 10 −30 esu). Both halogen groups withdraw electronic density towards themselves, and the ICT transfer is low compared to other above-discussed compounds. For compounds 2a-2f, β o values are in a range from 2.00 × 10 −30 to 3.54 × 10 −30 esu. Compounds (2a-2f) have a lower β o compared to compounds (1a-2f), but the overall results are similar 1a-2f series. The reason is the acetamide group where resonance occurs between the nitrogen and carbonyl groups, and it also acts as an electron-donating group, although the extent of donation is low. Strong electron-donating methoxy of 2b (3.54 × 10 −30 esu), moderate electron-donating pyridine, and thiophene groups of 2d (2.40 × 10 −30 esu) and 2f (2.37 × 10 −30 esu) enhance delocalization of the π bonding and electron-withdrawing group bromine on the other side accept the electronic density. The shifting of density is more in the case of methoxy group as compared to thiophene and pyridine. The sulfur of thiophene (2f) has the opposite effect due to the electronegativity of sulfur, so its β o is lower than the methoxy group. A low electron-donating group methyl in compounds 2a (2.02 × 10 −30 esu), 2e (2.01 × 10 −30 esu) has a lower tendency to undergo ICT and have comparatively low β o . The lowest β o value is obtained for 2c due to electron-withdrawing groups' attachment on the opposite side but due to the presence of moderately electron-donating group the ICT is more compared to 1c and a value of 2.00 × 10 −30 is seen.

UV-VIS Absorption Analysis
UV-VIS absorption analysis is performed by using the TD-DFT method to get an insight into electronic excitation that occurs from a lower energy state to a higher energy state [52]. UV-VIS spectra of the parent compound (1) and substituted compounds (1a-2f and 2a-2f) are given in Figure 5. Results of oscillating strength (f o ), excitation energies (∆E), and wavelength (λ) of all compounds (1,1a-2f and 2a-2f) are given in Table 6 The redshift is seen in all substituted compounds (1a-2f) in comparison to the parent compound (1). The pronounced effect is seen for 1a compound in series of 1a-2f compounds. The electronic excitation takes place from donor-acceptor groups. In compounds where on both sides of the parent skeleton electron-donating and withdrawing groups are attached, this excitation is more prominent. In compound 1a, the electron-withdrawing bromine group (-Br) is attached, while, on the other side, the electron-donating methyl group (-CH 3 ) is attached. A prominent redshift of 248 nm wavelength is observed for it (1a). Based on the electrondonating capability of other groups, the corresponding redshift is observed. The increasing trend of wavelength is as follows: 1a (248 nm) > 1b (226 nm) > 1d (221 nm) > 1f (217 nm) > 1e (213 nm) > 1c (212 nm). Table 6. Oscillating strength (f o ), excitation energies (∆E) and wavelength (λ max ) of parent and substituted compounds (1,1a-2f and 2a-2f). tached, this excitation is more prominent. In compound 1a, the electron-withdrawing bromine group (-Br) is attached, while, on the other side, the electron-donating methyl group (-CH3) is attached. A prominent redshift of 248 nm wavelength is observed for it (1a). Based on the electron-donating capability of other groups, the corresponding redshift is observed. The increasing trend of wavelength is as follows: 1a (248 nm) > 1b (226 nm) > 1d (221 nm) > 1f (217 nm) > 1e (213 nm) > 1c (212 nm). UV-VIS analysis for 2a-2f series of compounds is also performed where hydrogen is replaced by an amide group ( Figure 5). Again, a redshift is seen compared to parent compound 1, which indicates the electronic excitation; the bathochromic shift is seen where peaks move toward a higher wavelength. Maximum wavelength is obtained for com- Figure 5. Calculated UV-VIS absorption spectra of parent (1) and substituted compounds (1a-1f (above) and 2a-2f (below)), using GaussView 05 software. The width at half maximum used for the Gaussian broadening is 28.98 σ.

Compounds
UV-VIS analysis for 2a-2f series of compounds is also performed where hydrogen is replaced by an amide group ( Figure 5). Again, a redshift is seen compared to parent compound 1, which indicates the electronic excitation; the bathochromic shift is seen where peaks move toward a higher wavelength. Maximum wavelength is obtained for compound 2b, where the peak is located at 222 nm. The amide group increases the wavelength of the compounds that have a lower wavelength in the 2a-2f series, because of the presence of electron-withdrawing groups at both opposite terminals (compound 2e). The maximum wavelength for these compounds 2a-2f ranges from 214 to 223 nm. Amide is an electron-withdrawing group, and it facilitates the push-pull mechanism and delocalization of electrons. The observed electronic excitations enable these compounds (1a-1f and 2a-2f) as efficient NLO materials usage in a second-harmonic generation. These compounds show transparency below 200 nm, so they can also be used in UV laser technology. Su and co-workers also worked on the single group substitution of benzimidazole; they observed that electron-donating groups decrease the wavelength, whereas the electron-withdrawing group increases the observed wavelength [53]. Their calculated wavelength is lower than our reported maximum wavelength for compounds 1a and 1b.

Experimental
All the purchased chemicals (Sigma Aldrich, St. Louis, MO, USA) and solvents were purified before use by distillation, or, for extra purity, some of the solvents were also dried in the lab. All the reactions were completed under and open-air environment. The reaction progress was checked by using a TLC card. Organic solutions were evaporated from the reaction by using a rotatory evaporator (Buchi, R-210, Allschwil, Switzerland) and vacuum pump (Buchi, V-700, Flawil, Switzerland). NMR spectra were calculated in CDCl 3 and DMSO-d 6 by using Bruker ARX 400 and 125 MHz FT-NMR spectrometers (Billercia, MA, USA). Different aryl boronic acid, bases, and solvents were used for the synthesis of required compounds by following specific methods. After that, 1.88 mmol of base Et3N or TMEDA, 1.128 mmol of aryl boronic acid, and 2 mL of water were added and stirred again for two hours. The reaction progress was monitored by TLC. After completion, the reaction mixture was filtered and concentrated by rotary evaporated to get residue [29]. The residue was column chromatographed to get the desired product, which was further analyzed by using NMR spectroscopic techniques.

Synthesis of N-(5-Bromo-1H-benzo[d]imidazol-2-yl) Acetamide (2)
The mixture of 5-bromo-2-aminobenzimidazole (0.1 g, 0.47 mmol) in 5 mL of acetic anhydride was heated at 40 • C for 4 h. The precipitates of the reaction mixture were cooled to 0 • C, then filtered off and washed with DCM and water. Finally, they were dried to give the final product with white crystals [28]. Then, the reaction mixture was stirred at room temperature, for 36 h; TLC examination was carried out to monitor the reaction completion. Upon completion, the reaction mixture was filtered and washed with ethyl acetate (5 mL). Then, the solvent was evaporated by rotary evaporation to obtain the residue [15,30]. The residue was column chromatographed by utilizing a mixture of hexane and ethyl acetate (60/40) to obtain the wanted product with a good % yield, which was further analyzed by using NMR spectroscopic techniques. The NMR spectra are given in Supplementary Material (Figures S1-S18).

Conclusions
Consequently, our investigation concentrated on the cross-coupling of 5-bromo-2aminobenzimidazole and aryl boronic acids under various reaction conditions. The impact of different bases and solvents was observed, and we noted that the reaction is not potent with Et 3 N and pyridine; however, the reaction products obtained a good yield with TMEDA. Finally, all synthesized derivatives (1a-1f and 2a-2f) were computationally studied by applying DFT calculations. Frontier molecular orbital (FMO) analysis provides information about the electronic properties and reactivity of compounds. Out of 1a-1f, compound 1d, and out of 2a-2f, compound 2d presented the lowest HOMO-LUMO energy gap of 4.25 and 4.79 eV, respectively. The isodensity shift was more pronounced for these compounds as compared to others. Reactivity descriptor parameters were also calculated to determine the chemical reactivity relations of all synthesized derivatives. Compounds 1d and 2d have the lowest values of chemical hardness (2.13 and 2.39 eV), lower electronic chemical potential (−3.89 and −4.26 eV), and highest values of electrophilicity index (3.56 and 3.80 eV), respectively. Also, Compounds 2d has lowest electronic chemical potential (−4.26 eV) among all compounds. The compounds 1a-1f displayed better results for all properties when compared to compounds 2a-2f, because the acetamide group shows conjugation between nitrogen and carbonyl groups, where electronic density shifting is low. Moreover, polarizability (α o ) and first hyperpolarizability (β o ) parameters were used to measured NLO response. Overall, 1a has a high αo value of 228 au, and 2e exhibited the highest value of 247 au due to electron-withdrawing and donating groups attached at opposite terminals of the rings. The highest value of βo is shown by 1b (5.66 × 10 −30 au). UV-VIS absorption analysis is performed to understand electronic excitation in the designed compounds. These results suggested the use of these compounds in the field of optics and nonlinear optics.