Next Article in Journal
Electrochemical Measurement of Interfacial Distribution and Diffusion Coefficients of Electroactive Species for Ion-Exchange Membranes: Application to Br2/Br Redox Couple
Next Article in Special Issue
Structure, Thermal Properties and Proton Conductivity of the Sulfonated Polyphenylquinoxalines
Previous Article in Journal
Separation of Soluble Benzene from an Aqueous Solution by Pervaporation Using a Commercial Polydimethylsiloxane Membrane
Previous Article in Special Issue
Nanocomposite Membranes Based on Fluoropolymers for Electrochemical Energy Sources
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanism Understanding of Li-ion Separation Using A Perovskite-Based Membrane

by
Mahsa Golmohammadi
1,
Meysam Habibi
2,
Sima Rezvantalab
3,
Yasin Mehdizadeh Chellehbari
4,
Reza Maleki
5,* and
Amir Razmjou
6,7,*
1
Department of Polymer Engineering & Color Technology, Amirkabir University of Technology, Tehran 15916-34311, Iran
2
School of Chemical Engineering, College of Engineering, University of Tehran, Tehran 14176-14411, Iran
3
Renewable Energies Department, Faculty of Chemical Engineering, Urmia University of Technology, Urmia 57166-419, Iran
4
Department of Chemical and Petroleum Engineering, Sharif University of Technology, Tehran 11155-9161, Iran
5
Computational Biology and Chemistry Group (CBCG), Universal Scientific Education and Research Network (USERN), Tehran 14197-33141, Iran
6
School of Engineering, Edith Cowan University, 270 Joondalup Drive, Joondalup, Perth 6027, WA, Australia
7
UNESCO Centre for Membrane Science and Technology, School of Chemical Engineering, University of New South Wales, Sydney 2052, NSW, Australia
*
Authors to whom correspondence should be addressed.
Membranes 2022, 12(11), 1042; https://doi.org/10.3390/membranes12111042
Submission received: 19 September 2022 / Revised: 13 October 2022 / Accepted: 24 October 2022 / Published: 26 October 2022
(This article belongs to the Special Issue Advanced Membranes for Energy Storage and Conversion)

Abstract

:
Lithium ions play a crucial role in the energy storage industry. Finding suitable lithium-ion-conductive membranes is one of the important issues of energy storage studies. Hence, a perovskite-based membrane, Lithium Lanthanum Titanate (LLTO), was innovatively implemented in the presence and absence of solvents to precisely understand the mechanism of lithium ion separation. The ion-selective membrane’s mechanism and the perovskite-based membrane’s efficiency were investigated using Molecular Dynamic (MD) simulation. The results specified that the change in the ambient condition, pH, and temperature led to a shift in LLTO pore sizes. Based on the results, pH plays an undeniable role in facilitating lithium ion transmission through the membrane. It is noticeable that the hydrogen bond interaction between the ions and membrane led to an expanding pore size, from (1.07 Å) to (1.18–1.20 Å), successfully enriching lithium from seawater. However, this value in the absence of the solvent would have been 1.1 Å at 50 °C. It was found that increasing the temperature slightly impacted lithium extraction. The charge analysis exhibited that the trapping energies applied by the membrane to the first three ions (Li+, K+, and Na+) were more than the ions’ hydration energies. Therefore, Li+, K+, and Na+ were fully dehydrated, whereas Mg2+ was partially dehydrated and could not pass through the membrane. Evaluating the membrane window diameter, and the combined effect of the three key parameters (barrier energy, hydration energy, and binding energy) illustrates that the required energy to transport Li ions through the membrane is higher than that for other monovalent cations.

1. Introduction

Nowadays, Li is known as one of the most demanding materials with ever-increasing applications in different fields of research such as medicine [1,2], energy storage devices [3,4,5], air treatment products [6], Li-ion batteries [3,7,8], ceramics and glass [9,10,11], lubricating greases, polymer production, and alloys [12,13,14]. Li is the lightest alkali metal with a density of 0.534 g/cm3, which is electrochemically active with a high electrode potential of −3.05 V and has the highest specific heat capacity of any solid element. These properties make Li compounds very appealing in a variety of commercial applications. Based on the abovementioned applications, massive sources of Li compounds are expected to exist. Li resources are divided into two main categories: solids (e.g., mineral ores, recycled-waste Li-ion batteries, and electronic waste) and liquids (e.g., salt-lake brine, geothermal brine, and seawater). Seawater has the highest amount of Li resources (2.5 × 1011 tone), but the average Li concentration is low at 0.17 mg/L, which means it has a concentration of 0.1–0.2 ppm; thus, it is not of technical relevance in the spotlight. However, Li sources—deposits that contain Li [15]—are estimated to be equal to 21,000,000 metric tons [16], and will be exhausted by 2080 [17,18]. The extraction of Li can be performed via different methods such as precipitation [19,20], calcination methods [21], adsorption [22], electro-dialysis [23,24,25], electrolysis [26,27], and nanofiltration membranes [28,29,30,31,32,33]. Because of its simplicity and cost effectiveness, the adsorption process is a potential technology for Li recovery from seawater and salt-lake brine. Inorganic Li ions sieve with a high selectivity, capacity, and stability, and have obtained a lot of interest.
Among the aforementioned methods, Li extraction in the nanofiltration membranes has been conducted via a combination of the two main extraction mechanisms, including size sieving and charge repulsion effects [34]. These nanoporous membranes attract massive attention due to their desirable advantages, including their scale-up possibility [18,35], low energy consumption, facile operation, and environmental friendliness [36]. However, the decisive factor which makes nanofiltration membranes suitable for Li extraction is the adaptability of the membrane method with ocean solutions, which contain high concentrations of competing ions (e.g., Na+, Mg2+, Ca2+, K+, and Li+) [18]. The low Li+ concentration in solution but the high concentration of the interfering ions (e.g., Mg2+, Ca2+, Na+, and K+) is a significant obstacle for the highly efficient and environmentally friendly extraction of Li. In the nanofiltration membrane, Li permeation is usually measured by the extraction ratio of Li+ and Mg2+, which is known as Li+/Mg2+ selectivity [34]. For porous membranes, membrane porosity determines permeability, whereas membrane pore size determines selectivity. High permeability and selectivity require both high porosity and proper pore sizes. For instance, in some discussions, molecular dynamics simulations revealed that Li+ mobility was enhanced, whereas the K+ and Na+ mobility was reduced. Consequently, a higher selectivity of Li+ with K+ and Na+ was obtained. Meanwhile, for membrane selection, a better understanding of the trade-off between membrane permeability and Li selectivity is expected [37]. Therefore, it is essential to design Li-Selective Membranes (LSMs) to tackle the challenges in the path of Li extraction. In this regard, different porous structures were proposed and investigated computationally and experimentally for utilizing Li-selective membranes. For instance, Wang et al. studied the Li extraction efficiency of the A-site La0.33NbO3 perovskite structure. Based on their results, per La0.33NbO3 structure, near 0.67 Li were inserted electrochemically. Additionally, the charge transfer resistance in the La0.33NbO3 structure has an inverse relationship with the diffusion coefficient due to the Li insertion into the La0.33NbO3 perovskite structure [38]. In terms of computational investigation, the Li conductivity of Li3OCl0.75Br0.25 was predicted using an Ab Initio Molecular Dynamics (AIMD) simulation [39]. However, among the proposed structures, Lithium Lanthanum Titanate (LLTO) perovskite, Li3xLa(2/3−x)(1/3−2x)TiO3 (0.04 < x < 0.16), attracts huge attention [40]. The plural A-site vacancies in its structure make it suitable for the seperation of the Li ions. Additionally, the bulk and grain-boundary conductivities of the LLTO structure at room temperature equal 10−3 S/cm and 10−4–10−5 S/cm, respectively [41], where the grain boundaries influence the Li+ diffusion through the LLTO samples. In the LLTO structure, a substantial number of “A” site vacancies are created by the disordered arrangement of Li and La ions, which facilitates Li ion transportation. To migrate from one “A” site to the next available site, Li+ ions have to pass through a bottleneck of four surrounding oxygen atoms [42]. Ionic conduction is caused by Li ion migration [43]. Many attempts have been made and discussed to increase ionic conductivity in the literature. Ionic conductivity (σ), an important property of bulk materials, defined as the velocity multiplied by the number of ions that are sequentially moving, is linked to Li migration on a micrometer length scale [44]. The Li diffusion pathways in the LLTO structure have been proposed, and the vacancies at the La site are important for Li+ migration [45]. High Li+ conductivity and selectivity in the extraction membranes are directly related to the crystal structure; in this respect, LLTO is known as one of the solid-state Li-ion superconductors [18]. In terms of the crystal structure, the LLTO structure contains TiO6 octahedra, which are interconnected to form cubic cages that can be suitable places for accommodating the Li+ and La3+ ions. The diameter of the created cages, which the four neighboring TiO6 tetrahedra have created, is 1.07 Å, and Li should be passed through them. However, by comparing the LLTO window diameter (1.07 Å) with the Li ion size (1.18 Å), it is obvious that the LLTO pore diameter is larger. Because of this miniature difference, small distortion and enlargement in the windows should happen [18].
In the current research, Molecular Dynamic (MD) simulations were conducted on the mechanism of Li extraction via the LLTO membrane to study the mechanism of the change in LLTO pores as well as the effect of the most dominant parameters (e.g., pH and temperature) on the LLTO pores expansion and Li extraction through them. If it is preferred to analyze the factors that could be used to distinguish Li+ ions from chemically identical monovalent and divalent ions in seawater, theoretical models and experimental verification reveal that the spacing, surface charge, functional group, and nanopore morphology are all significant for Li+ selectivity [46,47]. For this fact, in Li extraction from saline brines, the absolute separation efficiency depends on the relative selectivity of Li and other ions. The selective separation of Li+ and other ions requires molecular scale control of chemical/physical homogeneous nanomaterials. The lack of practical analytical methods within the nanoscale spatial and temporal resolution makes it difficult to fully comprehend the Li+ ion separation mechanism and pore geometry/chemistry interaction with Li+ ions. Although molecular simulations can provide insight into solute transport, the simulation capacity and membrane chemical heterogeneity preclude a complete depiction of ion transport through a realistic long pore. The results confirmed the undeniable role of pH and temperature on the LLTO pores size increment and facilitation of Li ions transmission through the membrane. Additionally, differences in hydration are important in explaining the extremely precise selectivity in cations. In the case of partial or complete dehydration during the process of transport, a more detailed approach is required to describe the interaction of water molecules with the ion as well as the interaction between the hydrated ion and the pore. The pore size directly influences ion transport, and this mechanism states that the energy barrier for ion transfer through small pores increases; hence, if water remains in these pores during the equilibrium, ion transport is the most hindered [46,47].

2. Simulation Methods

2.1. Study Plan

A simulation of the LLTO membrane with specific applications of Li extraction was performed in five different simulations, as outlined in Table 1.
Simulation 1, 2, and 3: Boxes in simulations 1 and 2 contain two LLTO structures and one Li ion. While in simulation 3, two LLTO membranes were placed in the middle of the simulation box surrounded by 500 ions (125 Li+, 125 K+, 125 Na+, and 125 Mg2+). In these three simulations, we aimed to investigate the effect of temperature and pH on the Li hydration, LLTO selectivity, and mechanisms of the LLTO pore changes. Finally, the interaction energies between the ions and membranes were scrutinized in more detail.
Simulation 4 and 5: In these cases, neither water nor other ions entered the simulation box; however, there were only two LLTO membranes in the vacuum condition. The objective of comparing the LLTO structure without water was to observe the impact of temperature on the expansion of membrane cavities. Moreover, the obtained results were also explained by utilizing molecular and atomic analyses, including interaction energies, hydration radii, etc.

2.2. Material and Methods

