Next Article in Journal
Novel Gel Formulations as Topical Carriers for the Essential Oil of Bidens tripartita for the Treatment of Candidiasis
Next Article in Special Issue
Optimization of the Synthesis of Glycerol Derived Monoethers from Glycidol by Means of Heterogeneous Acid Catalysis
Previous Article in Journal
Sorbitol as a Chain Extender of Polyurethane Prepolymers to Prepare Self-Healable and Robust Polyhydroxyurethane Elastomers
Previous Article in Special Issue
Application of Ionic Liquid-Based Ultrasonic-Assisted Extraction of Flavonoids from Bamboo Leaves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lewis Acidity and Basicity of Mixed Chlorometallate Ionic Liquids: Investigations from Surface Analysis and Fukui Function

College of Chemistry and Chemical Engineering, Inner Mongolia University, Hohhot 010021, China
*
Author to whom correspondence should be addressed.
Molecules 2018, 23(10), 2516; https://doi.org/10.3390/molecules23102516
Submission received: 29 August 2018 / Revised: 23 September 2018 / Accepted: 29 September 2018 / Published: 30 September 2018
(This article belongs to the Special Issue Green and Sustainable Solvents)

Abstract

:
Mixed chlorometallate ionic liquids (ILs) have been regarded as potential solvents, catalysts, and reagents for many organic processes. The acidity and basicity of these ILs were correlated with theoretically estimated parameters such as electrostatic surface potential maxima and minima, average local surface ionization energy, and Fukui and dual descriptor functions. The introduction of metal chloride into the anions would influence the acidity/basicity of ILs by withdrawing the electron density from the cationic counterpart. For the [C4mim]-based ILs with the mixed-metal anions, the acidity tends to attenuate while the basicity becomes stronger, as compared to the corresponding chloroaluminate ILs. However, the acidity of [(C2H5)3NH]-based ILs with the mixed-metal anions are greater than that of the net chloroaluminate ILs. The Fukui function values showed that most of the mixed chlorometallate ILs belong to bifunctional distribution. The mixed chlorometallate ILs both have electrophilic and nucleophilic sites, which would be beneficial for their applications.

1. Introduction

Room-temperature ionic liquids (ILs) have received much attention in the last few decades [1]. Ionic liquids exhibit interesting physical properties, such as low vapor pressure, a large liquidus range, high thermal stability, and favorable solvation behavior [2]. A key feature of ionic liquids is the large number of ion combinations that are possible, which suggests the feasibility of designing a suitable liquid for a specific task [3]. Due to the versatility in their physical and chemical properties, ionic liquids have been widely studied in many fields [4]. These fields usually include the production of biodiesel [5], the production of tensides/detergents by the Friedel-Crafts reaction [6,7,8], polymer chemistry [9], liquid-liquid extraction [10,11], and even enzymatic reactions [12]. On the other hand, although a great number of ionic liquids theoretically exist [13], only hundreds of them are commercially available and used as catalysts [14]. Among these ionic liquids, chlorometallate ionic liquids hold an important place in the catalysis field [15].
Primary applications of chloroaluminate ionic liquids are industrial Friedel-Crafts alkylation and acylation, and oligomerization and isomerization reactions of olefins and aromatics [16]. For example, a pilot-scale oligomerization of the olefins process at BP and an alkylation of the benzene process of Akzo-Nobel all used acidic chloroaluminate ionic liquids [(CH3)3NH]Cl-AlCl3 as catalysts [17]. Zhang et al. reported on the isomerization of n-pentane with [(C2H5)3NH]Cl-AlCl3 [18]. Huang et al. [19] and Wang et al. [20] studied the isomerization of dicyclopentadienes by [Hpy]Cl-AlCl3. In particular, the alkylation reaction of isobutane and butenes catalyzed by chloroaluminate ionic liquids, such as [C4mim]Cl-AlCl3 [21], [(C2H5)3NH]Cl-AlCl3 [22], and amide-AlCl3-based IL [23], have been studied intensively for years [15]. Compared with chloroaluminate ILs, the chlorometallate systems that only contain Fe-, Cu-, Zn-, In-, or Sn-based anions usually present lower Lewis acidity [24]. Thus, there are only comparatively few examples in their performance as acid catalysts. Cui et al. [25] have found that the alkylation of benzene with p-formaldehyde could be catalyzed by a series of acidic chlorometallate ionic liquids, such as [(C2H5)3NH]Cl-FeCl3, [(C2H5)3NH]Cl-ZnCl2, and [(C2H5)3NH]Cl-SnCl2. Gunaratne et al. also reported that phenols and catechols could be alkylated in weakly acidic chloroindate ionic liquids. Particularly, this acidic system could control the product selectivity and prevent polyalkylation that would occur using stronger acids [26]. In addition, [C4mim][FeCl4] has been used as a mild Lewis acid to catalyze the benzylation of arenes [27].
Adding metal salts to chlorometallate ionic liquids can either lead to changes in chlorometallate anion activity, or to the coordination of the added metal ions by chlorometallate anions [28]. This feature might be one of the key advantages for some reactions. For example, we used the [(C2H5)3NH]Cl-AlCl3-CuCl system for isobutane alkylation, and high quality gasoline was obtained [29,30]. Zinurov et al. [31] found that using the acidic [(CH3)3NH]Cl-AlCl3 ionic liquid doped with copper salts as catalysts could control the route of the isomerization of n-pentane. It is also found that adding titanium tetrachloride to the slightly acidic [C4mim]Cl-AlCl3 ionic liquid could lead to the formation of high yields of branched polymers [32]. Additionally, when the mixed chlorometallate [(C2H5)3NH]Cl-FeCl3-CuCl was used, the conversion and selectivity of the isobutene oligomerization could be enhanced significantly [33]. The best known example of catalysis in mixed chlorometallate ionic liquids may be the Institut Francais Du Petrole (IFP) Difasol process for olefin dimerization, which uses the [C4mim][AlCl4]-EtAlCl2-NiCl2 ionic liquid as a catalyst [34,35].
In order to improve the efficiency of IL-based processes, it is necessary to obtain accurate knowledge of physicochemical properties of ILs. Several experimental studies on the acidity of chlorometallate ILs have been reported. Using aromatic hydrocarbons as indicators, the acidity of HCl in chloroaluminate ILs can be determined by ultraviolet-visible (UV-vis) spectroscopy [36]. Osteryoung and co-workers [37,38], and Swadzźba-Kwasśny and co-workers [39] measured the Gutmann acceptor number (AN) of chlorometallate ILs using the chemical shift in the 31P NMR spectroscopy of triethylphosphine oxide (TEPO) as a probe molecule. X-ray photoelectron spectroscopy (XPS) was also employed to study the Lewis acidity and hydrogen bond basicity of halometallate-based ILs [40]. In addition, Kou [41] and Hu [42] used pyridine or acetonitrile as probe molecules to determine the acidity of chloroaluminate ILs by infrared spectroscopy.
The strength of interaction of a Lewis acid-base pair depends on the size, shape, and relative energies of the acid and the base. Consequently, the strength of a Lewis acid depends on the base that it is interacting with [15]. However, the impossibility of establishing a universal scale of Lewis acidity/basicity does not prevent the determination of the quantitative behavior of Lewis acids/bases [43]. There are a number of established Lewis acidity scales by experimental studies like the above mentioned, and the strength is also quantified by computational calculations [44]. Contreras et al. [45] have introduced two models of the Lewis molecular acidity: the excess electronic chemical potential and the local charge capacity. Furthermore, regional electrophilic and nucleophilic Fukui functions were also proposed to explain the Lewis acidity and basicity of ILs [46]. According to these models, several quantum-chemical parameters were used to determine the acid-base properties of the ILs. For example, Parveen et al. have performed a precise surface and Fukui analysis for several chlorometallate ionic liquids. The acidity and basicity of these ILs were correlated with theoretically parameters [47]. Similar works have also been done by Wu and co-workers [48].
As mentioned above, the isobutane-butene alkylation reaction catalyzed by mixed chlorometallate ILs is an important chemical process. The improvement of product selectivity is often attributed to the introduction of metal salts to the chloroaluminate ILs. Generally, it is believed that adding metal salts would reduce the super acidity of the chloroaluminate ILs [49]. However, the acidity of these ILs is still ambiguous. In addition, why are mixed chlorometallate ILs are better Lewis acidic catalysts than chloroaluminate ILs? What are the roles of anion and cation on the total acidity of ILs? These problems are also important for the better application of the mixed chlorometallate ILs. In this work, the acidity and basicity of mixed chlorometallate ILs were investigated by the electrostatic analysis of the molecular surface. In order to quantitatively determine the Lewis acidity of a region, the maximum surface electrostatic potential (Vs,max) was employed. Additionally, the basicity of the ILs was estimated by the minimum surface electrostatic potential (Vs,min) and the lowest average local ionization energy on the surface ( I ¯ s,min). Moreover, the Fukui function and the dual descriptor function were also applied to assess the roles of anions and cations on the Lewis acidity/basicity of ILs.

