Next Article in Journal
Improved Dynamic Compressive and Electro-Thermal Properties of Hybrid Nanocomposite Visa Physical Modification
Next Article in Special Issue
A Review on the Progress of Optoelectronic Devices Based on TiO2 Thin Films and Nanomaterials
Previous Article in Journal
Simulation and Optimization of CNTs Cold Cathode Emission Grid Structure
Previous Article in Special Issue
Capacitive Behavior of Aqueous Electrical Double Layer Based on Dipole Dimer Water Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Surface Roughness on Partition and Selectivity of Ionic Liquids Mixture in Porous Electrode

National & Local Joint Engineering Research Center for Mineral Salt Deep Utilization, Key Laboratory for Palygorskite Science and Applied Technology of Jiangsu Province, Huaiyin Institute of Technology, Huai’an 223003, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(1), 51; https://doi.org/10.3390/nano13010051
Submission received: 1 December 2022 / Revised: 17 December 2022 / Accepted: 18 December 2022 / Published: 22 December 2022

Abstract

:
Understanding the influence of surface roughness on the adsorption of ions from an ionic liquids (ILs) mixture is essential for designing supercapacitors. The classical density functional theory (DFT) is applied to investigate the adsorption behavior of ILs mixtures in rough nanopores. The model parameters for each ion are determined by fitting experimental data of pure IL density. The results show that the smaller anions are densely accumulated near the rough surface and are the dominant species at a high positive potential. The exclusion of larger anions is enhanced by roughness at almost all potentials. At negative potential, the surface roughness promotes the adsorption of cations, and the partition coefficient increases with roughness. The partition coefficient of smaller anions is virtually independent of roughness. At positive potential, the surface roughness only promotes the adsorption of smaller anions and raises the partition coefficient. The partition coefficient of smaller anions is far greater than one. The selectivity of smaller anions for rough surfaces is very high and increases with roughness. The surface charge of a more uneven surface is significantly higher (about 30%) at a high potential.

1. Introduction

Ionic liquids are liquid salts consisting of organic cations and organic or inorganic anions. They can be in a liquid state below 100 °C. Because of the unique properties of ILs, they have been extensively used in many fields including catalysis, cellulose processing, batteries, supercapacitors, electrodeposition, and lubricant [1,2,3]. The knowledge of solid-IL interface is essential for the applications of ILs. The structure and adsorption of ILs at electrified interfaces significantly influence the performance of ILs in electrochemical processes [4,5].
ILs have been considered suitable electrolytes for supercapacitors that can store electrical energy through the reversible adsorption of ions. The ILs exhibit the electrical double layer (EDL) structure near the charged interface. Numerous studies have focused on the interfacial behavior of ILs near smooth surfaces. However, the solid surface of a porous electrode is usually not perfectly smooth but is somewhat rough. The surface roughness influences the ionic structure and affects the performance of ILs on electrified interfaces. Understanding the partition of ILs in rough nanopores is essential for their use as electrolytes in real applications. Recently some theoretical and simulation work has shown that electrode roughness can improve capacitive performance [6,7,8,9,10]. The surface roughness of porous electrodes has been experimentally shown to enhance the specific capacitance. The molecular dynamics simulations have indicated that the differential capacitance versus surface potential curves at rough surfaces are more complex.
Mixtures of ILs have gained considerable attention in tuning the properties of ILs. Many IL mixtures have been used in gas absorption, chromatography, liquid-liquid extraction, supercapacitors, and batteries [11,12,13,14,15,16,17,18]. For example, mixing ILs has been demonstrated experimentally to improve the gas absorption performance [13]. IL mixtures have been successfully utilized for dissolution and in situ catalytic degradations of cellulose in a batch reactor [19]. In electrochemistry, some mixing ILs can lead to lower viscosity and conductivity [15]. Recently, for applications of ILs in supercapacitors, the mixture of 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([C2mim][Tf2N]) and 1-ethyl-3-methylimidazolium tetrafluoroborate ([C2mim][BF4]) have been shown to provide symmetric charge storage and a maximum capacitance, and expand the operational potential window [20,21]. The mixture of tetramethylammonium tetrafluoroborate ([TMA][BF4]) and [C2mim][BF4] can provide simultaneous enhanced power and energy densities [22]. However, the fundamental understanding of mixed ILs in rough nanopores of electrode materials remains limited.
The classical density functional theory (DFT) is a successful and efficient method for studying interfacial properties of ILs at electrified surfaces by taking into account excluded volume effects and electrostatic correlations. The DFT using the fundamental measure theory and mean spherical approximation can reasonably model the EDL structure and capacitance of ILs [21,23,24]. For example, Jiang et al. predicted the capacitance oscillation for [C2mim][Tf2N] in nanopores [25]. Lian et al. studied the effects of curvature and pore size on the differential and integral capacitances [26]. Liu et al. studied the influence of polar solvents on the performance of IL supercapacitor [27]. Qing et al. studied the charging of asymmetric ILs supercapacitor using a time-dependent DFT [28]. The effects of shape asymmetry and formation of ion pair on ILs EDL have also been studied [29,30,31,32,33,34]. A key point in DFT is to find the accurate approximations of the Helmholtz free energy functional for the inhomogeneous system. A popular class of Helmholtz energy functionals is based on the statistical associating fluid theory (SAFT) equation of state (EoS). The DFT model based on a SAFT EoS can be used to study the bulk and inhomogeneous systems in a consistent way.
In this work, we employ the PC-SAFT-based DFT to study the adsorption and selectivity of mixed ILs in charged rough nanopores, which has rarely been considered in previous studies. We mainly focus on the impact of surface roughness. The Helmholtz energy functional of the DFT model accounts for the hard-sphere, the dispersive, and the electrostatic correlation interactions. A brief overview of the DFT model and a description of the confined system are presented in Section 2. In Section 3, the molecular parameters for each ion are fitted to the experimental density data of pure ILs. The DFT model is applied to explore the impact of roughness on ionic structure, adsorption, and selectivity of the ILs mixtures. Finally, conclusions are given in Section 4.