For the LLTO design, the LLTO’s CIF file was obtained from the Crystallography Open Database. Then, water molecules and ions were inserted into the simulation boxes with different concentrations and numbers using GROMACS software. The OPLSA force field and the tip3p water model were used to simulate water in these systems. The resulting structure was parameterized using the Polypargen or OBGMX sites and the GROMACS x2top command. Atom charges were calculated using CP2K software via ESP force fields in the Gaussian approach. The simulation box and LLTO structure were designed using the gmx editconf command and their sizes were equaled to 15 × 6.7 × 6.7 and 1.2 × 6.5 × 6.5 nm3, respectively.
The migration of the ions from the ionic part to the ion-free part in this system is caused by the concentration gradients, which act as a driving force; thus, the ions pass through the LLTO membrane. When the concentration changes, the bulk molecules suffer changes caused by the competition of three main forces: Coulomb interactions, due to the presence of an electrical charge; a steric interactions force, produced by the overlap of the molecules; and the elastic forces, due to entropic effects. These forces will be discussed in detail in the simulation 2 explanation. Considering that the membrane is a perovskite, it can suffer expansion according to its physicochemical properties. When increasing the concentration, an osmotic pressure appears and tries to separate them (steric interaction), and as a consequence, an elastic force appears because of the change in the conformational entropy. The electric repulsion produces an expansion of the pore size and also contributes to the overlapping between molecules. Under this behavior, the steric interaction and the elastic force control the pore size.
The simulation times were also considered 100 ns. The LINCS algorithm was also used to constrain the hydrogen bonds. In NVT and NPT simulations, Berendsen thermostat and barostat algorithms were used to stabilize the temperature at adjusted temperatures and pressure at 1 bar, respectively, with a simulation time of 600 ps. Additionally, the Nose–Hoover thermostat and Parrinello–Rahman barostat were applied in the MD simulation to stabilize the temperature and pressure during the 100 ns simulation, respectively. The simulation box was first minimized with a steep integrator. Then, it was equilibrated with the algorithms for 0.5 nanoseconds in the NVT and NPT ensembles. Finally, the simulations were performed for 100 ns. Coulomb energy relations and Lennard-Jones energy were used to calculate electrostatic energy and van der Waals (vdW) energy. The basis for calculating their energy was based on the average distance between the charges of the atoms and the geometric mean distance of the charges of the atoms from each other, respectively. The electric charge corresponds to each atom. K is also a constant of the Coulomb energy equation [48,49]. The average energy at each time was also obtained from the simulation. All the hydrogen-containing bonds were restrained using the LINCS algorithm [50].

3. Results and Discussions

