Synthesis Characterization and X-ray Structure of 2-(2,6-Dichlorophenylamino)-2-imidazoline Tetraphenylborate: Computational Study

The title compound tetraphenylborate salt of clonidine (Catapres®), 2-(2,6-dichlorophenylamino)-2-imidazoline tetraphenylborate (3), was prepared in 76 % yield by the reaction of 2-(2,6-dichlorophenylamino)-2-imidazoline hydrochloride (clonidine hydrochloride) (1) with sodium tetraphenylborate (2) in deionized water through anion exchange reaction at ambient temperature. The structure of the title borate salt was characterized by UV, thermal analysis, mass and NMR analyses. White crystals of (3) suitable for an X-ray structural analysis were obtained by slow growing from acetonitrile. The molecular structure of the titled compound (3) was crystallized in the acetonitrile, P21/c, a = 9.151 (3) Å, b = 12.522 (3) Å, c = 25.493 (6) Å, β = 105.161 (13)◦ V = 2819.5 (13) Å3, Z = 4. A DFT quantum chemistry calculation method was employed to investigate the interaction mechanism of clonidine with tetraphenylborate. The stable configurations of the complexes of clonidine with tetraphenylborate with electrostatic interactions were obtained. Finally, the interaction strength and type of the complexes were studied through the reduced density gradient (RDG) function. This study provides new theoretical insight into the interaction mechanism and a guide for screening and designing the optimal clonidine and tetraphenylborate reacting to form the complex.


Introduction
Clonidine is an imidazoline-derivative hypotensive agent and is a centrally acting alpha 2-adrenergic agonist. It crosses the blood-brain barrier and acts in the hypothalamus to induce a decrease in blood pressure. It may also be administered as an epidural infusion as an adjunct treatment in the management of severe cancer pain that is not relieved by opiate analgesics alone [1]. Clonidine may be used for differential diagnosis of pheochromocytoma in hypertensive patients. Other uses for clonidine include prophylaxis of vascular migraine headaches, treatment of severe dysmenorrhea, management of vasomotor symptoms associated with menopause, rapid detoxification in the management of opiate withdrawal, treatment of alcohol withdrawal used in conjunction with benzodiazepines, management of nicotine dependence, topical use to reduce intraocular pressure in the treatment of open-angle and secondary glaucoma and hemorrhagic glaucoma associated with hypertension, and in the treatment of attention-deficit hyperactivity disorder (ADHD) [2]. Clonidine also exhibits some peripheral activity. The chemical structure is 2-(2,6-dichlorophenylamino)-2-imidazoline hydrochloride [3].
The determination of clonidine by the potentiometric PVC membrane sensor was reported, using clonidine-tetraphenyl borate as electroactive material [4]. The proposed

X-ray Analysis
The titled compound of 3 was obtained as single crystals by slow evaporation from the acetonitrile solution of the pure compound at room temperature. Data were collected on "a Bruker APEX-II D8 Venture area diffractometer, equipped with graphite monochromatic Mo Kα radiation, λ = 0.71073 Å at 100 (2) K". Cell refinement and data reduction were carried out by Bruker SAINT. SHELXT [10,11] was used to solve the structure. The final refinement was carried out by full-matrix least-squares techniques with anisotropic thermal data for Appl. Sci. 2022, 12, 3568 3 of 18 nonhydrogen atoms on F. CCDC 1453705 contains the supplementary crystallographic data for this compound and can be obtained free of charge from the Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/data_request/cif (accessed on 20 January 2022).

Computational Study 2.4.1. Geometry Study
In this study, we examine the geometry of clonidine-tetraphenylborate as compared with data obtained by an X-ray diffraction experiment. We used the result as reference to benchmark the computed results for geometric parameters from a density functional theory (DFT) [12]. Starting with the experimental geometry, the geometry of clonidinetetraphenylborate was extensively optimized using the B3LYP functionals with the combined 6-31G(d,p). On the potential energy surface, computation of the Hessian matrix verified that all structures are minimal. Equation (1) was used to calculate the mean absolute error (MAE) between the computed and experimental data [12].
where x iT represents the chosen theoretical data, x iE represents the related experimental data, and n represents the number of terms in the sum.

