Adsorption of choline benzoate ionic liquid on graphene , silicene , germanene and boron-nitride nanosheets : a DFT perspective

The adsorption of choline benzoate ([CH][BE]) ionic liquid (IL) on the surface of different hexagonal nanosheets has been studied using Density Functional Theory (DFT) methods. For this, the interaction mechanism, binding energies and electronic structure of [CH][BE] ionic liquid on four types of nanosheets, i.e., graphene, silicene, germanene and boron-nitride, were estimated and compared. The adsorption of [CH][BE] ionic liquid on different nanosheets is mainly featured by van der Waals forces, leading to strong benzoate ion–surface p-stacking. Likewise, there is also an important charge transfer from the anion to the sheet. The electronic structure analysis shows that Siand Ge-based sheets lead to the largest changes in the HOMO and LUMO levels of choline benzoate. This paper provides new insights into the capability of DFT methods to provide useful information about the adsorption of ionic liquids on nanosheets and how ionic liquid features could be tuned through the adsorption on the suitable nanosheet.


Introduction
Graphene could be defined as a single layer of carbon atoms with sp 2 hybridization arranged in a two dimensional (2D) honeycomb structure. 1,2 Since its discovery in 2004, 3 a new area of research involving materials science, physics, chemistry and biotechnology has emerged. 2,4 Graphene is attractive because it shows many extraordinary characteristics, which are a direct consequence of its unique atomic structure, such as intrinsic zero-gap semiconductor character, high carrier mobility, high surface area, excellent thermal conductivity or biocompatibility. 2,4,5 Because of its remarkable properties, graphene-based applications have been reported in a wide range of areas, such as optoelectronic devices, energy generation, chemical sensors, or hybrid materials. 1,2,4 Likewise, the potential applications of graphene can be extended and task-specifically tuned through the functionalization (or chemical adsorption) of molecules and nanoparticles. 6 The mechanism of adsorption of relevant molecules, and also the changes in graphene properties upon molecular adsorption have been widely studied using both experimental and computational approaches. [6][7][8] A wide variety of molecules could be absorbed on the graphene surface, and a very relevant field of study is the interaction of ionic liquids (ILs) with carbon nanostructures and their behavior with regard to the carbon nanostructure, because of its importance both for basic science and for applied purposes. [8][9][10][11][12] ILs are emerging as an attractive alternative to conventional organic solvents due to their special chemical and physical properties: negligible vapor pressures, high thermal and chemical stability, non-flammability or good solvent capacity for a wide range of organic, inorganic, polymeric and organometallic compounds, and the possibility of designing task-specific ILs through the suitable combination of cations and anions. 13 Therefore, ILs have been considered for applications in different technological fields such as lubrication, solvent extraction, catalytic processes or electrochemical applications. 14 In addition, the behavior of ILs with regard to graphene surfaces has also been studied, using both experimental and theoretical approaches, which have led to propose relevant applications in several technological fields, such as full cells, supercapacitors, solar cells or storage devices. 8,10,12,15 As a matter of fact, both theoretical and experimental studies have been conducted to elucidate the mechanism of interaction between ILs and graphene. 8,10,11,[15][16][17] Nevertheless, most available studies are limited to a reduced number of classic types of ILs, such as imidazole based ionic liquids.
Beyond the graphene sheets, other hexagonal 2D nanosheets such as silicene, germanene and boron-nitride analogues have also attracted increasing attention. 1,18 Silicene and germanene are the graphene counterparts based on silicon and germanium, which are also elements of group IV. These materials also possess a honeycomb structure similar to graphene, and thus, they also show the most outstanding properties of graphene such as high carrier mobility, ferromagnetism or zero band gap. 1,18 Similarly, Si and Ge based nanosheets are expected to offer an alternative for the enhancement of the performance and scalability of the traditional silicon-based devices. 1,19 Nonetheless, some remarkable differences should be noted between C and Si/Ge based nanosheets. The buckled honeycomb structure of silicene and germanene leads to a significantly high chemical reactivity than graphene, and thus, to a much stronger trend for adsorption of atoms and molecules. 1,19,20 In silicene/germanene sheets, silicon/germanium shows sp 3 hybridization, which is the most favorable configuration in comparison with the sp 2 or the mixed sp 2 -sp 3 orbitals. This sp 3 hybridization leads to common covalent Si-Si/Ge-Ge bonds in a low buckled structure. 1,19,21 2D heterostructures, such as boron nitride (BN), are composed of an equal number of alternating boron and nitrogen atoms with sp 2 hybridization in a honeycomb arrangement, with a similar structure to graphene. BN sheets are intrinsically insulators (or wide gap semiconductors), stable up to 1000 K and more resistant to oxidation than graphene. 22 All these properties justify the great interest on Si, Ge and BN nanosheets. As for graphene, properties of different Si, Ge and BN sheets can be tuned through the adsorption of different molecules. This has led to several studies in the last few years dealing with the surface functionalization of Si, 1,18,19,[23][24][25]19,20 or BN 1,22,[26][27][28][29] nanosheets. These studies point out that a new hybrid system based on hexagonal nanosheets with concrete features could be designed through the adsorption of adequate molecules on the surface. In spite of the wide interest on 2D surface functionalization through molecular adsorption, there is still scarce information about the adsorption of different molecules on silicene, germanene and BN sheets. Experimental techniques, such as atomic force microscopy, scanning tunneling microscopy or X-ray spectroscopy, as well as theoretical methods like ab initio quantum chemistry methods (mainly Density Functional Theory, DFT), molecular dynamics (MD) and Monte Carlo (MC) simulations have been applied to study the ionic-liquid-graphene interface. 8,10,11,[15][16][17] As previously noted, those studies are mainly limited to classic ILs. Few studies on the adsorption of ILs on BN nanosheets have been also reported, 29,30 while studies on the adsorption of IL on silicene or germanene surfaces have been not published.
Although DFT simulations are limited to relatively simple theoretical models, they have proven their ability in providing useful information about the design of functionalized nanosheets. [6][7][8] In addition, DFT methods provide a deeper knowledge of the interaction of ILs at the surface, explaining their electronic structure and shedding light on the interaction mechanism. Therefore, a DFT study on the adsorption of ILs on C, Si, Ge and BN based hexagonal nanosheets is reported in this work. Concretely, the adsorption of choline benzoate ([CH][BE], Fig. 1) ionic liquid was studied here. Choline ([CH]) based ionic liquids are a new family of ionic liquids 31 with suitable properties such as null toxicity, high biodegradability and low cost. Likewise, the combination of a choline cation with a benzoate anion ([BE]) leads to an ionic liquid composed completely of biomaterials, 32 which can be produced at very low cost, 33 and with null environmental impact. Thus, [CH][BE] ionic liquid has been selected in this work as a first example in the study of the adsorption of ILs on different nanosheets.