The elucidation of ion transport through membranes is an important issue to improve membrane separation. Steric interaction, electrostatic, dielectric, vdW, frictional, and viscous interactions are the various mechanisms of transporting ions through the membranes. The steric effect poses a high energy barrier, particularly for rigid pores with dimensions smaller than the bare size of the ion. Electrostatic interactions and the intrinsic charges at the pore entrance are further affected by selectivity, which favors counter-ions diffusion but excludes co-ions. For a solute with a hydration diameter (dh) larger than the pore size (dp), the necessary dehydration (or partial dehydration) sets an energy barrier as well [46]. A dehydration reaction is a chemical reaction between two compounds where one of them is water that involves the loss of water from the reacting molecule or ion. The dielectric effect regulates the energy barrier during ion permeability once the ions become partially dehydrated and penetrate the pore, resulting in the selectivity of similar monovalent ions across membranes. After penetrating the pore, pore surface friction and viscous drag forces play an essential role in elucidating its mechanism, so that the affinity of the pore wall with the electrostatic and vdW forces affects solute hindering [46]. Designing a highly selective membrane requires a rational incorporation of the above mechanisms.
The current study addresses two aspects of ion transport through the LLTO membrane: the impact of pH and temperature in two main parts. The first part comprises three individual simulations to investigate the pH of the solution. Figure 1a–c depict visual highlights observed in the transportation of the single solvated Li+ ion through the LLTO membrane at pH 7 and 5.5.
Based on the results of the first simulation, we observed no difference in the LLTO membrane at room temperature and pH 7. Apart from the electrical interactions, the Grotthuss-like mechanism [51] describes the separation mechanism for the neutral media. The ions are conducted via the two-surface charge-governed transport mechanism and Grotthuss-like mechanism by jumping from one membrane side to another while moving in the direction of the electric field. Additionally, mass transfer is carried from the bulk to the boundary before transporting inside the membrane. Combining this theory with an appropriate concentration driving force model can determine the performance of membrane and selectivity. Based on the results, the membrane’s pore size confronted no variation and remained at 1.07 Å during the simulation (Figure 1b).
In the second simulation, the pH of the water solution was altered to 5.5 and the temperature of the simulation box was adjusted to 50 °C. Figure 1c displays the dehydration and permeation of Li+ ions through the LLTO membrane, which has expanded pores in the acidic pH. The snapshot shows the peeling of the water shell off the Li+ ion while passing through the enlarged pore. As mentioned earlier, the pore size of the LLTO membrane at pH 5.5 expanded to 1.18–1.2 Å. The miniature expansion (0.03–0.13 Å) in the size of the LLTO pores can be attributed to the change in both the temperature and pH value. The pore size is a function of the temperature and hydrogen bond interaction between ions and the membrane.
In the second simulation, we sought to analyze the impact of charge (acidic condition). As is well known, Li+ cations attract oxygen ions in the solutions while repelling hydrogen cations because they have the same charge. The lithium ions were dehydrated completely due to the simultaneous existence of repulsive and attractive forces between lithium ions and the LLTO membrane. Figure 1d shows that the interaction energies of Li+-H+ was repulsive, while the Li+-O2− interaction energies are shown in Figure 1e as being attractive forces. These simultaneous repulsive–attractive forces pull the competing ions back and forth, which can be deemed one of the leading reasons for the ions’ dehydration. In addition to the lithium ion dehydration, these repulsive–attractive interactions caused the LLTO membrane’s pore size to expand up to 1.18 to 1.2 Å, which can bring efficient lithium selectivity for the LLTO membrane.
We investigated the interaction of Li+ with the LLTO membrane after analyzing the interactions of O2− and H+ ions with Li+. In this regard, the surface charge of the LLTO membrane was calculated (Figure 1f); because of the negative charge of the LLTO membrane, positive lithium ions were attracted to it and were kept as close to the membrane as possible. The Li ions were thus affected by the forces of attraction from the LLTO membrane side. The Lithium ions moved closer to the membrane during this period and became dehydrated when positioned between the attraction–repulsion forces. As a result of the lithium ion dehydration, the possibility of lithium transferring across the membrane will be provided. Notably, increased temperature and acidic conditions expand the LLTO membrane’s pore size and make it easier for lithium ions to pass through the LLTO membrane.
In order to conduct further research, average hydrogen bond interactions between Li+ and the LLTO membrane along with the three simulations were quantified (Figure 2a) [52]. As can be seen, among the simulations 1, 2, and 3, the number of hydrogen bonds in the latter simulation increased the H bonds to 54 during the simulation, which can be attributed to both the pH and the existence of other ions as well as lithium in the simulation box. The slight difference in the number of hydrogen bonds between the second simulation (average H bond = 49) and third simulation can be related to the number of competing ions. The solution pH also plays a critical role in local and long-range ion−water interactions. It is thus reasonable to expect the solution pH to influence the structure and size of solvated ions. Hydrogen bonding interaction changed membrane pore parameters and covered the membrane surface. The presence of more ions in a water solution causes the formation of more hydrogen bonds in the system, increasing the total number of hydrogen bonds.
After filtration, the lithium ions formed almost no monohydrate (H2O) Li+. The low proportion of lithium monohydrates may be attributed to the abnormally low proportion of monohydrates in the bulk solution because of the high charge density of Li+ and high hydration energy. Regardless of the differences between these alkali metal ions, dehydration appears to be ubiquitous when the pore size is proportional to the hydrated size of the ion.
To provide deeper insight into the governed interaction energies, the first part of the simulations was compared, as seen in Figure 2b. The relative similarity in interaction energies for simulations 2 and 3 can be attributed to their similar conditions (50 °C and pH 5.5); however, due to the existence of more ions in the third simulation, it possesses a slightly higher difference in interaction energies compared with the second simulation. By comparing simulations 1 and 2, it can be firmly confirmed that the stronger interactions in simulation 2 were due to the higher temperature and acidic condition. Subsequently, the results approve the conditioned high Li+ selectivity in the third simulation. In other words, acidic conditions and higher temperatures provide a better situation to achieve higher Li selectivity over other existing ions (Figure 2c).
In the third simulation, two LLTO membrane structures were surrounded by 500 ions (125 Li+, 125 K+, 125 Na+, and 125 Mg2+) in the solution condition of pH 5.5 at 50 °C. Figure 2d reveals the effect of the presence of other ions in the separation of Li+ ions to visually represent the ratio of ions that passed the barrier.
The interaction energies between the Li+ and LLTO membranes in the presence of other cations (simulation 3) were investigated (Figure 2e). The negative sign of the interaction energies approves the existence of the attractive forces among the Li+, Na+, K+, and Mg2+ ions and the LLTO structure. As can be noticed, the vdW interactions had relatively minor contributions to the total energies, and its negligible value stemmed from the small radii of the competing ions. Thus, the constituent part of the total energies was formed by the electrostatic interaction energies, consisting of the ion–ion, water–ion, and water–water interactions energies. Electrostatic energies are more determinative in the membrane selectivity, in which more differences in the electrostatic interaction energies lead to higher selectivity and better membrane performance. The total energy of Li+, Na+, K+, and Mg2+ were −790 kJ/mol, −850 kJ/mol, −940 kJ/mol, and −1040 kJ/mol, respectively. By comparing the ion’s total interaction energies, it is obvious that they were higher than the hydration energies of Li+, Na+, K+, and Mg2+, which are equal to −515, −405, −321, and −1922 kJ/mol [53]. Suppose the trapping energy or total energy is stronger than the hydration energy of each ion. In that case, they can overcome the ion’s hydration energies, and as a consequence, the ions will be dehydrated. The obtained results showed that only Na+, K+, and Li+ could be fully dehydrated among present ions, and Mg2+ ions were partially dehydrated.
The fully dehydrated diameter of Li+ (0.136 nm) was only slightly larger than that of Mg2+ (0.13 nm), but was comparably smaller than that of Na+ (0.19 nm) and K+ (0.266 nm). The Na+ and K+ dehydrated diameters were larger than the LLTO membrane pore size, and the Mg2+ ions were partially dehydrated and consequently could not pass through the membrane (Figure 3). Li+ had a whole situation to pass through the membrane. The results exhibited the remarkable Li selectivity of the LLTO membrane against other ions. Nanofiltration under the various conditions indicated that dehydration was not merely driven by the complete exclusion of larger Mg2+ hydrates from the membrane; rather, a partial dehydration-driven transformation was occurring. In other words, the large Mg2+ hydrates were reducing their effective size by shedding part of their solvation shell to partition into and permeate across the nanometer membrane pores. In this case, the synergic effect of the temperature increment and pH decrement can help the selectivity of the LLTO membrane by widening its pores and effective ion dehydration. The difficulties in separating Mg2+ from Li+ lie in their similar ionic radii, either in bare or hydration states (Table 2). The fully dehydrated radius of Li+ (0.068 nm) was only slightly larger than that of Mg2+ (0.065 nm), but the hydration energy of Mg2+ was comparably larger than the trapping energy. Therefore, the Mg2+ ions could not pass through the LLTO membrane. Additionally, the hydration number of Li+ is five. Still, the lifetime or exchange rate of the water molecules in the hydration layer was 5 × 10−9 s, while Mg2+ had a hydration number of six with a significantly much slower exchange rate of 10−6 s. Because the extraction process took place at the aqueous/organic interface, the ions needed to overcome an energy barrier to cross the interface [46].
In the fourth and fifth simulations, the LLTO structure was placed in the simulation box without any ions or solvent. In these cases, only temperature changes from 25 °C to 50 °C occurred in simulations 4 and 5, respectively. As it is well known, temperature is an important factor governing mass transfer in membrane separation processes and has a significant impact on NF membrane performance. By comparing the results of these two simulations, the slight effect of the temperature change on the pore expansion was firmly confirmed (Figure 3). At 25 °C, the pore size remained at 1.07 nm (Figure 3a), while at 50 °C, the pore size in the LLTO membrane increased to 1.1 Å (Figure 3b). It can be seen that the pore size changed depending on temperature.