Interaction Energies (IE)
The energy difference between the compound and the sum of the energies of separated fragments was used to compute the interaction energies (IE) of these molecules using Equation (2): The whole-molecule refers to the clonidine-tetraphenylborate system, while tetraphenylborate and clonidine are the separated molecules. This computation includes all local minimum structures, including zero-point and BSSE corrections.

Non-Covalent Interaction (NCI) Index
The non-covalent interaction (NCI) index [13], which maps the non-covalent interaction zone in real space qualitatively, enables qualitative visualization of non-covalent contact. The approach uses two scalar fields to map the bonding properties, namely electron density (q) and reduced density gradient (RDG, s). When s and q are combined, real space is roughly divided into two bonding regions: noncovalent interactions have low s and low q, whereas covalent interactions have high s and high q. H-bonding, van der Waals, and steric interaction can all be distinguished from each other using NCI plots. The non-covalent interaction zone is defined for low q and low s. In the negative low q region, the low s value is enticing; in the positive low q region, the value is repellent. The Multiwfn [14] application was used to conduct the NCI analysis [14].

Quantum Theory of Atoms in Molecules (QTAIM)
By employing Bader's QTAIM [15], it is possible to investigate the topology of electron density in the compound. Covalent contacts often have a substantial electron density (>0.2 a.u) at the bond critical point (BCP), and a large and negative Laplacian (∇ 2 ρ(r)), indicating that the interaction is strong and stable. As opposed to these interactions, when it comes to ionic, van der Waals, or hydrogen bond interactions, ρ(r) is modest (less than 0.10 au) and ∇ 2 ρ(r) is positive. Due to the fact that their features are expressed as characteristics of a real electron density of a system under investigation, QTAIM is considered one of the most suitable methodologies for analyzing various intra-and intermolecular interactions. Using the Multiwfn [14] program, we conducted a QTAIM analysis at the B3LYP/6-311G(d,p) level of theoretical analysis.