Theoretical methodology
DFT simulations were carried out using the generalized gradient approximation of Perdew, Burke and Ernzerhof (PBE) 34 as implemented in the SIESTA 3.2 package, 35 along with norm-conserving Troullier-Martins pseudopotentials 36 and numerical double-z polarized (DZP) basis sets. All calculations were done using an energy mesh cut-off of 400 Ry and a k-point mesh of 8 Â 8 Â 8 in the Monkhorst-Pack scheme. 37 Model nanosheets were described to have a 2D honeycomb structure of 128 atoms under periodic boundary conditions (PBC) applied in the three space directions. PBC in the direction normal to the nanosheet surface was 30 Å, long enough to avoid interaction between image layers. Structural relaxation was done by conjugate gradients, with convergence criteria that forces acting on all atoms do not exceed 0.03 eV Å À1 .
The PBE functional has been successfully applied to study interaction between different molecules and C, 38,39 Si, 24,40 Ge 41 or BN 22,28 nanosheets. Long range dispersion interactions are expected to be important for an adequate description of ILnanosheet systems under study. In this sense, it is well known that the most common GGA (generalized gradient approximation) functionals have shortcomings of the adequate description of long range dispersion interactions. 42 Some reported results have shown that dispersion corrections are not required in some cases for which molecules develop van der Waals interactions. 7,39 On the other hand, many studies have pointed out the need for the adequate description of dispersive interactions. 7,20,25,26,43 Therefore, dispersion corrections, according to Grimme's scheme, 44 over the PBE functional (PBE-D2) were also employed to obtain information about the effect of dispersion corrections on calculated properties.
[CH][BE]-nanosheet models were built by (randomly) placing the geometry optimized ionic liquid onto the nanosheet. In this starting geometry, the IL was placed to allow a p-stacking interaction between the phenyl motif and the surface with a distance of around 3.5 Å (typical p-stacking distance). These starting geometries were optimized. Then, the relative disposition (through longitudinal displacements) over the surface as well as the rotational angle through the main bond between both ions (see d 1 , Fig. 1) holding the above mentioned p-stacking between the phenyl motif and the surface, have been assessed. These analyses were carried out based on previously optimized geometries through single point calculations. From previous single points, choline benzoate ionic liquid was placed at the most stable disposition over each nanosheet and these [CH][BE]nanosheet models were optimized, which led to the optimized structures discussed here.
Binding energies, BEs, between [CH][BE] and the corresponding nanosheets were estimated as where E S , E IL , and E IL-S stand for the total energy of the sheet, ionic pair and total energy of the IL-nanosheet, respectively. BE corresponding to the interaction between ions is defined as where E cat and E ani are the total energy of the cation and the anion, respectively. Charge transfer (CT) between graphene and IL has been studied through Hirshfeld 45 and ChelpG 46 atomic charges. The ChelpG scheme has proven to be adequate for describing atomic charges in ILs, 47,48 while the Hirshfeld model is an alternative definition of atomic charges based on partitioning of the electron density. 45 Interactions between the IL and the nanosheets were featured through both a topological analysis of the electron density, according Bader's 49 theory (Atoms in Molecules, AIM), and the analysis of the reduced density gradient (RGD) at low densities. 50 According to AIM theory, 49 there are four kinds of critical points, but because of the characteristics of the studied systems, and for improving and clarifying data analysis, attention has been paid to bond critical points (BCPs), which improves the criteria for considering the presence of intermolecular interactions, through the computed electronic density (r) and its Laplacian (r 2 r). Nevertheless, some information regarding the adsorption process could also be obtained from ring and cage critical points (RCP and CCP, respectively). RGD analysis is able to find non-covalent interactions based on the peaks that appear at low densities. Thus, the visualization of RGD iso-surfaces for these peaks allows the localization of weak interactions. The strength and nature of the interactions are determined through the sign of the second density Hessian eigenvalue. 50 AIM and RGD analyses were carried out using the MultiWFN code. 51 Partial density of states (PDOS) was extracted using the program developed by A. Postnikov. 52 Regarding [CH][BE] on top of the graphene sheet, intermolecular interactions d 1 -d 4 are not dramatically affected during the adsorption process. For example, the average bond variation after the adsorption process is 0.026 Å/0.070 Å using the PBE/PBE-D2 functional. These small bond variations point out that the main features of the ionic liquid are also present for the ionic liquid on top of the nanosheet. A similar trend was found for atomic charges, which are similar to those estimated in the absence of the graphene surface. At the PBE-D2/DZP level and by using the ChelpG scheme, a cation/anion loses its positive/negative charge (0.07/0.13 e À ), while the graphene surface has a charge of À0.06 e À . The negative sign of graphene charge indicates a small charge transfer from the ionic liquid (concretely from the anion) to the surface. The remaining anionic lost charge stands for the increase noted in the cationic charge. Thus the adsorption of IL on graphene leads to a charge transfer between ions of 0.34 e À according the ChelpG scheme at the PBE/DZP level. In general terms, similar conclusions are obtained for the interaction between both ions in presence of boron-nitride sheets. The Hirshfeld model predicts a higher charge transfer between the anion and the BN-based surface (0.26 e À ).