4. Conclusions

In the case of ion separation, developing membranes with nanopores tailored to match the ion solvation structure could significantly enhance ion−ion selectivity and lead to increased energy and process efficiencies. The novelty in this study lies in systematically evaluating the separation and dehydration mechanism based on the desired driving force during pore transport for small ions with Molecular Dynamics (MD). The LLTO membrane, which was employed as a high-performance membrane in the Li extraction application, was simulated in different environmental conditions (e.g., temperature and pH) via molecular dynamics. The results of this research indicated that three main factors have a significant positive effect on the performance of the Li extraction of LLTO membranes, including acidic pH conditions; repulsion–attraction forces, which led to the ions’ dehydration; and finally the efficient affiliation between the LLTO structure and Li ions, which led to a higher rate of Li permission to the LLTO membrane due to the pore size expansion from 1.07 Å to approx. 1.2 Å. Our simulations therefore suggest that ion dehydration is important to transport when the size of the nanopore is smaller than the hydrated ion, leading to water molecule shedding from the outermost solvation layer. Filtration by LLTO membranes induced ion dehydration for all the various alkali metal ions (Li+, Na+, and K+) in the bulk solution investigated. The shedding of the water molecules from the solvation shell allowed ions to enter pores that were smaller than their bulk hydrated size. With the current MD simulation, the dehydration of ions was found to mainly dominate the ion selectivity of the pores, and the difficulty of the ion dehydration could be reflected by the ion energy barrier. According to the results, monovalent cations can pass more easily through the nanopores than the divalent cations due to their smaller energy barriers, indicating that the pores are more effective in separating these two kinds of cations. As a final point, three important conclusions should be mentioned: Firstly, energetic barriers were strongly dependent on pore size. Additionally, there were distinct regimes related to the required dehydration. Transport was strongly hindered when the size of the pore was smaller than the hydrated radius. Second, energy barriers were determined by ion type (and thus hydration properties), and ion selectivity varied depending on pore size. Finally, and perhaps most crucially, dehydration was the primary method of ion transfer in the pores. Even when charge repulsion was taken into account, partial dehydration was the dominant predictor of the energy barriers for small, strongly hydrated ions whose hydrated radius was bigger than the pore size.

Author Contributions