Molecular Electrostatic Potential (MESP)
The molecular electrostatic potential is a well-known method for predicting chemical reactivity [16][17][18][19][20]. The application of MESP was pioneered for the understanding of intermolecular interaction [21,22], which was a breakthrough. In the field of MESP, the most negative value of the MESP (V min ) is characterized as the lone pair or p electron density or "heap," while the most positive value of the MESP (V max ) is characterized as the electron deficient region or "hole" [23,24]. The noncovalent contact in these systems is caused by the interaction of this hole and heap. The MESP method was used to estimate the structures and interaction energies of clonidine and tetraphenylborate interactions. The position of the hole and heap in clonidine and tetraphenylborate were also characterized using these MESP methods, which allowed us to gain a better understanding of the structure of these noncovalent systems. MESP calculations were carried out with the help of the Multiwfn program [14].
The UV spectra of the titled compound 3 are shown in Figure 1 which show a shif the maximum absorption compared with compound 1 and 2.
The titled compound 3 was characterized by differential scanning calorimetry (DS which exhibits a significant exothermic peak at 160 °C. The differential scanning calor etry thermogram of the titled compound 3 is substantially depicted in Figure 2. The titled compound 3 was characterized by differential scanning calorimetry (DSC), which exhibits a significant exothermic peak at 160 • C. The differential scanning calorimetry thermogram of the titled compound 3 is substantially depicted in Figure 2.  The titled compound was heated from 35 °C to 650 °C at 10 °C/min. Thermal g metric analysis (TGA) of the titled compound (3) show two inflection points at abou °C and 300 °C compared with the start compounds (1) and (2), which show the infle at 325 °C and 425 °C, respectively. Figure 3 show the thermal gravimetric analysis o titled compound 3. The titled compound was heated from 35 • C to 650 • C at 10 • C/min. Thermal gravimetric analysis (TGA) of the titled compound (3) show two inflection points at about 200 • C and 300 • C compared with the start compounds (1) and (2), which show the inflection at 325 • C and 425 • C, respectively. Figure 3 show the thermal gravimetric analysis of the titled compound 3. The titled compound was heated from 35 °C to 650 °C at 10 °C/min. Thermal gravimetric analysis (TGA) of the titled compound (3) show two inflection points at about 200 °C and 300 °C compared with the start compounds (1) and (2), which show the inflection at 325 °C and 425 °C, respectively. Figure 3 show the thermal gravimetric analysis of the titled compound 3.

X-ray Crystallography
The title compound of C 33 H 30 BCl 2 N 3, the crystallographic data and refinement information are summarized in Table 1. The selected bond lengths and bond angles are listed in Table 2. The title compound crystallizes with one cation and anion in the asymmetric unit as shown in Figure 4, that is a 2-(2,6-dichlorophenylamino)-4,5-dihydro-1Himidazol-3-ium cation and a tetraphenylborate anion. As typical for tetraphenylborate, the tetraphenylborate anion is in a tetrahedral geometry around the B atom [C-B-C angles of 108.80 (19) • -110.13 (18) • ]. All the bond lengths and angles are in normal ranges [11]. The crystal structure is stabilized by C2-H2A···Cl2 hydrogen bonds along the c axis direction ( Figure 5, Table 3). Appl

X-ray Crystallography
The title compound of C33H30BCl2N3, the crystallographic data and refinement information are summarized in Table 1. The selected bond lengths and bond angles are listed in Table 2. The title compound crystallizes with one cation and anion in the asymmetric unit as shown in Figure 4, that is a 2-(2,6-dichlorophenylamino)-4,5-dihydro-1H-imidazol-3-ium cation and a tetraphenylborate anion. As typical for tetraphenylborate, the tetraphenylborate anion is in a tetrahedral geometry around the B atom [C-B-C angles of 108.80 (19)°-110.13 (18)°]. All the bond lengths and angles are in normal ranges [11]. The crystal structure is stabilized by C2-H2A•••Cl2 hydrogen bonds along the c axis direction ( Figure 5, Table 3).

X-ray Crystallography
The title compound of C33H30BCl2N3, the crystallographic data and refinement information are summarized in Table 1. The selected bond lengths and bond angles are listed in Table 2. The title compound crystallizes with one cation and anion in the asymmetric unit as shown in Figure 4, that is a 2-(2,6-dichlorophenylamino)-4,5-dihydro-1H-imidazol-3-ium cation and a tetraphenylborate anion. As typical for tetraphenylborate, the tetraphenylborate anion is in a tetrahedral geometry around the B atom [C-B-C angles of 108.80 (19)°-110.13 (18)°]. All the bond lengths and angles are in normal ranges [11]. The crystal structure is stabilized by C2-H2A•••Cl2 hydrogen bonds along the c axis direction ( Figure 5, Table 3).

Geometry Data
The mean absolute error, for the (12) bond distances and the (28) bond angles are listed in Table 4. The MAE values range from 0.022 to 0.002 for the complex bond distances. The complex bond angles have MAE values below 2 • . B3LYP was utilized to obtain these MAE values.

Interaction Energies (IE)
To obtain the most stable configurations and binding energies of structures, the DFT approach at the B3LYP level using the 6-31G+(d,p) basis set was widely used. In this research, clonidine and tetraphenylborate were investigated as a complex. The complexes (molar ratio 1:1) in the gas phase have binding energy (∆E b ), BSSE energy, and changes of Gibbs free energy (∆G) of, −107.22, 0.007642, and −87.43 kcal/mol, respectively. The clonidine-tetraphenylborate complex is the largest absolute value of ∆E b , and the BSSE energy indicates that it is the most stable clonidine-tetraphenylborate complex yet discovered. It also has a negative ∆G value, indicating that it is simple to generate spontaneously and that the result is consistent with the experimental method.

Non-Covalent Interaction (NCI) Index
On the basis of the RDG, the NCI between clonidine and tetraphenylborate was graphically established. NCI-RDG analysis can provide information on weak NCI in real space and depict different regions where NCI appear. Thus, using simple color codes, the approach can differentiate hydrogen bonding, van der Waals, and repulsive steric interactions [25,26]. discovered. It also has a negative G value, indicating that it is simple to generate spontaneously and that the result is consistent with the experimental method.

Non-Covalent Interaction (NCI) Index
On the basis of the RDG, the NCI between clonidine and tetraphenylborate was graphically established. NCI-RDG analysis can provide information on weak NCI in real space and depict different regions where NCI appear. Thus, using simple color codes, the approach can differentiate hydrogen bonding, van der Waals, and repulsive steric interactions [25,26].
High repulsion is indicated by red regions, strong attraction is indicated by blue regions, and electrostatic interactions are indicated by green regions. The sign of λ2 (the second largest eigenvalue of the Hessian electron density matrix) merges to create sign(λ2)ρ. Hydrogen bond interactions are distinguished from other interactions using the RDG function and sign(λ2)ρ. The values of sign(λ2)ρ range from −0.05 to 0.05 a.u., as illustrated in Figure 6a; −0.015 a.u. is the highest negative value in green, while 0.011 a.u. is the least positive value in reddish green, both of which are represented by the spikes on the left. The clonidine-tetraphenylborate complex's stability may be due to van der Waals and steric interactions [27].
The green and red regions of the RDG isosurface map (Figure 6a,b) suggest the presence of van der Waals and steric effects, respectively, inside the clonidine-tetraphenylborate complex. As demonstrated in Figure 6, the van der Waals and steric interactions are responsible for molecule binding, which is more favorable to system stabilization. This result agrees with the AIM analysis described in Section 2.4.  High repulsion is indicated by red regions, strong attraction is indicated by blue regions, and electrostatic interactions are indicated by green regions. The sign of λ 2 (the second largest eigenvalue of the Hessian electron density matrix) merges to create sign(λ 2 )ρ. Hydrogen bond interactions are distinguished from other interactions using the RDG function and sign(λ 2 )ρ. The values of sign(λ 2 )ρ range from −0.05 to 0.05 a.u., as illustrated in Figure 6a; −0.015 a.u. is the highest negative value in green, while 0.011 a.u. is the least positive value in reddish green, both of which are represented by the spikes on the left. The clonidine-tetraphenylborate complex's stability may be due to van der Waals and steric interactions [27].
The green and red regions of the RDG isosurface map (Figure 6a,b) suggest the presence of van der Waals and steric effects, respectively, inside the clonidine-tetraphenylborate complex. As demonstrated in Figure 6, the van der Waals and steric interactions are responsible for molecule binding, which is more favorable to system stabilization. This result agrees with the AIM analysis described in Section 2.4.  (55), there are six BCPs (BCP147, BCP140, BCP131, BCP108, BCP84, and BCP77) between clonidine and tetraphenylborate. The total electron density ρ(r), Laplacian electron density ∇ 2 ρ(r), and electron energy density H(r), which is composed of the electron potential energy density V(r), characterize the nature of chemical bonds and molecular reactivity. Table 5 lists the values of the AIM topological parameters for the complex of clonidine and tetraphenylborate.   As a result, the higher the (r) value at the BCP, the greater the concentration of electronic charge on the surface at this point and the stronger the considered contact.
It is a lot smaller than in covalent bonding, where the value of ρ(r) is about ∼10 −1 au. For hydrogen bonding, the value of ρ(r) is ∼10 −2 and for van der Waals interactions, the value of ρ(r) is ∼10 −3 au.
The absolute potential energy density |V(r)|/G(r), the (r) value at the BCP of interacting atoms and the value of its Laplacian ∇ 2 ρ(r), the total energy density H(r), and the ratio of the absolute potential energy density to the kinetic energy density |V(r)|/G(r) are the most commonly used terms to describe the nature and strength of bonding interactions. In order to explore topological characteristics, such as the Laplacian of the electron density (∇ 2 ρ(r)), potential energy (V(r)), and electron density (ρ(r)), among others. This approach resulted in the formation of a density gradient that originates from a location between two atoms, which is referred to as the critical bond point (BCP). The examination of this critical point BCP reveals the type and nature of the interaction between the clonidine and tetraphenylborate ions as a result of their interaction with their active sites. It should be remembered that the total electronic energy (H (r)) is equal to the sum of G(r) and V(r)) and was estimated from the BCP using Equation (9).
The interaction energy can be determined using the equation provided by Equation (10).
The parameters of the QTAIM analysis are shown in Table 4. ∇ 2 ρ(r) >0 is present in all interactions, indicating that they are non-covalent. This indicates the presence of intermolecular interactions. For the complex, the values of ∇ 2 ρ(r) and H(r) are positive, and the absolute potential energy density |V(r)|/G(r) is less than 1, indicating that the interactions are van der Waals interactions and weak electrostatic H-bonds. The strengths of interactions between the clonidine and tetraphenylborate are represented by the values of the electron density ρ(r). The values of ρ(r) are 0.00567 (u.a), 0.0151 (u.a), 0.00831 (u.a), 0.0213 (u.a), 0.0101 (u.a), and 0.00668 (u.a) for BCP147, BCP140, BCP131, BCP108, BCP84 and BCP77, respectively. Therefore, the values of ρ(r) indicate that N-H···C and N-H···Cl interacts more than C-H···C and C-H···Cl with the active site of the two ions. In addition, the interaction on N-H··· is more efficient and intensive because of the high values of ρ(r).