2. Results and Discussion

2.1. Acidity and Basicity of the Chlorometallate ILs

In this work, the most commonly used chloroaluminate ILs, i.e., [(C2H5)3NH]Cl-AlCl3 (triethylamine hydrochloride chloroaluminate IL), and [C4mim]Cl-AlCl3 (1-butyl-3-methylimidazolium chloroaluminate IL were selected for study. Some non-chloroaluminate ILs, [(C2H5)3NH]Cl-ZnCl2, and [C4mim]Cl-CuCl, were also employed for comparison. Figures S1–S16 shows the optimized geometries of these ILs. To some extent, the acidity of an IL depends on its ability to form electrostatic interactions with basic species. The surface electrostatic potential (Vs,max) may be one useful tool for investigating and predicting the acidic region of ILs at the microscopic level. The larger magnitude of Vs,max usually indicates that the ions have a stronger interaction or acidity. Some quantum-chemical parameters are listed in Table 1.
For the same cation, the acidity of chlorometallate ILs (Vs,max) showed a decrease in the order of Al > Zn > Cu, which is in agreement with the electronegativity trend of the three metals [37]. Different cations also significantly influence the Lewis acidity of the ILs with the same anion (e.g., [(C2H5)3NH]+[AlCl4] and [C4mim]+[AlCl4]). When the [C4mim]+ cation was presented, the value of Vs,max could reach to 48.275 kcal/mol. Additionally, it is well known that the Lewis acidity of chloroaluminate ILs can be controlled by the mole fraction of aluminum chloride used for the IL preparation. If the mole ratio of aluminum chloride to organic salt is greater than 2:1, the chloroaluminate IL would exhibit super-acidity and [Al2Cl7] would be the dominated anion [50]. The surface analysis parameters support the above conclusion. The Vs,max value of [C4mim]+[Al2Cl7] is much higher than that of [C4mim]+[AlCl4]: 52.596 vs. 48.275. Similarly, [(C2H5)3NH]+[Al2Cl7] has a larger Vs,max value than that of [(C2H5)3NH]+[AlCl4]. The non-chloroaluminate ILs [(C2H5)3NH]+[ZnCl3] and [(C2H5)3NH]+[Zn2Cl5] have the same trend. Therefore, in accordance with the experimental data, the Lewis acidity of ILs will become stronger with the increase of the mole fraction of metal chloride.
Figure 1 graphically illustrates the molecular surface properties of the ion pairs for the studied ILs. As Figure 1 shows, the locations of Vs,max for the ILs were all around the cations, which indicated that the acidity of the ILs resulted from the cation. This result was in conformity with the conclusions of other studies [4].
The magnitude of Vs,min can indicate the electrostatically reactive abilities of the molecules, while I ¯ s,min can show the electron-transfer ability. The two parameters often indicate the regional basicity of the ILs. Generally, a larger absolute value of Vs,min and a smaller I ¯ s,min usually indicate that the IL will exhibit a stronger basicity than the others. For the chlorometallate ILs, it is found that the sites of Vs,min and I ¯ s,min are all around the Cl atom in the anions (Figure 1), which means that the basicity of ILs is mainly focused on the anions. In general, the much larger size of the anions always leads to a more dispersed distribution of the negative charge. Therefore, the higher order complexation anions (e.g., [Al2Cl7] and [Zn2Cl5]) would have less basicity, consistent with the results from the Vs,msx analysis. Additionally, the metal chloride could markedly affect the basicity of the chlorometallate ILs, as Table 1 shows.

2.2. Acidity and Basicity of the Mixed Chlorometallate ILs

Although chloroaluminate ILs have been widely used in solvent and catalysis processing, only a few surface analyses have focused on the mixed chlorometallate ILs (e.g., [(C2H5)3NH]Cl-AlCl3-CuCl, [C4mim]Cl-AlCl3-AgCl). Herein, a series of typical mixed chlorometallate ILs were chosen for investigation. These ILs are composed of the chloroaluminate anion and several other chlorometallate anions. The acidities of the ILs were compared with the conventional chloroaluminate ILs. Additionally, some ionic solvents for olefin absorption (e.g., benzene-CuAlCl4 and ether-CuClAl4 [51,52]) were also employed for comparison.
In general, the acidities of the chloroaluminate ILs were stronger than those of the corresponding mixed metal ILs; e.g., [C4mim]Cl-AlCl3 > [C4mim]Cl-AlCl3-CuCl (Table 1 and Table 2). Many studies believe that the Lewis acidity of chloroaluminate IL would be reduced by introducing another metal chloride [39,53]. Additionally, the ion pairs [C4mim]+[CuAlCl5] or [C4mimCu]+[AlCl4] could be formed in [C4mim]Cl-AlCl3-CuCl. As Table 2 shown, a surface analysis found that [C4mim]+[CuAlCl5] and [C4mimCu]+[AlCl4] all led to smaller Vs,max values. The location of Vs,max in [C4mim]+[AlCl4] was around the H atom of an imidazolium ring, similar to the Vs,max sites in copper−containing ILs (Figure 1 and Figure 2). Owing to the smaller values of Vs,max, copper-containing IL [C4mim]Cl-AlCl3-CuCl exhibited weaker acidity than that of [C4mim]Cl-AlCl3. These results are consistent with the previous conclusions.
From the determination of the Gutmann acceptor number, we learned that the acidity of [(C2H5)3NH]Cl-AlCl3-CuCl is larger than that of [(C2H5)3NH]Cl-AlCl3 [54]. Some researchers have considered that [(C2H5)3NCu]+ and [AlCl4] are the dominant ions of [(C2H5)3NH]Cl-AlCl3-CuCl. In other words, CuCl coordinates with the [(C2H5)3NH]+ cation to form [(C2H5)3NCu]+ and to release HCl [55]. However, the Vs,max value of [(C2H5)3NCu]+[AlCl4] is much smaller than that of [(C2H5)3NH]+[AlCl4] (Table 1), which is not consistent with the experimental results. By contrast, the Vs,max of [(C2H5)3NH]+[CuAlCl5] is larger than that of the corresponding net chloroaluminate IL, indicating that [(C2H5)3NH]+[CuAlCl5] exhibits a stronger acidity than that of [(C2H5)3NH]+[AlCl4]. For comparison, the benzene–CuAlCl4 solvent was also employed for the surface analysis, and the Vs,max of benzene-CuAlCl4 was higher than that of benzene-[AlCl4]. All of the results suggest that a large number of [(C2H5)3NH]+ and [CuAlCl5] should be presented in the ionic liquid [(C2H5)3NH]Cl-AlCl3-CuCl.
Vs,min and I ¯ s,min can determine the ability of ILs to act as H-bond acceptors. No matter which type of cation is, the sites of Vs,min for the studied ILs are all around the Cl atom of the anion. Nonetheless, the Vs,min of the chlorometallate ILs was more negative than those of the net chloroaluminate ILs. Therefore, the mixed chlorometallate ILs would be stronger H-bond acceptors with respect to the conventional chloroaluminate ILs.
As Figure 1 and Figure 2 shown, the I ¯ s,min site of the chloroaluminate ILs were located around the chloride atom, whereas the I ¯ s,min locations of the mixed chlorometallate ILs were all around the transition metal atoms (e.g., Cu, Ag) of the ions. A low value of I ¯ s suggests that the electron at this position was not tightly bounded, and the site with the lowest I ¯ s (i.e., I ¯ s,min) on the surface is usually recognized as the location where an electron is mostly inclined to escape from the molecule. Therefore, the electron of the transition metal atom may leave the molecular surface more easily, leading to I ¯ s,min sites around the transition metal atom in mixed chlorometallate ILs. In addition, for triethylamine hydrochloride-based ILs, less negative interaction energy of [(C2H5)3NCu]+[AlCl4] (Table 2) indicates that the cation–anion interaction of [(C2H5)3NCu]+[AlCl4] is weaker than that of [(C2H5)3NH]+[CuAlCl5]. Thus, the [(C2H5)3NCu]+ cations are more inclined to leave the [AlCl4] anions and become unstable.
Generally, the Vs,max values have a strong dependence on the cation-anion charge transfer (QCT). To further investigate the Lewis acidity, the atomic charges and charge transfer properties of the ILs under study were calculated, and the data are listed in Table 3. Generally, chlorometallate-based ILs were synthesized by direct addition of the required metal chloride MClx to C4mim-Cl in an appropriate molar ratio and under an inert atmosphere. Thus, the charge transfer QCT should take place from the ion pair unit C4mim-Cl, in which the Cl atom of C4mim-Cl is close to the imidazole ring, and it forms a H-bond with the C2-H moiety, to the MClx unit. Thus, the charge transfer in the IL can be described as QCT(IP→MClx). As Figure 3 and the data collected in Table 3 show, the charge transfer QCT of [C4mim]+[AlCl4] indeed occurs from the ion pair unit C4mim-Cl to the AlCl3 unit. Moreover, a larger amount of charge transfer corresponds to the stronger acidity of the ILs. Relative to [C4mim]+[CuAlCl5], QCT(IP→MClx) appears to be greater in magnitude for the [C4mim]+[AlCl4] that exhibits stronger acidity: 407 me vs. 345 me. Similarly, [C4mim]+[Zn2Cl5] shows stronger acidity than that of [C4mim]+[ZnCl3]. These results are consistent with the conclusions drawn from the Vs,max values.

2.3. Analysis of Fukui Function for Mixed Chlorometallate ILs

The Fukui function (FF) is an important concept in conceptual density functional theory, which is defined as Equation (1) [56]:
f ( r ) = ( ρ ( r ) N ) + = ρ N + 1 ( r ) ρ N ( r )
where N is the number of electrons in the present system. Cerda-Monje et al. have proposed that the condensed to atom electrophilic FF f atom + can be considered as the natural distributor of Lewis molecular acidity (Equation (2)), and the nucleophilic FF f atom is the corresponding distributor of Lewis molecular basicity (Equation (3)) [45,46]. Thus, the condensed Fukui function of the molecule might be readily calculated by using the charges of the atoms:
f atom + = ρ N + 1 ( r ) ρ N ( r ) = q N atom q N + 1 atom
f atom = ρ N ( r ) ρ N - 1 ( r ) = q N - 1 atom q N atom
If we integrate over the corresponding basins Ω k ± (local area of a molecule), the basins FF could be produced as Equation (4):
N Ω ± = Ω atom ± ρ ( r ) d r = atom Ω atom ± f atom ±
The basins Ω atom ± represent the maximum molecular regions of the electrophilic FF ( Ω atom + ) and nucleophilic FF ( Ω atom ), respectively. Thus, the region with the larger N Ω + is the favorable site to accept the electron charge, while the region with the greater N Ω is the location that prefers to donate the electron charge. According to the suggestion of Cerda-Monje et al., the concept of “normal distribution’’, ‘‘bifunctional distribution’’, and ‘‘border-line distribution’’ was employed to describe the studied ILs [46]. The “normal distribution” indicates that the Lewis molecular acidity ( N Ω + ) is mainly centered on the cation fragment, and the Lewis molecular basicity ( N Ω ) is mostly centered on the anion. In our cases, the N Ω + of the cation was usually greater than 0.9 eV, and N Ω was less than 0.2 eV. Meanwhile, the values of the anion moiety should be N Ω + < 0.2 eV and N Ω > 0.9 eV, respectively. The “bifunctional distribution” usually indicates that both the electrophilic and nucleophilic reactive sites are located at the same region of the ILs. When N Ω + and N Ω are great than 0.5 eV and 0.8 eV, respectively, the anion will exhibit “bifunctional distribution”. For the cation, the “bifunctional distribution” means that N Ω + > 0.8 eV and N Ω > 0.5 eV. The remaining values belong to the “borderline distribution”, which will not clearly indicate the location that prefers to accept or donate the electron charge.
As the data collected in Table 4 shown, most net chlorometallate ILs (e.g., [(C2H5)3NH]+[AlCl4], [(C2H5)3NH]+[ZnCl3]) exhibited a normal distribution, indicating that the Lewis acidity and basicity were located on the cation and the anion, respectively. However, the IL [(C2H5)3NH]+[CuAlCl5] showed a bifunctional distribution, implying that the location to accept and to donate the electronic charges were all situated on the anions [CuAlCl5]. As Figure 4 shows, the presence of the Cu atom enhances the ability of electrophilic attack of the anion. In previous work, we found that the [(C2H5)3NH]+[CuAlCl5] ionic liquid could produce better alkylation gasoline than that of net [(C2H5)3NH]+[AlCl4] [29,30]. According to the experimental results, the authors proposed a possible mechanism and raised a reasonable explanation for the possible mechanism, that is, the higher catalytic selectivity is attributed to the interaction between iC4+ carbonium ion and the transition metal atom [57]. As demonstrated in this work, the IL [(C2H5)3NH]+[CuAlCl5] exhibited a bifunctional distribution ( N Ω + > 0.5 eV and N Ω > 0.8 eV for the anion), indicating that the anion not only can provide Lewis acidity, but it can also absorb a considerable number of iC4+ carbonium ion.
A dual descriptor is another useful function that is used to reveal the reactive sites [58]. Based on the Fukui function, the dual descriptor can be easily calculated by f atom + and f atom (Equation (5)).
f atom = f atom + f atom
If the distribution of f atom around a site A is more positive than another site B, then one can say that A is a more favorable site for nucleophilic attack than B, and in the meantime, B is a more preferential site for electrophilic attack than A. As can be seen from Table S2, the Cu atom of [(C2H5)3NH]+[CuAlCl5] gave a more negative f atom value than that of any other atom; meanwhile, the Al atom showed a maximum f atom value. These results again indicate that the mixed chlorometallate anion [CuAlCl5] had both favorable sites for nucleophilic and electrophilic properties with respect to the net chloroaluminate anion [AlCl4]. The ionic liquid [(C2H5)3NH]+[AgAlCl5] also exhibited bifunctional distribution. The anion N Ω ± values were 0.879 and 0.932, respectively.
Fukui analysis further indicates that many mixed chlorometallate ILs, such as [C4mim]Cl-AlCl3-CuCl, [C4mim]Cl-AlCl3-AgCl, [(C2H5)3NH]Cl-AlCl3-ZnCl2, benzene-CuAlCl4, and ether-AgAlCl4 all can show bifunctional behavior. For imidazolium-based ionic liquids, the data of Fukui function shows that N Ω + > 0.5 eV and N Ω > 0.8 eV for [CuAlCl5] anion, and N Ω + > 0.8 eV and N Ω > 0.5 eV for [C4mimAg]+ cation. The results strongly implies that both nucleophilic and electrophilic sites are located on the same ion of [C4mim]Cl-AlCl3-MCl (M=Cu, Ag) ILs. Additionally, ∆f of the transition metal atoms in [(C2H5)3NH]Cl-AlCl3-CuCl, [C4mim]Cl-AlCl3-CuCl, [(C2H5)3NH]Cl-AlCl3-AgCl, ether-AgAlCl4, and [(C2H5)3NH]Cl-AlCl3-ZnCl2 ILs all gave the most negative values (see Tables S2, S3, S8, S9 and S12)), suggesting that the corresponding atoms were the most favorable sites for electrophilic (iC4+) attack. The transition metal mixed ILs could provide a higher catalytic selectivity of C4 alkylation as compared to the traditional catalysts (e.g., H2SO4) and the conventional chlorometallate ILs (e.g., [C4mim]Cl-AlCl3) [57]. Other than as a Lewis acid, it is believed that the mixed chlorometallate ILs would attract more iC4+ to the surrounding transition metal atoms. Thus, better catalytic selectivity can be obtained due to an increase of the I/O (iC4+ to olefins) ratio around the anion.

