1
|
Liu L, Guan D, Lu Y, Sun M, Liu Y, Zhao J, Yang L. A Molecular Dynamics Study on Xe/Kr Separation Mechanisms Using Crystal Growth Method. ACS OMEGA 2024; 9:25822-25831. [PMID: 38911791 PMCID: PMC11191100 DOI: 10.1021/acsomega.4c00108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 03/31/2024] [Accepted: 05/23/2024] [Indexed: 06/25/2024]
Abstract
The separation of xenon/krypton gas mixtures is a valuable but challenging endeavor in the gas industry due to their similar physical characteristics and closely sized molecules. To address this, we investigated the effectiveness of the hydrate-based gas separation method for mixed Xe-Kr gas via molecular dynamics (MD) simulations. The formation process of hydrates facilitates the encapsulation of guest molecules within hydrate cages, offering a potential strategy for gas separation. Higher temperatures and pressures are advantageous for accelerating the hydrate growth rate. The final occupancy of guest molecules and empty cages within 512, 51264, and all hydrate cages were thoroughly examined. An increase in the pressure and temperature enhanced the occupancy rates of Xe in both 512 and 51264 cages, whereas elevated pressure alone improved the occupancy of Kr in 51264 cages. However, the impact of temperature and pressure on Kr occupancy within 512 cages was found to be minimal. Elevated temperature and pressure resulted in a reduced occupancy of empty cages. Predominantly, 51264 cages were occupied by Xe, whereas Kr showed a propensity to occupy the 512 cages. With increasing simulated pressure, the final occupancy of Xe molecules in all cages rose from 0.37 to 0.41 for simulations at 260 K, while the final occupancy of empty cages decreased from 0.24 to 0.2.
Collapse
Affiliation(s)
- Liangliang Liu
- Shenyang
Aircraft Design Institute Shenyang 110042, China
| | - Dawei Guan
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Yi Lu
- Shenyang
Aircraft Design Institute Shenyang 110042, China
| | - Mingrui Sun
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Yu Liu
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Jiafei Zhao
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Lei Yang
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| |
Collapse
|
2
|
Bassani CL, Engel M, Sum AK. Mesomorphology of clathrate hydrates from molecular ordering. J Chem Phys 2024; 160:190901. [PMID: 38767264 DOI: 10.1063/5.0200516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2024] [Accepted: 03/13/2024] [Indexed: 05/22/2024] Open
Abstract
Clathrate hydrates are crystals formed by guest molecules that stabilize cages of hydrogen-bonded water molecules. Whereas thermodynamic equilibrium is well described via the van der Waals and Platteeuw approach, the increasing concerns with global warming and energy transition require extending the knowledge to non-equilibrium conditions in multiphase, sheared systems, in a multiscale framework. Potential macro-applications concern the storage of carbon dioxide in the form of clathrates, and the reduction of hydrate inhibition additives currently required in hydrocarbon production. We evidence porous mesomorphologies as key to bridging the molecular scales to macro-applications of low solubility guests. We discuss the coupling of molecular ordering with the mesoscales, including (i) the emergence of porous patterns as a combined factor from the walk over the free energy landscape and 3D competitive nucleation and growth and (ii) the role of molecular attachment rates in crystallization-diffusion models that allow predicting the timescale of pore sealing. This is a perspective study that discusses the use of discrete models (molecular dynamics) to build continuum models (phase field models, crystallization laws, and transport phenomena) to predict multiscale manifestations at a feasible computational cost. Several advances in correlated fields (ice, polymers, alloys, and nanoparticles) are discussed in the scenario of clathrate hydrates, as well as the challenges and necessary developments to push the field forward.
Collapse
Affiliation(s)
- Carlos L Bassani
- Institute for Multiscale Simulation, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - Michael Engel
- Institute for Multiscale Simulation, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - Amadeu K Sum
- Phases to Flow Laboratory, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| |
Collapse
|
3
|
Abstract
Clathrate hydrates have diverse crystal structures, and among them, the three (sI, sII, and sH) most prevalent ones cover nearly all known structures, while the norm is to consider other structures only when specific guest molecules are present. Here we report the observation of a hidden clathrate structure: the tetragonal structure (TS-I) in commonly formed gas hydrates, as evidenced from molecular dynamics simulations. We show that when two (or more) sI crystal grains with different growth directions come into contact or when the growth of a sI crystal encounters geometrical frustration, the TS-I results as a cocrystal. We give evidence that TS-I may also play an important role in the combination and/or transition between sI and sII. These results imply that this previously neglected structure may be commonly present whenever sI or sII is formed. This hidden structure must be identified, experimentally and in simulations; confining the possible structures may hinder an in-depth understanding of clathrate hydrates.
Collapse
Affiliation(s)
- Yong Chen
- Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou 510640, P. R. China
- Phases to Flow Laboratory, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, United States
| | - Satoshi Takeya
- National Metrology Institute of Japan (NMIJ), National Institute of Advanced Industrial Science and Technology (AIST), Central 5, 1-1-1 Higashi, Tsukuba 305-8565, Japan
| | - Amadeu K Sum
- Phases to Flow Laboratory, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, United States
| |
Collapse
|
4
|
Liu Y, Sun J, Chen C, Li W, Qin Y, Wang Y. Molecular insights into gas hydrate formation in the presence of graphene oxide solid surfaces. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119309] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
5
|
Lu Y, Lv X, Li Q, Yang L, Zhang L, Zhao J, Song Y. Molecular behavior of hybrid gas hydrate nucleation: separation of soluble H 2S from mixed gas. Phys Chem Chem Phys 2022; 24:9509-9520. [PMID: 35388810 DOI: 10.1039/d1cp05302g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Soluble H2S widely exists in natural gas or oil potentially corroding oil/gas pipelines. Furthermore, it can affect the hydrate formation condition, resulting in pipeline blockage; the nucleation mechanism from mixed gas including H2S is still largely unclear. Molecular dynamics simulations were performed to reveal the effects of different initial mixed H2S/CH4 compositions on the hydrate nucleation and growth process. The geometric details of the nanobubbles and gas composition in the nanobubbles were analyzed; the size of the nanobubbles was found to decrease from 3.4 nm to 1.4 nm. With the increase in the initial H2S proportion, the diameter of the nanobubbles decreased; more guest molecules were dissolved in the water, which improved the initial concentration of guest molecules in the water. A multi-site nucleation process was observed, and separate hydrate clusters could grow independently until the simulation box limited their growth due to high local H2S concentration as a potential nucleation location. When the initial proportion of mixed gas approaches, H2S preferred to occupy and stabilize the incipient cage. Moreover, 512, 4151062, and 51262 cages accounted for approximately 95% of the first hydrate cage. Nucleation rates were shown to increase from 4.62 × 1024 to 9.438 × 1026 nuclei cm-3 s-1. The present high subcooling and H2S concentration provided a high driving force to promote mixed hydrate nucleation and growth. The proportion of cages occupied by H2S increased with increasing initial H2S proportion, but the largest enrichment factor of 1.38 occurred at 10% initial H2S/CH4 mixed gas.
Collapse
Affiliation(s)
- Yi Lu
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, 116024, China.
| | - Xin Lv
- State Key Laboratory of Natural Gas Hydrate, Beijing, 100028, China
| | - Qingping Li
- State Key Laboratory of Natural Gas Hydrate, Beijing, 100028, China
| | - Lei Yang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, 116024, China.
| | - Lunxiang Zhang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, 116024, China.
| | - Jiafei Zhao
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, 116024, China.
| | - Yongchen Song
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, 116024, China.
| |
Collapse
|
6
|
Arjun A, Bolhuis PG. Homogenous nucleation rate of CO 2 hydrates using transition interface sampling. J Chem Phys 2021; 154:164507. [PMID: 33940852 DOI: 10.1063/5.0044883] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Carbon dioxide and water can form solid clathrate structures in which water cages encapsulate the gas molecules. Such hydrates have sparked much interest due to their possible application in CO2 sequestration. How the solid structure forms exactly from the liquid phase via a homogenous nucleation process is still poorly understood. This nucleation event is rare on the molecular timescale even under moderate undercooling or supersaturation conditions because of the large free energy barrier toward crystallization, rendering a brute force simulation of hydrate nucleation unfeasible for moderate undercooling or supersaturation. Here, we perform transition interface sampling simulations to quantify the homogenous nucleation rate for CO2 hydrate formation using accurate atomistic force fields at 500 bars for three different temperatures between 260 and 273 K. Collecting more than 100 000 pathways comprising roughly two milliseconds of simulation time, we computed a nucleation rate in the amorphous phase of ∼1021 nuclei s-1 cm-3 for a temperature of 260 K and a rate of ∼1012 nuclei s-1 cm-3 for a temperature of 265 K. For a temperature of 273 K, we find that the hydrate forms an sI crystalline phase with a rate of order of ∼101 nuclei s-1 cm-3. We compare these rates to classical nucleation theory estimates as well as experiments, and to nucleation rate estimates for methane hydrates and discuss possible causes of the observed differences. Our findings shed light on the kinetics of this important clathrate and should assist in future hydrate formation investigation.
Collapse
Affiliation(s)
- A Arjun
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| | - Peter G Bolhuis
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| |
Collapse
|
7
|
Arjun A, Bolhuis PG. Molecular Understanding of Homogeneous Nucleation of CO 2 Hydrates Using Transition Path Sampling. J Phys Chem B 2021; 125:338-349. [PMID: 33379869 PMCID: PMC7816195 DOI: 10.1021/acs.jpcb.0c09915] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2020] [Revised: 12/11/2020] [Indexed: 11/29/2022]
Abstract
Carbon dioxide hydrate is a solid built from hydrogen-bond stabilized water cages that encapsulate individual CO2 molecules. As potential candidates for reducing greenhouse gases, hydrates have attracted attention from both the industry and scientific community. Under high pressure and low temperature, hydrates are formed spontaneously from a mixture of CO2 and water via nucleation and growth. Yet, for moderate undercooling, i.e., moderate supersaturation, studying hydrate formation with molecular simulations is very challenging due to the high nucleation barriers involved. We investigate the homogeneous nucleation mechanism of CO2 hydrate as a function of temperature using transition path sampling (TPS), which generates ensembles of unbiased dynamical trajectories across the high barrier between the liquid and solid states. The resulting path ensembles reveal that at high driving force (low temperature), amorphous structures are predominantly formed, with 4151062 cages being the most abundant. With increasing temperature, the nucleation mechanism changes, and 51262 becomes the most abundant cage type, giving rise to the crystalline sI structure. Reaction coordinate analysis can reveal the most important collective variable involved in the mechanism. With increasing temperature, we observe a shift from a single feature (size of the nucleus) to a 2-dimensional (size and cage type) variable as the salient ingredient of the reaction coordinate, and then back to only the nucleus size. This finding is in line with the underlying shift from an amorphous to a crystalline nucleation channel. Modeling such complex phase transformations using transition path sampling gives unbiased insight into the molecular mechanisms toward different polymorphs, and how these are determined by thermodynamics and kinetics. This study will be beneficial for researchers aiming to produce such hydrates with different polymorphic forms.
Collapse
Affiliation(s)
- A. Arjun
- van ’t Hoff Institute
for Molecular Sciences, University of Amsterdam, PO Box 94157, 1090 GD Amsterdam, The Netherlands
| | - P. G. Bolhuis
- van ’t Hoff Institute
for Molecular Sciences, University of Amsterdam, PO Box 94157, 1090 GD Amsterdam, The Netherlands
| |
Collapse
|