Next Article in Journal
Recent Advances in Molecular Dynamics Simulations of Tau Fibrils and Oligomers
Next Article in Special Issue
Stimuli-Responsive Track-Etched Membranes for Separation of Water–Oil Emulsions
Previous Article in Journal
Zero-Material Cost Production of Soil-Coated Fabrics with Underwater Superoleophobicity for Antifouling Oil/Water Separation
Previous Article in Special Issue
Performance of Hypersaline Brine Desalination Using Spiral Wound Membrane: A Parametric Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Efficient CO2 Separation at High Pressure by Grain-Boundary-Controlled CHA Zeolite Membrane Investigated by Non-Equilibrium Molecular Dynamics

1
Graduate School of Engineering, Kogakuin University, Tokyo 192-0015, Japan
2
Department of Environmental Chemistry and Chemical Engineering, School of Advanced Engineering, Kogakuin University, Tokyo 192-0015, Japan
*
Author to whom correspondence should be addressed.
Membranes 2023, 13(3), 278; https://doi.org/10.3390/membranes13030278
Submission received: 10 February 2023 / Accepted: 23 February 2023 / Published: 26 February 2023
(This article belongs to the Special Issue Membrane Materials and Processes for Liquid and Gas Separation)

Abstract

:
The CO2 permeability and selectivity of CHA-type zeolite membranes in the separation of a CO2/CH4 mixture gas at high pressure were evaluated using non-equilibrium molecular dynamics (NEMD). It was found that in a perfectly crystalline, defect-free CHA membrane, the adsorption of CH4, which diffuses slowly in the pores, hinders CO2 permeation. Therefore, an increase in the amount of CH4 adsorbed at high pressure decreases the CO2 permeability and significantly reduces the CO2 selectivity of the CHA membrane. CHA membranes with grain boundaries parallel to the permeation direction were found to show higher CO2 selectivity than perfectly crystalline CHA membranes at high pressure, as the blocking effect of CH4 on CO2 permeation occurring within the grain boundary is not significant. This paper is the first to show that the CO2 permeability of CHA membranes with controlled grain boundaries can exceed the intrinsic performance of fully crystalline zeolite membranes at high pressure.

1. Introduction

Greenhouse gases are considered to be the main cause of global warming and the reduction of CO2 emissions into the atmosphere is necessary; CCS (carbon dioxide capture and storage) and CCUS (carbon dioxide capture, utilization, and storage) are attracting worldwide attention, and CO2 separation is a fundamental technology for them. Natural gas refining is a well-known CO2 separation process [1]. Associated gas produced from gas fields contains CO2, so a process to purify it from methane is essential for the utilization of natural gas. The gas concentration of the associated gas varies from producing area to producing area, and the CO2 concentration is known to vary greatly from 10 to 90% [2].
Chemical absorption is the current process in natural gas plants but is expensive for the regeneration of the absorbent when separating high concentrations of CO2. The membrane separation is a process driven by the chemical potential gradient (concentration and pressure). Membrane separations are simple devices and can be operated continuously so that the pressure of gas emanating from gas fields can be used as the driving force, such as in natural gas plants, which are the target of separation. Polymer membranes have already been commercialized in the field of natural gas, but there are issues with CO2 separation performance and membrane durability [3]. Furthermore, they show a high affinity for aromatic hydrocarbons, which are present in trace amounts in natural gas, making their practical application difficult depending on the compositional conditions of the natural gas produced. Zeolite membranes, on the other hand, are promising as CO2 separation membranes due to their excellent properties such as molecular sieving effect, selective adsorption, and high mechanical strength derived from their regular pore structure. Recently, zeolite membranes with an 8-membered ring pore size (3.8 × 3.8 Å), such as chabazite (CHA) and Deca-dodecasil 3R (DDR), have been successfully produced and have been reported to show high separation performance with separation factors exceeding 100 for CO2/CH4 mixture separation [4,5]. This is due to the fact that the kinetic diameter of the CO2 molecule (3.3 Å) is smaller than the pore size and that of the CH4 molecule (3.8 Å) is equal to the pore size of CHA zeolite.
High-silica CHA-type zeolite membranes exhibit high gas permeability due to their three-dimensional pore structure and large pore volume [6,7,8]. However, it has been reported that the CO2 permeability of CHA-type zeolite membranes decreases in high-pressure CO2/CH4 mixture gas systems [9,10]. Therefore, improved gas permeability is required for the application of zeolite membranes in high-pressure CO2 separation, such as in natural gas purification.
Grain boundaries formed during the formation process of zeolite membranes have been considered a problem to be solved for practical application as a factor that reduces the molecular sieving properties. Zeolite membranes have a polycrystalline structure in which zeolite crystals are densely deposited on a porous support, as mechanical strength and densification become problems without support [11]. In recent years, membrane formation methods that can suppress the formation of grain boundaries and repair methods after formation have been proposed [12,13,14]. Recent studies have attempted to quantitatively assess the influence of grain boundaries in Silicalite-1 (MFI) membranes by combining fluorescence confocal optical microscopy (FCOM) and image processing techniques, and have also shown that cracks (crystal cracks) are the biggest factor in reducing the molecular sieving effect of zeolite membranes [15]. The most important factor that reduces the molecular sieving effect of zeolite films is cracking [16].
However, it is difficult to completely remove grain boundaries from zeolite membranes, and the inherent performance of zeolites cannot be exceeded. The mechanism of gas separation at grain boundaries is also not well understood. We have reported the investigation of gas diffusion at grain boundaries inside MFI using the equilibrium molecular dynamics (EMD) simulation [17,18], but the effect of grain boundaries in mixed gas systems by EMD is difficult to determine due to the effect of the chemical potential gradient not being considered directly. On the other hand, the non-equilibrium molecular dynamics (NEMD) simulation can directly simulate gas permeation in membranes. For the direct simulation of membrane permeation, the dual-control volume grand canonical MD (DCV-GCMD) and simple NEMD have been reported. [19,20] We have used NEMD simulations to study the permeation mechanism of zeolite membranes with grain boundaries. In the process, it was found that zeolite membranes with slit-like controlled grain boundaries exhibit a higher separation performance than zeolite-specific CO2 separation performance [21]. In particular, CHA membranes controlled to have grain boundaries parallel to the direction of gas permeation were found to exhibit high CO2 permeability.
In this study, a perfectly crystalline model and a membrane model with controlled grain boundaries were prepared for a high-silica CHA-type zeolite membrane, and the effect of controlled grain boundaries on the CO2 permeability was discussed by simulating CO2/CH4 mixture gas permeation using the NEMD method at different pressure conditions. The NEMD results show that the presence of CH4 with slow diffusion within the zeolite crystals decreases the CO2 permeability, suggesting that the CO2 aggregation effect at the grain boundaries can be used to improve the CO2 permeability of the CHA membrane.

2. Calculation Method

For the CO2/CH4 permeation simulation in CHA-type zeolite membranes, our previously reported NEMD calculation scheme was used, with atomic interaction parameters optimized to reproduce the experimental adsorption isotherm data [21]. In this study, the parameters of PCFF were applied for van der Waals between all atoms, and the atomic charges of silica, oxygen in crystal, and oxygen at the terminal were 1.244, −0.622, and −0.311, respectively. Figure 1 shows the unit cells of the perfectly crystalline model and the polycrystalline model in the NEMD. The gas molecules were generated on the generation region with a given velocity according to the Boltzmann distribution toward the membrane model, and permeated molecules were removed from the unit cell on the deletion region. This method makes it possible to control the number of gases impinging on the membrane on the feed and permeation sides.
The framework used a CHA-type zeolite structure with SiO2/Al2O3 = ∞, which contains no Al atoms at all, a perfectly crystalline model without defect, and a polycrystalline model consisting of surfaces with different crystal orientations fabricated. High-silica CHA membranes without Al are hydrophobic, which makes them very advantageous zeolite membranes for natural gas processing containing vapor [6]. The polycrystalline model was assumed to be a membrane with grain boundaries perpendicular to the surface. The grain boundaries in polycrystalline models were considered to be a slit pore with a width of 0.6 nm. The membrane thickness of all models was around 0.45 nm. The periodic boundary condition was used in two Cartesian coordinates along the membrane surface in the NEMD. To maintain the crystalline nature of the CHA zeolite, the coordinates of the zeolite’s constituent atoms, Si and O, were fixed during the NEMD. The oxygen atoms exposed to the surface were terminated with a hydrogen atom; on the other hand, the bared O atom on the grain boundary surface was not terminated with a hydrogen atom. Although the surface of the grain boundary of the CHA membrane would be terminated by the -OH group, because the test calculation of NEMD using the all OH terminated model showed no significant difference in the calculated membrane performance, the unterminated boundary surfaces were considered to be defined by the grain boundary size clearly.
The equimolar CO2/CH4 mixture was considered as the feed gas. Permeation simulations by NEMD were carried out by varying the time interval f for gas molecule generation, thereby varying the partial pressure conditions on the feed side (Equation (1)) [22].
f = P S 2 π m k B T  
where S is the membrane surface area, m is the molecular weight of the gas molecule, k B is Boltzmann constant, and T is the temperature. All NEMD simulations were carried out at 298 K and the temperature of the system was controlled by the velocity scaling method. Table 1 summarizes the appearance intervals f used in this study. Gas molecules that were reflected at the membrane surface and far enough away from the membrane were deleted as not involved in permeation. The gas molecules permeating the membrane were removed just before they reached the permeate side and the pressure on the permeate side in all NEMD simulations was zero. The permeate flux J is calculated by counting the number of molecules, n, permeated during time Δt at steady state and dividing by the surface area, S, of the membrane (Equation (2)).
J = n Δ t · S

3. Results and Discussion

3.1. Perfectly Crystalline Model

NEMD simulations of the CO2/CH4 mixture gas at different pressure conditions were performed on a perfectly crystalline membrane model of CHA zeolite with (1 0 0) crystalline planes oriented on the surface. Figure 2 shows snapshots of the NEMD results on the perfectly crystalline model at the feed pressures of 0.5 MPa and 8.0 MPa. Note that gas molecules are not shown in the snapshot as they are removed when they permeate the bottom of the zeolite membrane. Under the pressure condition of 0.5 MPa, CO2 molecules were observed to gradually permeate from the membrane surface, occupying the pores, irrespective of the orientation of the crystal planes. On the other hand, the number of CH4 molecules in the membrane was few. This is because of the molecular sieving effect of the membrane. The direction of the CO2 molecules in the membrane was parallel to the permeation direction and they were observed to be arranged in a row. This mode of diffusion is known as “single-file diffuse” [23]. This sequence is considered to indicate the direction of permeation of the CO2 molecules; the direction of the 8-membered ring pores in the CHA zeolite membrane. As the atomic interactions between the gas molecules and zeolite used in this paper were based on parameters that could reproduce experimentally reported adsorption isotherms [24], it can be assumed that CO2 molecules permeate according to the single-file diffusion within the zeolite crystals resulting in a high CO2 selectivity.
As shown in Figure 2b, the NEMD at the feed side pressure of 8.0 MPa shows that the number of CO2 molecules in the membrane is lower than at 0.5 MPa. Figure 3 shows the change over time in the number of permeated gas molecules at 0.5 MPa and 8.0 MPa. The number of permeated CH4 molecules at 0.5 MPa was very small and it takes a longer calculation time (>100 ns) to calculate a precise flux. This may be due to the fact that in the defect-free, perfectly crystalline model, the kinetic molecular diameter of CH4 is almost equal to the pore diameter of the CHA zeolite, which requires a long calculation time for diffusion through the zeolite [25]. The number of CO2 molecules permeated increased almost monotonically after 20 ns, irrespective of the pressure conditions. At 8.0 MPa, the number of CH4 molecules permeated also increased monotonically after 20 ns. This is probably because of an increase in the frequency of CH4 impinging on the membrane surface at high pressure and an increase in the number of CH4 molecules entering the membrane pore. In other words, in a defect-free, perfectly crystalline CHA membrane, the permeation of CH4 molecules is significantly reduced at low pressure, indicating that although CO2, which is smaller than the pore size, is selectively permeated, this advantage may be lost at higher pressures.
Figure 4 shows the permeability and separation factor from the NEMD with the feed side pressures of 0.5, 2.0, 4.0, 6.0, and 8.0 MPa: the CO2 permeability decreases at higher pressures, while the CH4 permeability increases slightly. The results indicate that the CO2 selectivity decreases in the perfectly crystalline CHA membrane at high pressure. This is consistent with the trend in the experimental data of Kida et al. [9]. Analysis of the concentration distribution of gas molecules in the perfectly crystalline CHA membrane model for different pressure conditions on the feed side (Figure 5) shows that the CO2 concentration near the membrane surface on the feed side decreases with increasing pressure. On the other hand, the CH4 concentration increased at higher pressures, suggesting that the CH4 molecules adsorbed on the membrane surface inhibit the CO2 permeation.
As shown in Figure 5, the molecular concentration of CH4 in the membrane tended to increase at high pressure and was almost the same above 4.0 MPa. On the other hand, the molecular concentration of CO2 in the membrane decreased at high pressure. This may be due to the fact that the permeation of CO2 molecules is inhibited by CH4 molecules, which diffuse more slowly in the pores, resulting in blocking the CO2 diffusion at higher pressures. Krishina et al. reported that the diffusion of gas molecules in zeolite membranes was strongly influenced by slower diffusing molecules [26]. Figure 6 shows the results of the analysis of the gas permeation velocity in the membrane with respect to the direction of permeation: at 0.5 MPa, CO2 molecules are permeated faster, while at higher pressures they are permeated more slowly. The higher the pressure, the slower the permeation velocity. On the other hand, the opposite trend was observed for CH4. At low pressure, the CH4 permeation is inhibited because CO2 molecules occupy the pores, but as the number of CH4 adsorbed in the membrane increases, the permeation of CH4 molecules, which diffuse more slowly, is considered to be the rate-limiting factor. A permeation test by Li et al. revealed that CO2 inhibits CH4 adsorption, and CH4 slows the CO2 diffusion in the SAPO-34 membrane that has pores that are similar to CHA [27]. In other words, it is clear that the contribution of diffusion within the zeolite crystals is greater than adsorption at the membrane surface to the reduction in the CO2 permeability in the CHA membrane at high pressure. Indeed, the adsorption selectivity of CHA-type zeolites in an equimolar CO2/CH4 gas mixture at 8.0 MPa is around eight [19,28], suggesting that the NEMD simulation results at 8.0 MPa show little CO2 selectivity due to diffusivity.

3.2. Polycrystalline Model

NEMD simulations of CO2 and CH4 mixtures were performed using a polycrystalline membrane model at 8.0 MPa. Similar to our previously reported results [21], selective aggregation of CO2 molecules was observed in the grain boundaries opened to the membrane surface. After 20 ns, when the number of permeated molecules becomes steady, CO2 molecules are selectively adsorbed at the grain boundaries rather than in the crystalline regions, as shown in Figure 7, and the permeability in the grain boundary membrane model is estimated to be a significant contribution of gas permeation.
Figure 8 shows the permeability and separation factors for the polycrystalline and perfectly crystalline membranes. The gas permeability in the polycrystalline model was more than ten times larger than in the perfectly crystalline membrane. The CO2 selectivity in the polycrystalline model at low pressure was smaller than in the perfectly crystalline model. This is due to the grain boundary being larger than the kinetic diameter of CH4, which reduced the molecular sieving effect, consistent with the trend obtained experimentally [8]. However, at pressures above 4.0 MPa in Figure 8c, the CO2 selectivity in the polycrystalline model was greater than that of the perfectly crystalline model. This is because the polycrystalline model has major CO2 permeation in the grain boundary region and is less susceptible to CO2 permeation inhibition by CH4, which occurs in the crystalline region. When separating CO2 from CH4 at high pressure using a high-density crystal membrane such as the perfectly crystalline model in this study, CO2 molecules have difficulty overtaking CH4 molecules in the CHA crystal. This is thought to be due to the one-dimensional diffusion of CO2 molecules between pores. On the other hand, the structure of the polycrystalline membrane modeled in this study is slit-like, which allows surface diffusion (two-dimensional diffusion) of CO2 molecules to overtake CH4 molecules. The illustration of mechanism of the effect by a grain boundary is shown in Figure 9. It is very difficult for CO2 molecules to overtake CH4 molecules on the permeation pathway in zeolite crystals because of kinetic diameters. On the other hand, in the presence of grain boundaries, there are permeation pathways from within the zeolite crystal to the grain boundary. This means that even if the CH4 molecules occupy the pathway within the crystal, the CO2 molecules can overtake the CH4 molecules by diffusing from the crystal to the grain boundary. Interestingly, the CO2 selectivity of the polycrystalline model at 8.0 MPa is greater than the adsorption selectivity (α = 8) in the CHA zeolite. At high-pressure conditions, the amount of gas adsorbed at the grain boundary is reported to be slightly greater than that of zeolite crystals [29]. Therefore, in CHA membranes with grain boundaries, CO2 molecules are not affected by the permeation-inhibiting effect of CH4 molecules within the crystals, and the polycrystalline membrane may have the ability to show higher separation performance than the perfectly crystalline model under high-pressure conditions. This means that fine-controlled grain boundaries may exceed the separation performance of zeolite crystals.