Reactivity Descriptors Global Reactivity Descriptors
Global reactivity descriptor for tetraphenylborate-clonidine complex parameters are shown in Table 6. Chemical reactivity changes with the structural configuration of molecules according to various characteristics. Clonidine has a lower chemical hardness (softness) value than the tetraphenylborate complex. Clonidine reacts as a cation, while tetraphenylborate reacts as an anion. As shown in Table 6, the tetraphenylborate structure has a greater electronegativity and electrophilicity index than clonidine. The electrophilicity of various chemical substances and reaction rates in biological systems were discovered to be related.

Local Reactivity Descriptors
The Fukui function (FF) is a technique for analyzing chemical reactions that offers information on the local site reactivity inside the molecule. These numbers refer to qualitative reactivity descriptors for various atoms in the molecule. The B3LYP/6-311 G (d, p) level of theory was used to create Fukui functions for electrophilic and nucleophilic attacks. The following equation was used to determine local Fukui functions (f + k , f − k ), and the dual descriptor, local softness values (s + k , s − k ), and local electrophilicity indices (ω + k , ω − k ) using Mulliken atomic charges of cationic and anionic states: for electrophilic attack : for radical attack : However, there are situations in which the results of both descriptions overlap, making it impossible to derive accurate conclusions. Rather than that, the dual descriptor Df(r) or DD can unambiguously characterize nucleophilic and electrophilic locations in a molecule. Thus, a graphical representation of the DD for both ions and complex is provided in Figure 8, clearly indicating the regions within the molecules where Df > 0 and Df < 0 occur, allowing for a better comprehension of these molecules' local chemical reactions. Other local reactivity descriptors are the local softness values and electrophilicity indices, which were calculated using the following equation.
where + and − signs show nucleophilic and electrophilic attacks, respectively. To compare the results obtained from the two ions and complex, the predictions for the specific reaction sites derived from the electrophilic and nucleophilic Fukui functions analyses for each ion were collated into two tables along with the local softness and electrophilicity indices. For completeness, the radical Fukui function f 0 k was provided, which can be thought of as an average of f + k , and f − k and denotes the most suitable locations for a radical attack. In addition to the tables, a graphical depiction of these descriptors is given in the form of a series of figures, such as the comparison of the results obtained for both ions and their complex under the accurate investigations.
By reviewing the numerical and graphical data shown in Tables 7 and 8 and Figure 8, it is possible to conclude that there is a good agreement between the numerical values and the graphical representations produced by both ions and the complex. When inspecting the numerical component, it is clear that the reactive sites predicted by the Fukui functions are accurate in all cases. This discovery is supported by considering the graphical representations of the Fukui functions and numerical results of the dual descriptor, which, despite being extended over the atoms within the molecules, are more compact than the other Fukui function descriptors. When tetraphenylborate and clonidine are utilized as anion and cation molecules, respectively, molecular electrostatic interaction between cation and anion molecules is thought to be one of the most fundamental driving forces for the formation of inclusion complexes. DFT calculations of the electrostatic potential for the cation and anion molecules were conducted to test the above idea. Figure 9A,G illustrates the outcomes. As a  Table 7. Fukui functions ( f + k , f − k ), local softness (s + k , s − k ), and local electrophilicity indices (ω + k , ω − k ) for selected atomic sites of tetraphenylborate, using Mulliken population analysis at B3LYP/6-311 G (d, p) level.  Table 8. Fukui functions ( f + k , f − k ), local softness (s + k , s − k ), and local electrophilicity indices (ω + k , ω − k ) for selected atomic sites of clonidine, using Mulliken population analysis at B3LYP/6-311 G (d, p) level. We provide Fukui functions, local softness values, and local electrophilicity indices for various atomic locations in clonidine and tetraphenylborate, as well as Figure 8. The greatest value for local nucleophilic reactivity descriptors f − k , s − k , and ω − k was at two phenyl groups in tetraphenylborate that were far away from the clonidine molecule in the complex, and that site is more susceptible to electrophilic assault. The largest value for local electrophilic reactivity descriptors f + k , s + k , and ω + k was shown in the clonidine molecule in a complex whose site is more sensitive to nucleophilic attack, as seen in Tables 7-9, which forms the bond between two compounds.