2. Model and Method

In this work, we consider mixtures of [C2mim][Tf2N] and [C2mim][BF4] confined in the rough slit-shaped pore of the porous electrode. To model the confined system with PC-SAFT-based DFT, the coarse-grained method is used to represent the ILs. Each IL-ion is modeled as a charged spherical particle of diameter σ. The detailed atomistic configuration of ions is not considered. The coarse-grained interaction potential between ions is assumed to be the sum of the soft repulsion of Chen and Kreglewski [35], Lennard-Jones dispersion and electrostatic interaction. The bulk properties of ILs can be calculated using the PC-SAFT-MSA EoS. The two pure IL-ion parameters in the model are the particle-size parameter σ and the dispersion-energy parameter ε. The Helmholtz free energy of PC-SAFT-MSA EoS is provided in previous publications [32,36].
Figure 1a schematically shows the mixture of ILs confined in the rough nanopore. We consider rough walls with arrays of identical square shape roughness elements of width we and height he. The rough walls are infinite in the y–z plane. The distance between the roughness elements is de.
The nanopore is considered to be an open system in contact with a bulk reservoir. The equilibrium properties of inhomogeneous ILs can be obtained from classical DFT. Therefore, it is convenient to use the grand potential Ω . Ω is a functional of fluid density, and it is written as
Ω = F + i d r ρ i ( r ) ( V i , e x t ( r ) μ i )
In Equation (1), F is the Helmholtz energy, ρ i is the density profile of species i, and μi is the bulk chemical potential. V i , e x t is the external potential composed of hard-wall and electrostatic contributions. The ionic profiles are uniform along the z-axis. Thus, we just calculate the two-dimensional ionic profile ρ ( x , y ) . In the y-direction, the external potential is periodic, and the ionic profiles have the period L y = w e + d e . We just need to calculate the ionic profiles in the finite range 0 y L y (see Figure 1b). Figure 1b shows a schematic for the computational domain. The line segments represent the solid wall, such as the line segment CD. The dotted line denotes the average profile of a rough surface. In order to get a better comparison between different surfaces, the pore volume accessible to the ions is a fixed value. The width of the smooth nanopore is equal to L x a v . The width of a rough nanopore L x is adjusted according to L x a v ( w e + d e ) = L x ( w e + d e ) 2 w e h e .
The Helmholtz energy functional F is a functional of density ρ i . F can be decomposed into an ideal term Fid and an excess term Fex
F = F i d + F e x
In Equation (2), Fex can be decomposed into four terms based on molecular interactions F e x = F h s + F d i s p + F c + F e l . Fc denotes the direct Coulomb contribution and is included in the mean electrostatic potential ϕ ( r ) . Fhs, Fdisp, and Fel denote hard-sphere repulsion, dispersive attraction, and electrostatic correlation contributions to the functional F. These contributions simplify to the corresponding contributions to Helmholtz energy of PC-SAFT-MSA EoS for a homogeneous fluid. Here we briefly introduce the functionals of PC-SAFT-MSA DFT. The modified fundamental measure theory (MFMT) is used to calculate the hard-sphere contribution Fhs [37]
β F h s = d r Φ h s [ n α ( r ) ]
Φ h s is the energy density from MFMT, it is a function of six weighted densities n α ( r ) proposed by Rosenfeld [38]. The MFMT can accurately describe inhomogeneous hard-sphere fluids.
Fdisp is the extension of the dispersive term of PC-SAFT based on the weighted density approximation [34].
β F d i s p = d r Φ d i s p [ ρ ¯ i ( r ) ]
where Φ d i s p is similar to the dispersive term of the PC-SAFT, it is a function of averaged densities ρ ¯ i ( r ) . More details of Fdisp are described in previous work [34].
Fel accounts for the electrostatic correlation contribution and is calculated by expanding the functional around the bulk fluid [39].
β F e l = β F e l [ ρ i , b ] i d r Δ C i 1 , e l ( r ) ( ρ i ( r ) ρ i , b ) 1 2 i j d r d r Δ C i j 2 , e l ( r , r ) ( ρ i ( r ) ρ i , b ) ( ρ j ( r ) ρ i , b )
Here, ρ i , b is the bulk density of species i. Δ C i 1 , e l represents the electrostatic contribution to the chemical potential. The expressions of Δ C i j 2 , e l from the mean spherical approximation are used [40].
In equilibrium, the grand potential is minimal. Thus minimization of Ω with respect to ρ i yields
δ Ω δ ρ i ( r ) = 0
Using Equation (1) and the expressions of F Equation (3) can be explicitly written as
ρ i ( r ) = exp ( β μ i δ β ( F h s + F d i s p + F e l ) δ ρ i ( r ) β Z i e ϕ ( r ) β V i , e x t ( r ) )
where β = 1 / k B T , kB is the Boltzmann constant, T is the temperature, Z i is the valence of component i, and e is the unit charge. Z i equals +1 for a cation and −1 for an anion. V i , e x t is the hard-wall part of V i , e x t , and the electrostatic part of V i , e x t is included in ϕ ( r ) .
For the calculation of the mixture, the Lorentz–Berthelot combining rules are used to obtain the size and energy parameter
{ σ i j = σ i + σ j 2 ϵ i j = ε i ε j
The binary interaction parameter between two components i and j was not applied.
We have used the fast Fourier transform to calculate the convolution integrals of functional derivatives in Equation (4). The local electrostatic potential ϕ ( r ) is related to the Poisson equation
2 ϕ ( r ) = e ε 0 ε r i Z i ρ i ( r )
ε 0 is the permittivity of the vacuum. ε r is the relative dielectric permittivity and is set at 12, this value is close to the experimental data [41]. It is assumed that surface polarization can be ignored.
The Poisson equation with appropriate boundary conditions is discretized using the finite difference method [42]. Equations (4) and (5) are solved self consistently using an iteration method for a specified surface potential ϕ s . The hybrid algorithm combining Picard iteration and Anderson acceleration is applied to speed up the calculation. The numerical details for DFT calculations can be found in the previous publication [34,43].
From the ionic profiles, we can calculate the total surface charge Q s (C/Å) of a solid surface based on the charge neutrality condition
Q s = 1 2 d x d y i e Z i ρ i ( x , y )
The average density ρ ¯ i of species i in the nanopore is calculated as
ρ ¯ i = 1 L x L y 2 w p h p d x d y ρ i ( x , y )
The partition coefficient is calculated as
Γ i = ρ ¯ i ρ i , b
The selectivity of the smaller anion [BF4] over the larger anion [Tf2N] is defined as S = Γ BF 4 / Γ Tf 2 N .

3. Results and Discussion

3.1. Bulk Density of Pure Ionic Liquids and Their Mixtures

In order to determine the model parameters for each ion, we first apply the PC-SAFT-MSA EoS to model the bulk density of pure ILs. The ionic parameters σ and ε are fitted to the experimental data of density of [C2mim][Tf2N] and [C2mim][BF4] [44,45,46]. The parameters are σ C 2 mim = 5.25   Å and ( ε / k B ) C 2 mim = 697.77   K , σ Tf 2 N = 6.44   Å and ( ε / k B ) Tf 2 N = 855.12   K , σ BF 4 = 4.54   Å and ( ε / k B ) BF 4 = 597.15   K . The absolute average relative deviation (AARD) is considered. The deviation is defined by AARD = 100 N i = 1 N | ρ c a l ρ e x p | / ρ e x p , ρ c a l and ρ e x p are the calculated density and experimental density, respectively. The AARD of 0.5% is obtained for pure ILs density. Figure 2a shows an example of model results at atmospheric pressure. The symbols denote experimentally measured density, and the lines indicate the model result. The model result is in good agreement with the experimental data.
Then, the ionic parameters are used to predict the density of mixtures of [C2mim][Tf2N]/[C2mim][BF4] without any binary parameter. Figure 2b shows an example of model results at a temperature of 293.1K and atmospheric pressure. The symbols are experimental data [47], and the line represents model prediction. The model prediction agrees with experimental data, and the AARD is about 2%. We have also calculated the relative deviation between the predicted density and experimental data for each composition, the result is shown in Figure S1 in the supplementary material. The maximum relative deviation is about 4%.

3.2. Ionic Liquids Mixtures in Rough Nanopores

In this section, the DFT model is used to study mixtures of ILs confined in the rough pore. The temperature and the pressure are fixed at 323 K and 1 bar. Two ILs mixtures were considered in this work with molar fraction of [C2mim][BF4] in the bulk phase XBF4 equals 0.2 and 0.5, respectively. The calculated number density of the mixture in the bulk phase is 4.97 nm−3 for XBF4 = 0.2 and 5.73 nm−3 for XBF4 = 0.5. For rough pores, the roughness width we is set at 1.4 nm, and the distance between roughness element de = we. The average size of the pore L x a v is set at 3 nm. The grid resolution of the two-dimensional DFT computation is 0.01 nm.
We have studied the adsorption of mixed ILs in nanopores with roughness height he ranging from 0 to 0.56 nm. The surface potential ϕ s varies from −1 V to 1 V. The Wenzel roughness factor is an important parameter to characterize a rough surface. The Wenzel roughness factor is the ratio of the true area to the apparent area of the solid surface [50]. Thus the roughness factor for the considered surface can be calculated as r e = 1 + 2 h e / ( w e + d e ) . For a smooth surface, re equals one. For a rough surface, re is greater than one and increases with he.
We first investigate the two-dimensional ion structure in rough nanopores. The rough walls induce oscillations in the density profiles in both normal (perpendicular to the substrate wall) and lateral (parallel to the substrate wall) directions. For example, Figure 3 and Figure 4 show two-dimensional density profiles for ions in a rough pore with the roughness factor r e = 1.4 , several normal density curves taken at the specified location, as well as the density profiles for ions in the smooth pore. More examples of normal and lateral density curves are shown in Figure S2. The molar fraction of [BF4] in the bulk ILs mixture is XBF4 = 0.2, and the surface potential ϕ s is 1 V. The smaller anions are densely packed near the rough surface with a significant lateral structure, and the density profile shows sharp peaks. The larger anions and cations are excluded from the surface.
In the x direction, the ions have a multilayer EDL structure. The smaller anions can get closer to the surface and occupy a smaller space, and they are more widely distributed in the pore. The contact density of smaller anions around the rough wall is higher than that around the smooth wall. The larger anions are mainly located near the middle of the pore, and almost none of them exist at the surface. Thus, at high surface potential, the smaller anions are more strongly adsorbed to the surface. In the y direction, at the bottom (line segment AB as shown in Figure 1b) and top (line segment CD) of the rough surface, there is considerable structure. As one can see in Figure 3c, there are three strong layers of smaller anions [BF4] at the bottom of the surface, the layer thickness is of the same order as the anion’s diameter. There are almost no cations and larger anions in this region. To gain a further understanding of the structure, we investigate the local excess adsorption, which is calculated as
N i ( x ) = d y ( ρ i ( x , y ) ρ i , b )
We have calculated the local excess adsorptions of three ions for two bulk concentrations, Figure 5a shows the result for bulk concentration XBF4 = 0.2, and the result for XBF4 = 0.5 is shown in Figure S3, all other parameters are the same as in Figure 3. There are four layers of excess ions near a surface. If we calculate the integral of ionic density in the layer we get the number of adsorbed ions. Figure 5b shows the number of ions adsorbed in each layer. The excess adsorption of smaller anions is higher and exhibits two peaks near the solid wall, and the small anions are the dominant species in the near-wall region. The first peak indicates the adsorption of small anions at the bottom of rough surface (line segments AB and EF), and the second peak indicates the adsorption around the top of rough surface (line segment CD). The smaller anions accumulate in the first and second layers and have high density, even at a low bulk concentration. The excess adsorption of smaller anions away from the wall is close to zero, and the local average density will be almost the same as the bulk value. When the bulk composition of smaller anions increases to 0.5, the value of excess adsorption peak for smaller ions around the wall shows little change. The number of larger anions in the first and second layers is close to zero and the total number of larger anions is much lower than smaller anions. Therefore, the larger anions tend to be excluded from the pore, and they are only enriched in the central region of the pore.
Then, we investigate the effect of surface roughness on the adsorption and selectivity of ions. Figure 6 shows the calculated partition coefficient of ions as a function of ϕ s , and the corresponding average density of each ion adsorbed inside the pore is shown in Figure S4, the bulk composition XBF4 is 0.2, and roughness factor r e = 1, 1.3 and 1.4. The partition coefficient is directly proportional to the average density as shown in Equation (8). At low absolute potential, the increasing surface roughness results in a decrease in the average density of cations and larger anions, so cations and larger anions tend to be excluded from the rough pores. ρ ¯ of larger anions [Tf2N] is significantly greater than that of smaller anions [BF4] because there are more anions [Tf2N] in the connected bulk phase. For the bulk composition XBF4 = 0.5, ρ ¯ of larger anions is only slightly greater than that of smaller anions.
At negative potential, the cations become enriched in the nanopore, and their average density monotonically increases with absolute potential. When the surfaces are sufficiently negatively charged, ρ ¯ of all anions decreases with increasing voltage. As the surface roughness increases, ρ ¯ of cations increases, and ρ ¯ of larger anions [Tf2N] decreases. Nevertheless, the average density of smaller anions [BF4] is almost independent of surface roughness. Overall, the adsorption of cations and the exclusion of larger anions are significantly enhanced by roughness.
At positive potential, the anions become enriched in the nanopore. As the surface roughness becomes higher, ρ ¯ of smaller anions monotonically increases, whereas the ρ ¯ of larger anions first increases, passes through a maximum, and then decreases. As the surface roughness increases, ρ ¯ of cations and larger anions [Tf2N] decreases, while ρ ¯ of smaller anions [BF4] significantly increases. Thus, the surface roughness induces exclusion toward larger anions in almost the entire range of potentials. The adsorption of smaller anions is promoted by roughness at positive potentials and is more energetically favorable than that of larger anions.
At low potential, ρ ¯ of all ions is lower than the bulk density, so all ions tend to be excluded from the nanopore. ρ ¯ of larger anions [Tf2N] is lower than the bulk density in the entire range of potentials. Thus, the partition coefficient of larger anions is less than one and decreases with increasing surface roughness.
At high positive potential, the high value of ρ ¯ of smaller anions indicates that smaller anions neutralize the surface charge favorably. ρ ¯ of smaller anions [BF4] is much higher than the bulk density. Thus, the partition coefficient of smaller anions is far greater than one and increases with increasing surface roughness. Table 1 and Table 2 show the partition coefficient and selectivity of smaller anions [BF4] for two bulk compositions, the potential ϕ s is 1 V. For the bulk composition XBF4 = 0.5, the partition coefficient of smaller anions at high positive potential is significantly lower (around 3).
At high negative potential, the partition coefficient of cations is greater than one, and is slightly lower for XBF4 = 0.5. The partition coefficient of smaller anions is virtually independent of roughness and slightly decreases with increasing absolute potential.
Figure 7 shows the selectivity of [BF4] over [Tf2N] as a function of surface potential for two bulk compositions. At low absolute potential, the selectivity S < 1, larger anions are preferred inside the pore. The surface roughness does not influence the selectivity.
When the surface is strongly negative-charged, the selectivity S > 1, and the smaller anions [BF4] are selected. A similar trend has been reported by Neal et al. for ILs mixture in smooth nanopores [51].
At moderate and high positive potential, adsorption of the smaller anions [BF4] is energetically more favorable, so the selectivity S > 1. The rough surfaces can enhance the adsorption of smaller anions, a significant selectivity for rough surfaces is observed, and the selectivity increases with increasing surface potential. The selectivity for rough surfaces with re = 1.4 is almost 1.6 times the value for smooth surfaces. The selectivity is lower for a higher bulk composition of [BF4].
Overall, at high potential, the electrosorption of small anions inside the pore is more energetically favorable. An increase in surface roughness leads to an increase in selectivity.
In the end, the surface charge Q s is investigated. For example, Figure 8 shows the surface charge versus ϕ s for XBF4 = 0.2. Q s monotonically increases with ϕ s . When the absolute value of ϕ s is high, the absolute value of Q s of a rough pore is considerably higher than that of a smooth pore. The absolute surface charge is higher for higher roughness. Because the volumes of different walls are alike, the rough pore tends to give a higher volume-specific integral capacitance.
Overall, the increase in surface roughness also induces a rise in the absolute surface charge. As the roughness increases, the reason for increasing absolute surface charge at negative potential is that the adsorption of cations and exclusion of larger anions are enhanced. The reason for increasing surface charge at positive potential is that the adsorption of smaller anions and exclusion of cations are enhanced.

4. Conclusions

In this work, we have investigated the adsorption and selectivity of mixed ILs in rough nanopores through a PC-SAFT-MSA-based DFT. The molecular parameters of ILs are fitted to the experimental data of pure liquid density. The molecular parameters are used to predict the bulk density of mixed ILs. The predicted density agrees with experimental data. For the confined system, the results show that at a high positive potential, the smaller anions are densely accumulated near the rough surface and are the dominant species. The larger anions tend to be excluded from the pore. The surface roughness enhances the exclusion of larger anions in the entire range of potentials, and the partition coefficient decreases with increasing roughness. At negative potential, the adsorption of cations is promoted by surface roughness, and the partition coefficient increases with increasing roughness. However, the surface roughness does not affect the adsorption of smaller anions; thus, the partition coefficient of smaller anions is independent of roughness. At positive potential, the surface roughness significantly enhances the adsorption of smaller anions. The partition coefficient of smaller anions is far greater than one and increases with surface roughness. The selectivity of smaller anions for a rough pore is higher than for a smooth pore and increases with roughness. At high surface potential, the absolute value of the surface charge of a rough pore is significantly larger than that of a smooth pore. Therefore, the rough electrode can hold more charge.
This work gives insight into the selective adsorption of mixed ILs in rough nanopores and helps to understand the charge storage on real mixed ILs/electrode interface. The results indicate that both composition and atomic scale roughness of surface are important factors that have to be considered for the optimization of ILs capacitors. In the present study, some other factors that may impact the partition and selectivity of mixed ILs are neglected, for example, the atomic configuration of ions, surface polarization, and the association between the cation and the anion. The IL is fully dissociated in the model. Therefore, the model needs to be further improved for some protic ILs that contain neutral species.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano13010051/s1, Figure S1: Relative deviation for density of mixture [C2mim][Tf2N]/[C2mim][BF4]. Relative deviation is calculated as ( ρ c a l ρ e x p ) / ρ e x p ; Figure S2: One-dimensional density curves of ions taken at the specified x and y positions. Other conditions are the same as in Figure 3; Figure S3: Local excess adsorption of ions in nanopore with rough surfaces, XBF4 = 0.5, other conditions are the same as in Figure 3; Figure S4: Average density of ions as a function of ϕ s , the bulk molar fraction XBF4 is 0.2.

Author Contributions

Conceptualization, G.S.; methodology, G.S.; validation, H.Y.; investigation, G.S. and Y.H.; resources, H.L. and K.H.; writing—original draft preparation, G.S.; writing—review and editing, X.Z. and F.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by National Natural Science Foundation of China (No. 21606096).

Data Availability Statement

The data that support the findings of this study are available within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Q.; Wang, Q.; Zhang, S.; Lu, X.; Zhang, X. Electrodeposition in Ionic Liquids. ChemPhysChem 2016, 17, 335–351. [Google Scholar] [CrossRef] [PubMed]
  2. Dong, K.; Liu, X.; Dong, H.; Zhang, X.; Zhang, S. Multiscale Studies on Ionic Liquids. Chem. Rev. 2017, 117, 6636–6695. [Google Scholar] [CrossRef] [PubMed]
  3. Li, L.; Zhang, Y.; Sun, Y.; Sun, S.; Shen, G.; Zhao, P.; Cui, J.; Qiao, H.; Wang, Y.; Zhou, H. Manufacturing pure cellulose films by recycling ionic liquids as plasticizers. Green Chem. 2020, 22, 3835–3841. [Google Scholar] [CrossRef]
  4. Fedorov, M.V.; Kornyshev, A.A. Ionic Liquids at Electrified Interfaces. Chem. Rev. 2014, 114, 2978–3036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Lian, Z.; Chao, H.; Wang, Z.-G. Effects of Confinement and Ion Adsorption in Ionic Liquid Supercapacitors with Nanoporous Electrodes. ACS Nano 2021, 15, 11724–11733. [Google Scholar] [CrossRef]
  6. Vatamanu, J.; Cao, L.; Borodin, O.; Bedrov, D.; Smith, G.D. On the Influence of Surface Topography on the Electric Double Layer Structure and Differential Capacitance of Graphite/Ionic Liquid Interfaces. J. Phys. Chem. Lett. 2011, 2, 2267–2272. [Google Scholar] [CrossRef]
  7. Xing, L.; Vatamanu, J.; Smith, G.D.; Bedrov, D. Nanopatterning of Electrode Surfaces as a Potential Route to Improve the Energy Density of Electric Double-Layer Capacitors: Insight from Molecular Simulations. J. Phys. Chem. Lett. 2012, 3, 1124–1129. [Google Scholar] [CrossRef]
  8. Vatamanu, J.; Borodin, O.; Bedrov, D.; Smith, G.D. Molecular Dynamics Simulation Study of the Interfacial Structure and Differential Capacitance of Alkylimidazolium Bis(trifluoromethanesulfonyl)imide [Cnmim][TFSI] Ionic Liquids at Graphite Electrodes. J. Phys. Chem. C 2012, 116, 7940–7951. [Google Scholar] [CrossRef]
  9. Vatamanu, J.; Vatamanu, M.; Bedrov, D. Non-Faradaic Energy Storage by Room Temperature Ionic Liquids in Nanoporous Electrodes. ACS Nano 2015, 9, 5999–6017. [Google Scholar] [CrossRef]
  10. Wei, J.; Li, Y.; Dai, D.; Zhang, F.; Zou, H.; Yang, X.; Ji, Y.; Li, B.; Wei, X. Surface Roughness: A Crucial Factor to Robust Electric Double Layer Capacitors. ACS Appl. Mater. Inter. 2020, 12, 5786–5792. [Google Scholar] [CrossRef]
  11. Li, W.; Liu, Y.; Wang, L.; Gao, G. Using Ionic Liquid Mixtures to Improve the SO2 Absorption Performance in Flue Gas. Energy Fuels 2017, 31, 1771–1777. [Google Scholar] [CrossRef]
  12. Pinto, A.M.; Rodríguez, H.; Colón, Y.J.; Arce, A., Jr.; Arce, A.; Soto, A. Absorption of Carbon Dioxide in Two Binary Mixtures of Ionic Liquids. Ind. Eng. Chem. Res. 2013, 52, 5975–5984. [Google Scholar] [CrossRef]
  13. Wang, M.; Zhang, L.; Gao, L.; Pi, K.; Zhang, J.; Zheng, C. Improvement of the CO2 Absorption Performance Using Ionic Liquid [NH2emim][BF4] and [emim][BF4]/[bmim][BF4] Mixtures. Energy Fuels 2013, 27, 461–466. [Google Scholar] [CrossRef]
  14. Baltazar, Q.Q.; Leininger, S.K.; Anderson, J.L. Binary ionic liquid mixtures as gas chromatography stationary phases for improving the separation selectivity of alcohols and aromatic compounds. J. Chromatogr. A 2008, 1182, 119–127. [Google Scholar] [CrossRef]
  15. Taige, M.; Hilbert, D.; Schubert, T.J.S. Mixtures of Ionic Liquids as Possible Electrolytes for Lithium Ion Batteries. Z. Phys. Chem. 2012, 226, 129–139. [Google Scholar] [CrossRef]
  16. García, S.; Larriba, M.; García, J.; Torrecilla, J.S.; Rodríguez, F. Liquid–liquid extraction of toluene from n-heptane using binary mixtures of N-butylpyridinium tetrafluoroborate and N-butylpyridinium bis(trifluoromethylsulfonyl)imide ionic liquids. Chem. Eng. J. 2012, 180, 210–215. [Google Scholar] [CrossRef]
  17. Costa, R.; Voroshylova, I.V.; Cordeiro, M.N.D.S.; Pereira, C.M.; Silva, A.F. Enhancement of differential double layer capacitance and charge accumulation by tuning the composition of ionic liquids mixtures. Electrochim. Acta 2018, 261, 214–220. [Google Scholar] [CrossRef]
  18. Osti, N.C.; Gallegos, A.; Dyatkin, B.; Wu, J.; Gogotsi, Y.; Mamontov, E. Mixed Ionic Liquid Improves Electrolyte Dynamics in Supercapacitors. J. Phys. Chem. C 2018, 122, 10476–10481. [Google Scholar] [CrossRef]
  19. Long, J.; Guo, B.; Li, X.; Jiang, Y.; Wang, F.; Tsang, S.C.; Wang, L.; Yu, K.M.K. One step catalytic conversion of cellulose to sustainable chemicals utilizing cooperative ionic liquid pairs. Green Chem. 2011, 13, 2334–2338. [Google Scholar] [CrossRef]
  20. Van Aken, K.L.; Beidaghi, M.; Gogotsi, Y. Formulation of Ionic-Liquid Electrolyte to Expand the Voltage Window of Supercapacitors. Angew. Chem. Int. Ed. 2015, 54, 4806–4809. [Google Scholar] [CrossRef]
  21. Lian, C.; Liu, K.; Van Aken, K.L.; Gogotsi, Y.; Wesolowski, D.J.; Liu, H.L.; Jiang, D.E.; Wu, J.Z. Enhancing the Capacitive Performance of Electric Double-Layer Capacitors with Ionic Liquid Mixtures. ACS Energy Lett. 2016, 1, 21–26. [Google Scholar] [CrossRef]
  22. Wang, X.; Mehandzhiyski, A.Y.; Arstad, B.; Van Aken, K.L.; Mathis, T.S.; Gallegos, A.; Tian, Z.; Ren, D.; Sheridan, E.; Grimes, B.A.; et al. Selective Charging Behavior in an Ionic Mixture Electrolyte-Supercapacitor System for Higher Energy and Power. J. Am. Chem. Soc. 2017, 139, 18681–18687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Lian, C.; Liu, H.; Wu, J. Ionic Liquid Mixture Expands the Potential Window and Capacitance of a Supercapacitor in Tandem. J. Phys. Chem. C 2018, 122, 18304–18310. [Google Scholar] [CrossRef]
  24. Ma, K.; Janssen, M.; Lian, C.; van Roij, R. Dynamic density functional theory for the charging of electric double layer capacitors. J. Chem. Phys. 2022, 156, 084101. [Google Scholar] [CrossRef] [PubMed]
  25. Jiang, D.-e.; Jin, Z.; Wu, J. Oscillation of Capacitance inside Nanopores. Nano Lett. 2011, 11, 5373–5377. [Google Scholar] [CrossRef]
  26. Lian, C.; Jiang, D.-e.; Liu, H.; Wu, J. A Generic Model for Electric Double Layers in Porous Electrodes. J. Phys. Chem. C 2016, 120, 8704–8710. [Google Scholar] [CrossRef]
  27. Liu, K.; Wu, J. Boosting the Performance of Ionic-Liquid-Based Supercapacitors with Polar Additives. J. Phys. Chem. C 2016, 120, 24041–24047. [Google Scholar] [CrossRef]
  28. Qing, L.; Jiang, J. Double-Edged Sword of Ion-Size Asymmetry in Energy Storage of Supercapacitors. J. Phys. Chem. Lett. 2022, 13, 1438–1445. [Google Scholar] [CrossRef]
  29. Ma, K.; Forsman, J.; Woodward, C.E. Influence of ion pairing in ionic liquids on electrical double layer structures and surface force using classical density functional approach. J. Chem. Phys. 2015, 142, 174704. [Google Scholar] [CrossRef]
  30. Liu, K.; Wu, J. Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids. J. Chem. Phys. 2020, 152, 054708. [Google Scholar] [CrossRef]
  31. Gallegos, A.; Lian, C.; Dyatkin, B.; Wu, J. Side-chain effects on the capacitive behaviour of ionic liquids in microporous electrodes. Mol. Phys. 2019, 117, 3603–3613. [Google Scholar] [CrossRef]
  32. Shen, G.; Sun, Y.; Wang, Y.; Lu, X.; Ji, X. Interfacial structure and differential capacitance of ionic liquid/graphite interface: A perturbed-chain SAFT density functional theory study. J. Mol. Liq. 2020, 310, 113199. [Google Scholar] [CrossRef]
  33. Yang, J.; Lian, C.; Liu, H. Chain length matters: Structural transition and capacitance of room temperature ionic liquids in nanoporous electrodes. Chem. Eng. Sci. 2020, 227, 115927. [Google Scholar] [CrossRef]
  34. Shen, G.; Zhang, D.; Hu, Y.; Zhang, X.; Zhou, F.; Qian, Y.; Lu, X.; Ji, X. Effect of surface roughness on partition of ionic liquids in nanopores by a perturbed-chain SAFT density functional theory. J. Chem. Phys. 2022, 157, 014701. [Google Scholar] [CrossRef]
  35. Chen, S.S.; Kreglewski, A. Applications of the Augmented van der Waals Theory of Fluids.: I. Pure Fluids. Ber. Bunsenges. Phys. Chem. 1977, 81, 1048–1052. [Google Scholar] [CrossRef]
  36. Gross, J.; Sadowski, G. Perturbed-Chain SAFT:  An Equation of State Based on a Perturbation Theory for Chain Molecules. Ind. Eng. Chem. Res. 2001, 40, 1244–1260. [Google Scholar] [CrossRef]
  37. Yu, Y.X.; Wu, J.Z. Structures of hard-sphere fluids from a modified fundamental-measure theory. J. Chem. Phys. 2002, 117, 10156–10164. [Google Scholar] [CrossRef] [Green Version]
  38. Rosenfeld, Y. Free-energy model for the inhomogeneous hard-sphere fluid mixture and density-functional theory of freezing. Phys. Rev. Lett. 1989, 63, 980. [Google Scholar] [CrossRef]
  39. Wu, J.; Jiang, T.; Jiang, D.-e.; Jin, Z.; Henderson, D. A classical density functional theory for interfacial layering of ionic liquids. Soft Matter 2011, 7, 11222. [Google Scholar] [CrossRef]
  40. Hiroike, K. Supplement to Blum’s theory for asymmetric electrolytes. Mol. Phys. 1977, 33, 1195–1198. [Google Scholar] [CrossRef]
  41. Weingaertner, H. Understanding Ionic Liquids at the Molecular Level: Facts, Problems, and Controversies. ChemInform 2008, 39. [Google Scholar] [CrossRef]
  42. LeVeque, R.J. Finite Difference Methods for Ordinary and Partial Differential Equations; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2007. [Google Scholar]
  43. Shen, G.; Sun, Y.; Zhang, X.; Gao, X.; Qian, Y.; Lu, X.; Ji, X. Partition and selectivity of electrolytes in cylindrical nanopores with heterogeneous surface charge. J. Mol. Liq. 2021, 340, 116839. [Google Scholar] [CrossRef]
  44. Gardas, R.L.; Freire, M.G.; Carvalho, P.J.; Marrucho, I.M.; Fonseca, I.M.A.; Ferreira, A.G.M.; Coutinho, J.A.P. PρT Measurements of Imidazolium-Based Ionic Liquids. J. Chem. Eng. Data 2007, 52, 1881–1888. [Google Scholar] [CrossRef]
  45. Sanmamed, Y.A.; González-Salgado, D.; Troncoso, J.; Romani, L.; Baylaucq, A.; Boned, C. Experimental methodology for precise determination of density of RTILs as a function of temperature and pressure using vibrating tube densimeters. J. Chem. Themodyn. 2010, 42, 553–563. [Google Scholar] [CrossRef]
  46. Taguchi, R.; Machida, H.; Sato, Y.; Smith, R.L. High-Pressure Densities of 1-Alkyl-3-methylimidazolium Hexafluorophosphates and 1-Alkyl-3-methylimidazolium Tetrafluoroborates at Temperatures from (313 to 473) K and at Pressures up to 200 MPa. J. Chem. Eng. Data 2009, 54, 22–27. [Google Scholar] [CrossRef]
  47. Finotello, A.; Bara, J.E.; Narayan, S.; Camper, D.; Noble, R.D. Ideal Gas Solubilities and Solubility Selectivities in a Binary Mixture of Room-Temperature Ionic Liquids. J. Phys. Chem. B 2008, 112, 2335–2339. [Google Scholar] [CrossRef]
  48. Yadav, A.; Guha, A.; Pandey, A.; Pal, M.; Trivedi, S.; Pandey, S. Densities and dynamic viscosities of ionic liquids having 1-butyl-3-methylimidazolium cation with different anions and bis(trifluoromethylsulfonyl)imide anion with different cations in the temperature range (283.15 to 363.15)K. J. Chem. Themodyn. 2018, 116, 67–75. [Google Scholar] [CrossRef]
  49. Neves, C.M.S.S.; Kurnia, K.A.; Coutinho, J.A.P.; Marrucho, I.M.; Lopes, J.N.C.; Freire, M.G.; Rebelo, L.P.N. Systematic Study of the Thermophysical Properties of Imidazolium-Based Ionic Liquids with Cyano-Functionalized Anions. J. Phys. Chem. B 2013, 117, 10271–10283. [Google Scholar] [CrossRef]
  50. Wenzel, R.N. Resistance of solid surfaces to wetting by water. Ind. Eng. Chem. 1936, 28, 988–994. [Google Scholar] [CrossRef]
  51. Neal, J.N.; Wesolowski, D.J.; Henderson, D.; Wu, J. Ion distribution and selectivity of ionic liquids in microporous electrodes. J. Chem. Phys. 2017, 146, 174701. [Google Scholar] [CrossRef]
Figure 1. (a) A schematic of the structure of IL-ions and a rough nanopore. (b) A schematic of the computational domain.
Figure 1. (a) A schematic of the structure of IL-ions and a rough nanopore. (b) A schematic of the computational domain.
Nanomaterials 13 00051 g001
Figure 2. (a) Bulk density of pure ILs [C2mim][Tf2N] and [C2mim][BF4] at T = 293.1K and atmospheric pressure. Symbols: experimental data [48,49]; lines: model results. (b) Density of mixture [C2mim][Tf2N]/[C2mim][BF4] as a function of mole fraction of [C2mim][BF4] (XBF4), T = 293.1K. Symbols: experimental data [47]; line: model prediction.
Figure 2. (a) Bulk density of pure ILs [C2mim][Tf2N] and [C2mim][BF4] at T = 293.1K and atmospheric pressure. Symbols: experimental data [48,49]; lines: model results. (b) Density of mixture [C2mim][Tf2N]/[C2mim][BF4] as a function of mole fraction of [C2mim][BF4] (XBF4), T = 293.1K. Symbols: experimental data [47]; line: model prediction.
Nanomaterials 13 00051 g002
Figure 3. Two-dimensional density profiles of [C2mim]+ (a), [Tf2N] (b), and [BF4] (c) in a nanopore at surface potential ϕ s = 1   V , XBF4 = 0.2, and r e = 1.4 that corresponds to he = 0.56 nm.
Figure 3. Two-dimensional density profiles of [C2mim]+ (a), [Tf2N] (b), and [BF4] (c) in a nanopore at surface potential ϕ s = 1   V , XBF4 = 0.2, and r e = 1.4 that corresponds to he = 0.56 nm.
Nanomaterials 13 00051 g003
Figure 4. (a) One-dimensional density curves of ions taken at the specified y position. Other conditions are the same as in Figure 3. (b) Density profiles for ions in a smooth pore at surface potential ϕ s = 1   V , and XBF4 = 0.2.
Figure 4. (a) One-dimensional density curves of ions taken at the specified y position. Other conditions are the same as in Figure 3. (b) Density profiles for ions in a smooth pore at surface potential ϕ s = 1   V , and XBF4 = 0.2.
Nanomaterials 13 00051 g004
Figure 5. (a) Local excess adsorption of ions in nanopore with rough surfaces. (b) Number of ions adsorbed in the four layers near the surface (i = 1 is the first layer, i = 2 is the second layer, i = 3 is the third layer, i = 4 is the fourth layer). Other conditions are the same as in Figure 3.
Figure 5. (a) Local excess adsorption of ions in nanopore with rough surfaces. (b) Number of ions adsorbed in the four layers near the surface (i = 1 is the first layer, i = 2 is the second layer, i = 3 is the third layer, i = 4 is the fourth layer). Other conditions are the same as in Figure 3.
Nanomaterials 13 00051 g005
Figure 6. Partition coefficient of ions as a function of ϕ s , the bulk molar fraction XBF4 is 0.2.
Figure 6. Partition coefficient of ions as a function of ϕ s , the bulk molar fraction XBF4 is 0.2.
Nanomaterials 13 00051 g006
Figure 7. Selectivity of [BF4] over [Tf2N] as a function of ϕ s . (a) XBF4 = 0.2 and (b) XBF4 = 0.5.
Figure 7. Selectivity of [BF4] over [Tf2N] as a function of ϕ s . (a) XBF4 = 0.2 and (b) XBF4 = 0.5.
Nanomaterials 13 00051 g007
Figure 8. Surface charge as a function of ϕ s , XBF4 = 0.2.
Figure 8. Surface charge as a function of ϕ s , XBF4 = 0.2.
Nanomaterials 13 00051 g008
Table 1. Partition coefficient and selectivity of [BF4], ϕ s = 1   V and XBF4 = 0.2.
Table 1. Partition coefficient and selectivity of [BF4], ϕ s = 1   V and XBF4 = 0.2.
re
11.31.4
ΓBF46.67.38
S10.712.715.7
Table 2. Partition coefficient and selectivity of [BF4], ϕ s = 1   V and XBF4 = 0.5.
Table 2. Partition coefficient and selectivity of [BF4], ϕ s = 1   V and XBF4 = 0.5.
re
11.31.4
ΓBF42.72.93.1
S4.44.95.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shen, G.; Yang, H.; Hu, Y.; Zhang, X.; Zhou, F.; Li, H.; Hong, K. Impact of Surface Roughness on Partition and Selectivity of Ionic Liquids Mixture in Porous Electrode. Nanomaterials 2023, 13, 51. https://doi.org/10.3390/nano13010051

AMA Style

Shen G, Yang H, Hu Y, Zhang X, Zhou F, Li H, Hong K. Impact of Surface Roughness on Partition and Selectivity of Ionic Liquids Mixture in Porous Electrode. Nanomaterials. 2023; 13(1):51. https://doi.org/10.3390/nano13010051

Chicago/Turabian Style

Shen, Gulou, Haoguang Yang, Yongke Hu, Xiaojie Zhang, Feng Zhou, Huaju Li, and Kun Hong. 2023. "Impact of Surface Roughness on Partition and Selectivity of Ionic Liquids Mixture in Porous Electrode" Nanomaterials 13, no. 1: 51. https://doi.org/10.3390/nano13010051

APA Style

Shen, G., Yang, H., Hu, Y., Zhang, X., Zhou, F., Li, H., & Hong, K. (2023). Impact of Surface Roughness on Partition and Selectivity of Ionic Liquids Mixture in Porous Electrode. Nanomaterials, 13(1), 51. https://doi.org/10.3390/nano13010051

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop