Next Article in Journal
Phosphorus Concentration in Water Affects the Biofilm Community and the Produced Amount of Extracellular Polymeric Substances in Reverse Osmosis Membrane Systems
Previous Article in Journal
Development of Pure Silica CHA Membranes for CO2 Separation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Membranes
Volume 11
Issue 12
10.3390/membranes11120927
membranes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient Helium Separation with Two-Dimensional Metal–Organic Framework Fe/Ni-PTC: A Theoretical Study

by
Jingyuan Wang
1,
Yixiang Li
1,
Yanmei Yang
2,
Yongqiang Li
1,
Mingwen Zhao
1,
Weifeng Li
1,
Jing Guan
1,* and
Yuanyuan Qu
1,*
1
School of Physics, Shandong University, Jinan 250100, China
2
College of Chemistry, Chemical Engineering and Materials Science, Shandong Normal University, Jinan 250014, China
*
Authors to whom correspondence should be addressed.
Membranes 2021, 11(12), 927; https://doi.org/10.3390/membranes11120927
Submission received: 1 November 2021 / Revised: 19 November 2021 / Accepted: 23 November 2021 / Published: 26 November 2021
Browse Figures
Review Reports Versions Notes

Abstract

:
Helium (He) is one of the indispensable and rare strategic materials for national defense and high-tech industries. However, daunting challenges have to be overcome for the supply shortage of He resources. Benefitted from the wide pore size distribution, sufficient intrinsic porosity, and high specific surface area, metal–organic framework (MOF) materials are prospective candidates for He purification in the membrane-based separation technology. In this work, through first-principles calculations and molecular dynamics (MD) simulations, we studied the permeability and filtration performance of He by the newly synthesized two-dimensional Fe-PTC MOF and its analogue Ni-PTC MOF. We found that both Fe-PTC and Ni-PTC have superior high performance for He separation. The selectivity of He over N2 was calculated to be ~1017 for Fe-PTC and ~1015 for Ni-PTC, respectively, both higher than most of the previously proposed 2D porous membranes. Meanwhile, high He permeance (10−4~10−3 mol s−1 m−2 Pa−1) can be obtained for the Fe/Ni-PTC MOF for temperatures ranging from 200 to 500 K. Therefore, the present study offers a highly prospective membrane for He separation, which has great potential in industrial application.

1. Introduction

As a rare strategic material, helium (He) is urgently needed in the medical, scientific research, and aerospace industries [1,2,3]. Albeit being the second most abundant element in the universe, He is still scarce in the atmosphere. Currently, the main source of He is the by-product in the exploitation of natural gas [4,5]. Purification of He with membrane-based separation technology has become a promising method, because of its low energy consumption, eco-friendly nature, simple operation, and other advantages [6,7,8,9,10].
The recent quick advances in the development of nanoporous two-dimensional (2D) structures have provided superior membrane candidates for He separation, in the merit of their intrinsic porosity, atomic thickness, and high specific surface area. For instance, graphene [9], silicene [11], and MoS2 [12] with drilling holes have been proposed for He separation and demonstrated to be successful. However, the requirement of precise treatment to these materials to introduce properly-sized hole structure becomes a technique hindrance, thus 2D structures with intrinsically distributed nanopores have attracted more attention which are highly expected to be efficient filtration membranes. For example, in 2015, the experimentally available g-C3N4 membrane has been demonstrated to possess high selectivity (107~1065 at 300 K) for separating He from impure gas molecules (H2, N2, CO, CH4, Ne, and Ar) in natural gas [13]. Then, Wang et al. theoretically predicted high He selectivity (102~1038 at 300 K) and permeance (1.0 × 103 GPU) of the monolayer CTF-0 membrane for He treatment [14]. Very recently, Liu et al. has shown that the nano-porous g-C2O membrane is a promising membrane for He separation with high selectivity (30~1028 at 300 K) and permeance (1.03 × 107 GPU) [15]. In addition, other intrinsically porous 2D materials, including C2N [16], IGP [17] and graphenylenee-1 [18], have also shown good performance for He separation. It should be noted that the trade-off between the selectivity and permeance is inherent in the system with porous structure, as large pores are favored by high permeance and small pores are advantageous to the selectivity of small molecules like He. Therefore, searching for a membrane with properly-sized pores for He separation is of great importance to achieve the goal of optimized balance of high permeance with sufficiently high selectivity.
Metal–organic frameworks (MOFs) are a new class of nanoporous materials, which is constituted of metal nodes connected by organic linkers [19,20]. Due to the high diversity of functional groups and ligands, the structure of MOFs can be tailored during the synthesis process, resulting in a large variety of structures with diverse geometries, pore sizes, and functional groups [21]. To date, a great deal of MOFs has been reported with excellent properties, such as high pore volume, large surface area, wide range of pore size from micro to meso scale, and reasonable thermal and mechanical stability. Based on these properties, application studies of MOFs for gas storage and separation [22,23,24,25,26], drug delivery [27], and catalysis [25,28] are widely explored. In particular, because of the rich porous structure, many MOFs have been proposed to be effective membrane materials in gas separation. Ranjan et al. [29] showed that microporous metal–organic framework (MMOF) membranes exhibit industrially acceptable selectivity of 23 for H2/N2 separation at 363 K. Cao et al. [30] demonstrated that the Cu−BTC membrane can be used for He separation with moderate selectivities and permeance for He over CO2, N2, and CH4, comparable to previously reported MOF membranes. In 2018, a highly crystalline ZIF-68 membrane was synthesized by Kasik et al., which offered a high He permeance of 2.6 × 10−7 mol s−1 m−2 Pa−1 [31]. The above-mentioned reports are all based on three-dimensional (3D) MOF structures, and relative application studies of 2D MOFs in gas separation are rare.
In 2018, Dong et al. successfully synthesized the black polycrystalline 2D Fe-PTC MOF by solvothermal method, which exhibits a high electrical conductivity of 10 S cm−1 (300 K) and a ferromagnetism (<20 K), implying its potential as a ferromagnetic semiconductor for application in spintronics [32]. It is noticed that the pore size of the Fe-PTC MOF is about 5.65 Å [33], similar to that of IGP (5.47 Å) [17] which has good performance on He separation. In addition, by replacing the transition metal atoms in the Fe-PTC MOF, a group of analogues can be obtained with similar pore sizes [33]. Given their atomic thickness and proper pore sizes, the He purification of this group of 2D MOF is highly expected.
In this work, using combined first-principles calculations and MD simulations, we theoretically demonstrated that the 2D Fe-PTC MOF and its analogues Ni-PTC MOF can efficiently separate He from natural gas and inert gases. The He permeance of Fe/Ni-PTC MOF reached ~6 × 10−4 mol s−1 m−2 Pa−1. It is worth mentioning that this permeability is five orders of magnitude higher than the industrially acceptable value (6.9 × 10−9 mol s−1 m−2 Pa−1), indicating its potential in industrial application. Moreover, the He/N2 selectivity is estimated to be 1015~1017, higher than most of the previously proposed 2D porous membranes. Hence, our results provide a new type of competitive He separation membranes with promising industrial applications.

2. Computational Methods

The Vienna ab initio simulation program (VASP) in the framework of density functional theory (DFT) is used to perform all first-principles calculations [34,35]. The projected-augmented wave (PAW) potential [36] was employed to describe the electron–ion interaction. The Perdew–Burke–Ernzerhof (PBE) functional under the generalized gradient approximation (GGA) was used to describe the electron–electron interaction [37]. The weak van der Waals interaction between the gas molecules and the Fe/Ni-PTC membrane was corrected by the Grimme method (DFT−D2) [38]. The energy cutoff was set to be 500 eV and iterative convergences of energy and force were 10−5 eV and 10−2 eV/Å. In order to prevent the interaction between adjacent images, a 20 Å vacuum zone was applied in the z direction. A k-point mesh of 5 × 5 × 1 was chosen for sampling the Brillouin zone for the unit cell [39]. The convergence of the system energy with respect to energy cutoff and the k-point mesh has been tested and shown in Table S1. For the pristine membrane, fully structural optimization was performed. For the adsorption interaction between gases and membranes, the cell vectors were fixed as negligible influence (membrane areas or adsorption energies) was observed upon gas adsorption (Table S2). Meanwhile, the z-coordinates of the membranes were frozen to maintain a planar surface, while the x- and y-coordinates of the membrane and the coordinates of the gas molecules were fully relaxed. The cohesive energies of these two MOF membranes and the adsorption energies of gas molecules adsorbed on a membrane are all obtained from single point energy calculation. The climbing image nudged elastic band (cNEB) method was adopted for searching the minimum energy pathway of the gas molecule passing through the 2D MOF membrane [40]. In the transition state calculation, we used convergence criteria of energy and force of 10−5 eV and 0.02 eV/Å. Moreover, 5 images were used along the minimum energy pathway, and the z-coordinates of the membrane were kept frozen to maintain a planar surface.
The molecular dynamics (MD) simulations were performed within the GROMACS package using the universal force field (UFF) [41,42]. The three partial charge models for CO2 and N2 molecules were used to account for the quadrupole distinction of the CO2 and N2 molecules (Table S3) [43,44]. The atomic charges of CH4 and CO were adopted from previous work [17], and are shown in Table S4. The atomic charges of Fe/Ni-PTC membranes were obtained from the ESP fitting based on Merz–Kollman scheme by Gaussian 09 at the B3LYP/6-31g(d) level [45]. The detailed charge information for the Fe/Ni-PTC is shown in Figure S1 in the ESI. In the molecular dynamics simulation system, seven prototype gas molecules, including He, Ne, CO2, N2, Ar, CO, and CH4 (50 for each molecule), were randomly placed into a simulation box sandwiched by a graphene sheet and a Fe/Ni-PTC membrane. The graphene sheet acts as a piston in the system, which was applied with a constant force of 115.1 kJ mol−1 nm−1 (equivalent to a constant pressure of 100 Bar), while the coordinates of the Fe/Ni-PTC membrane were kept frozen. The volume of the simulation box is 4.068 nm × 4.698 nm × 90 nm, and 3D periodic boundary conditions were applied.
The potential of mean force (PMF) [46] for He and N2 gases are calculated by means of the umbrella sampling method. Sixteen sampling windows with 0.1 nm intervals were taken on one side of the MOF membrane, and the total sampling length was 1.5 nm. For each sampling point, the position of the selected gas molecule was restrained using a harmonic potential along the reaction coordinate (z direction). The system was firstly equilibrated for 1 ns, and then the PMF was generated from 20 ns simulated force data in the canonical ensemble using the weighted histogram analysis method [47].

3. Results and Discussion

3.1. The Structural Stability of Fe/Ni-PTC

The metal−PTC 2D MOF has several members depending on the specific types of transition metal in the nodes. In this study, we took Fe-PTC and Ni-PTC as representative models to study their filtration performance for He separation. The optimized structures at DFT level are shown in Figure 1, where the lattice constants for Fe-PTC and Ni-PTC are determined to be a = b = 1.356 nm with the angle between them α = 120°, in good agreement with previous studies [33]. The diameters of the pore are 5.65 Å and 5.67 Å for Fe-PTC and Ni-PTC, respectively. For the Fe-PTC, the lengths of C-C bond, C-S bond, and Fe-S bond are 1.450 Å, 1.729 Å, and 2.094 Å, while for the Ni-PTC, that of the C-C, C-S, and Ni-S bonds are 1.449 Å, 1.710 Å, and 2.128 Å, which are all consistent with previously reported values [32,33]. The cohesive energy is used to examine the structural stability of the Fe/Ni-PTC monolayer, which represents the energy required to break the monolayer down into individual atoms, as defined by [48]:
E c o h = ( n 1 E c + n 2 E s + n 3 E F e / N i − E m e m ) / ( n 1 + n 2 + n 3 )
where Ec, Es, EFe/Ni, and Emem denote the energies of a single C atom, a single S atom, a single Fe/Ni atom, and the total energy of the Fe/Ni-PTC monolayer, respectively; n1, n2, and n3 denote the total numbers of the C atoms, the S atoms, and the Fe/Ni atoms in the unit cell. The cohesive energy of the monolayer Fe-PTC was computed to be 7.10 eV/atom while that of Ni-PTC was 7.21 eV/atom, slightly higher than that of Fe-PTC, indicating its high energetic stability and the possibility of experimental synthesis. Moreover, a previous study has shown that both Fe-PTC and Ni-PTC are all thermally stable [33], which is crucial for gas separation.

3.2. DFT Energy Barriers

In this study, several prototype gas molecules (He, Ne, CO2, CO, N2, Ar, and CH4) from the by-product during natural gas exploitation [49] were selected to assess the He filtration performance of Fe-PTC and Ni-PTC membranes. Prior to exploring the energy barrier of an individual gas molecule penetrating the Fe/Ni-PTC membrane, the adsorption ability of Fe/Ni-PTC monolayers to each gas molecule has been investigated. The adsorption energy between the Fe/Ni-PTC membrane and a gas molecule can be expressed by the following formula:
E a d = E g a s + m e m − E g a s − E m e m
where Egas represents the energy of an individual gas molecule, Emem represents the energy of a pristine membrane, both of which were obtained from the single point energy calculations computed from the geometry-optimized structures of the isolated components, and Egas+mem represents the total energies of a gas molecule adsorbed on the membrane at the most stable adsorption state. The adsorption energy and adsorption height of the most stable state for each molecule on the Fe-PTC and Ni-PTC membranes are summarized in Table 1 and Table 2. It can be seen that the adsorption energy ranges from −2 to −17 kJ mol−1 and the adsorption height is in the range from 2.8 to 3.0 Å, indicating that the gas molecules are physically adsorbed on the Fe/Ni-PTC membrane. Among these molecules, the He molecule has the weakest interaction with the Fe/Ni-PTC (−2.60/−2.31 kJ mol−1), which may facilitate the diffusion of He through the Fe/Ni-PTC membrane. The illustration of the adsorption state of each gas molecule on the Fe/Ni-PTC membrane can be found in Figure S2, where guest molecules uniformly prefer to adsorb on the pore area. Thus, the transition metal cations (Fe or Ni) do not play an important role on stabilization of the guest molecules, because no direct interactions between them were formed.
With the most stable adsorption state, the translocating pathway of an individual gas molecule through the Fe/Ni-PTC membrane was searched by the cNEB approach [40]. The energy barrier (Eb) of a gas molecule passing through the membrane can be defined as [50]:
E b = E T S − E S S
where ETS denotes the energy of a gas molecule and the Fe/Ni-PTC membrane at the transition state (namely, the energy saddle point) and ESS represents that at the most stable adsorption state. The energy profiles of all seven gas molecules passing through the Fe/Ni-PTC membrane are plotted in Figure 1c,d respectively, and the energy barriers are summarized in Table 1 and Table 2. It is obvious that the energy barriers for seven gases are increased following the order of He < Ne < CO2 < N2 < Ar < CO < CH4, where He has the lowest energy barriers to pass through the Fe-PTC (14.64 kJ mol−1) and Ni-PTC (13.00 kJ mol−1) membranes, respectively. The electron density difference maps for the transition state for both gas-Fe-PTC and gas-Ni-PTC systems in Figure S3 also show that the coupling between He and the membranes are the weakest, while that between CH4 and the membranes are the strongest, because the most significant electron transfer happened between CH4 and the membranes. Moreover, for electron density difference maps at the states near the transition states in Figure S4, dipole-shaped interactions can be clearly observed for CO2 and N2, which are induced by the inherent electric field of the membranes. These results reveal that He is more inclined to be separated from other gas molecules by the Fe/Ni-PTC membrane.

3.3. Electron Density Isosurfaces

To better understand the different barriers of Fe/Ni-PTC to seven gas molecules, we then calculated the electron density of gas molecules interacting with the Fe-PTC monolayer, because electron density overlap between them results in repulsive interaction between the gas molecules and the membrane [13]. As shown in Figure 2, with isovalue of 0.01 e Å−3, there is no obvious electron density overlap between the He molecule and the Fe-PTC membrane, which is well consistent with the lowest energy barrier for He diffusing through the Fe-PTC membrane. For other molecules, clear electron density overlap is observed, resulting in higher energy barriers. The electron overlap between CH4 and Fe-PTC membrane is the most significant, thus the highest energy barrier was obtained for CH4 translocating the Fe-PTC membrane. The electron density isosurfaces of the gas molecule interacting with the Ni-PTC monolayer are depicted in Figure S5, where a similar trend has been observed.

3.4. DFT Selectivity of He

Based on the energy barriers of gas molecule translocating the Fe/Ni-PTC membrane, the selectivity of the He molecule to other gases can be evaluated via the Arrhenius equation [51]:
S H e / g a s = r H e r g a s = A H e e − E H e / R T A g a s e − E g a s / R T
where r, A, and E represent the diffusion rate, the diffusion prefactor, and the diffusion energy barrier summarized in Table 1. The diffusion prefactor of all gases in this work was assumed to be 1011 s−1, following previous studies [52,53]. The diffusion rates of all the gas molecules and the selectivity of He over other gas molecules as a function of temperature are shown in Figure 3. With the increase of temperature, the diffusion rate of gas molecules uniformly rises, which can be attributed to the increased kinetic energy of gas molecules in the high-temperature region. In addition, the selectivity of He over other molecules of the Fe/Ni-PTC membrane decreases with the increase of temperature.

3.5. MD Simulations on He Filtration

For industrial He separation, the operation is normally carried out at room temperature [54]. Therefore, MD simulations were conducted to verify the separation performance of Fe/Ni-PTC membranes at 300 K. In the simulation model shown in Figure 4a, 50 He, 50 Ne, 50 Ar, 50 CO, 50 CO2, 50 N2, and 50 CH4 were mixed and randomly distributed between a graphene piston and a Fe/Ni-PTC membrane. Three parallel trajectories of 90 ns have been generated and the results are plotted in Figure 4b,c. It can be seen that during the 90 ns simulation, consecutive events of He passing through the Fe/Ni-PTC membrane were monitored and after 90 ns, nearly all of the He molecules had penetrated through the Fe/Ni-PTC membrane. In contrast, other molecules remained between the two membranes as indicated by the green lines. The insets show the He flow rates perpendicular to the membranes, which decrease with time as increasingly more He molecules have escaped away. Therefore, the results from the MD simulation demonstrated that Fe/Ni-PTC is highly selective to He over other gas molecules. The representative trajectory snapshots are shown in Figures S6 and S7 in the ESI.
In addition to selectivity, He permeance is another indicative factor to judge the separation efficiency of a membrane. Thus, we further calculated the He permeance of Fe/Ni-PTC membrane at temperatures ranging from 200–500 K, following the method reported in a previous work [55] which can be found in the ESI. The results of He permeance are listed in Table 3. It can be seen that as the temperature increases, the permeance gradually increased from 10−4 to 10−3 mol s−1 m−2 Pa−1, which is attributed to the increased kinetic energy of He particles at high temperatures. At room temperature of 300 K, the permeance of He reached ~ 6 × 10−4 mol s−1 m−2 Pa−1 for both Fe-PTC and Ni-PTC membrane, which is five orders of magnitude higher than the industrially acceptable value (6.79 × 10−9 mol s−1 m−2 Pa−1) [56], implying the plausibility of the Fe/Ni-PTC membrane for the industrial applications.

3.6. Comparison between Membranes

After the pre-treatment of natural gas, the major component in the remaining gas stream is N2. Thus, the selectivity of the Fe/Ni-PTC for He over N2 is of great significance in practical applications. We compared the He permeance and the He/N2 selectivity of both Fe-PTC and Ni-PTC membrane with previously proposed nanoporous membranes to evaluate the He separation efficiency of the two MOFs. As shown in Table 4, it is obvious that although the Fe/Ni-PTC has relatively lower He permeance, the He/N2 selectivity of Fe/Ni-PTC MOF is significantly higher than those of g-C2O, C2N and 6N-PG membranes. While compared to the CTF-0 membrane, the selectivity of Fe/Ni-PTC MOF is lower, but the permeance is three orders of magnitude higher. It should be noted that compared to the IGP membrane, with almost the same permeance, the He selectivity of Fe/Ni-PTC is higher by five orders of magnitude, indicating its superior performance for He harvest.

3.7. Free Energy Barrier Analysis

The above results have demonstrated that the Fe/Ni-PTC membrane have excellent selectivity and permeance for He separation, both in the high selectivity and permeance. To get a deeper insight into the improved He/N2 selectivity and He permeance at room temperature, we have estimated the free energy barriers by calculating the PMF of He and N2 passing through the Fe/Ni-PTC membrane and the results are summarized in Figure 5. As can be seen, the free energy barrier values for He penetrating Fe/Ni-PTC membrane are 5.62 and 5.18 kBT, respectively, which are much lower than those of N2 (73.8 and 66.2 kBT), resulting in both the high He/N2 selectivity and high He permeance. It is worth noting that compared with the energy barriers based on DFT calculation for N2 (117.61/100.95 kJ mol−1, equivalent to 47.2/40.5 kBT for Fe-PTC and Ni-PTC respectively), the free energy barriers obtained from PMF analysis are increased. This is because the temperature-related conformational entropy for rod-like molecules such as N2 is explicitly considered.

4. Conclusions

To conclude, we performed combined first-principles calculations and MD simulations to theoretically investigate the feasibility of the 2D Fe/Ni-PTC MOF materials for He separation from natural gas and other noble gases. High permeance (~6 × 10−4 mol s−1 m−2 Pa−1) and He/N2 selectivity (1015~1017) of Fe/Ni-PTC MOF are predicted, superior to most of the previously proposed 2D porous membranes. Moreover, quantitative free energy analysis indicates that the free energy barriers for He passing through the Fe/Ni-PTC MOF are only ~5 kBT, which is significantly lower than that of N2 (~70 kBT), resulting in the high selectivity for He separation at room temperature. Overall, our results have demonstrated that the newly synthesized 2D Fe-PTC and its analogue Ni-PTC are promising membranes for He separation, which might have great prospects in industrial production.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/membranes11120927/s1, Figure S1: Atomic charges, Figure S2: The schematic diagram of the adsorption states, Figures S3 and S4: The electron density difference maps, Figure S5: Electron-density isosurfaces, Figures S6 and S7: Snapshots, Figure S8: Model for calculating permeance, Figure S9: Number of permeated He molecules versus simulation time. Table S1: The convergence of the energy, Table S2: The vector changes and the energy differences, Table S3: Force field parameters, Table S4: Atomic charges. Supplementary Method: Method for calculating permeance.

Author Contributions

Conceptualization, Y.Q.; methodology, J.W.; software, Y.L. (Yixiang Li); validation, Y.L. (Yixiang Li) and J.W.; formal analysis, J.W.; investigation, J.W.; data curation, J.W.; writing—original draft preparation, J.W., Y.Q. and W.L.; writing—review & editing, Y.Q. and J.G.; supervision, Y.Q., J.G. and W.L.; project administration, Y.Y., Y.L. (Yongqiang Li) and J.G.; funding acquisition, Y.Q., W.L. and M.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number Nos.12074215, 11874238, the Natural Science Foundation of Shandong Province of China (No. ZR2020JQ04) and the Taishan Scholar Project of Shandong Province.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

We would like to acknowledge all support for this work by the Shandong University, the National Natural Science Foundation of China and the Taishan Scholar Project of Shandong Province.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaplan, K.H. Helium shortage hampers research and industry. Phys. Today 2007, 60, 31. [Google Scholar] [CrossRef]
  2. Cho, A. Helium-3 shortage could put freeze on low-temperature research. Science 2009, 326, 778–779. [Google Scholar] [CrossRef]
  3. Nuttall, W.J.; Clarke, R.H.; Glowacki, B.A. Resources: Stop squandering helium. Nature 2012, 485, 573–575. [Google Scholar] [CrossRef] [PubMed]
  4. Smith, Z.P.; Tiwari, R.R.; Dose, M.E.; Gleason, K.L.; Murphy, T.M.; Sanders, D.F.; Gunawan, G.; Robeson, L.M.; Paul, D.R.; Freeman, B.D. Influence of Diffusivity and Sorption on Helium and Hydrogen Separations in Hydrocarbon, Silicon, and Fluorocarbon-Based Polymers. Macromolecules 2014, 47, 3170–3184. [Google Scholar] [CrossRef]
  5. Yuan, W.; Chen, J.; Shi, G. Nanoporous graphene materials. Mater. Today 2014, 17, 77–85. [Google Scholar] [CrossRef]
  6. Hosseini, S.S.; Teoh, M.M.; Chung, T.S. Hydrogen separation and purification in membranes of miscible polymer blends with interpenetration networks. Polymer 2008, 49, 1594–1603. [Google Scholar] [CrossRef]
  7. Hu, W.; Wu, X.; Li, Z.; Yang, J. Porous silicene as a hydrogen purification membrane. Phys. Chem. Chem. Phys 2013, 15, 5753–5757. [Google Scholar] [CrossRef]
  8. Koenig, S.P.; Wang, L.; Pellegrino, J.; Bunch, J.S. Selective molecular sieving through porous graphene. Nat. Nanotechnol. 2012, 7, 728–732. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Bunch, J.S.; Verbridge, S.S.; Alden, J.S.; van der Zande, A.M.; Parpia, J.M.; Craighead, H.G.; McEuen, P.L. Impermeable atomic membranes from graphene sheets. Nano. Lett. 2008, 8, 2458–2462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Fischbein, M.D.; Drndic′, M. Electron beam nanosculpting of suspended graphene sheets. Appl. Phys. Lett. 2008, 93, 113107. [Google Scholar] [CrossRef] [Green Version]
  11. Hu, W.; Wu, X.; Li, Z.; Yang, J. Helium separation via porous silicene based ultimate membrane. Nanoscale 2013, 5, 90629066. [Google Scholar] [CrossRef]
  12. Kedi, Y.; Shengxi, H.; Xiaofei, C.; Xinwei, W.; Jing, K.; Yan, C.; Jianming, X. Generating Sub-nanometer Pores in Single-Layer MoS2 by Heavy-Ion Bombardment for Gas Separation: A Theoretical Perspective. ACS Appl. Mater. Interfaces 2018, 10, 28909–28917. [Google Scholar]
  13. Li, F.; Qu, Y.; Zhao, M. Efficient helium separation of graphitic carbon nitride membrane. Carbon 2015, 95, 51–57. [Google Scholar] [CrossRef]
  14. Wang, Y.; Li, J.; Yang, Q.; Zhong, C. Two-Dimensional Covalent Triazine Framework Membrane for Helium Separation and Hydrogen Purification. ACS Appl. Mater. Interfaces 2016, 8, 8694–8701. [Google Scholar] [CrossRef]
  15. Liu, X.; Chang, X.; Zhu, L.; Li, X. High-efficiency helium separation through g-C2O membrane: A theoretical study. Comput. Mater. Sci. 2019, 157, 1–5. [Google Scholar] [CrossRef]
  16. Zhu, L.; Xue, Q.Z.; Li, X.F.; Wu, T.T.; Jin, Y.K.; Xing, W. C2N: An excellent two-dimensional monolayer membrane for He separation. J. Mater. Chem. 2015, 3, 21351–21356. [Google Scholar] [CrossRef]
  17. Wang, L.; Li, F.; Wang, J.; Li, Y.; Li, W.; Yang, Y.; Zhao, M.; Qu, Y. High-efficiency helium separation through an inorganic graphenylene membrane: A theoretical study. Phys. Chem. Chem. Phys. 2020, 22, 9789–9795. [Google Scholar] [CrossRef]
  18. Rezaee, P.; Naeij, H.R. Graphenylenee1 membrane: An excellent candidate for hydrogen purification and helium separation. Carbon 2020, 157, 779–787. [Google Scholar] [CrossRef] [Green Version]
  19. Meek, S.T.; Greathouse, J.A.; Allendorf, M.D. Metal−organic frameworks: A rapidly growing class of versatile nanoporous materials. Adv. Mater. 2001, 23, 249. [Google Scholar] [CrossRef] [PubMed]
  20. Eddaoudi, M.; Li, H.; Yaghi, O. Highly porous and stable metal-organic frameworks: Structure design and sorption properties. J. Am. Chem. Soc. 2000, 122, 1391–1397. [Google Scholar] [CrossRef]
  21. Stock, N.; Biswas, S. Synthesis of metal-organic frameworks (MOFs): Routes to various MOF topologies, morphologies, and composites. Chem. Rev. 2011, 112, 933–969. [Google Scholar] [CrossRef] [PubMed]
  22. Li, J.-R.; Kuppler, R.J.; Zhou, H.-C. Selective Gas Adsorption and Separation in Metal-Organic Frameworks. Chem. Soc. Rev. 2009, 38, 1477–1504. [Google Scholar] [CrossRef] [PubMed]
  23. Li, J.; Sculley, J.; Zhou, H. Metal-organic frameworks for separations. Chem. Rev. 2012, 112, 869–932. [Google Scholar] [CrossRef]
  24. Adil, K.; Belmabkhout, Y.; Pillai, R.S.; Cadiau, A.; Bhatt, P.M.; Assen, A.H.; Maurin, G.; Eddaoudi, M. Gas/Vapour Separation Using Ultra-Microporous Metal−Organic Frameworks: Insights into the Structure/Separation Relationship. Chem. Soc. Rev. 2017, 46, 3402–3430. [Google Scholar] [CrossRef]
  25. Basdogan, Y.; Keskin, S. Simulation and modelling of MOFs for hydrogen storage. CrystEngComm 2015, 17, 261–275. [Google Scholar] [CrossRef] [Green Version]
  26. Murray, L.J.; Dinca, M.; Long, J.R. Hydrogen storage in metal-organic frameworks. Chem. Soc. Rev. 2009, 38, 1294–1314. [Google Scholar] [CrossRef]
  27. Horcajada, P.; Serre, C.; Maurin, G.; Ramsahye, N.A.; Balas, F.; Vallet-Regi, M.; Sebban, M.; Taulelle, F.; Ferey, G. Flexible porous metal-organic frameworks for a controlled drug delivery. J. Am. Chem. Soc. 2008, 130, 6774–6780. [Google Scholar] [CrossRef] [PubMed]
  28. Zhu, L.; Liu, X.-Q.; Jiang, H.-L.; Sun, L.-B. Metal-Organic Frameworks for Heterogeneous Basic Catalysis. Chem. Rev. 2017, 117, 8129–8176. [Google Scholar] [CrossRef]
  29. Ranjan, R.; Tsapatsis, M. Microporous metal organic framework membrane on porous support using the seeded growth method. Chem. Mater. 2009, 21, 4920–4924. [Google Scholar] [CrossRef]
  30. Cao, F.; Zhang, C.; Xiao, Y.; Huang, H.; Zhang, W.; Liu, D.; Zhong, C.; Yang, Q.; Yang, Z.; Lu, X. Helium recovery by a Cu-BTC metal–organic-framework membrane. Ind. Eng. Chem. Res. 2012, 51, 11274–11278. [Google Scholar] [CrossRef]
  31. Kasik, A.; Dong, X.; Lin, Y. Synthesis and stability of zeolitic imidazolate framework 68 membranes. Micropor. Mesopor. Mater. 2015, 204, 99–105. [Google Scholar] [CrossRef]
  32. Dong, R.; Zhang, Z.; Tranca, D.C.; Zhou, S.; Wang, M.; Adler, P.; Liao, Z.; Liu, F.; Sun, Y.; Shi, W.; et al. A coronene-based semiconducting two-dimensional metal-organic framework with ferromagnetic behavior. Nat. Commun. 2018, 9, 2637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Wang, J.; Wang, J.; Song, X.; Qi, S.; Zhao, M. Multifunctional electrocatalytic activity of coronene-based two-dimensional metal-organic frameworks: TM-PTC. Appl. Surf. Sci. 2020, 511, 145393. [Google Scholar] [CrossRef]
  34. Kresse, G.; Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  35. Kresse, G.; Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  36. Blochl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [Green Version]
  37. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  38. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [Green Version]
  39. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  40. Kolsbjerg, E.L.; Groves, M.N.; Hammer, B. An automated nudged elastic band method. J. Chem. Phys 2016, 145, 094107. [Google Scholar] [CrossRef]
  41. Pronk, S.; Pall, S.; Schulz, R.; Larsson, P.; Bjelkmar, P.; Apostolov, R.; Shirts, M.R.; Smith, J.C.; Kasson, P.M.; van der Spoel, D.; et al. GROMACS 4.5: A high-throughput and highly parallel opensource molecular simulation toolkit. Bioinformatics 2013, 29, 845–854. [Google Scholar] [CrossRef]
  42. Garberoglio, G.; Comput, J. OBGMX: A web-based generator of GROMACS topologies for molecular and periodic systems using the universal force field. J. Comput. Chem. 2012, 33, 2204–2208. [Google Scholar] [CrossRef]
  43. Harris, J.G.; Yung, K.H. Carbon Dioxides Liquid-Vapor Coexistence Curve And Critical Properties As Predicted By A Simple Molecular-Model. J. Phys. Chem. 1995, 99, 12021–12024. [Google Scholar] [CrossRef]
  44. Qu, Y.; Li, F.; Zhao, M. Theoretical Design of Highly Efficient CO2/N2 Separation Membranes Based on Electric Quadrupole Distinction. J. Phys. Chem. C 2017, 121, 17925–17931. [Google Scholar] [CrossRef]
  45. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, Revision A02. Gaussian. Inc. Wallingford. CT 2009, 7. [Google Scholar] [CrossRef] [Green Version]
  46. Torrie, G.M.; Valleau, J.P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187–199. [Google Scholar] [CrossRef]
  47. Kumar, S.; Rosenberg, J.M.; Bouzida, D.; Swendsen, R.H.; Kollman, P.A. The weighted histogram analysis method for free-energy calculations on biomolecules.1. The method. J. Comput. Chem. 1992, 13, 1011–1021. [Google Scholar] [CrossRef]
  48. Li, X.; Guo, T.; Zhu, L.; Ling, C.; Xue, Q.; Xing, W. Charge-modulated CO2 capture of C3N nanosheet: Insights from DFT calculations. Chem. Eng. J. 2018, 338, 92–98. [Google Scholar] [CrossRef]
  49. Adewole, J.K.; Ahmad, A.L.; Ismail, S.; Leo, C.P. Current challenges in membrane separation of CO2 from natural gas: A review. Int. J. Greenh. Gas. Control 2013, 17, 46–65. [Google Scholar] [CrossRef]
  50. Chang, X.; Xue, Q.; He, D.; Zhu, L.; Li, X.; Tao, B. 585 divacancy-defective germanene as a hydrogen separation membrane: A DFT study. Int. J. Hydrog. Energy 2017, 42, 24189–24196. [Google Scholar] [CrossRef]
  51. Arrhenius, S. Über die Reaktionsgeschwindigkeit bei der inversion von Rohrzucker durch Sauren. Phys. Chem. 1889, 4, 226–248. [Google Scholar]
  52. Blankenburg, S.; Bieri, M.; Fasel, R.; Müllen, K.; Pignedoli, C.A.; Passerone, D. Porous graphene as an atmospheric nanofilter. Small 2010, 6, 2266–2271. [Google Scholar] [CrossRef]
  53. Zhang, Y.; Shi, Q.; Liu, Y.; Wang, Y.; Meng, Z.; Xiao, C.; Deng, K.; Rao, D.; Lu, R. Hexagonal Boron Nitride with Designed Nanopores as a High-Efficiency Membrane for Separating Gaseous Hydrogen from Methane. J. Phys. Chem. C 2015, 119, 19826–19831. [Google Scholar] [CrossRef]
  54. Sullivan, D.M.; Bruening, M.L. Ultrathin, Gas-Selective Polyimide Membranes Prepared from Multilayer Polyelectrolyte Films. Chem. Mater. 2003, 15, 281–287. [Google Scholar] [CrossRef]
  55. Xu, J.; Zhou, S.; Sang, P.; Li, J.; Zhao, L. Inorganic graphenylene as a promising novel boron nitrogen membrane for hydrogen purification: A computational study. J. Mater. Sci. 2017, 52, 10285–10293. [Google Scholar] [CrossRef]
  56. Zhu, Z. Permeance should be used to characterize the productivity of a polymeric gas separation membrane. J. Membr. Sci. 2006, 281, 754–756. [Google Scholar] [CrossRef]
  57. Wei, S.; Zhou, S.; Wu, Z.; Wang, M.; Wang, Z.; Guo, W.; Lu, X. Mechanistic insights into porous graphene membranes for helium separation and hydrogen purification. Appl. Surf. Sci. 2018, 441, 631–638. [Google Scholar] [CrossRef]
Membranes 11 00927 g001 550
Figure 1. The top views of the optimized unit cell of the (a) Fe-PTC and (b) Ni-PTC monolayer, where a → and b → represent the cell vectors and α represents the angle between them. The energy profiles for seven prototype gas molecules translocating the (c) Fe-PTC and (d) Ni-PTC membranes.
Figure 1. The top views of the optimized unit cell of the (a) Fe-PTC and (b) Ni-PTC monolayer, where a → and b → represent the cell vectors and α represents the angle between them. The energy profiles for seven prototype gas molecules translocating the (c) Fe-PTC and (d) Ni-PTC membranes.
Membranes 11 00927 g001
Membranes 11 00927 g002 550
Figure 2. Electron density isosurfaces for the gas molecules: (a) He, (b) Ne, (c) CO2, (d) N2, (e) Ar and (f) CH4 translocating the Fe-PTC membrane (isovalue of 0.01 Å−3).
Figure 2. Electron density isosurfaces for the gas molecules: (a) He, (b) Ne, (c) CO2, (d) N2, (e) Ar and (f) CH4 translocating the Fe-PTC membrane (isovalue of 0.01 Å−3).
Membranes 11 00927 g002
Membranes 11 00927 g003 550
Figure 3. (a) The diffusion rates and (b) the selectivity for gases to pass through the Fe-PTC membrane as a function of temperature; (c) The diffusion rates and (d) the selectivity for gases to pass through the Ni-PTC membrane as a function of temperature. Different colors denote different gas molecules.
Figure 3. (a) The diffusion rates and (b) the selectivity for gases to pass through the Fe-PTC membrane as a function of temperature; (c) The diffusion rates and (d) the selectivity for gases to pass through the Ni-PTC membrane as a function of temperature. Different colors denote different gas molecules.
Membranes 11 00927 g003
Membranes 11 00927 g004 550
Figure 4. (a) The snapshot of the simulation model; pure He permeation by (b) Fe-PTC and (c) Ni-PTC membranes at 300 K, where individual results obtained from three independent trajectories are represented by red, blue, and black dots. The dashed lines are the numerical fitted results. The insets show the gas flow rates calculated from the numerical fitted results.
Figure 4. (a) The snapshot of the simulation model; pure He permeation by (b) Fe-PTC and (c) Ni-PTC membranes at 300 K, where individual results obtained from three independent trajectories are represented by red, blue, and black dots. The dashed lines are the numerical fitted results. The insets show the gas flow rates calculated from the numerical fitted results.
Membranes 11 00927 g004
Membranes 11 00927 g005 550
Figure 5. The PMF of He and N2 passes through the (a) Fe-PTC and (b) Ni-PTC membranes by umbrella sampling.
Figure 5. The PMF of He and N2 passes through the (a) Fe-PTC and (b) Ni-PTC membranes by umbrella sampling.
Membranes 11 00927 g005
Table
Table 1. The adsorption height (H0), the corresponding adsorption energy (Ead) between the gas molecule and the Fe-PTC monolayer, and the energy barriers (Eb) of the seven gas molecules translocating the Fe-PTC monolayer.
Table 1. The adsorption height (H0), the corresponding adsorption energy (Ead) between the gas molecule and the Fe-PTC monolayer, and the energy barriers (Eb) of the seven gas molecules translocating the Fe-PTC monolayer.
Gas MoleculeH0 (Å)Ead (kJ mol−1)Eb (kJ mol−1)
He3.000−2.6014.64
Ne2.901−4.4340.26
CO23.337−9.2560.87
N22.836−13.20117.61
Ar2.987−10.98141.78
CO3.000−13.97152.96
CH42.950−14.74185.42
Table
Table 2. The adsorption height (H0), the corresponding adsorption energy (Ead) between the gas molecule and the Ni-PTC monolayer, and the energy barriers (Eb) of the seven gas molecules translocating the Ni-PTC monolayer.
Table 2. The adsorption height (H0), the corresponding adsorption energy (Ead) between the gas molecule and the Ni-PTC monolayer, and the energy barriers (Eb) of the seven gas molecules translocating the Ni-PTC monolayer.
Gas MoleculeH0 (Å)Ead (kJ mol−1)Eb (kJ mol−1)
He2.997−2.3113.00
Ne3.000−4.4334.19
CO22.929−6.8462.51
N22.918−13.00100.95
Ar3.000−10.79125.61
CO2.999−13.87143.91
CH42.683−16.28171.36
Table
Table 3. He permeance (mol s−1 m−2 Pa−1) at the temperature range of 200–500 K.
Table 3. He permeance (mol s−1 m−2 Pa−1) at the temperature range of 200–500 K.
MembraneTemperature (K)Permeance
Fe-PTC2002.9 × 10−4
3006.3 × 10−4
4008.9 × 10−4
5001.3 × 10−3
Ni-PTC2003.6 × 10−4
3006.6 × 10−4
4009.2 × 10−4
5001.5 × 10−3
Table
Table 4. Comparison results of He permeance (GPU) and the He/N2 selectivity between Fe/Ni-PTC MOF and previously proposed porous membranes at 300 K. (1 GPU = 3.3 × 10−10 mol s−1 m−2 Pa−1).
Table 4. Comparison results of He permeance (GPU) and the He/N2 selectivity between Fe/Ni-PTC MOF and previously proposed porous membranes at 300 K. (1 GPU = 3.3 × 10−10 mol s−1 m−2 Pa−1).
MembranePermeanceSelectivity
Fe-PTC 11.9 × 1069.1 × 1017
Ni-PTC 12.1 × 1062.2 × 1015
g-C2O 21.0 × 1071.5 × 106
IGP 32.0 × 1061.0 × 1012
CTF-0 41.0 × 1032.0 × 1027
C2N 51.0 × 1073.0 × 1012
6N-PG 66.9 × 1076.0 × 108
1 This study. 2 Ref. [15]. 3 Ref. [17]. 4 Ref. [14]. 5 Ref. [16]. 6 Ref. [57].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, J.; Li, Y.; Yang, Y.; Li, Y.; Zhao, M.; Li, W.; Guan, J.; Qu, Y. Efficient Helium Separation with Two-Dimensional Metal–Organic Framework Fe/Ni-PTC: A Theoretical Study. Membranes 2021, 11, 927. https://doi.org/10.3390/membranes11120927

AMA Style

Wang J, Li Y, Yang Y, Li Y, Zhao M, Li W, Guan J, Qu Y. Efficient Helium Separation with Two-Dimensional Metal–Organic Framework Fe/Ni-PTC: A Theoretical Study. Membranes. 2021; 11(12):927. https://doi.org/10.3390/membranes11120927

Chicago/Turabian Style

Wang, Jingyuan, Yixiang Li, Yanmei Yang, Yongqiang Li, Mingwen Zhao, Weifeng Li, Jing Guan, and Yuanyuan Qu. 2021. "Efficient Helium Separation with Two-Dimensional Metal–Organic Framework Fe/Ni-PTC: A Theoretical Study" Membranes 11, no. 12: 927. https://doi.org/10.3390/membranes11120927

APA Style

Wang, J., Li, Y., Yang, Y., Li, Y., Zhao, M., Li, W., Guan, J., & Qu, Y. (2021). Efficient Helium Separation with Two-Dimensional Metal–Organic Framework Fe/Ni-PTC: A Theoretical Study. Membranes, 11(12), 927. https://doi.org/10.3390/membranes11120927

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