4. Conclusions

To investigate the CO2 separation performance of CHA zeolite membranes at high pressure, permeation simulations of a CO2/CH4 mixture gas up to 8.0 MPa were carried out using the NEMD. The NEMD results showed that CO2 permeation was blocked by the adsorption of CH4, which diffuses slowly inside the CHA zeolite membrane, and that this blocking effect reduced the CO2 selectivity at high pressure. On the other hand, in the CHA membrane with grain boundaries opened to the membrane surface, gases selectively permeated through the grain boundaries, thus suppressing the blocking effect caused by CH4 in the zeolite crystals. It was also suggested that the active utilization of fine-controlled grain boundaries can maintain CO2 selectivity and improve the permeation flux. This paper shows that grain boundary CHA films exhibit higher gas permeability than fully crystallized ones. These results demonstrate that precise control of the grain boundary structure can separate CO2 more efficiently than conventional zeolite membranes under high pressure.

Author Contributions

Methodology, F.H.; data curation, F.H.; writing—original draft preparation, F.H.; visualization, F.H.; investigation, F.H.; supervision, M.M. and H.T.; writing—reviewing and editing, M.M. and H.T.; conceptualization, H.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rubin, E.S.; Davison, J.E.; Herzog, H.J. The cost of CO2 capture and storage. Int. J. Greenh. Gas Control. 2015, 40, 378–400. [Google Scholar] [CrossRef]
  2. Baker, R.W.; Lokhandwala, K. Natural Gas Processing with Membranes: An Overview. Ind. Eng. Chem. Res. 2008, 47, 2109–2121. [Google Scholar] [CrossRef]
  3. Koros, W.J.; Mahajan, R. Pushing the limits on possibilities for large scale gas separation:  which strategies? J. Membr. Sci. 2000, 175, 181–196. [Google Scholar] [CrossRef]
  4. Kong, X.; Qiu, H.; Meng, D.; Tang, X.; Yang, S.; Guo, W.; Zhang, Y.; Kong, L.; Zhang, Y.; Zhang, Z. Reproducible synthesis of all-silica CHA zeolite membranes in a homogeneous mother liquor. Sep. Purif. Tech. 2021, 274, 119104. [Google Scholar] [CrossRef]
  5. Wang, L.; Zhang, C.; Gao, X.; Peng, L.; Jiang, J.; Gu, X. Preparation of defect-free DDR zeolite membranes by eliminating template with ozone at low temperature. J. Membr. Sci. 2017, 539, 152–160. [Google Scholar] [CrossRef]
  6. Kida, K.; Maeta, Y.; Yogo, K. Preparation and gas permeation properties on pure silica CHA-type zeolite membranes. J. Membr. Sci. 2017, 522, 363–370. [Google Scholar] [CrossRef] [Green Version]
  7. Karakiliç, P.; Wang, X.; Kapteijn, F.; Nijmeijer, A.; Winnubst, L. Defect-free high-silica CHA zeolite membranes with high selectivity for light gas separation. J. Membr. Sci. 2019, 586, 34–43. [Google Scholar] [CrossRef]
  8. Araki, S.; Okubo, Y.; Maekawa, K.; Imasaka, S.; Yamamoto, H. Preparation of a high-silica chabazite-type zeolite membrane with high CO2 permeability using tetraethylammonium hydroxide. J. Membr. Sci. 2020, 613, 118480. [Google Scholar] [CrossRef]
  9. Kida, K.; Maeta, Y.; Yogo, K. Pure silica CHA-type zeolite membranes for dry and humidified CO2/CH4 mixtures separation. Sep. Purif. Tech. 2018, 197, 116–121. [Google Scholar] [CrossRef]
  10. Zhou, J.; Gao, F.; Sun, K.; Jin, X.; Zhang, Y.; L, B.; Zhou, R. Green Synthesis of Highly CO2-Selective CHA Zeolite Membranes in All-Silica and Fluoride-Free Solution for CO2/CH4 Separations. Energy Fuels 2020, 34, 11307–11314. [Google Scholar] [CrossRef]
  11. Xomeritakis, G.; Gouzinis, A.; Nair, S.; Okubo, T.; He, M.; Overney, R.M.; Tsapatsis, M. Growth, microstructure, and permeation properties of supported zeolite (MFI) films and membranes prepared by secondary growth. Chem. Eng. Sci. 1999, 54, 3521–3531. [Google Scholar] [CrossRef]
  12. Korelskiy, D.; Ye, P.; Nabavi, M.S.; Hedlund, J. Selective blocking of grain boundary defects in high-flux zeolite membranes by coking. J. Mater. Chem. A 2017, 5, 7295–7299. [Google Scholar] [CrossRef] [Green Version]
  13. Choi, J.; Jeong, H.; Snyder, M.A.; Stoenger, J.A.; Masel, R.I.; Tsapatisis, M. Grain Boundary Defect Elimination in a Zeolite Membrane by Rapid Thermal Processing. Science 2009, 325, 590–593. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, J.; Jang, E.; Hong, S.; Kim, D.; Kim, E.; Ricther, H.; Simon, A.; Choi, N.; Korelskiy, D.; Fouladvand, S.; et al. Microstructural control of a SSZ-13 zeolite film via rapid thermal processing. J. Membr. Sci. 2019, 591, 117342. [Google Scholar] [CrossRef]
  15. Hong, S.; Kim, D.; Richter, H.; Moon, J.; Choi, N.; Nam, J.; Choi, J. Quantitative elucidation of the elusive role of defects in polycrystalline MFI zeolite membranes on xylene separation performance. J. Membr. Sci. 2019, 569, 91–103. [Google Scholar] [CrossRef]
  16. Park, S.; Lee, M.; Hong, S.; Jeong, Y.; Kim, D.; Choi, N.; Nam, J.; Baik, H.; Choi, J. Low-temperature ozone treatment for p-xylene perm-selective MFI type zeolite membranes: Unprecedented revelation of performance-negating cracks larger than 10 nm in polycrystalline membrane structures. J. Membr. Sci. 2023, 668, 121212. [Google Scholar] [CrossRef]
  17. Takaba, H.; Yamamoto, A.; Hayamizu, K.; Oumi, Y.; Sano, T.; Akiba, E.; Nakao, S. Dependence of the diffusion coefficients of methane in silicalite on diffusion distance as investigated by 1H PFG NMR. Chem. Phys. Lett. 2004, 393, 87–91. [Google Scholar] [CrossRef]
  18. Takaba, H.; Yamamoto, A.; Hayamizu, K.; Nakao, S. Gas Diffusion in Polycrystalline Silicalite Membranes Investigated by 1H Pulse Field-Gradient NMR. J. Phys. Chem. B 2005, 109, 13871–13876. [Google Scholar] [CrossRef] [PubMed]
  19. Arya, G.; Chang, H.; Maginn, E.J. Acritical comparison of equilibrium, non-equilibrium and boundary-driven molecular dynamics techniques for studying transport in microporous materials. J. Chem. Phys. 2001, 115, 8112–8124. [Google Scholar] [CrossRef]
  20. Yoshioka, T.; Tsuru, T.; Asaeda, M. Molecular dynamics studies on gas permeation properties through microporous silica membranes. Sep. Purif. Tech. 2001, 25, 441–449. [Google Scholar] [CrossRef]
  21. Hirosawa, F.; Miyagawa, M.; Takaba, H. Selectivity enhancement by the presence of grain boundary in chabazite zeolite membranes investigated by non-equilibrium molecular dynamics. J. Membr. Sci. 2021, 632, 119348. [Google Scholar] [CrossRef]
  22. Yoshioka, T.; Asaeda, M.; Tsuru, T. A molecular dynamics simulation of pressure-driven gas permeation in a micropore potential field on silica membranes. J. Membr. Sci. 2007, 293, 81–93. [Google Scholar] [CrossRef] [Green Version]
  23. Nagumo, R.; Takaba, H.; Nakao, S. A Methodology to Estimate Transport Diffusivities in ‘Single-File’ Permeation through Zeolite Membranes Using Molecular Simulations. J. Chem. Eng. Jpn. 2007, 40, 1045–1055. [Google Scholar] [CrossRef]
  24. Miyamoto, M.; Fujioka, Y.; Yogo, K. Pure silica CHA type zeolite for CO2 separation using pressure swing adsorption at high pressure. J. Mater. Chem. 2012, 22, 20186–20189. [Google Scholar] [CrossRef]
  25. Krishna, R.; van Baten, J.M. Insights into diffusion of gases in zeolites gained from molecular dynamics simulations. Micropor. Mesopor. Mater. 2008, 109, 91–108. [Google Scholar] [CrossRef]
  26. Krishna, R.; van Baten, J.M. Maxwell–Stefan modeling of slowing-down effects in mixed gas permeation across porous membranes. J. Membr. Sci. 2011, 383, 289–300. [Google Scholar] [CrossRef]
  27. Li, S.; Martinek, J.G.; Falconer, J.L.; Noble, R.D.; Gardner, T.Q. High-Pressure CO2/CH4 Separation Using SAPO-34 Membranes. Ind. Eng. Chem. Res. 2005, 44, 3220–3228. [Google Scholar] [CrossRef]
  28. Krishna, R.; van Baten, J.M. Using the spreading pressure to inter-relate the characteristics of unary, binary and ternary mixture permeation across microporous membranes. J. Membr. Sci. 2022, 643, 120049. [Google Scholar] [CrossRef]
  29. Takaba, H.; Yamamoto, A.; Nakao, S. Modeling of methane permeation through a defective region in MFI-type zeolite membranes. Desalination 2006, 192, 82–90. [Google Scholar] [CrossRef]
Figure 1. The cross view of CHA zeolite membrane on NEMD; (a) is a perfectly crystalline model and (b) is a polycrystalline model. The atoms by CPK model represent the silica atoms in yellow, oxygen atoms in red. The layer at the top and the layer at the bottom shows the gas-generation region and the gas-deletion region, respectively.
Figure 1. The cross view of CHA zeolite membrane on NEMD; (a) is a perfectly crystalline model and (b) is a polycrystalline model. The atoms by CPK model represent the silica atoms in yellow, oxygen atoms in red. The layer at the top and the layer at the bottom shows the gas-generation region and the gas-deletion region, respectively.
Membranes 13 00278 g001
Figure 2. The snapshots of the unit cell at different simulated times were obtained from CO2/CH4 mixture gas permeation simulations using NEMD. CO2 and CH4 molecules by the CPK model represented the carbon atoms in gray, oxygen atoms in red, and hydrogen atoms in white, whereas the CHA zeolite was represented by the stick model. (a) 0.5 MPa, (b) 8.0 MPa, (1) 5 ns, (2) 10 ns, (3) 15 ns.
Figure 2. The snapshots of the unit cell at different simulated times were obtained from CO2/CH4 mixture gas permeation simulations using NEMD. CO2 and CH4 molecules by the CPK model represented the carbon atoms in gray, oxygen atoms in red, and hydrogen atoms in white, whereas the CHA zeolite was represented by the stick model. (a) 0.5 MPa, (b) 8.0 MPa, (1) 5 ns, (2) 10 ns, (3) 15 ns.
Membranes 13 00278 g002
Figure 3. Changes in the number of permeated gas molecules through the perfectly crystalline model of CHA zeolite membrane. (a) 0.5 MPa, (b) 8.0 MPa.
Figure 3. Changes in the number of permeated gas molecules through the perfectly crystalline model of CHA zeolite membrane. (a) 0.5 MPa, (b) 8.0 MPa.
Membranes 13 00278 g003
Figure 4. NEMD results on the perfectly crystalline model of zeolite membrane at different pressure conditions. (a) Permeability, (b) separation factor.
Figure 4. NEMD results on the perfectly crystalline model of zeolite membrane at different pressure conditions. (a) Permeability, (b) separation factor.
Membranes 13 00278 g004
Figure 5. Concentration profile of gas molecules in the perfectly crystalline model on NEMD simulation at different pressure conditions. (a) Schematic of membrane model, (b) concentration profile of CO2, and (c) concentration profile of CH4.
Figure 5. Concentration profile of gas molecules in the perfectly crystalline model on NEMD simulation at different pressure conditions. (a) Schematic of membrane model, (b) concentration profile of CO2, and (c) concentration profile of CH4.
Membranes 13 00278 g005
Figure 6. The velocity of gas molecules along permeation direction in a perfectly crystalline model of CHA zeolite membrane on NEMD simulation at different pressure conditions.
Figure 6. The velocity of gas molecules along permeation direction in a perfectly crystalline model of CHA zeolite membrane on NEMD simulation at different pressure conditions.
Membranes 13 00278 g006
Figure 7. The snapshots of the cross view of the membrane in the unit cell were obtained from the NEMD for CO2/CH4 mixture gas permeation at high pressure. CO2 and CH4 molecules by the CPK model represented the carbon atoms in gray, oxygen atoms in red, and hydrogen atoms in white, whereas the CHA zeolite was represented by the stick model. (a) 5 ns, (b) 10 ns, and (c) 15 ns.
Figure 7. The snapshots of the cross view of the membrane in the unit cell were obtained from the NEMD for CO2/CH4 mixture gas permeation at high pressure. CO2 and CH4 molecules by the CPK model represented the carbon atoms in gray, oxygen atoms in red, and hydrogen atoms in white, whereas the CHA zeolite was represented by the stick model. (a) 5 ns, (b) 10 ns, and (c) 15 ns.
Membranes 13 00278 g007
Figure 8. Comparison of NEMD results on the polycrystalline and perfectly crystalline membrane models at different pressure conditions. (a) Permeability, (b) separation factor, (c) enlarged (b) at high pressure.
Figure 8. Comparison of NEMD results on the polycrystalline and perfectly crystalline membrane models at different pressure conditions. (a) Permeability, (b) separation factor, (c) enlarged (b) at high pressure.
Membranes 13 00278 g008
Figure 9. The illustration of the mechanism of effect in permeability by grain boundaries at low and high pressure. (a) Perfectly crystalline, (b) polycrystalline. (1) Schematic of membrane model, (2) low pressure, (3) high pressure, The symbols of the circle, the triangle, and the cross mean the high CO2 selectivity, low CO2 selectivity, and no selectivity, respectively.
Figure 9. The illustration of the mechanism of effect in permeability by grain boundaries at low and high pressure. (a) Perfectly crystalline, (b) polycrystalline. (1) Schematic of membrane model, (2) low pressure, (3) high pressure, The symbols of the circle, the triangle, and the cross mean the high CO2 selectivity, low CO2 selectivity, and no selectivity, respectively.
Membranes 13 00278 g009aMembranes 13 00278 g009b
Table 1. Calculations of the time interval f for generating a gas molecule on each membrane model.
Table 1. Calculations of the time interval f for generating a gas molecule on each membrane model.
f (fs)
0.5 MPa2.0 MPa4.0 MPa6.0 MPa8.0 MPa
Perfectly crystallineCO221,5245381269017941345
CH412,978324516221082811
PolycrystallineCO2914622871143762572
CH455151379689460345
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hirosawa, F.; Miyagawa, M.; Takaba, H. High Efficient CO2 Separation at High Pressure by Grain-Boundary-Controlled CHA Zeolite Membrane Investigated by Non-Equilibrium Molecular Dynamics. Membranes 2023, 13, 278. https://doi.org/10.3390/membranes13030278

AMA Style

Hirosawa F, Miyagawa M, Takaba H. High Efficient CO2 Separation at High Pressure by Grain-Boundary-Controlled CHA Zeolite Membrane Investigated by Non-Equilibrium Molecular Dynamics. Membranes. 2023; 13(3):278. https://doi.org/10.3390/membranes13030278

Chicago/Turabian Style

Hirosawa, Fumiya, Masaya Miyagawa, and Hiromitsu Takaba. 2023. "High Efficient CO2 Separation at High Pressure by Grain-Boundary-Controlled CHA Zeolite Membrane Investigated by Non-Equilibrium Molecular Dynamics" Membranes 13, no. 3: 278. https://doi.org/10.3390/membranes13030278

APA Style

Hirosawa, F., Miyagawa, M., & Takaba, H. (2023). High Efficient CO2 Separation at High Pressure by Grain-Boundary-Controlled CHA Zeolite Membrane Investigated by Non-Equilibrium Molecular Dynamics. Membranes, 13(3), 278. https://doi.org/10.3390/membranes13030278

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