Molecular Electrostatic Potential (MESP)
When tetraphenylborate and clonidine are utilized as anion and cation molecules, respectively, molecular electrostatic interaction between cation and anion molecules is thought to be one of the most fundamental driving forces for the formation of inclusion complexes. DFT calculations of the electrostatic potential for the cation and anion molecules were conducted to test the above idea. Figure 9A,G illustrates the outcomes. As a map of electrostatic potential on an isoelectronic surface, the positive electrostatic potentials of amino groups are centered around the hydrogen atoms, shown in Figure 9A,C, while the negative electrostatic potentials are localized on nitrogen and chlorine atoms in clonidine molecules, shown in Figure 9B,C. The negative potential is concentrated on the carbon class boron section of the tetraphenylborate molecule. The anion of the tetraphenylborate molecule has no positive potential ( Figure 9E), but the positive charge is found in the para carbon of the phenyl groups in the neutral atom as shown in Figure 9F. It is worth noting that the positive potentials of clonidine's amino hydrogen match the negative potentials on the carbon of phenyl groups that connect with the boron atom on the tetraphenylborate molecule, and the negative potential of the tetraphenylborate molecule matches the positive potentials observed in the clonidine molecule in the complexes shown in Figure 7G. The charge distribution in clonidine and tetraphenylborate appears to be a primary driving mechanism for the creation of electrostatic attraction complexes. Table 9. Fukui functions ( f + k , f − k ), local softness (s + k , s − k ), and local electrophilicity indices (ω + k , ω − k ) for selected atomic sites of complex, using Mulliken population analysis at B3LYP/6-311 G (d, p) level.

