2.1. Thermodynamics of Ion Selectivity
In thermodynamics, Gibbs free energy (Δ
G) can be utilized to investigate whether a process is likely to occur. As the solutes are embedded into water, this includes the solute–solute, the solute–water, and the water–water interactions. The total Gibbs energy may be expressed as
In fact, this may be applied to investigate the ion selectivity of channels as the ions pass through the filters.
For the selectivity filter of the KcsA channel, the ion selectivity represents a competition of K
+ and Na
+ for a binding site in the filter. The selectivity for K
+ over Na
+ is due to the difference in free energies of the ions in the pore, which departs from the corresponding difference in the bulk solution. It follows that the ion selectivity of binding is fundamentally governed by differences in relative free energies [
31], and the ion selectivity can be stated from the point of view of thermodynamic equilibrium as
Thermodynamically, it can be understood to allow the preferential permeation of K+ over Na+ by providing an unfavorable environment for Na+ in the presence of K+.
Regarding the KcsA K
+ channel, there are mainly O-atom quaternion rings with three sizes: 3.12/3.33 Å, 3.89 Å, and 4.22 Å. When the ions permeate the nanopore channels, one unambiguous observation based on the crystallographic X-ray structures of K
+ channels is that the K
+ ions must become essentially dehydrated to pass through the narrow selectivity filter [
6] and directly interact with the filter. Therefore, as the ions pass through the nanopore channels, this undoubtedly leads to the changes in water structure.
As the ion permeates the nanopore, the ion selectivity of the nanopore may be closely related to not only the changes in hydrogen bondings in water, but also the direct interactions between the ion and channel. Therefore, the ion selectivity of the nanopore may be related to various interactions, such as those due to the solute–solute, solute–water, and water–water interactions. In fact, before the ions are affected by the interactions between ions and the channel, it is necessary for the ions to approach the filter. This may be due to the rearrangement of water molecules as ions approach the nanopore, which is related to ΔGWater-water and ΔGSolute-water. Therefore, it is important to investigate the structure of liquid water, and the effects of dissolved solute on the water structure.
Numerous experimental and theoretical works have been conducted to investigate the structure of liquid water. To date, a variety of structural models for water have been proposed, broadly categorized into two groups: (a) mixture models and (b) distorted hydrogen bond or continuum models [
52]. The mixture model postulates the concurrent existence of two distinct types of structures. In the latter model, this implies that water consists of a random, three-dimensional network of hydrogen bonds with a wide distribution of O-H∙∙∙O hydrogen bond angles and distances. However, unlike in the mixture models, the water networks cannot be “broken” or separated into distinct molecular species.
The OH vibrations are sensitive to hydrogen bondings and are employed to study the structure of water. Based on the dependence of OH vibrations on water clusters ((H
2O)
n) [
42,
43,
44], it is found that, as three-dimensional hydrogen bondings emerge (
n ≥ 6), OH vibrations are predominantly influenced by the local hydrogen bondings (within the first shell) of water molecules, and the effects of hydrogen bonding beyond the first shell on OH vibrations are weak. Hence, different OH vibrations may reasonably be attributed to the OH vibrations engaged into various local hydrogen-bonded networks of a water molecule.
Under ambient conditions, the Raman OH stretching bands of water can be fitted into five sub-bands. These may be attributed to OH vibrations involved in various local hydrogen bonds, including DDAA (double donor-double acceptor, tetrahedral hydrogen bonding), DDA (double donor-single acceptor), DAA (single donor-double acceptor), and DA (single donor-single acceptor) hydrogen bondings, respectively. In ambient water, a water molecule interacts with neighboring water molecules (the first shell) through various local hydrogen-bonded networks [
44]. It is different from the mixture and continuum structural models of water. Additionally, under ambient conditions, the local hydrogen-bonded networks of a water molecule may be influenced by factors such as temperature, pressure, dissolved salt, and confined environments.
When a solute is embedded in water, an interface appears between the solute and water, which undoubtedly affects the structure of liquid water. As the OH vibration is primarily influenced by the local hydrogen bondings of a water molecule, the dissolved solute mainly affects the structure of interfacial water (the topmost water layer at the solute/water interface). In theory, vibrational sum frequency generation (SFG) spectroscopy is a technique selectively applicable to interfaces and is widely used to examine the structure of the air/water interface. Based on the SFG study [
45] on the air/water interface, it is observed that, in contrast to bulk water, no DDAA (tetrahedral) hydrogen bondings are present in interfacial water.
The solute mainly affects the hydrogen bondings of interfacial water. The formation of the solute/water interface is related to the loss of DDAA hydrogen bonding in interfacial water layer. Once the ratio of the interfacial water layer to the volume is established, it can be employed to calculate the Gibbs free energy between solute and water (interfacial water). Therefore, ΔG
Solute-water can reasonably be expressed as
where R
Interfacial water represents the ratio of interfacial to bulk water, ∆G
DDAA is the Gibbs energy of DDAA hydrogen bondings, and n
HB means the average hydrogen bondings of a water molecule. For DDAA hydrogen bonidngs, n
HB is equal to 2.
The hydration free energy refers to the change in Gibbs energy when an ion (or molecule) is transferred from a vacuum (or the gas phase) to a solvent. Thermodynamically, the process of inserting a hard sphere into water is equivalent to the formation of an empty spherical cavity in water. After the solute is regarded as an ideal sphere, the hydration free energy may be expressed as [
47]
where r
H2O is the average radius of an H
2O molecule, and R represents the radius of the solute. At ambient conditions, the Gibbs energy of DDAA hydrogen bondings is determined to be −2.66 kJ/mol [
47].
In principle, the lower the hydration free energy, the more thermodynamically stable the system becomes. The hydration free energy is the sum of Gibbs energy of pure water (ΔG
Water-water) and that of interfacial water (ΔG
Solute-water) related to the dissolved solute. Therefore, the structural transition is expected when ΔG
Water-water equals ΔG
Solute-water (
Figure 1).
where R
c represents the critical radius of solute. At 293 K and 0.1 MPa, R
c is calculated to be 6.5 Å. With increasing the solute size (or concentrations), it is divided into the initial and hydrophobic solvation processes [
47,
48]. Additionally, the Gibbs free energy of interfacial water (ΔG
Solute-water) is inversely proportional to the radius of solute (or the ratio of surface area to volume of solute). As the solute size (or solute concentrations) increases, different dissolved behaviors can be expected for various processes (
Figure 1) [
47,
48].
In the initial solvation process (∆GSolute-water < ∆GWater-water, both of which are negative), the hydration free energy is dominated by the Gibbs energy of interfacial water (ΔGSolute-water). To achieve greater thermodynamic stability, the solutes tend to be dispersed in water to maximize ΔGSolute-water. Consequently, water molecules are expected to exist between solutes. In other words, it appears that there are “repulsive” forces between the solutes as they are brought closer together. In the hydrophobic solvation process (∆GSolute-water > ∆GWater-water), the hydration free energy is primarily governed by the Gibbs energy of water (ΔGWater-water). To maximize the Gibbs energy of water, the dissolved solutes are expected to be aggregated to minimize the interfacial water. It seems that there are “attractive” forces between the dissolved solutes.
Hydrophobic effects are generally defined as the tendency of non-polar molecules (or molecular surfaces) to aggregate in an aqueous solution. In our recent investigations [
47,
48] on hydrophobic interactions, they were reasonably characterized as the tendency to minimize the ratio of the surface area to the volume of the solutes in order to maximize the hydrogen bondings of water. Of course, this is due to the dissolved solute mainly affecting the structure of interfacial water, and the hydrogen bondings of interfacial water are weaker than bulk water. Additionally, the origin of hydrophobic interactions arises from the structural competition between interfacial and bulk water [
47,
48].
The Gibbs free energy between solute and water (ΔG
Solute-water) is closely related to the water molecules of the interfacial water layer. When two solutes are embedded into water, the interfacial water of two solutes may be inversely proportional to the distance between them. Therefore, ΔG
Solute-water is reasonably expressed as
where R
Solute-solute means the distance between solutes. As two identical spheres are dissolved in ambient water, the R
c is 3.25 Å [
48]. As the distance between solutes is larger or less than R
c, hydrophobic interactions between solutes may be regarded as the “attractive” or the “repulsive” forces [
48].
To investigate the ion selectivity of nanopores, the test solute is a sphere, and the target filter is the membrane embedded with a nanopore (
Figure 2). As the distance between the ion and filter is decreased, the hydration free energy may reasonably be expressed as
where ΔG
Ion-filter is the Gibbs energy of interfacial water arising from the ion and filter, which is inversely proportional to the distance between them, ∝1/d
Ion-filter, where d
Ion-filter is the distance between the ion and filter (
Figure 2).
In fact, as the nanopore is filled with a circle with corresponding radius R, the filter becomes the membrane with the same size as the filter (
Figure 2). Therefore, as the ion approaches the filter, the Gibbs free energy (ΔG
Ion-filter) related to the ion and filter equals the difference in Gibbs energy between the ion and membrane (ΔG
Ion-membrane), and that between the ion and filled circle with radius R (ΔG
Ion-circle) (
Figure 2), which may be expressed as
As the ion enters the nanopore, the interfacial water molecules related to the ion and nanopore are changed into bulk water, which may be related to the ΔGIon-circle.
In comparison with the membrane, the test ion is treated as a point after the size of the ion is ignored. Additionally, the ion is restrained to perpendicularly move to the center of the membrane. Therefore, the Gibbs energy between the ion and membrane (ΔG
Ion-membrane) may mathematically be determined as (
Figure 2)
where w is the radius of membrane, and d means the distance between the centers of the ion and membrane. In addition, ΔG
Ion-membrane is higher than ΔG
Water-water (ΔG
Ion-membrane >ΔG
Water-water), and the hydrophobic interactions between the test ion and membrane substrate can reasonably be regarded as the “attractive” force [
48].
Additionally, Equation (6) may also be used to determine the Gibbs free energy between the ion and the circle (ΔG
Ion-circle). In this work, the nanopore (or circle) is simply regarded as a sphere with the radius R
Pore. Of course, ΔG
Ion-circle is also dependent on the distance between them. This may be reasonably expressed as
where d
Ion-circle is the distance between the centers of the ion and the circle (or nanopore), R
Pore is the pore radius, r
Ion is the ion size, and sep is the separation between the ion and the circle (
Figure 2). Therefore, when the ion approaches the nanopore, the Gibbs energy of interfacial water may be related to the pore radius (R
Pore), ion size (r
Ion), and the separation between them. However, different from the ΔG
Ion-membrane, the hydrophobic interactions between the ion and the circle (ΔG
Ion-circle) may be the “attractive” or “repulsive” forces, which may be closely related to the ion size and pore radius.
As the ion and filter are embedded into water, hydration free energy is the sum of the Gibbs energy of water and interfacial water, and the latter is composed of the interactions between the ion and membrane substrate (ΔG
Ion-membrane), and those between the ion and the circle with radius R (ΔG
Ion-circle). From Equation (10), it can be seen that a decrease in the d
Ion-circle may increase the ΔG
Ion-circle. In our recent works on hydrophobic interactions [
47,
48], dc (or R
c) was expected to be
From the equation, dc is expected as the ion and channel come into contact (or with the separation between them being zero, sep = 0). Therefore, as an ion enters the channel, the dc may be closely related to the ion size (rIon) and the radius of the channel (RPore). Regarding an ion with a specific size rIon, the corresponding pore radius (RPore,c) may be expected for the ion, RPore,c for rIon,c = rIon, or vice versa. Based on the discussion on hydrophobic interactions, with reference to dc, hydrophobic interactions between the ion and filter may be either “repulsive” or “attractive” forces when the ion approaches the filter.
Regarding the specific ion with radius r
Ion, as the radius of the pore is larger than R
Pore,c (R
Pore > R
Pore,c and r
Ion = r
Ion,c), it is associated with the hydrophobic process. Therefore, as the ion approaches the filter, an “attractive” force may be expected between them, and no energy barrier is expected (
Figure 3). In other words, the ion is expected to penetrate the pore of the filter, similar to the ion diffusion in aqueous solutions.
Additionally, as the radius of the nanopore is less than R
Pore,c (R
Pore < R
Pore,c and r
Ion = r
Ion,c), it is related to the initial solvation process. In the initial solvation process, the “repulsive” force between solutes is expected as they pushed together in water. In other words, as the ion approaches the nanopore, the energy barrier may be expected between the ion and the filter (
Figure 3).
From the above, as the ion approaches the nanopore, the energy barrier between the ion and the channel may be related to the difference in Gibbs energy between ΔG
Ion-circle and ΔG
Ion-circle (d = dc). This may be expressed as
Therefore, when an ion approaches the nanopore, the energy barrier between them is related to the ion size (rIon) and pore radius (RPore).
Additionally, regarding the specific pore with the radius being R
Pore, the energy barrier may also be expected as the ion size is less than r
Ion,c, R
Pore = R
Pore,c and r
Ion < r
Ion,c. For the filter with specific pore radius (R
Pore = R
Pore,c), it can be derived that the energy barrier may be dependent on the ion size (r
Ion). In comparison with r
Ion,c, the lower energy barrier may be expected for the larger ion, and the higher energy barrier is related to the smaller ion (
Figure 4). This may be utilized to understand the ion selectivity of the nanopore.
When the ions permeate the KcsA K
+ channel, they become essentially dehydrated to pass through the narrow selectivity filter. Therefore, as the ion enters the nanopore, direct interactions, such as van der Waals interactions, are expected between the ion and nanopore. This means that the energy barrier may also be affected by the intermolecular interactions between the ion and nanopore. In other words, the ion selectivity of the nanopore may be modulated through changing the molecular polarity of the channel. In fact, this can be demonstrated with many theoretical simulations [
31,
32].
As ion size is less than rIon,c (rIon < rIon,c for RPore = RPore,c), an energy barrier is expected between the ion and filter, and the ion may be suppressed to enter the nanopore of the filter. In other words, the external driving force, such as electrical potential or pressure, is necessarily added to the ion to surmount the energy barrier and pass through the nanopore. Of course, this is in accordance with the experimental studies on ion channel of the filter, which means that the electrical transmembrane potential is necessary inside and outside of the membrane to drive the ion through the channel.
From the above, the energy barrier between the ion and nanopore is related to ion size and pore radius. Regarding the nanopore with specific radius RPore (RPore = RPore,c), a lower energy barrier is expected for large ions, and a higher energy barrier is related to small ions. Therefore, compared to small ions, the weaker external force may be enough to make large ions pass through the nanopore. This means that the ion selectivity of nanopores may also be related to the strength of external forces.
When the ion approaches the nanopore, “repulsive” or “attractive” forces may be expected between the ion and nanopore. From the above, dc is related to not only the ion size (rIon) but also the radius of nanopore (RPore). Regarding the nanopore with specific RPore (RPore = RPore,c), the corresponding critical ion size (rIon,c) may be expected. As ion size is less than rIon,c, the energy barrier between the ion and nanopore is expected, rIon < rIon,c and RPore = RPore,c. With reference to rIon,c, the lower energy barrier may be expected for large ions. Additionally, as the ion passes through the filter, the direct interactions between the ion and channel, such as van der Waals interactions, become important, thereby influencing the energy barrier. In addition, the external driving force is necessarily added to the ion so that it can pass through the filter. This may be utilized to understand the ion selectivity of nanopores.
2.2. MD Simulations
From the above discussion, the ion selectivity of nanopores is closely related to the energy barrier as the ion enters the nanopore of the filter. Furthermore, the energy barrier between the ions and nanopore is related to both ion size and the radius of the pore. To understand the dependence of ion selectivity on the ion size, different ions (K+ and Na+) are forced to pass through the specific nanopore with a radius of 3.7 Å. In addition, to investigate the effects of nanopore pore radius on ion selectivity, MD simulations are also conducted on the systems in which a K+ ion enters the nanopores with different radii, such as 1.8 Å, 3.7 Å, and 5.6 Å.
In this work, the PMFs are calculated through ABF methods as a K
+ ion perpendicularly approaching the nanopore with a radius of 3.7 Å (
Figure 5). With a decrease in the distance between the K
+ ions and nanopores, various minima can be found on the calculated PMFs (
Figure 5). The minimum at 7.5 Å is found in the PMF, corresponding to a double water molecular layer located between the ion and nanopore (
Figure 6a). Another minimum is located at 4.7 Å, attributed to the solvent-separated PMF, which means that only a water molecular layer can be found between the K
+ ion and the graphene (
Figure 6b).
From the calculated PMFs at 300 K, with decreasing distance between the K
+ ion and graphene embedded with nanopores, energy barriers can be found between the neighboring minima (
Figure 5). Of course, they may be related to the expulsion of a single water layer in the confined volume as the solutes moved closer together. Additionally, an obvious energy barrier can be found on the calculated PMFs as the distance between the K
+ ions and nanopores is less than 5.0 Å, in which only one water layer can be found between them (
Figure 6b), especially the distance less than 2.4 Å where K
+ lies at the graphene surface (
Figure 6c). Due to the obvious energy barrier between K
+ ion and graphene, the ion is suppressed (or forbidden) from entering the nanopore of graphene (
Figure 6d).
In addition, as the ion enters the nanopore, no water molecules can be found between the K
+ ion and the nanopore of graphene, and the ion comes into contact with the graphene. This indicates that the energy barrier between the K
+ ion and the nanopore may be closely related to the expulsion of interfacial water when the ion enters the nanopore. In other words, the above energy barrier may be related to the water molecules changing from interfacial to bulk water as the K
+ ion approaches the nanopore of graphene. In combination with our recent studies [
47,
48], it may be concluded that this is closely related to hydrophobic interactions between the ion and nanopore.
Based on the above discussion on hydrophobic interactions, the dissolved K+ ion and graphene may be dispersed in water in the initial solvation process. To be more thermodynamically stable, water molecules can be expected between the K+ ion and graphene. In other words, it seems that there is a “repulsive” force as the K+ ion enters the filter, in which no water molecules can be found between them. Therefore, when the K+ ion enters the nanopore, the energy barrier may be due to the expulsion of interfacial water between the K+ ion and the graphene, and the K+ ion comes into contact with graphene in solutions. Of course, this is also in agreement with the above thermodynamic analysis on the energy barrier.
In our Raman spectroscopic studies [
42,
43,
44], the Raman OH vibrations are mainly dependent on the local hydrogen-bonded network of a water molecule as three-dimensional hydrogen bondings appear in the water. Therefore, in ambient conditions, the dissolved solute mainly affects the structure of interfacial water. Of course, this is demonstrated by other studies on the structure and dynamics of water around ions using neutron and X-ray diffraction [
53], X-ray absorption spectroscopy [
54], femtosecond time-resolved infrared (fs-IR) vibrational spectroscopy [
55,
56,
57], and optical Kerr effect spectroscopy [
58]. These mean that the effects of ions on water may be largely limited to the first solvation shell.
The dissolved solutes primarily affect the structure of interfacial water. When they are dissolved in water, they can be divided into interfacial and bulk water. To investigate the structural reorganization between interfacial and bulk water, they are calculated as the K
+ ion permeates the nanopore (
Figure 7). With the decrease in the distance between the ion and nanopore, there is a decrease in the molecular number for interfacial water (
Figure 7a), but an increase in the number for bulk water (
Figure 7b). When the K
+ ion enters the nanopore, the interfacial water associated with the ion and nanopore may be expelled into bulk water. Additionally, this is also related to the appearance of the energy barrier between the ion and nanopore.
To investigate the changes in hydrogen bondings when the K
+ ion is forced to enter the nanopore of graphene, the hydrogen bondings of water are also determined. In this study, the geometrical definition of hydrogen bonding is utilized to determine the hydrogen bondings in water [
59]. According to the geometrical definition, hydrogen bonding is considered to exist between two neighboring water molecules if the oxygen–oxygen distance (r
OO) and ∠OOH angle between two water molecules are less than 3.5 Å and 30°, respectively.
From the MD simulations, the hydrogen bonding number (average number of hydrogen bonds per water molecule, n
HB) can be calculated (
Figure 8). In comparison with bulk water, due to the truncations of hydrogen bonding at the solute/water interface, the hydrogen bonding of interfacial water is less than that of bulk water (
Figure 8). Therefore, the origin of hydrophobic interactions is due to the competition between the hydrogen bondings of bulk and interfacial water.
Regarding the nanopore with specific radius R
Pore, the energy barrier between the ion and nanopore is expected as the ion size is less than r
Ion,c, r
Ion < r
Ion,c and R
Pore = R
Pore,c. Additionally, the above energy barrier may be dependent on the ion size. To understand the dependence of the energy barrier on the ion size, MD simulations are carried out on different ions (K
+ and Na
+) passing through the specific nanopore with a radius of 3.7 Å. From the calculated PMFs, the obvious barrier can also be found as the distance between the Na
+ ion and the filter is less than 3.1 Å (
Figure 5). However, compared to the Na
+ ion, a lower energy barrier can be expected as the K
+ ion enters the nanopore.
From the above discussion, the origin of the energy barrier between the ion and nanopore is closely related to the expulsion of the interfacial water layer between them as the ion enters the nanopore. This is in accordance with the fact that the ion will lose the hydration shell and become essentially dehydrated as the ion enters the pore. The dissolved ion mainly affects the structure of the interfacial water, and the coordination number within the first hydration shell is related to the ion size. Regarding Na+, the number of water molecules in the first hydration layer is 5.4, which is less than the corresponding number of K+ (8.0). This may be applied to understand the dependence of the energy barrier on ion size.
Due to the energy barrier, when the ion enters the nanopore, the ion may be suppressed or forbidden from entering the channel. Therefore, the external driving force, such as the transmembrane electrical voltage or pressure, is necessarily added to the ion so that it may surmount the energy barrier and pass through the filter. Regarding the channel with a specific radius pore RPore, the energy barrier is dependent on the ion size. Under the same electrical voltage, it can be expected that the nanopore preferentially allows the specific ion with the lower barrier to pass through it. In comparison with small ions, a lower energy barrier is expected for larger ions. Therefore, compared to the smaller Na+ ions, the nanopore prefers the larger K+ ions to pass through the channel. This may be utilized to understand the ion selectivity of nanopores, such as KcsA for K+ over Na+.
In this study, a few MD simulations were conducted to investigate the effects of electrical voltage on ion selectivity (
Supplementary Materials). In the simulations, various electrical voltages were added to the aqueous solution, which was composed of 10 KCl and 10 NaCl molecules. It is necessary to provide enough voltage (2.0 kcal·mol
−1·Å
−1·e
−1) so that the K
+ ions may firstly pass through the filter with a pore radius of 3.7 Å (
Figure 9). With increasing voltage, more K
+ ions may pass through the filter. However, some Na
+ ions were found to penetrate the nanopore at high voltages. As the voltage is 3.0 kcal·mol
−1·Å
−1·e
−1, the filter gave the suitable ion selectivity of K
+/(K
+ + Na
+) (
Figure 9). Regarding the specific ion, it is important to design the corresponding nanopore for the ion. Additionally, it is also necessary to choose the reasonable voltage to drive the specific ion to pass through the corresponding nanopore of the filter.
In addition, the energy barrier between the ion and nanopore may be affected by not only the ion size but also the pore radius of the nanopore (R
Pore < R
Pore,c and r
Ion = r
Ion,c). In this study, to investigate the dependence of ion selectivity on the radius of the nanopore, MD simulations were also conducted on the systems in which the K
+ ion approaches the nanopore with different radii, namely 1.8 Å, 3.7 Å, and 5.6 Å. From the calculated PMFs (
Figure 10), the energy barrier can be found as the K
+ ion enters the nanopore with the different radii. However, with an increase in the radius of the nanopore, this decreases the energy barrier between the K
+ ion and nanopore. This is in agreement with the above discussion on the effects of pore radius on the energy barrier.
From the calculated PMFs (
Figure 10), a decrease in pore radius of the nanopore may lead to the increase in an energy barrier as the K
+ ion enters the nanopore. Therefore, the strong transmembrane electrical voltage or driving pressure is necessarily added to the ion so that it can surmount the energy barrier and pass through the small nanopore. In fact, due to the strong energy barrier when K
+ approaches the nanopore with a pore radius of 1.8 Å, the ion may be forbidden from entering the small nanopore even under high electrical voltage. Of course, the filter with this small pore radius is not suitably applied to separate the ion from others.
Additionally, an increase in the pore radius of the nanopore leads to a decrease in the energy barrier as K
+ permeates the nanopore (
Figure 10). In fact, no obvious energy barrier is found as K
+ enters the nanopore with a pore radius of 5.6 Å. Therefore, the K
+ ion may be expected to approach and freely enter the nanopore. In other words, this indicates that, with an increase in the pore radius, the hydrophobic interactions between the K
+ ion and the nanopore may change from “repulsive” to “attractive” forces. Additionally, this also means that the nanopore with a pore radius of 5.6 Å may lose the ion selectivity for K
+.
From the above discussion, it can be stated that the ion selectivity of the nanopore is associated with the energy barrier between the ion and the nanopore, arising from the transformation of water molecules from interfacial to bulk water as the ion enters the nanopore. Therefore, if the interfacial water of the ion is still preserved when it passes through the channel, it may be expected that no energy barrier can be found. In fact, it can be found that the interfacial water of the K
+ ion can still be found as K
+ permeates the nanopore with a pore radius of 5.6 Å (
Figure 11). In addition, this also indicates that the effects of the dissolved ions on water structure may be mainly limited to the first hydration shell of the ion.
The ion selectivity of the nanopore is due to the difference in the energy barrier as various ions pass through the channel. Thermodynamically, the nanopore may prefer the specific ion with a lower barrier to permeate the filter. From this study, it can be derived that the energy barrier is closely related to the initial solvation process in which the Gibbs free energy of interfacial water is lower than that of bulk water. As the ions pass through the nanopore, they may become dehydrated. Therefore, the energy barrier is due to the water molecules transforming from interfacial to bulk water, which is dependent on not only the ion size but also the pore radius of the nanopore. Additionally, as the ion enters the channel, the energy barrier may also be affected by the direct interactions between the ion and nanopore. Due to the energy barrier between the ions and nanopore, the driving forces, such as the electrical voltage, may be necessarily added to the ions so that they may surmount the energy barrier and pass through the nanopore (
Figure 12).