Choline benzoate on nanosheets
The adsorption of choline benzoate on silicene leads to a proton transfer between both ions (Fig. 3). As can be seen in Table 3, Table 1 Main molecular parameters for the most stables structures optimized at PBE/DZP and PBE-D2/DZP levels for [CH][BE] ionic liquid. AIM parameters related with intermolecular interactions (electronic density, r, and its lLplacian, r 2 r) as well as the total charge over choline (q + ) and benzoate (q À ) ions computed according to ChelpG/Hirshfeld schemes are also collected. See Fig. 1  while there is a decrease in the inter-ionic CT. The PBE-D2 functional also predicts a proton transfer between both ions. However, there is a strengthening of d 0 interaction (which decreases/increase its length/electronic density) in comparison with the optimized structure at the PBE level. Both ChelpG and Hirshfeld models yield smaller inter-ionic CT (of around 0.30 e À ) than isolated IL, while CT to the silicene surface increases up to 0.65 e À .
The optimized structure of [CH][BE] on the Ge based surface at the PBE/DZP level is shown in Fig. 4. As can be seen in Table 4, d 1 suffers a shortening (0.261 Å) with an increase in electronic density due to the adsorption process, which results in the strengthening of the interaction between both ions. In addition, the relative disposition of choline hinders the interaction d 4 (the corresponding the BCP was not found). According to the Hirshfeld model, cation/anion charge becomes more [Be]-Be) system along the structural parameters related to intermolecular interactions (left) and RGD iso-surfaces (right), whose green colour indicates van der Waals interactions. Atom colour code: (green) carbon, (red) oxygen, (blue) nitrogen, (white) hydrogen and (yellow) germanium. Red and purple points stand for RCP and CCP, respectively, related to intermolecular interactions. BCPs were omitted for simplicity.  Table 4).
In short, the most dramatic changes during the adsorption process are the proton transfer in the IL-Si system (also in IL-Ge at the PBE-D2/DZP level), and CT computed according to the Hirshfeld model. Isolated choline benzoate, IL-graphene and IL-BN systems yield an inter-ionic CT of B0.50 e À , while this CT is equal to 0.28 e À /0.19 e À for the IL on the silicene/ germanene surface (at the PBE/DZP level). In addition, due to the adsorption process, the new arrangement between both ions leads to a new cage critical point (labeled as CCP1, Fig. 2

Adsorption mechanism: [CH][BE]-nanosheet interactions
As concerns the changes in the structure of studied nanosheets during the adsorption process, the Dd X-X parameter (see Tables 2-5) has been defined as the difference value of X-X after and before the adsorption process (X-X = C-C, Si-Si, Ge-Ge or B-N bonds).
The studied IL-nanosheet systems and both theoretical levels (without dispersion corrections or dispersion corrected functional), Dd X-X o 0.009 Å, point out that the nanosheet structure does not suffer important change upon the ionic liquid adsorption. In addition, graphene and boron-nitride sheets retain their planarity after the adsorption process, while silicene and germanene sheets also maintain their buckled structure.
Results reported in Fig. 2-5 show that the phenyl (benzoate) motif adopts a configuration parallel to the surface. The average interplanar distances between the phenyl moiety and different sheets (labeled as d eq , Tables 2-5) were found to be of roughly 3.65 Å for IL-graphene and IL-silicene systems and B3.74 Å for IL-germanene and IL-boron nitride at the PBE/DZP level. For [CH][BE] on the graphene sheet, d eq computed at the PBE/DZP level lies close to the typical p-stacking distance (3.5-4.0 Å). However, the PBE-D2/DZP level brings an approach (B0.50 Å) between both planes.
The adsorption of the phenyl motif on graphene sets up three interactions with the surface (d 5 -d 7 ), whose average distance is 3.625 Å (at PBE/DZP level). Besides, the carboxylate group also develops a link with the surface (d 8 ). The dihedral angle between phenyl and carboxylate motifs (t 1 ) is 0.0 degrees, which hinders the interaction between the carboxylate group [BE]-graphene system, AIM parameters related with intermolecular interactions (electronic density, r, and its laplacian, r 2 r) as well as the total charge over choline (q + ), benzoate (q À ) ions and sheet (q s ) computed according to ChelpG/ Hirshfeld schemes are also collected. See Fig. 2 represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. b P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. c P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. d Difference value of C-C bonds corresponding to nanosheets after and before of the adsorption process. [BE]-silicene system, AIM parameters related to intermolecular interactions (electronic density, r, and its Laplacian, r 2 r) as well as the total charge over choline (q + ), benzoate (q À ) ions and sheet (q s ) computed according to ChelpG/ Hirshfeld schemes are also collected. See Fig. 3 represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. b P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. c P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. d Difference value of Si-Si bonds corresponding to nanosheets after and before the adsorption process.
(which keeps most of the negative charge) and the surface. Hence, the CT (0.02 e À according to the Hirshfeld scheme at the PBE/DZP level) from the anion up to graphene sheets is close to zero. Silicene/germanene surfaces establish 3/2 intermolecular bonds with the phenyl motifs (labeled as d 5 -d 7 /d 5 , d 6 ), whose distances are larger than those found for the IL-graphene system (Fig. 3, 4 and Tables 3, 4). Nonetheless, all of them yield larger electronic density values, which would point out to stronger interactions between the phenyl motif and the surface. Due to larger diameters of those rings forming the Si-Ge surface, bonds between the surface and the phenyl are not perpendicular to the surface plane. In this case, t 1 = 15.3 degrees/16.7 degrees, which allows a shortening of the distance between the carboxylate group and the surface. Two/one interactions were found between the carboxylate group and the Si/Ge sheet (d 8 , d 9 /d 7 ). The CT (0.54 e À /0.68 e À according to the Hirshfeld scheme at the PBE/DZP level) between the anion and the germanene surface agrees with the shortest distance between the carboxyl group and the surface (2.318 Å at the PBE/DZP level). Likewise, the benzoate anion is able to develop five interactions with the boron-nitride surface (Fig. 5), four of them through the phenyl motif (d 5 -d 8 ) and the other one through the COO À group (d 9 ). For [CH][BE] on top of the BN sheet, t 1 is close to the planarity. This dihedral angle obstructs an approach (d 9 = 3.900 Å at PBE/DZP level) between the COO À group of benzoate and the BN surface, leading to small CT from the anion to the BN sheet (0.11 e À according to the Hirshfeld scheme at the PBE/DZP level).
[CH]-surface interactions are always carried out through two intermolecular bonds between the oxygen atom and the methylene group on the surface. Moreover, the germanene sheet allows another additional interaction with the methylene group (see Fig. 4). The shortest distance between [CH] and the surface is found for the IL-silicene system (d 10 = 1.841 Å at the PBE/DZP level, Table 3). The presence of choline benzoate on the silicene surface leads to a deformation of the silicene sheet, which allows a shorter length for d 10 (d 10 = 1.841 Å for the IL-Si system at PBE/DZP). A different behavior is found with regard to the interaction strength (based on the distance and electronic density values) between the oxygen atom and the methylene group on the surface (d 10 = 2.762 Å/d 11 = 2.867 Å/d 9 = 3.032 Å and d 10 = 3.056 Å/d 11 = 2.66 Å for IL-C/IL-Si/IL-Ge/IL-BN systems at PBE/DZP).
In addition to the information related with intermolecular IL-sheet interactions (labeled as d 5 -d 11 ) and the AIM properties for their corresponding BCPs, main AIM parameters for RCPs Table 4 Main molecular parameters for the most stables structures optimized at PBE/DZP and PBE-D2/DZP levels. For the [CH][BE]-germanene system, AIM parameters related to intermolecular interactions (electronic density, r, and its Laplacian, r 2 r) as well as the total charge over choline (q + ), benzoate (q À ) ions and the sheet (q s ) computed according to Hirshfeld schemes are also collected. See Fig. 4  a P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. b P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. c P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related to intermolecular interactions between the IL and the sheet. d ChelpG atomic charges were not computed since Breneman radii (which is needed to compute ChelpG charges) for the Ge atom is not available in the literature. e Difference value of Ge-Ge bonds corresponding to nanosheets after and before the adsorption process. [BE]-boronnitride system, AIM parameters related to intermolecular interactions (electronic density, r, and its laplacian, r 2 r) as well as total charge over choline (q + ), benzoate (q À ) ions and sheet (q s ) computed according to ChelpG/Hirshfeld schemes are also collected. See Fig. 5 represents the sum of r or r 2 r for those BCP/RCP/CCP related with intermolecular interactions between the IL and the sheet. b P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related with intermolecular interactions between the IL and the sheet. c P (BCP)/ P (RCP)/ P (CCP) represents the sum of r or r 2 r for those BCP/RCP/CCP related with intermolecular interactions between the IL and the sheet. d Difference value of B-N bonds corresponding to nanosheets after and before the adsorption process. and CCPs (red and purple points in Fig. 2-5) have also been collected. With the aim of looking for a relationship between the interaction strength for IL-nanosheet systems and AIM parameters, the total value of electronic density of the different calculated critical points is also gathered in Tables 2-5. Some relationships were found between AIM parameters for intermolecular interaction and binding energies. For example, Si and Ge based surfaces provide the largest value for the total electronic density for the BCPs related to the interactions between the ionic liquid and the surface ( P (BCP)). Silicene and germanene sheets allow stronger interactions than graphene or boron-nitride surfaces.
As seen for d eq , the PBE-D2 functional estimates, in general, shorter intermolecular distances between the IL and the surfaces, which leads to an increase of some factors related with the intermolecular interactions, such as interaction strength (based on electronic density value) and IL-surface charge transfer.
The nature of the interaction between [CH][BE] ionic liquid and different nanosheets was elucidated based on RGD isosurfaces. For [CH][BE]-nanosheet systems, the RGD iso-surfaces are displayed in Fig. 2-5. The green color for the regions between the IL and the nanosheets points out that van der Waals interactions are the main forces responsible for the IL adsorption on nanosheets regardless of the considered nanosheet. The largest iso-surface between the phenyl motif and the BN surface is in concordance with the parallel disposition between the phenyl motif and the BN sheet. In short, anion-surface interactions should mainly contribute to the IL-surface binding energies.
As mentioned in previous sections, the available literature dealing with the adsorption of IL on nanosheets is very scarce (mainly limited to imidazole based ILs on graphene). Recently, Herrera et al. previously studied the adsorption of imidazole cation on the amino acid based anion graphene surface, 8 while Shakourian-Fard et al. analyzed IL adsorption on BN nanosheets for classical ILs, 29 both of them through DFT methods. The first employed functional was PBE-D2/DZP, while a meta-GGA functional was selected later. Both approximations should provide a good description for the dispersion forces. 44,53 However, despite the IL (choline benzoate, aminoacid based or classical ILs), the surface (graphene or boron-nitride) or the applied methods (DFT-D2 or meta-GGA), similar structural features were inferred, i.e., the IL-surface distance of around 3.0 Å was found; IL-surface interactions were mainly featured by p-stacking interactions between aromatic ions (the imidazole cation for previously reported studies or the benzoate anion in this paper) and the surface.

Binding energies
Main intermolecular interactions between both ions as well as between the ionic liquid and the nanosheets were analyzed in previous sections. In addition, binding energies can be used as a measurement of the interaction strength. The binding energy of the interaction between both ions (BE IL , Fig. 6) was 84.31 kcal mol À1 at the PBE/DZP level. This high value is mainly due to the important contribution from the Coulombic forces between both ions. Fig. 6 also shows BE IL of the IL on top of the nanosheets, which was calculated from the previously optimized IL geometry on top of each nanosheets. As seen in Fig. 6, BE IL decreases from IL on top of graphene up to IL on top of germanene, while [CH][BE] on top of BN yields higher BE IL . IL adsorption on graphene, silicene, or germanene leads to a BE IL of up to 16.34 kcal mol À1 , while choline benzoate on top of boron-nitride yields a BE IL of 83.94 kcal mol À1 . Hence, the adsorption of choline benzoate on the studied homo-nanosheets has important drawbacks on BE IL . Albeit there is an important proton transfer from choline up to benzoate for IL-silicene systems; this fact would not have important effects on the interaction between both ions on the silicene surface. BE IL at the PBE-D2/DZP level (which is somewhat higher, mainly due to smaller intermolecular distances) follows a similar pattern during the adsorption process. Coulombic forces between ions play an important role in the interionic interactions. In fact, [CH][BE] on silicene and germanene provides the lowest interionic CT and the lowest BE IL . As said, intermolecular interactions (such as hydrogen bonds) are also an important contribution to BE IL . In this sense, Fig. 6 also displays dispersion energy contribution (BE dis,IL ) to the total binding energy (quantified according to Grimme's approach 44 for the PBE-D2 functional). BE dis,IL , which is very close to the difference between the binding energy between both (dispersion) corrected and non-corrected functionals, only provides a small contribution to the total BE IL , Fig. 6. Fig. 7 shows binding energies of the interaction between the ionic liquid and different nanosheets (BE IL-S ). At the PBE/DZP level BE IL-S follows as: graphene C boron-nitride o germanene o silicene. The adsorption of [CH][BE] on top of graphene and boron-nitride yields BE IL-S B 12.24 kcal mol À1 , while the values increase up to 19.41 kcal mol À1 and 34.49 kcal mol À1 for Ge and Si based sheets, respectively. As previously noted, interactions between choline benzoate and silicene lead to surface distortion, which allows a short intermolecular distance between the cation and the surface (d 8 , see Fig. 3). However, this effect is not displayed in the IL-germanene system. At the PBE-D2/DZP level, computed BE IL-S values follow a similar pattern, however both IL-silicene and IL-germanene systems yield similar BE IL-S (C49.65 kcal mol À1 ). Fig. 7 also draws the contribution from BE dis,IL-S to the total BE IL-S estimated using the PBE-D2 functional. Note that BE dis,IL-S C 28.67 kcal mol À1 , except for IL-germanene systems, whose BE dis,IL-S = 39.43 kcal mol À1 . As seen in Tables 2-5, most intermolecular IL-germanene distances are, in general, shorter than those obtained for the remaining systems, which would lead to an enhancement of the dispersion contribution. These high BE dis,IL-S values prove that dispersion interactions are the main driving forces between the selected IL and the surfaces. Those results agree with iso-surface features described previously. In general, the highest values obtained for IL-silicene and IL-germanene systems would be due to their buckled honeycomb structure, which allow a higher chemical reactivity than graphene, leading to a much stronger adsorption of atoms and molecules. 1,19,20 Previously reported studies estimated that BE IL-S for different ionic liquid on top of graphene 8 or boron-nitride 29 sheets were roughly 63.00 kcal mol À1 (at the PBE-D2/DZP level) and 19.74 kcal mol À1 (PBE-D3/TZP), respectively. In these studies, nanosheets mainly tends to interact with the cation. However in this work, the larger contribution to the BE IL-S comes from anion-sheet interactions. ionic liquids, as well as the contributions from both ions are also reported in Fig. 8. There are several studies dealing with the electronic structure of pristine nanosheets, 1,18,19,40,54 and the results reported in Fig. 8 are in agreement with the available data, which show zero gap semiconductor characteristics for all cases, except for boron-nitride nanosheets which show the insulator character (wide gap). The adsorption of choline benzoate ionic liquid on graphene, silicene or germanene does not have a dramatic effect on the nanosheet electronic structure, Fig. 8, with only a slight increase the energies of both valence/conductivity bands (B0.20 eV). For IL on silicene and germanene, the highest/ lowest occupied/unoccupied orbitals corresponding to the ionic liquid (i.e. the HOMO and the LUMO over the ionic liquid motif) Fig. 7 Binding energies for the interaction between the IL and the nanosheets (BE IL-S ) at PBE/DZP (black) and PBE-D2/DZP (red + blue) theoretical levels, wherein blue bar stands for BE dis contributions. are far from the valence and conduction bands of the ILnanosheet systems. The same features are also noted for the LUMO over IL for IL-graphene systems, while the HOMO is closer in energy to the valence band. Therefore, IL adsorption does not change its zero band gap character. With regards to the boron-nitride nanosheet, the presence of IL leads to a decrease of its band gap, from 3.93 eV up to 3.45 eV (at the PBE/DZP level). After the adsorption process, the highest occupied level is due to the IL, which is higher in energy that the boron-nitride valence band. In addition, interaction with the IL leads to conductivity band destabilization, which yields energies near to the LUMO.

Electronic structure
From PDOS plotted in Fig. 8 information about cation and anion contributions to their corresponding HOMOs and LUMOs is inferred. The same information is easily obtained from the molecular orbital contours shown in Fig. 9. The LUMO of the isolated ionic liquid is mainly localized over the benzoate anion. The same behavior is also found for [CH][BE] on top of the selected nanosheets. The adsorption on graphene/boron-nitride nanosheets leads to a decrease in the LUMO energy of 0.36 eV/ 0.48 eV, while the adsorption on silicene/germanene leads to a decrease in the LUMO energy of C1.71 eV (at the PBE/DZP level). Similar patterns are obtained for the HOMO energies, i.e. the adsorption on graphene, boron-nitride or silicene/germanene leads to a decrease in a HOMO energies of 0.32 eV, 0.44 eV or 2.32 eV. Hence, isolated ionic liquid as well as ionic liquid on top of graphene or boron-nitride surfaces yield a HOMO-LUMO difference of around 3.46 eV; while the larger HOMO energy change for choline benzoate on top of silicene or germanene provides a HOMO-LUMO difference of C4.07 eV (at PBE/DZP level). However, the contribution from each ion to the HOMO varies as a nanosheet function. The HOMO is located over the choline motif, although it also shows a very important contribution from the benzoate anion, for isolated [CH][BE] and choline benzoate on top of the graphene/boron-nitride sheet. For these systems, interionic CT is C0.32 e À /0.51 e À (according ChelpG/Hirshfeld models), while CT from the benzoate to the sheet is close to zero (for both ChelpG/Hirshfeld schemes) at the PBE/DZP theoretical level. For the ionic liquid on top of silicene or germanene sheets, interionic CT is C0.28 e À or 0.19 e À , respectively, according to Hirshfeld populations. Both IL-silicene and IL-germanene systems yield an important CT from the ionic liquid (concretely from the anion) to the sheet, equal to 0.54 e À and 0.68 e À , respectively. These high CT values lead to the HOMO being mainly located on the benzoate anion. In addition, these HOMO contours also show some similarities with the HOMOÀ1 orbital for the pristine ionic liquid. ILnanosheet CTs according to the ChelpG model are close to zero. There should be a relationship between the IL-sheet CT transfer according to the Hirshfeld method and the HOMO energy: the greatest/smallest CT is obtained for those nanosheets which are able to provide the lowest/highest HOMO (from the ionic liquid) energy. Therefore, the Hirshfeld model could be more appropriate to study the adsorption of IL on top of nanosheets.

Conclusions
The adsorption of choline benzoate ionic liquid on the surface of 2D nanosheets, graphene, silicene, germanene and boronnitride was studied using the DFT approach. The interaction mechanisms, binding energies and electronic structures were calculated using the PBE functional as well as PBE plus dispersion corrections (PBE-D2). The latter provided useful information on the contributions from dispersion forces to the interaction between both ions (which is mainly due to Coulombic forces) as well as the interaction between the IL and the nanosheet, which is mainly due to dispersion forces. Charge transfer from the anion to the surface also plays an important role. There is also a relationship between interionic and ionic liquid-nanosheet interactions. Therefore, the larger interionic binding energies are obtained for choline benzoate on top of carbon and boron nitride based sheets, which provides the lowest IL-sheet interaction energies. The analysis results showed that the considered sheets interact with the anion through p-stacking, although there is also a remarkable charge transfer in IL-silicene and IL-germanene. The electronic structure analysis exposed that the nanosheet electronic structure is not affected during the IL adsorption process. However, HOMO and LUMO molecular orbitals located on the IL suffer changes in their energies upon adsorption. Thus, the largest changes are noted for choline benzoate on top of silicene and germanene sheets. For both systems, there is also a larger charge transfer from the IL to the sheet, which leads to changes in the molecular fragment contributions to the HOMO contours. In this work, ChelpG and Hirshfeld models have been used to compute charge transfer processes, with results showing that only the Hirshfeld method would provide coherent values for IL-sheet charge transfer processes and their relationship with binding energies or with the molecular orbital shapes. In short, IL-graphene and IL-BN would provide similar features (such as interaction mechanism, CT, binding energies or electronic structure), while IL-silicene and IL-germanene would bring greater changes on assessed properties.
The reported results provide useful information on the ionic liquid adsorption on different nanosheets and their interaction mechanisms. They could also provide insights for the tuning of ionic liquid properties through the adsorption on the adequate nanosheet.