Conclusions
In conclusion, the title compound 3 was prepared efficiently by the reaction of clonidine hydrochloride 1 with sodium tetraphenylborate 2 in deionized water at ambient temperature. The structure of the suggested compound (3) was confirmed by UV, thermal analysis, NMR, and mass and was established under the basis of its X-ray single-crystal analysis. This paper provides new theoretical insight into the selection and design of optimal clonidine with tetraphenylborate complex: at first, DFT calculations were used to obtain the optimized geometries of the clonidine with tetraphenylborate complex, and it was found to be matching with experimental geometries, which was provided from the X-ray single crystal. The theoretical results reported in this manuscript stress the importance of the interplay between noncovalent interactions involving different ions pair systems that can lead to stable complex. These effects were studied using genuine nonadditivity energies, the Bader's AIM theory and both reactivity descriptors (global and local), which offer information on the local site reactivity inside the molecule. Finally, a visualized analysis of interaction nature was explored by the RDG method in combination with interaction energy calculations to choose the optimal clonidine with tetraphenylborate complex. Generally speaking, the electronegativity of the anions in clonidine and tetraphenylborate determines the strength of the electrostatic bonds between clonidine and tetraphenylborate. It is concluded that tetraphenylborate with strong

Conclusions
In conclusion, the title compound 3 was prepared efficiently by the reaction of clonidine hydrochloride 1 with sodium tetraphenylborate 2 in deionized water at ambient temperature. The structure of the suggested compound (3) was confirmed by UV, thermal analysis, NMR, and mass and was established under the basis of its X-ray single-crystal analysis. This paper provides new theoretical insight into the selection and design of optimal clonidine with tetraphenylborate complex: at first, DFT calculations were used to obtain the optimized geometries of the clonidine with tetraphenylborate complex, and it was found to be matching with experimental geometries, which was provided from the X-ray single crystal. The theoretical results reported in this manuscript stress the importance of the interplay between noncovalent interactions involving different ions pair systems that can lead to stable complex. These effects were studied using genuine non-additivity energies, the Bader's AIM theory and both reactivity descriptors (global and local), which offer information on the local site reactivity inside the molecule. Finally, a visualized analysis of interaction nature was explored by the RDG method in combination with interaction energy calculations to choose the optimal clonidine with tetraphenylborate complex. Generally speaking, the electronegativity of the anions in clonidine and tetraphenylborate determines the strength of the electrostatic bonds between clonidine and tetraphenylborate. It is concluded that tetraphenylborate with strong electronegativity will be considered for interaction as an anion (donor), while the clonidine with law electronegativity will be considered for interaction as a cation (acceptor).