3. Methods

Recently, an assessment of density functional theory and Møller-Plesset perturbation theory for ILs has been performed [59]. Li et al. have recommended that the density functional of the Minnesota family of the M0X type with a diffusion function basis set (aug-cc-pvdz) would give reliable results for IL calculations [60]. Similarly, Zahn et al. noted that dispersion-corrected density functionals, such as B3LYP/6-31++G**-gd3bj (DFT-D3), could also lead to reliable results for the IL calculations [61]. In this work, we selected the M06-2X/aug-cc-pVDZ level of theory to optimize the geometries of the ILs under study [62]. This level is known to afford good geometry at low computational cost, which is particularly appropriate to be used for the calculation of cation-anion interactions [63]. To calculate the Ag-containing ILs, the basis set aug-cc-pVDZ-PP (Ag, 28 core electrons) that was obtained from the EMSL Basis Set Exchange was employed, and the Dunning’s correlation-consistent basis set aug-cc-pVDZ was applied for the remaining atoms [64]. Optimized geometries of the studied ILs were confirmed as minima on the potential energy surface by the absence of imaginary vibrational frequencies. The Gaussian 09 package [65] was employed to carry out all of these computations.
In this work, the scheme of restrained electrostatic potentials (RESP) employing the ChelpG point selection algorithms and the density matrix partitioning scheme Natural Population Analysis (NPA) were employed to calculate the atomic partial charges in the ionic liquids. Rigby and Izgoridina [66] have recently assessed different atomic partial charge schemes for polarization and charge transfer effects in ionic liquids. They suggest that RESP schemes are preferred when producing atomic partial charges, as RESP can show an appreciable degree of charge transfer. Among the RESP schemes, ChelpG is the least systematic for ion pairs of ionic liquids, and it should be used. Therefore, the RESP-ChelpG model was selected, and the atomic charges were finally calculated by the Multiwfn 3.6 (dev) program [67]. Because of largely basis set independent, the NPA scheme was also selected for the calculation. NPA atomic charges were obtained by the natural bond orbital (NBO) program [68] using M06-2X/aug-cc-pVDZ(-PP) for the ion pairs. Though the charges from the two schemes were slightly different in magnitude, the conclusions drawn are almost no difference.
Electrostatic potential, V(r), has been widely used for prediction of nucleophilic and electrophilic sites, as well as the molecular recognition mode, which is defined as Equation (6):
V ( r ) = A Z A | r - R A | ρ ( r ) | r - r | d r
where RA and ZA denote the position vector and the nuclear charge of atom A, respectively. If the pseudo-potential was used, ZA would be the number of explicitly expressed electrons. The analyses of electrostatic potential are commonly performed on molecular Van der Waals (vdW) surfaces. In our analyses, we took the 0.001 a.u. isosurface of the electron density as the vdW surface. Although the definition of such a surface is arbitrary, the above definition may well reflect specific electron structural features of a molecule. The strength and orientation of weak interactions, including regional acidity, hydrogen bonding, and halogen bonding, can be well predicted and explained by analysing the magnitudes and positions of minima (Vs,min) and maxima (Vs,max) on the surface [67]. The surface average local ionization energy I ¯ ( r ) can be written as Equation (7):
I ¯ ( r ) = i ρ i ( r ) | ε i | ρ ( r )
where ρ(r) and εi are the electron density function and orbital energy of the ith molecular orbital, respectively. Typical DFT functionals are suitable for computing I ¯ ( r ) . A lower value of I ¯ ( r ) indicates that the electrons at this point are more weakly bounded. It is proven that the minima of the average local ionization energy ( I ¯ s,min) on the vdW surface are good indicators to reveal the basicity of ILs [69].
The Multiwfn 3.6 program was employed to determine the values of Vs,min, Vs,max, and I ¯ s,min, respectively. Using the optimized geometries, wave functions, and electrostatic potential cube files were generated from Gaussian 09. The graphs of the electrostatic potential and the average local ionization energy were depicted by the visualization tool, VMD [70]. The regional Fukui function and dual descriptor, which are two well-known local descriptors for electron gain and donation, were also calculated by the Multiwfn 3.6 (dev) program.
Cation–anion interaction energies of the ILs were calculated according to Equation (8). Furthermore, the basis set superposition error (BSSE) was corrected using the counterpoise method when calculating the interaction energies [71]:
EIL(kcal/mol) = 627.51 × [Eion-pair (au) − Ecation (au) − Eanion (au)]

4. Conclusions

Using a surface analysis, Fukui function, and a dual descriptor function, the acidity and basicity of the mixed chlorometallate ILs were investigated. The quantum-chemical parameter Vs,max was adopted to quantitatively describe the Lewis acidity of ILs, while Vs,min, and I ¯ s,min were proposed to interpret their basicity. Many mixed metal-containing ILs exhibit weaker Lewis acidity with respect to the corresponding chloroaluminate ILs. However, due to the presence of the coordinated anions MAlCl5 (M=Cu, Ag), [(C2H5)3NH]Cl-AlCl3-MCl has higher acid strength than that of [(C2H5)3NH]Cl-AlCl3. Mixed chlorometallate ILs have more complicated ion species, while the MAlCl5 ion and MAlCl4 complex would abound in the ILs. The analyses of the Fukui function and the dual descriptor function show that most of the mixed chlorometallate ILs belong to the bifunctional distribution. Both the electrophilic and the nucleophilic sites are located in the same region of ILs, which have benefits in improving the product selectivity of the isobutane–butene alkylation process. These results can provide a deep understanding of the nature of the chlorometallate ILs. Moreover, they might give useful instruction for preparing novel functional ILs.

Supplementary Materials

The following are available online, Figures S1–S16: structures of optimized ionic liquids, Tables S1–S12: calculation data of Fukui function and dual description, Tables S13–S28: Cartesian coordinates of ILs.