Conceptualization, M.G., M.H., S.R., Y.M.C. and R.M.; Formal analysis, M.G. and M.H.; Methodology, M.H.; Software, M.G.; Supervision, A.R.; Writing—original draft, M.H. and S.R.; Writing—review & editing, M.G., Y.M.C. and R.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ajayi, O.S.; Oluyamo, S.S.; Ofiwe, C.U.; Aborisade, C.A. Zinc Oxide Nanoparticle in Lithium Triborate Microparticle System: Visibility for Application in Dosimetry. Mater. Today Proc. 2021, 38, 879–886. [Google Scholar] [CrossRef]
  2. Saray, A.A.; Kaviani, P.; Shahbazi-Gahrouei, D. Dosimetric Characteristics of Lithium Triborate (LiB3O5) Nanophosphor for Medical Applications. Radiat. Meas. 2021, 140, 106502. [Google Scholar] [CrossRef]
  3. Chen, T.; Jin, Y.; Lv, H.; Yang, A.; Liu, M.; Chen, B.; Xie, Y.; Chen, Q. Applications of Lithium-Ion Batteries in Grid-Scale Energy Storage Systems. Trans. Tianjin Univ. 2020, 26, 208–217. [Google Scholar] [CrossRef] [Green Version]
  4. Marín, P.E.; Milian, Y.; Ushak, S.; Cabeza, L.F.; Grágeda, M.; Shire, G.S.F. Lithium Compounds for Thermochemical Energy Storage: A State-of-the-Art Review and Future Trends. Renew. Sustain. Energy Rev. 2021, 149, 111381. [Google Scholar] [CrossRef]
  5. Cabeza, L.F.; Gutierrez, A.; Barreneche, C.; Ushak, S.; Fernández, Á.G.; Inés Fernádez, A.; Grágeda, M. Lithium in Thermal Energy Storage: A State-of-the-Art Review. Renew. Sustain. Energy Rev. 2015, 42, 1106–1112. [Google Scholar] [CrossRef] [Green Version]
  6. Wu, J.; Yuan, L.; Li, Z.; Xie, X.; Huang, Y. Air-Stable Means More: Designing Air-Defendable Lithium Metals for Safe and Stable Batteries. Mater. Horizons 2020, 7, 2619–2634. [Google Scholar] [CrossRef]
  7. Menazea, A.A.; Ismail, A.M.; Elashmawi, I.S. The Role of Li4Ti5O12 Nanoparticles on Enhancement the Performance of PVDF/PVK Blend for Lithium-Ion Batteries. J. Mater. Res. Technol. 2020, 9, 5689–5698. [Google Scholar] [CrossRef]
  8. Sun, Y.; Guan, P.; Liu, Y.; Xu, H.; Li, S.; Chu, D. Recent Progress in Lithium Lanthanum Titanate Electrolyte towards All Solid-State Lithium Ion Secondary Battery. Crit. Rev. Solid State Mater. Sci. 2019, 44, 265–282. [Google Scholar] [CrossRef]
  9. Siyal, S.H.; Javed, M.S.; Jatoi, A.H.; Lan, J.L.; Yu, Y.; Saleem, M.; Yang, X. In Situ Curing Technology for Dual Ceramic Composed by Organic–Inorganic Functional Polymer Gel Electrolyte for Dendritic-Free and Robust Lithium–Metal Batteries. Adv. Mater. Interfaces 2020, 7, 2000830. [Google Scholar] [CrossRef]
  10. Zhu, Y.; Gonzalez-Rosillo, J.C.; Balaish, M.; Hood, Z.D.; Kim, K.J.; Rupp, J.L.M. Lithium-Film Ceramics for Solid-State Lithionic Devices. Nat. Rev. Mater. 2020, 6, 313–331. [Google Scholar] [CrossRef]
  11. Hartmann, M.; Pfaffinger, M.; Stampfl, J.; Stampfl, J. The Role of Solvents in Lithography-Based Ceramic Manufacturing of Lithium Disilicate Cell-Compatible Hydrogels for Two-Photon Polymerization View Project Hydrogels for Two Photon Polymerization View Project The Role of Solvents in Lithography-Based Ceramic Manufacturing of Lithium Disilicate. Materials 2021, 14, 1045. [Google Scholar] [CrossRef]
  12. Mateus, R.; Costa, M.B.; Alves, L.C.; Guedes, M.; Alves, E.; Ferro, A.C. Lithium Dilution in Li-Sn Alloys. Nucl. Mater. Energy 2020, 25, 100783. [Google Scholar] [CrossRef]
  13. Liang, C.; Gao, M.; Pan, H.; Liu, Y.; Yan, M. Lithium Alloys and Metal Oxides as High-Capacity Anode Materials for Lithium-Ion Batteries. J. Alloys Compd. 2013, 575, 246–256. [Google Scholar] [CrossRef]
  14. Nithyadharseni, P.; Reddy, M.V.; Nalini, B.; Kalpana, M.; Chowdari, B.V.R. Sn-Based Intermetallic Alloy Anode Materials for the Application of Lithium Ion Batteries. Electrochim. Acta 2015, 161, 261–268. [Google Scholar] [CrossRef]
  15. Stringfellow, W.T.; Dobson, P.F. Technology for the Recovery of Lithium from Geothermal Brines. Energies 2021, 14, 6805. [Google Scholar] [CrossRef]
  16. U.S. Geological Survey. Mineral Commodity Summaries 2021; Geological Survey: Reston, VA, USA, 2021; 200p. [Google Scholar] [CrossRef]
  17. Yang, S.; Zhang, F.; Ding, H.; He, P.; Zhou, H. Lithium Metal Extraction from Seawater. Joule 2018, 2, 1648–1651. [Google Scholar] [CrossRef] [Green Version]
  18. Li, Z.; Li, C.; Liu, X.; Cao, L.; Li, P.; Wei, R.; Li, X.; Guo, D.; Huang, K.W.; Lai, Z. Continuous Electrical Pumping Membrane Process for Seawater Lithium Mining. Energy Environ. Sci. 2021, 14, 3152–3159. [Google Scholar] [CrossRef]
  19. Lai, X.; Xiong, P.; Zhong, H. Extraction of Lithium from Brines with High Mg/Li Ratio by the Crystallization-Precipitation Method. Hydrometallurgy 2020, 192, 105252. [Google Scholar] [CrossRef]
  20. Zhang, Y.; Hu, Y.; Wang, L.; Sun, W. Systematic Review of Lithium Extraction from Salt-Lake Brines via Precipitation Approaches. Miner. Eng. 2019, 139, 105868. [Google Scholar] [CrossRef]
  21. Kuai, Y.; Yao, W.; Ma, H.; Liu, M.; Gao, Y.; Guo, R. Recovery Lithium and Potassium from Lepidolite via Potash Calcination-Leaching Process. Miner. Eng. 2021, 160, 106643. [Google Scholar] [CrossRef]
  22. Wang, S.; Chen, X.; Zhang, Y.; Zhang, Y.; Zheng, S. Lithium Adsorption from Brine by Iron-Doped Titanium Lithium Ion Sieves. Particuology 2018, 41, 40–47. [Google Scholar] [CrossRef]
  23. Iizuka, A.; Yamashita, Y.; Nagasawa, H.; Yamasaki, A.; Yanagisawa, Y. Separation of Lithium and Cobalt from Waste Lithium-Ion Batteries via Bipolar Membrane Electrodialysis Coupled with Chelation. Sep. Purif. Technol. 2013, 113, 33–41. [Google Scholar] [CrossRef]
  24. Ji, Z.Y.; Chen, Q.B.; Yuan, J.S.; Liu, J.; Zhao, Y.Y.; Feng, W.X. Preliminary Study on Recovering Lithium from High Mg2+/Li+ Ratio Brines by Electrodialysis. Sep. Purif. Technol. 2017, 172, 168–177. [Google Scholar] [CrossRef]
  25. Song, Y.; Zhao, Z. Recovery of Lithium from Spent Lithium-Ion Batteries Using Precipitation and Electrodialysis Techniques. Sep. Purif. Technol. 2018, 206, 335–342. [Google Scholar] [CrossRef]
  26. Lv, W.; Ruan, D.; Zheng, X.; Li, L.; Cao, H.; Wang, Z.; Zhang, Y.; Sun, Z. One-Step Recovery of Valuable Metals from Spent Lithium-Ion Batteries and Synthesis of Persulfate through Paired Electrolysis. Chem. Eng. J. 2021, 421, 129908. [Google Scholar] [CrossRef]
  27. Liu, X.; Chen, X.; He, L.; Zhao, Z. Study on Extraction of Lithium from Salt Lake Brine by Membrane Electrolysis. Desalination 2015, 376, 35–40. [Google Scholar] [CrossRef]
  28. Li, Q.; Liu, H.; He, B.; Shi, W.; Ji, Y.; Cui, Z.; Yan, F.; Mohammad, Y.; Li, J. Ultrahigh-Efficient Separation of Mg2+/Li+ Using an in-Situ Reconstructed Positively Charged Nanofiltration Membrane under an Electric Field. J. Memb. Sci. 2022, 641, 119880. [Google Scholar] [CrossRef]
  29. Luo, H.; Peng, H.; Zhao, Q. High Flux Mg2+/Li+ Nanofiltration Membranes Prepared by Surface Modification of Polyethylenimine Thin Film Composite Membranes. Appl. Surf. Sci. 2022, 579, 152161. [Google Scholar] [CrossRef]
  30. Yang, Z.; Fang, W.; Wang, Z.; Zhang, R.; Zhu, Y.; Jin, J. Dual-Skin Layer Nanofiltration Membranes for Highly Selective Li+/Mg2+ Separation. J. Memb. Sci. 2021, 620, 118862. [Google Scholar] [CrossRef]
  31. Feng, Y.; Peng, H.; Zhao, Q. Fabrication of High Performance Mg2+/Li+ Nanofiltration Membranes by Surface Grafting of Quaternized Bipyridine. Sep. Purif. Technol. 2022, 280, 119848. [Google Scholar] [CrossRef]
  32. He, R.; Dong, C.; Xu, S.; Liu, C.; Zhao, S.; He, T. Unprecedented Mg2+/Li+ Separation Using Layer-by-Layer Based Nanofiltration Hollow Fiber Membranes. Desalination 2022, 525, 115492. [Google Scholar] [CrossRef]
  33. Wu, M.-B.; Ye, H.; Zhu, Z.-Y.; Chen, G.-T.; Ma, L.-L.; Liu, S.-C.; Liu, L.; Yao, J.; Xu, Z.-K. Positively-Charged Nanofiltration Membranes Constructed via Gas/Liquid Interfacial Polymerization for Mg2+/Li+ Separation. J. Memb. Sci. 2022, 644, 119942. [Google Scholar] [CrossRef]
  34. Peng, H.; Zhao, Q. A Nano-Heterogeneous Membrane for Efficient Separation of Lithium from High Magnesium/Lithium Ratio Brine. Adv. Funct. Mater. 2021, 31, 1–7. [Google Scholar] [CrossRef]
  35. Bashiri, A.; Nikzad, A.; Maleki, R.; Asadnia, M.; Razmjou, A. Rare Earth Elements Recovery Using Selective Membranes via Extraction and Rejection. Membranes 2022, 12, 80. [Google Scholar] [CrossRef] [PubMed]
  36. Xu, F.; Dai, L.; Wu, Y.; Xu, Z. Li+/Mg2+ Separation by Membrane Separation: The Role of the Compensatory Effect. J. Memb. Sci. 2021, 636, 119542. [Google Scholar] [CrossRef]
  37. Li, X.; Mo, Y.; Qing, W.; Shao, S.; Tang, C.Y.; Li, J. Membrane-Based Technologies for Lithium Recovery from Water Lithium Resources: A Review. J. Memb. Sci. 2019, 591, 117317. [Google Scholar] [CrossRef]
  38. Wang, G.X.; Yao, P.; Bradhurst, D.H.; Dou, S.X.; Liu, H.K. Structure Characterisation and Lithium Insertion in La0.33NbO3 Perovskite. Solid State Ionics 1999, 124, 37–43. [Google Scholar] [CrossRef]
  39. Deng, Z.; Radhakrishnan, B.; Ong, S.P. Rational Composition Optimization of the Lithium-Rich Li3OCl1−XBrx Anti-Perovskite Superionic Conductors. Chem. Mater. 2015, 27, 3749–3755. [Google Scholar] [CrossRef]
  40. Inaguma, Y.; Liquan, C.; Itoh, M.; Nakamura, T.; Uchida, T.; Ikuta, H.; Wakihara, M. High Ionic Conductivity in Lithium Lanthanum Titanate. Solid State Commun. 1993, 86, 689–693. [Google Scholar] [CrossRef]
  41. Chen, C.H.; Amine, K. Ionic Conductivity, Lithium Insertion and Extraction of Lanthanum Lithium Titanate. Solid State Ionics 2001, 144, 51–57. [Google Scholar] [CrossRef]
  42. Abhilash, K.P.; Selvin, P.C.; Nalini, B.; Nithyadharseni, P.; Pillai, B.C. Investigations on Pure and Ag Doped Lithium Lanthanum Titanate (LLTO) Nanocrystalline Ceramic Electrolytes for Rechargeable Lithium-Ion Batteries. Ceram. Int. 2013, 39, 947–952. [Google Scholar] [CrossRef]
  43. Ruiz, A.I.; López, M.L.; Veiga, M.L.; Pico, C. Electrical Properties of La1.33−xLi3xTi2O6 (0.1 < x < 0.3). Solid State Ionics 1998, 112, 291–297. [Google Scholar] [CrossRef]
  44. Hayamizu, K.; Seki, S.; Haishi, T. 7Li NMR Diffusion Studies in Micrometre-Space for Perovskite-Type Li0.33La0.55TiO3 (LLTO) Influenced by Grain Boundaries. Solid State Ionics 2018, 326, 37–47. [Google Scholar] [CrossRef]
  45. Yashima, M.; Itoh, M.; Inaguma, Y.; Morii, Y. Crystal Structure and Diffusion Path in the Fast Lithium-Ion Conductor La0.62Li0.16TiO3. J. Am. Chem. Soc. 2005, 127, 3491–3495. [Google Scholar] [CrossRef] [PubMed]
  46. Xu, S.; Song, J.; Bi, Q.; Chen, Q.; Zhang, W.M.; Qian, Z.; Zhang, L.; Xu, S.; Tang, N.; He, T. Extraction of Lithium from Chinese Salt-Lake Brines by Membranes: Design and Practice. J. Memb. Sci. 2021, 635, 119441. [Google Scholar] [CrossRef]
  47. Razmjou, A.; Hosseini, E.; Cha-Umpong, W.; Korayem, A.H.; Asadnia, M.; Moazzam, P.; Orooji, Y.; Karimi-Maleh, H.; Chen, V. Effect of Chemistry and Geometry of GO Nanochannels on the Li Ion Selectivity and Recovery. Desalination 2020, 496, 114729. [Google Scholar] [CrossRef]
  48. Kumar, G.; Mishra, R.R.; Verma, A. Introduction to Molecular Dynamics Simulations. In Lecture Notes in Applied and Computational Mechanics; Springer: Singapore, 2022; Volume 99, pp. 1–19. ISBN 3000126414. [Google Scholar]
  49. Mehdizadeh Chellehbari, Y.; Sayyad Amin, J.; Zendehboudi, S. How Does a Microfluidic Platform Tune the Morphological Properties of Polybenzimidazole Nanoparticles? J. Phys. Chem. B 2022, 126, 308–326. [Google Scholar] [CrossRef]
  50. Hess, B.; Bekker, H.; Berendsen, H.J.C.; Fraaije, J.G.E.M. LINCS: A Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 14631472. [Google Scholar] [CrossRef]
  51. Zakertabrizi, M.; Hosseini, E.; Korayem, A.H.; Razmjou, A.; Fane, A.G.; Chen, V. Insight from Perfectly Selective and Ultrafast Proton Transport through Anhydrous Asymmetrical Graphene Oxide Membranes under Grotthuss Mechanism. J. Memb. Sci. 2021, 618, 118735. [Google Scholar] [CrossRef]
  52. Hydrogen Bond—GROMACS 2022.3 Documentation. Available online: https://manual.gromacs.org/current/reference-manual/analysis/hydrogen-bonds.html (accessed on 11 October 2022).
  53. Smith, D.W. Ionic Hydration Enthalpies. J. Chem. Educ. 1977, 54, 540–542. [Google Scholar] [CrossRef]
Figure 1. Ion transportation through the LLTO membrane: (a) initial conformation of solvated Li+ ion in the simulation box containing two LLTO membranes; (b,c) dehydration and permeation of Li+ ion at pH 7 and 5.5, respectively. Additional insets provide details on the effect of pH and temperature on the pore size; (d,e) energy analysis of Li+-H+ and Li+-O2− interactions with the LLTO membrane, respectively; and (f) LLTO membrane’s surface charge.
Figure 1. Ion transportation through the LLTO membrane: (a) initial conformation of solvated Li+ ion in the simulation box containing two LLTO membranes; (b,c) dehydration and permeation of Li+ ion at pH 7 and 5.5, respectively. Additional insets provide details on the effect of pH and temperature on the pore size; (d,e) energy analysis of Li+-H+ and Li+-O2− interactions with the LLTO membrane, respectively; and (f) LLTO membrane’s surface charge.
Membranes 12 01042 g001
Figure 2. Li transportation through LLTO membrane in the presence of other cations: (a) average H bond between Li+ ion and LLTO membrane for the first part of simulations; (b) interaction energy between Li+ ion and LLTO membrane for the first part of simulations. Obviously, Li+ ion has the strongest interactions with the LLTO membrane in the 3rd simulation due to the pH and temperature conditions; (c) Gibbs free energy for three simulations reveals that the 3rd simulation displays the most stable interaction between Li+ ion and LLTO membrane; (d) transportation of ions through the membranes in the presence of 500 ions. Closer snapshots of both sides clearly show the high selectivity of the membrane toward Li+ ion in comparison with other ions (K, Mg, and Na); and (e) interaction energy of cations with the LLTO membrane in the 3rd simulation.
Figure 2. Li transportation through LLTO membrane in the presence of other cations: (a) average H bond between Li+ ion and LLTO membrane for the first part of simulations; (b) interaction energy between Li+ ion and LLTO membrane for the first part of simulations. Obviously, Li+ ion has the strongest interactions with the LLTO membrane in the 3rd simulation due to the pH and temperature conditions; (c) Gibbs free energy for three simulations reveals that the 3rd simulation displays the most stable interaction between Li+ ion and LLTO membrane; (d) transportation of ions through the membranes in the presence of 500 ions. Closer snapshots of both sides clearly show the high selectivity of the membrane toward Li+ ion in comparison with other ions (K, Mg, and Na); and (e) interaction energy of cations with the LLTO membrane in the 3rd simulation.
Membranes 12 01042 g002
Figure 3. The impact of temperature on the pore size of LLTO membrane: (a) T = 298 K, (b) T = 323 K.
Figure 3. The impact of temperature on the pore size of LLTO membrane: (a) T = 298 K, (b) T = 323 K.
Membranes 12 01042 g003
Table 1. Characteristics of all simulations.
Table 1. Characteristics of all simulations.
PartSimulation No.Temp. (K)pHSolventNum. Of Ions/Ion Type
1st12987Water1 Li+
23235.5Water1 Li+
33235.5Water125 Li+, 125 K+, 125 Na+, and 125 Mg2+
2nd4298--0
5323--0
Table 2. Characteristics of the common cations and anions in the salt-lake brine.
Table 2. Characteristics of the common cations and anions in the salt-lake brine.
IonBare Ion RadiusHydrated RadiusHydration NumberDiffusion CoefficientHydration Free Energy
[nm][nm]±1(10−9 m2 s−1 )[kJ. mol−1]
H3O+-0.2839.31-
Li+0.0680.3851.03474
Na+0.0950.3641.33364
K+0.1330.3331.96295
Mg+0.0650.4360.711828
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Golmohammadi, M.; Habibi, M.; Rezvantalab, S.; Mehdizadeh Chellehbari, Y.; Maleki, R.; Razmjou, A. Mechanism Understanding of Li-ion Separation Using A Perovskite-Based Membrane. Membranes 2022, 12, 1042. https://doi.org/10.3390/membranes12111042

AMA Style

Golmohammadi M, Habibi M, Rezvantalab S, Mehdizadeh Chellehbari Y, Maleki R, Razmjou A. Mechanism Understanding of Li-ion Separation Using A Perovskite-Based Membrane. Membranes. 2022; 12(11):1042. https://doi.org/10.3390/membranes12111042

Chicago/Turabian Style

Golmohammadi, Mahsa, Meysam Habibi, Sima Rezvantalab, Yasin Mehdizadeh Chellehbari, Reza Maleki, and Amir Razmjou. 2022. "Mechanism Understanding of Li-ion Separation Using A Perovskite-Based Membrane" Membranes 12, no. 11: 1042. https://doi.org/10.3390/membranes12111042

APA Style

Golmohammadi, M., Habibi, M., Rezvantalab, S., Mehdizadeh Chellehbari, Y., Maleki, R., & Razmjou, A. (2022). Mechanism Understanding of Li-ion Separation Using A Perovskite-Based Membrane. Membranes, 12(11), 1042. https://doi.org/10.3390/membranes12111042

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