Author Contributions

Methodology, Y.L.; Investigation, J.W. All authors have read and approved the final version.

Funding

This research and the APC were funded by the National Natural Science Foundation of China (grant numbers: 21766021 and 21266015).

Acknowledgments

The images of surface analysis were made with VMD software support. VMD is developed with NIH support by the Theoretical and Computational Biophysics group at the Beckman Institute, University of Illinois at Urbana-Champaign.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Radai, Z.; Kiss, N.Z.; Keglevich, G. An overview of the applications of ionic liquids as catalysts and additives in organic chemical reactions. Curr. Org. Chem. 2018, 22, 533–556. [Google Scholar] [CrossRef]
  2. Wang, B.; Qin, L.; Mu, T.; Xue, Z.; Gao, G. Are ionic liquids chemically stable? Chem. Rev. 2017, 117, 7113–7131. [Google Scholar] [CrossRef] [PubMed]
  3. Yue, C.; Fang, D.; Liu, L.; Yi, T.F. Synthesis and application of task-specific ionic liquids used as catalysts and/or solvents in organic unit reactions. J. Mol. Liq. 2011, 163, 99–121. [Google Scholar] [CrossRef]
  4. Amarasekara, A.S. Acidic ionic liquids. Chem. Rev. 2016, 116, 6133–6183. [Google Scholar] [CrossRef] [PubMed]
  5. Andreani, L.; Rocha, J.D. Use of ionic liquids in biodiesel production: A review. Braz. J. Chem. Eng. 2012, 29, 1–13. [Google Scholar] [CrossRef]
  6. Xin, H.; Wu, Q.; Han, M.; Wang, D.; Jin, Y. Alkylation of benzene with 1-dodecene in ionic liquids [Rmim]+Al2Cl6X(R = butyl, octyl and dodecyl; X = chlorine, bromine and iodine). Appl. Catal. A Gen. 2005, 292, 354–361. [Google Scholar] [CrossRef]
  7. Yi, Y.; Ding, F.; Li, W.; Ge, M.; Jin, G.; Gao, J. Alkylization of benzene with 1-dodecene catalyzed by Brønsted acidic ionic liquids. Electrochemistry 2009, 77, 591–593. [Google Scholar] [CrossRef]
  8. He, L.; Tao, G.H.; Liu, W.S.; Xiong, W.; Wang, T.; Kou, Y. One-pot synthesis of Lewis acidic ionic liquids for alkylation reaction. Chin. Chem. Lett. 2006, 17, 321–324. [Google Scholar]
  9. Mecerreyes, D. Applications of Ionic Liquids in Polymer Science and Technology; Springer: Berlin/Heidelberg, Germany, 2015; ISBN 978-3-662-44902-8. [Google Scholar]
  10. Bösmann, A.; Datsevich, L.; Jess, A.; Lauter, A.; Schmitz, C.; Wasserscheid, P. Deep desulfurization of diesel fuel by extraction with ionic liquids. Chem. Commun. 2001, 23, 2494–2495. [Google Scholar] [CrossRef]
  11. Huddleston, J.G.; Willauer, H.D.; Swatloski, R.P.; Visser, A.E.; Rogers, R.D. Room temperature ionic liquids as novel media for ‘clean’ liquid-iquid extraction. Chem. Commun. 1998, 1765–1766. [Google Scholar] [CrossRef]
  12. Itoh, T. Ionic liquids as tool to improve enzymatic organic synthesis. Chem. Rev. 2017, 117, 10567–10607. [Google Scholar] [CrossRef] [PubMed]
  13. Olivier-Bourbigou, H.; Magna, L.; Morvan, D. Ionic liquids and catalysis: Recent progress from knowledge to applications. Appl. Catal. A Gen. 2010, 373, 1–56. [Google Scholar] [CrossRef]
  14. Zhang, Q.; Zhang, S.; Deng, Y. Recent advances in ionic liquid catalysis. Green Chem. 2011, 13, 2619–2637. [Google Scholar] [CrossRef]
  15. Brown, L.C.; Hogg, J.M.; Swadźba-Kwaśny, M. Lewis Acidic ionic liquids. Top. Curr. Chem. 2017, 375, 1–40. [Google Scholar] [CrossRef]
  16. Estager, J.; Holbrey, J.D.; Swadźba-Kwaśny, M. Halometallate ionic liquids—Revisited. Chem. Soc. Rev. 2014, 43, 847–886. [Google Scholar] [CrossRef] [PubMed]
  17. Maase, M. Ionic Liquids in Synthesis; Wasserscheid, P., Welton, T.T., Eds.; Wiley-VCH Verlags GmbH & Co., KGaA: Weinheim, Germany, 2008; Volume 1, ISBN 3527312390. [Google Scholar]
  18. Zhang, R.; Meng, X.; Liu, Z.; Meng, J.; Xu, C. Isomerization of n-pentane catalyzed by acidic chloroaluminate ionic liquids. Ind. Eng. Chem. Res. 2008, 47, 8205–8210. [Google Scholar] [CrossRef]
  19. Huang, M.Y.; Wu, J.C.; Shieu, F.S.; Lin, J.J. Preparation of high energy fuel JP-10 by acidity-adjustable chloroaluminate ionic liquid catalyst. Fuel 2011, 90, 1012–1017. [Google Scholar] [CrossRef]
  20. Wang, L.; Zou, J.J.; Zhang, X.; Wang, L. Rearrangement of tetrahydrotricyclopentadiene using acidic ionic liquid: Synthesis of diamondoid fuel. Energy Fuel. 2011, 25, 1342–1347. [Google Scholar] [CrossRef]
  21. Chauvin, Y.; Hirschauer, A.; Olivier, H. Alkylation of isobutane with 2-butene using 1-butyl-3-methylimidazolium chloride-aluminium chloride molten salts as catalysts. J. Mol. Catal. 1994, 92, 155–165. [Google Scholar] [CrossRef]
  22. Zhang, J.; Huang, C.P.; Chen, B.H.; Ren, P.J.; Pu, M. Isobutane/2-butene alkylation catalyzed by chloroaluminate ionic liquids in the presence of aromatic additives. J. Catal. 2007, 249, 261–268. [Google Scholar] [CrossRef]
  23. Hu, P.; Wang, Y.; Meng, X.; Zhang, R.; Liu, H.; Xu, C.; Liu, Z. Isobutane alkylation with 2-butene catalyzed by amide-AlCl3-based ionic liquid analogues. Fuel 2017, 189, 203–209. [Google Scholar] [CrossRef]
  24. Pitula, S.; Mudring, A.V. Synthesis, structure, and physico-optical properties of manganate(II)-based ionic liquids. Chem. A Eur. J. 2010, 16, 3355–3365. [Google Scholar] [CrossRef] [PubMed]
  25. Cui, S.; Lu, B.; Xiao, X.; Han, Y.; Cai, Q. Synthesis of diphenylmethane from benzene and paraformaldehyde catalyzed by acidic ionic liquids. Catal. Lett. 2007, 119, 277–282. [Google Scholar] [CrossRef]
  26. Gunaratne, H.Q.N.; Lotz, T.J.; Seddon, K.R. Chloroindate(III) ionic liquids as catalysts for alkylation of phenols and catechol with alkenes. New J. Chem. 2010, 34, 1821–1824. [Google Scholar] [CrossRef]
  27. Gao, J.; Wang, J.Q.; Song, Q.W.; He, L.N. Iron(III)-based ionic liquid-catalyzed regioselective benzylation of arenes and heteroarenes. Green Chem. 2011, 13, 1182–1186. [Google Scholar] [CrossRef]
  28. Parshall, G.W. Catalysis in molten salt media. J. Am. Chem. Soc. 1972, 94, 8716–8719. [Google Scholar] [CrossRef]
  29. Liu, Y.; Hu, R.; Xu, C.; Su, H. Alkylation of isobutene with 2-butene using composite ionic liquid catalysts. Appl. Catal. A Gen. 2008, 346, 189–193. [Google Scholar] [CrossRef]
  30. Liu, Y.; Wang, L.; Li, R.; Hu, R. Reaction mechanism of ionic liquid catalyzed alkylation: Alkylation of 2-butene with deuterated isobutene. J. Mol. Catal. A Chem. 2016, 421, 29–36. [Google Scholar] [CrossRef]
  31. Zinurov, D.R.; Zinurov, R.R.; Akhmed’yanova, R.A.; Liakumovich, A.G. Skeletal isomerization of n-pentane in the presence of an AlCl3-based ionic liquid. Pet. Chem. 2010, 50, 376–380. [Google Scholar] [CrossRef]
  32. Stenzel, O.; Brüll, R.; Wahner, U.M.; Sanderson, R.D. Oligomerization of ole ns in a chloroaluminate ionic liquid. J. Mol. Catal. 2003, 192, 217–222. [Google Scholar] [CrossRef]
  33. Yang, S.; Liu, Z.; Meng, X.; Xu, C. Oligomerization of isobutene catalyzed by iron(III) chloride ionic liquids. Energy Fuel. 2009, 23, 70–73. [Google Scholar] [CrossRef]
  34. Chauvin, Y.; Gilbert, B.; Guibard, I. Catalytic dimerization of alkenes by nickel complexes in organochloroaluminate molten salts. J. Chem. Soc. Chem. Commun. 1990, 23, 1715–1716. [Google Scholar] [CrossRef]
  35. Gilbert, B.; Olivier-Bourbigou, H.; Favre, F. Chloroaluminate Ionic Liquids: From their Structural Properties to their Applications in Process Intensification. Oil Gas. Sci. Technol.-Rev. de l’IFP 2007, 62, 745–759. [Google Scholar] [CrossRef]
  36. Smith, G.P.; Dworkin, A.S.; Pagni, R.M.; Zingg, S.P. Bronsted superacidity of HC1 in a Liquid Chloroaluminate. A1C13-1-Ethyl-3-methyl-1H-imidazolium Chloride. J. Am. Chem. Soc. 1989, 111, 525–530. [Google Scholar] [CrossRef]
  37. Zawodzinski, T.A.; Osteryoung, R.A. Donor-Acceptor properties of ambient-temperature chloroaluminate melts. Inorg. Chem. 1989, 28, 1710–1715. [Google Scholar] [CrossRef]
  38. Koronaios, P.; King, D.; Osteryoung, R.A. Acidity of Neutral buffered 1-ethyl-3-methylimidazolium Chloride-AlCl3 Ambient-temperature molten salts. Inorg. Chem. 1998, 37, 2028–2032. [Google Scholar] [CrossRef]
  39. Estager, J.; Oliferenko, A.A.; Seddon, K.R.; Swadźba-Kwaśny, M. Chlorometallate(iii) ionic liquids as Lewis acidic catalysts—A quantitative study of acceptor properties. Dalton Trans. 2010, 39, 11375. [Google Scholar] [CrossRef] [PubMed]
  40. Taylor, A.W.; Men, S.; Clarke, C.J.; Licence, P. Acidity and basicity of halometallate-based ionic liquids from X-ray photoelectron spectroscopy. RSC Adv. 2013, 3, 9436. [Google Scholar] [CrossRef]
  41. Yang, Y.-L.; Kou, Y. Determination of the Lewis acidity of ionic liquids by means of an IR spectroscopic probe. Chem. Commun. 2004, 226–227. [Google Scholar] [CrossRef] [PubMed]
  42. Hu, P.; Zhang, R.; Meng, X.; Liu, H.; Xu, C.; Liu, Z. Structural and spectroscopic characterizations of amide-AlCl3-based ionic liquid analogues. Inorg. Chem. 2016, 55, 2374–2380. [Google Scholar] [CrossRef] [PubMed]
  43. Laurence, C.; Graton, J.; Gal, J.F. An overview of Lewis basicity and affinity scales. J. Chem. Educ. 2011, 88, 1651–1657. [Google Scholar] [CrossRef]
  44. Izgorodina, E.I.; Seeger, Z.L.; Scarborough, D.L.A.; Tan, S.Y.S. Quantum chemical methods for the prediction of energetic, physical, and spectroscopic properties of ionic liquids. Chem. Rev. 2017, 117, 6696–6754. [Google Scholar] [CrossRef] [PubMed]
  45. Contreras, R.; Aizman, A.; Tapia, R.A.; Cerda-Monje, A. Lewis molecular acidity of ionic liquids from empirical energy-density models. J. Phys. Chem. B 2013, 117, 1911–1920. [Google Scholar] [CrossRef] [PubMed]
  46. Cerda-Monje, A.; Ormazábal-Toledo, R.; Cárdenas, C.; Fuentealba, P.; Contreras, R. Regional electrophilic and nucleophilic Fukui functions efficiently highlight the Lewis acidic/basic regions in ionic liquids. J. Phys. Chem. B 2014, 118, 3696–3701. [Google Scholar] [CrossRef] [PubMed]
  47. Parveen, F.; Patra, T.; Upadhyayula, S. A structure-activity relationship study using DFT analysis of Bronsted-Lewis acidic ionic liquids and synergistic effect of dual acidity in one-pot conversion of glucose to value-added chemicals. New J. Chem. 2018, 42, 1423–1430. [Google Scholar] [CrossRef]
  48. Wu, W.; Lu, Y.; Ding, H.; Peng, C.; Liu, H. The acidity/basicity of metal-containing ionic liquids: Insights from surface analysis and the Fukui function. Phys. Chem. Chem. Phys. 2015, 17, 1339–1346. [Google Scholar] [CrossRef] [PubMed]
  49. Zhang, X.; Zhang, R.; Liu, H.; Meng, X.; Xu, C.; Liu, Z.; Klusener, P.A.A. quantitative characterization of Lewis acidity and activity of chloroaluminate ionic liquids. Ind. Eng. Chem. Res. 2016, 55, 11878–11886. [Google Scholar] [CrossRef]
  50. Wilkes, J.S.; Frye, J.S.; Reynolds, G.F. Aluminum-27 and carbon-13 NMR studies of aluminum chloride-dialkylimidazolium chloride molten salts. Inorg. Chem. 1983, 22, 3870–3872. [Google Scholar] [CrossRef]
  51. Turner, R.W.; Amma, E.L. Crystal and molecular structure of metal ion-aromatic complexes. I. The cuprous ion-benzene complex, C6H6·CuAlCl4. J. Am. Chem. Soc. 1963, 85, 4046–4047. [Google Scholar] [CrossRef]
  52. Roebuck, A.; Evering, B. Isobutane-olefin alkylation with inhibited aluminum chloride catalysts. Ind. Eng. Chem. Prod. Res. Dev. 1970, 9, 76–82. [Google Scholar] [CrossRef]
  53. Mantz, R.A.; Trulove, P.C.; Carlin, R.T.; Theim, T.L.; Osteryoung, R.A. Gutmann acceptor properties of LiCl, NaCl, and KCl buffered ambient-temperature chloroaluminate ionic liquids. Inorg. Chem. 1997, 36, 1227–1232. [Google Scholar] [CrossRef] [PubMed]
  54. Zheng, A.; Liu, S.; Deng, F. 31P NMR chemical shifts of phosphorus probes as reliable and practical acidity scales for solid and liquid catalysts. Chem. Rev. 2017, 117, 12475–12531. [Google Scholar] [CrossRef] [PubMed]
  55. Bui, T.L.T.; Korth, W.; Jess, A. Influence of acidity of modified chloroaluminate based ionic liquid catalysts on alkylation of iso-butene with butene-2. Catal. Commun. 2012, 25, 118–124. [Google Scholar] [CrossRef]
  56. Parr, R.G.; Parr, R.G.; Yang, W.; Yang, W. Density functional approach to the frontier-electron theory of chemical reactivity. J. Am. Chem. Soc. 1984, 106, 4049–4050. [Google Scholar] [CrossRef]
  57. Liu, Y.; Li, R.; Sun, H.; Hu, R. Effects of catalyst composition on the ionic liquid catalyzed isobutane/2-butene alkylation. J. Mol. Catal. A Chem. 2015, 398, 133–139. [Google Scholar] [CrossRef]
  58. Morell, C.; Grand, A.; Toro-Labbé, A. New dual descriptor for chemical reactivity. J. Phys. Chem. A 2005, 109, 205–212. [Google Scholar] [CrossRef] [PubMed]
  59. Zahn, S.; MacFarlane, D.R.; Izgorodina, E.I. Assessment of Kohn-Sham density functional theory and Moller-Plesset perturbation theory for ionic liquids. Phys. Chem. Chem. Phys. 2013, 15, 13664–13675. [Google Scholar] [CrossRef] [PubMed]
  60. Li, H.P.; Chang, Y.H.; Zhu, W.S.; Jiang, W.; Zhang, M.; Xia, J.X.; Yin, S.; Li, H.M. A DFT study of the extractive desulfurization mechanism by BMIM+AlCl4 ionic liquid. J. Phys. Chem. B 2015, 119, 5995–6009. [Google Scholar] [CrossRef] [PubMed]
  61. Grimme, S.; Hujo, W.; Kirchner, B. Performance of dispersion-corrected density functional theory for the interactions in ionic liquids. Phys. Chem. Chem. Phys. 2012, 14, 4875–4883. [Google Scholar] [CrossRef] [PubMed]
  62. Zhao, Y.; Truhlar, D.G. A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions. J. Chem. Phys. 2006, 125, 194101–194118. [Google Scholar] [CrossRef] [PubMed]
  63. Simeon, T.M.; Ratner, M.A.; Schatz, G.C. Nature of noncovalent interactions in catenane supramolecular complexes: Calibrating the MM3 force field with ab initio, DFT, and SAPT methods. J. Phys. Chem. A 2013, 117, 7918–7927. [Google Scholar] [CrossRef] [PubMed]
  64. Kendall, R.A.; Dunning, T.H.J.; Harrison, R.J. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 1992, 96, 6796–6806. [Google Scholar] [CrossRef]
  65. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09 (Revision E.01); Gaussian, Inc.: Wallingford, CT, USA, 2013. [Google Scholar]
  66. Rigby, J.; Izgorodina, E.I. Assessment of atomic partial charge schemes for polarisation and charge transfer effects in ionic liquids. Phys. Chem. Chem. Phys. 2013, 15, 1632–1646. [Google Scholar] [CrossRef] [PubMed]
  67. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef] [PubMed]
  68. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint. Chem. Rev. 1988, 88, 899–926. [Google Scholar] [CrossRef]
  69. Politzer, P.; Murray, J.S.; Bulat, F.A. Average local ionization energy: A review. J. Mol. Model. 2010, 16, 1731–1742. [Google Scholar] [CrossRef] [PubMed]
  70. Visual Molecular Dynamics. Available online: http://www.ks.uiuc.edu/Research/vmd/ (accessed on 1 July 2018).
  71. Boys, S.F.; Bernardi, F. The calculations of small molecular interaction by the difference of separate total energies. some procedures with reduced error. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
Sample Availability: Samples of the compounds are available from the authors.
Figure 1. Molecular surface properties of chloroaluminate ILs: (a) the electrostatic potential with a 0.001 a.u. contour of the electron density and location of the Vs,max (red arrow) and Vs,min (blue arrow) on the ILs; (b) the average local ionization energy with a 0.001 a.u. contour of the electron density and I ¯ s,min on the ILs (blue arrow).
Figure 1. Molecular surface properties of chloroaluminate ILs: (a) the electrostatic potential with a 0.001 a.u. contour of the electron density and location of the Vs,max (red arrow) and Vs,min (blue arrow) on the ILs; (b) the average local ionization energy with a 0.001 a.u. contour of the electron density and I ¯ s,min on the ILs (blue arrow).
Molecules 23 02516 g001
Figure 2. Molecular surface properties of mixed chlorometallate ILs: (a) electrostatic potential with a 0.001 a.u. contour of the electron density and the location of Vs,max (red arrow) and Vs,min (blue arrow) on the ILs; (b) average local ionization energy with a 0.001 a.u. contour of the electron density and I ¯ s,min on the ILs (blue arrow).
Figure 2. Molecular surface properties of mixed chlorometallate ILs: (a) electrostatic potential with a 0.001 a.u. contour of the electron density and the location of Vs,max (red arrow) and Vs,min (blue arrow) on the ILs; (b) average local ionization energy with a 0.001 a.u. contour of the electron density and I ¯ s,min on the ILs (blue arrow).
Molecules 23 02516 g002
Figure 3. Direction of the charge transfer when AlCl3 approaches the ion pair unit (C4mim-Cl).
Figure 3. Direction of the charge transfer when AlCl3 approaches the ion pair unit (C4mim-Cl).
Molecules 23 02516 g003
Figure 4. Electron donor Fukui function f atom electron-density isosurface plot at 0.004 a.u: (a) [(C2H5)3NH]+[CuAlCl5]; (b) [C4mim]+[CuAlCl5]; (c) [C4mimAg]+[AlCl4]; (d) [(C2H5)3NCu]+[AlCl4]; the blue color around Cu or Ag indicates that the transition metal atom is a favorite site for accepting the electrophilic reagent.
Figure 4. Electron donor Fukui function f atom electron-density isosurface plot at 0.004 a.u: (a) [(C2H5)3NH]+[CuAlCl5]; (b) [C4mim]+[CuAlCl5]; (c) [C4mimAg]+[AlCl4]; (d) [(C2H5)3NCu]+[AlCl4]; the blue color around Cu or Ag indicates that the transition metal atom is a favorite site for accepting the electrophilic reagent.
Molecules 23 02516 g004
Table 1. Cation-anion interaction energies and surface parameters of the Room-temperature ionic liquids (ILs) (kcal/mol).
Table 1. Cation-anion interaction energies and surface parameters of the Room-temperature ionic liquids (ILs) (kcal/mol).
[Cation][Anion]Vs,maxVs,min I ¯ s,min1E
[(C2H5)3NH][AlCl4]42.430−33.70610.224−90.361
[(C2H5)3NH][Al2Cl7]52.596−30.14510.480−96.637
[(C2H5)3NH][CuCl2]36.104−41.1208.730−96.637
[(C2H5)3NH][ZnCl3]41.400−39.7439.423−92.871
[C4mim][Cl]40.53−59.295.81−90.34
[C4mim][AlCl4]48.275−36.01010.048−84.086
[C4mim][Al2Cl7]52.596−30.14510.480−77.184
[C4mim][CuCl2]42.003−45.6178.337−84.086
[C4mim][Cu2Cl3]43.950−40.1949.324−56.476
[C4mim][ZnCl3]42.585−33.9189.710−85.969
[C4mim][Zn2Cl5]44.895−37.3929.834−75.765
1 The unit of I ¯ s,min is eV.
Table 2. Cation-anion interaction energies and the surface parameters of the ILs (kcal/mol).
Table 2. Cation-anion interaction energies and the surface parameters of the ILs (kcal/mol).
[Cation][Anion]Vs,maxVs,min I ¯ s,min1E
[(C2H5)3NH][CuAlCl5]43.261−37.119.261−88.479
[(C2H5)3NH][AgAlCl5]43.201−37.3369.243−86.854
[(C2H5)3NCu][AlCl4]27.965−23.68210.330−82.992
[(C2H5)3NCu][Al2Cl7]28.230−20.2299.286−81.109
[(C2H5)3NAg][Al2Cl7]26.725−15.9959.692−80.579
[(C2H5)3NH][ZnAlCl6]45.011−35.5729.677−85.969
[C4mim][CuAlCl5]43.445−36.6039.218−84.086
[C4mim][AgAlCl5]45.866−37.4709.257−83.952
[benzene][CuAlCl4]35.085−20.59610.891−26.982
[benzene][AlCl4]33.680−28.7609.724−28.987
[ether][CuAlCl4]35.305−23.54010.871−15.060
[ether][AlCl4]40.184−24.1257.182−7.530
1 The unit of I ¯ s,min is eV.
Table 3. The magnitude of the charge transfer from the unit (C4mim-Cl) to the metal chloride [QCT(IP→MClx)], and the restrained electrostatic potentials (RESP) charges for the ILs under study 1.
Table 3. The magnitude of the charge transfer from the unit (C4mim-Cl) to the metal chloride [QCT(IP→MClx)], and the restrained electrostatic potentials (RESP) charges for the ILs under study 1.
[Cation][Anion]QCTN1C2H(C2)N3C4H(C4)C5
[C4mim][Cl]−0.0450.0080.2730.043−0.1420.225−0.226
[C4mim][AlCl4]−0.407(−0.419)0.136−0.1270.2390.165−0.2230.261−0.219
[C4mim][CuAlCl5]−0.345(−0.378)0.171−0.0730.2440.149−0.2480.305−0.256
[C4mim][ZnCl3]−0.297(−0.267)0.093−0.1150.2340.098−0.0780.253−0.187
[C4mim][Zn2Cl5]−0.334(−0.368)0.087−0.1140.2110.123−0.0950.256−0.198
1 All values are given in a.u. The atomic labels are shown in Figure 3. The values in parentheses are the amount of charge transfer computed with the Natural Population Analysis (NPA) scheme.
Table 4. Electrophilic ( N Ω + ) and nucleophilic ( N Ω ) Fukui function values over the cations and anions for the studied ILs (eV).
Table 4. Electrophilic ( N Ω + ) and nucleophilic ( N Ω ) Fukui function values over the cations and anions for the studied ILs (eV).
[Cation][Anion]CationAnionMin. ∆fatomMax. ∆fatomAcidity/Basicity
Distribution
N Ω + N Ω N Ω + N Ω
[(C2H5)3NH][AlCl4]0.9300.0810.0630.918ClHNormal
[(C2H5)3NH][ZnCl3]0.9170.0520.9480.065ClHNormal
[(C2H5)3NH][CuAlCl5]0.2710.0810.7020.918CuAlBifunctional
[(C2H5)3NCu][AlCl4]0.8620.8150.1310.185CuAlBifunctional
[(C2H5)3NH][AgAlCl5]0.0750.0680.8790.932AgClBifunctional
[C4mim][AlCl4]0.9140.0840.0830.916ClCNormal
[C4mim][CuCl2]0.9250.0730.0920.908CuCNormal
[C4mim][CuAlCl5]0.2730.0790.7190.921CuAlBifunctional
[C4mimAg][AlCl4]0.8160.8040.1170.196AgHBifunctional
[(C2H5)3NH][ZnAlCl6]0.6780.0570.3130.943ZnAlBifunctional 1
[Benzene][CuAlCl4]0.1660.1090.9130.924ClCuBifunctional
[Ether][CuAlCl4]0.4230.1560.5750.844CuAlBifunctional
1 [(C2H5)3NH][ZnAlCl6] is only close to the standard of bifunctional distribution for anion: N Ω + > 0.5 eV and N Ω > 0.8 eV.

Share and Cite

MDPI and ACS Style

Liu, Y.; Wang, J. Lewis Acidity and Basicity of Mixed Chlorometallate Ionic Liquids: Investigations from Surface Analysis and Fukui Function. Molecules 2018, 23, 2516. https://doi.org/10.3390/molecules23102516

AMA Style

Liu Y, Wang J. Lewis Acidity and Basicity of Mixed Chlorometallate Ionic Liquids: Investigations from Surface Analysis and Fukui Function. Molecules. 2018; 23(10):2516. https://doi.org/10.3390/molecules23102516

Chicago/Turabian Style

Liu, Ying, and Juanfang Wang. 2018. "Lewis Acidity and Basicity of Mixed Chlorometallate Ionic Liquids: Investigations from Surface Analysis and Fukui Function" Molecules 23, no. 10: 2516. https://doi.org/10.3390/molecules23102516

APA Style

Liu, Y., & Wang, J. (2018). Lewis Acidity and Basicity of Mixed Chlorometallate Ionic Liquids: Investigations from Surface Analysis and Fukui Function. Molecules, 23(10), 2516. https://doi.org/10.3390/molecules23102516

Article Metrics

Back to TopTop