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
|
Elías-Domínguez A, Alvarado JFJ, Pérez-Villaseñor F, Ortíz-Arroyo A, Castro-Agüero Á, López-Medina F, Medina-Velázquez DY. Computer Simulation of Three-Phase Equilibria for Some Water/ n-Alkane Binary Systems. J Phys Chem B 2024; 128:5044-5054. [PMID: 38727627 DOI: 10.1021/acs.jpcb.4c00992] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2024]
Abstract
In this work, the vapor-liquid-liquid equilibrium (VLLE) of the water/n-pentane, water/n-hexane, water/n-octane, and water/n-decane binary systems is calculated by computer simulation using the NVT-Gibbs ensemble (in the version of three simulation boxes) combined with the configurational bias Monte Carlo method. The combination of both methods, the molecular potential models used, and the simulation details allowed us to calculate the triphasic equilibrium properties of the systems studied: the densities of the three phases in equilibrium, their compositions, and potential energies. In previous works, these simulations were not carried out at a temperature range nor water/n-alkanes systems simulated in this work, probably because they are highly nonideal systems; so, to the best of our knowledge, this is the first time that this phenomenon is studied in detail. The results from VLLE simulations of the water/n-pentane system for temperatures from 343.2 to 435 K, the water/n-hexane system for temperatures from 373.11 to 473.15 K, the water/n-octane system for temperatures from 310.9 to 500 K, and for the water/n-decane system for temperatures from 374.15 to 525 K are reported here. The temperature range was selected in concordance with the experimental data available for an adequate study of the VLLE simulation results. The subcritical densities (vapor and liquid rich in n-alkane phases) at various temperatures fit well with the scaling law and the law of rectilinear diameters, allowing the estimation of upper critical end point temperature and density of the VLLE. The simulation results show a good prediction with experimental data reports in the literature.
Collapse
Affiliation(s)
- Arturo Elías-Domínguez
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Juan F J Alvarado
- Departamento de Ingeniería Química, Instituto Tecnológico de Celaya, Av. Tecnológico y A. García Cubas S/N, Celaya, Guanajuato CP 38010, México
| | - Fernando Pérez-Villaseñor
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Arturo Ortíz-Arroyo
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Ángel Castro-Agüero
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Friné López-Medina
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| | - Dulce Y Medina-Velázquez
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Ángel Solana S/N, San Luis Apizaquito, Apizaco, Tlaxcala CP 90341, México
| |
Collapse
|
3
|
Algaba J, Blazquez S, Feria E, Míguez JM, Conde MM, Blas FJ. Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate. J Chem Phys 2024; 160:164722. [PMID: 38687000 DOI: 10.1063/5.0201306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/03/2024] [Indexed: 05/02/2024] Open
Abstract
In this work, the effects of finite size on the determination of the three-phase coexistence temperature (T3) of the carbon dioxide (CO2) hydrate have been studied by molecular dynamic simulations and using the direct coexistence technique. According to this technique, the three phases involved (hydrate-aqueous solution-liquid CO2) are placed together in the same simulation box. By varying the number of molecules of each phase, it is possible to analyze the effect of simulation size and stoichiometry on the T3 determination. In this work, we have determined the T3 value at 8 different pressures (from 100 to 6000 bar) and using 6 different simulation boxes with different numbers of molecules and sizes. In two of these configurations, the ratio of the number of water and CO2 molecules in the aqueous solution and the liquid CO2 phase is the same as in the hydrate (stoichiometric configuration). In both stoichiometric configurations, the formation of a liquid drop of CO2 in the aqueous phase is observed. This drop, which has a cylindrical geometry, increases the amount of CO2 available in the aqueous solution and can in some cases lead to the crystallization of the hydrate at temperatures above T3, overestimating the T3 value obtained from direct coexistence simulations. The simulation results obtained for the CO2 hydrate confirm the sensitivity of T3 depending on the size and composition of the system, explaining the discrepancies observed in the original work by Míguez et al. [J. Chem Phys. 142, 124505 (2015)]. Non-stoichiometric configurations with larger unit cells show a convergence of T3 values, suggesting that finite-size effects for these system sizes, regardless of drop formation, can be safely neglected. The results obtained in this work highlight that the choice of a correct initial configuration is essential to accurately estimate the three-phase coexistence temperature of hydrates by direct coexistence simulations.
Collapse
Affiliation(s)
- J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - S Blazquez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Feria
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - J M Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y del Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| |
Collapse
|
4
|
Blazquez S, Algaba J, Míguez JM, Vega C, Blas FJ, Conde MM. Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate. J Chem Phys 2024; 160:164721. [PMID: 38686998 DOI: 10.1063/5.0201295] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/01/2024] [Indexed: 05/02/2024] Open
Abstract
Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field.
Collapse
Affiliation(s)
- S Blazquez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - J M Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y del Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| |
Collapse
|
5
|
Algaba J, Blazquez S, Míguez JM, Conde MM, Blas FJ. Three-phase equilibria of hydrates from computer simulation. III. Effect of dispersive interactions in the methane and carbon dioxide hydrates. J Chem Phys 2024; 160:164723. [PMID: 38686999 DOI: 10.1063/5.0201309] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/02/2024] [Indexed: 05/02/2024] Open
Abstract
In this work, the effect of the range of dispersive interactions in determining the three-phase coexistence line of the CO2 and CH4 hydrates has been studied. In particular, the temperature (T3) at which solid hydrate, water, and liquid CO2/gas CH4 coexist has been determined through molecular dynamics simulations using different cutoff values (from 0.9 to 1.6 nm) for dispersive interactions. The T3 of both hydrates has been determined using the direct coexistence simulation technique. Following this method, the three phases in equilibrium are put together in the same simulation box, the pressure is fixed, and simulations are performed at different temperatures T. If the hydrate melts, then T > T3. Conversely, if the hydrate grows, then T < T3. The effect of the cutoff distance on the dissociation temperature has been analyzed at three different pressures for CO2 hydrate: 100, 400, and 1000 bar. Then, we have changed the guest and studied the effect of the cutoff distance on the dissociation temperature of the CH4 hydrate at 400 bar. Moreover, the effect of long-range corrections for dispersive interactions has been analyzed by running simulations with homo- and inhomogeneous corrections and a cutoff value of 0.9 nm. The results obtained in this work highlight that the cutoff distance for the dispersive interactions affects the stability conditions of these hydrates. This effect is enhanced when the pressure is decreased, displacing the T3 about 2-4 K depending on the system and the pressure.
Collapse
Affiliation(s)
- J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J M Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y del Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| |
Collapse
|
6
|
Fernández-Fernández ÁM, Bárcena Á, Conde MM, Pérez-Sánchez G, Pérez-Rodríguez M, Piñeiro MM. Modeling oceanic sedimentary methane hydrate growth through molecular dynamics simulation. J Chem Phys 2024; 160:144107. [PMID: 38591679 DOI: 10.1063/5.0203116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Accepted: 03/23/2024] [Indexed: 04/10/2024] Open
Abstract
The crystallization process of methane hydrates in a confined geometry resembling seabed porous silica sedimentary conditions has been studied using molecular dynamics simulations. With this objective in mind, a fully atomistic quartz silica slit pore has been designed, and the temperature stability of a methane hydrate crystalline seed in the presence of water and guest molecule methane has been analyzed. NaCl ion pairs have been added in different concentrations, simulating salinity conditions up to values higher than average oceanic conditions. The structure obtained when the hydrate crystallizes inside the pore is discussed, paying special attention to the presence of ionic doping inside the hydrate and the subsequent induced structural distortion. The shift in the hydrate stability conditions due to the increasing water salinity is discussed and compared with the case of unconfined hydrate, concluding that the influence of the confinement geometry and pore hydrophilicity produces a larger deviation in the confined hydrate phase equilibria.
Collapse
Affiliation(s)
| | - Álvaro Bárcena
- Dpto. de Física Aplicada, Univ. de Vigo, Vigo 36310, Spain
| | - María M Conde
- Dpto. de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, Madrid 28006, Spain
| | | | | | | |
Collapse
|
7
|
Cao P, Wu J, Ning F. Mechanical properties of amorphous CO 2 hydrates: insights from molecular simulations. Phys Chem Chem Phys 2024; 26:9388-9398. [PMID: 38444360 DOI: 10.1039/d4cp00203b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Understanding physicochemical properties of amorphous gas hydrate systems is of great significance to reveal structural stabilities of polycrystalline gas hydrate systems. Furthermore, amorphous gas hydrates can occur ordinarily in the nucleation events of gas hydrate systems. Herein, the mechanical properties of amorphous carbon dioxide hydrates are examined by means of all-atom classical molecular dynamic simulations. Our molecular simulation results reveal that mechanical strengths of amorphous carbon dioxide hydrates are evidently governed by temperatures, confining pressures, and ratios of water to carbon dioxide molecules. Notably, under compressive loads, amorphous carbon dioxide hydrates firstly exhibit monotonic strain hardening, followed by an interesting distinct phenomenon characterized by a steady flow stress at further large deformation strains. Furthermore, structural evolutions of amorphous carbon dioxide hydrates are analyzed on the basis of the N-Hbond DOP order parameter. These important findings can not only contribute to our understanding of the structural stabilities of amorphous gas hydrate systems, but also help to develop fundamental understandings about grain boundaries of gas hydrate systems.
Collapse
Affiliation(s)
- Pinqiang Cao
- School of Resource and Environmental Engineering, Wuhan University of Science and Technology, Wuhan, Hubei 430081, China.
| | - Jianyang Wu
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Fujian Provincial Key Laboratory for Soft Functional Materials Research, Xiamen University, Xiamen 361005, China.
| | - Fulong Ning
- Faculty of Engineering, China University of Geosciences, Wuhan, Hubei 430074, China.
| |
Collapse
|
8
|
Cai X, Worley J, Phan A, Salvalaglio M, Koh C, Striolo A. Understanding the effect of moderate concentration SDS on CO 2 hydrates growth in the presence of THF. J Colloid Interface Sci 2024; 658:1-11. [PMID: 38091793 DOI: 10.1016/j.jcis.2023.11.136] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 11/17/2023] [Accepted: 11/21/2023] [Indexed: 01/12/2024]
Abstract
Hypothesis Additives like Tetrahydrofuran (THF) and Sodium Dodecylsulfate (SDS) improve Carbon Dioxide (CO2) hydrates thermal stability and growth rate when used separately. It has been hypothesised that combining them could improve the kinetics of growth and the thermodynamic stability of CO2 hydrates. Simulations and Experiments We exploit atomistic molecular dynamics simulations to investigate the combined impact of THF and SDS under different temperatures and concentrations. The simulation insights are verified experimentally using pendant drop tensiometry conducted at ambient pressures and high-pressure differential scanning calorimetry. Findings Our simulations revealed that the combination of both additives is synergistic at low temperatures but antagonistic at temperatures above 274.1 K due to the aggregation of SDS molecules induced by THF molecules. These aggregates effectively remove THF and CO2 from the hydrate-liquid interface, thereby reducing the driving force for hydrates growth. Experiments revealed that the critical micelle concentration of SDS in water decreases by 20% upon the addition of THF. Further experiments in the presence of THF showed that only small amounts of SDS are sufficient to increase the CO2 storage efficiency by over 40% compared to results obtained without promoters. Overall, our results provide microscopic insights into the mechanisms of THF and SDS promoters on CO2 hydrates, useful for determining the optimal conditions for hydrate growth.
Collapse
Affiliation(s)
- Xinrui Cai
- Thomas Young Centre and Department of Chemical Engineering, University College London, Torrington Place, London, WC1E 7JE, United Kingdom
| | - Joshua Worley
- Department of Chemical and Biological Engineering, Colorado School of Mines, Golden, CO 80401, United States
| | - Anh Phan
- School of Chemistry and Chemical Engineering, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford, Surrey GU2 7XH, United Kingdom
| | - Matteo Salvalaglio
- Thomas Young Centre and Department of Chemical Engineering, University College London, Torrington Place, London, WC1E 7JE, United Kingdom
| | - Carolyn Koh
- Department of Chemical and Biological Engineering, Colorado School of Mines, Golden, CO 80401, United States
| | - Alberto Striolo
- Thomas Young Centre and Department of Chemical Engineering, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; School of Sustainable Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, OK 73019, United States.
| |
Collapse
|
9
|
Tanaka H, Matsumoto M, Yagasaki T. Cage occupancies of CH4, CO2, and Xe hydrates: Mean field theory and grandcanonical Monte Carlo simulations. J Chem Phys 2024; 160:044502. [PMID: 38258922 DOI: 10.1063/5.0188679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 01/02/2024] [Indexed: 01/24/2024] Open
Abstract
We propose a statistical mechanical theory for the thermodynamic stability of clathrate hydrates, considering the influence of the guest-guest interaction on the occupancies of the cages. A mean field approximation is developed to examine the magnitude of the influence. Our new method works remarkably well, which is manifested by two sorts of grandcanonical Monte Carlo (GCMC) simulations. One is full GCMC, and the other is designed in the present study for clathrate hydrates, called lattice-GCMC, in which each guest can be adsorbed at one of the centers of the cage. In the latter simulation, only the guest-guest interaction is explicitly treated, incorporating the host-guest interaction into the free energy of the cage occupation without other guests. Critical phenomena for guest species, such as large density fluctuations, are observed when the temperature is low or the guest-guest interaction is strong.
Collapse
Affiliation(s)
- Hideki Tanaka
- Toyota Physical and Chemical Research Institute, Nagakute 480-1192, Japan
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Osaka 560-8531, Japan
| |
Collapse
|
10
|
Tanaka H, Matsumoto M, Yagasaki T. Efficiency and energy balance for substitution of CH4 in clathrate hydrates with CO2 under multiple-phase coexisting conditions. J Chem Phys 2023; 159:194504. [PMID: 37987520 DOI: 10.1063/5.0179655] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2023] [Accepted: 10/30/2023] [Indexed: 11/22/2023] Open
Abstract
Many experimental and theoretical studies on CH4-CO2 hydrates have been performed aiming at the extraction of CH4 as a relatively clean energy resource and concurrent sequestration of CO2. However, vague or insufficient characterization of the environmental conditions prevents us from a comprehensive understanding of even equilibrium properties of CH4-CO2 hydrates for this substitution. We propose possible reaction schemes for the substitution, paying special attention to the coexisting phases, the aqueous and/or the fluid, where CO2 is supplied from and CH4 is transferred to. We address the two schemes for the substitution operating in three-phase and two-phase coexistence. Advantages and efficiencies of extracting CH4 in the individual scheme are estimated from the chemical potentials of all the components in all the phases involved in the substitution on the basis of a statistical mechanical theory developed recently. It is found that although substitution is feasible in the three-phase coexistence, its working window in temperature-pressure space is much narrower compared to the two-phase coexistence condition. Despite that the substitution normally generates only a small amount of heat, a large endothermic substitution is suggested in the medium pressure range, caused by the vaporization of liquid CO2 due to mixing with a small amount of the released CH4. This study provides the first theoretical framework toward the practical use of hydrates replacing CH4 with CO2 and serves as a basis for quantitative planning.
Collapse
Affiliation(s)
- Hideki Tanaka
- Toyota Physical and Chemical Research Institute, Nagakute 480-1192, Japan
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Osaka 560-8531, Japan
| |
Collapse
|
11
|
Hao X, Li C, Meng Q, Sun J, Huang L, Bu Q, Li C. Molecular Dynamics Simulation of the Three-Phase Equilibrium Line of CO 2 Hydrate with OPC Water Model. ACS OMEGA 2023; 8:39847-39854. [PMID: 37901483 PMCID: PMC10601413 DOI: 10.1021/acsomega.3c05673] [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: 08/16/2023] [Accepted: 09/26/2023] [Indexed: 10/31/2023]
Abstract
The three-phase coexistence line of the CO2 hydrate was determined using molecular dynamics (MD) simulations. By using the classical and modified Lorentz-Berthelot (LB) parameters, the simulations were carried out at 10 different pressures from 3 to 500 MPa. For the OPC water model, simulations with the classic and the modified LB parameters both showed negative deviations from the experimental values. For the TIP4P/Ice water model, good agreement with experimental equilibrium data can be achieved when the LB parameter is adjusted based on the solubility of CO2 in water. Our results also show that the influence of the water model on the equilibrium prediction is much larger than the CO2 model. Current simulations indicated that the H2O-H2O and H2O-CO2 cross-interactions' parameters might contribute equally to the accurate prediction of T3. According to our simulations, the prediction of T3 values showed relatively higher accuracy while using the combination of TIP4P/Ice water and EPM2 CO2 with modified LB parameter. Furthermore, varied χ values are recommended for accurate T3 estimation over a wide pressure range. The knowledge obtained in this study will be helpful for further accurate MD simulation of the process of CO2/CH4 replacement.
Collapse
Affiliation(s)
- Xiluo Hao
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Chengfeng Li
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Qingguo Meng
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Jianye Sun
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Li Huang
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Qingtao Bu
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Congying Li
- Center
of Deep Sea Research, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| |
Collapse
|
12
|
Tanaka H, Matsumoto M, Yagasaki T. On the phase behaviors of CH4-CO2 binary clathrate hydrates: Two-phase and three-phase coexistences. J Chem Phys 2023; 158:2895252. [PMID: 37290087 DOI: 10.1063/5.0155143] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Accepted: 05/23/2023] [Indexed: 06/10/2023] Open
Abstract
We develop a statistical mechanical theory on clathrate hydrates in order to explore the phase behaviors of clathrate hydrates containing two kinds of guest species and apply it to CH4-CO2 binary hydrates. The two boundaries separating water and hydrate and hydrate and guest fluid mixtures are estimated, which are extended to the lower temperature and the higher pressure region far distant from the three-phase coexisting conditions. The chemical potentials of individual guest components can be calculated from free energies of cage occupations, which are available from intermolecular interactions between host water and guest molecules. This allows us to derive all thermodynamic properties pertinent to the phase behaviors in the whole space of thermodynamic variables of temperature, pressure, and guest compositions. It is found that the phase boundaries of CH4-CO2 binary hydrates with water and with fluid mixtures locate between simple CH4 and CO2 hydrates, but the composition ratios of CH4 guests in hydrates are disproportional to those in fluid mixtures. Such differences arise from the affinities of each guest species to the large and small cages of CS-I hydrates and significantly affect occupation of each cage type, which results in a deviation of the guest composition in hydrates from that in fluid on the two-phase equilibrium conditions. The present method provides a basis for the evaluation of the efficiency of the guest CH4 replacement to CO2 at the thermodynamic limit.
Collapse
Affiliation(s)
- Hideki Tanaka
- Toyota Physical and Chemical Research Institute, Nagakute 480-1192, Japan
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Osaka 560-8531, Japan
| |
Collapse
|
13
|
Phan A, Stamatakis M, Koh CA, Striolo A. Microscopic insights on clathrate hydrate growth from non-equilibrium molecular dynamics simulations. J Colloid Interface Sci 2023; 649:185-193. [PMID: 37348338 DOI: 10.1016/j.jcis.2023.06.032] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2023] [Revised: 05/03/2023] [Accepted: 06/06/2023] [Indexed: 06/24/2023]
Abstract
Clathrate hydrates form and grow at interfaces. Understanding the relevant molecular processes is crucial for developing hydrate-based technologies. Many computational studies focus on hydrate growth within the aqueous phase using the 'direct coexistence method', which is limited in its ability to investigate hydrate film growth at hydrocarbon-water interfaces. To overcome this shortcoming, a new simulation setup is presented here, which allows us to study the growth of a methane hydrate nucleus in a system where oil-water, hydrate-water, and hydrate-oil interfaces are all simultaneously present, thereby mimicking experimental setups. Using this setup, hydrate growth is studied here under the influence of two additives, a polyvinylcaprolactam oligomer and sodium dodecyl sulfate, at varying concentrations. Our results confirm that hydrate films grow along the oil-water interface, in general agreement with visual experimental observations; growth, albeit slower, also occurs at the hydrate-water interface, the interface most often interrogated via simulations. The results obtained demonstrate that the additives present within curved interfaces control the solubility of methane in the aqueous phase, which correlates with hydrate growth rate. Building on our simulation insights, we suggest that by combining data for the potential of mean force profile for methane transport across the oil-water interface and for the average free energy required to perturb a flat interface, it is possible to predict the performance of additives used to control hydrate growth. These insights could be helpful to achieve optimal methane storage in hydrates, one of many applications which are attracting significant fundamental and applied interests.
Collapse
Affiliation(s)
- Anh Phan
- School of Chemistry and Chemical Engineering, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford, Surrey GU2 7XH, UK.
| | - Michail Stamatakis
- Department of Chemical Engineering, University College London, London WC1E 7JE, UK
| | - Carolyn A Koh
- Center for Hydrate Research, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, CO 80401, United States
| | - Alberto Striolo
- Department of Chemical Engineering, University College London, London WC1E 7JE, UK; School of Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, OK 73019, United States.
| |
Collapse
|
14
|
Romero-Guzmán C, Zerón IM, Algaba J, Mendiboure B, Míguez JM, Blas FJ. Effect of pressure on the carbon dioxide hydrate-water interfacial free energy along its dissociation line. J Chem Phys 2023; 158:2890475. [PMID: 37184014 DOI: 10.1063/5.0139699] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 04/28/2023] [Indexed: 05/16/2023] Open
Abstract
We investigate the effect of pressure on the carbon dioxide (CO2) hydrate-water interfacial free energy along its dissociation line using advanced computer simulation techniques. In previous works, we have determined the interfacial energy of the hydrate at 400 bars using the TIP4P/Ice and TraPPE molecular models for water and CO2, respectively, in combination with two different extensions of the Mold Integration technique [J. Colloid Interface Sci. 623, 354 (2022) and J. Chem. Phys. 157, 134709 (2022)]. Results obtained from computer simulation, 29(2) and 30(2) mJ/m2, are found to be in excellent agreement with the only two measurements that exist in the literature, 28(6) mJ/m2 determined by Uchida et al. [J. Phys. Chem. B 106, 8202 (2002)] and 30(3) mJ/m2 determined by Anderson et al. [J. Phys. Chem. B 107, 3507 (2002)]. Since the experiments do not allow to obtain the variation of the interfacial energy along the dissociation line of the hydrate, we extend our previous studies to quantify the effect of pressure on the interfacial energy at different pressures. Our results suggest that there exists a correlation between the interfacial free energy values and the pressure, i.e., it decreases with the pressure between 100 and 1000 bars. We expect that the combination of reliable molecular models and advanced simulation techniques could help to improve our knowledge of the thermodynamic parameters that control the interfacial free energy of hydrates from a molecular perspective.
Collapse
Affiliation(s)
- Cristóbal Romero-Guzmán
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Iván M Zerón
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Jesús Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Bruno Mendiboure
- Laboratoire des Fluides Complexes et Leurs Réservoirs, UMR5150, Université de Pau et des Pays de l'Adour, B.P. 1155, Pau Cedex 64014, France
| | - José Manuel Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Felipe J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| |
Collapse
|
15
|
Algaba J, Zerón IM, Míguez JM, Grabowska J, Blazquez S, Sanz E, Vega C, Blas FJ. Solubility of carbon dioxide in water: Some useful results for hydrate nucleation. J Chem Phys 2023; 158:2889490. [PMID: 37158326 DOI: 10.1063/5.0146618] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/18/2023] [Indexed: 05/10/2023] Open
Abstract
In this paper, the solubility of carbon dioxide (CO2) in water along the isobar of 400 bar is determined by computer simulations using the well-known TIP4P/Ice force field for water and the TraPPE model for CO2. In particular, the solubility of CO2 in water when in contact with the CO2 liquid phase and the solubility of CO2 in water when in contact with the hydrate have been determined. The solubility of CO2 in a liquid-liquid system decreases as the temperature increases. The solubility of CO2 in a hydrate-liquid system increases with temperature. The two curves intersect at a certain temperature that determines the dissociation temperature of the hydrate at 400 bar (T3). We compare the predictions with T3 obtained using the direct coexistence technique in a previous work. The results of both methods agree, and we suggest 290(2) K as the value of T3 for this system using the same cutoff distance for dispersive interactions. We also propose a novel and alternative route to evaluate the change in chemical potential for the formation of hydrates along the isobar. The new approach is based on the use of the solubility curve of CO2 when the aqueous solution is in contact with the hydrate phase. It considers rigorously the non-ideality of the aqueous solution of CO2, providing reliable values for the driving force for nucleation of hydrates in good agreement with other thermodynamic routes used. It is shown that the driving force for hydrate nucleation at 400 bar is larger for the methane hydrate than for the carbon dioxide hydrate when compared at the same supercooling. We have also analyzed and discussed the effect of the cutoff distance of dispersive interactions and the occupancy of CO2 on the driving force for nucleation of the hydrate.
Collapse
Affiliation(s)
- Jesús Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Iván M Zerón
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - José Manuel Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Joanna Grabowska
- Department of Physical Chemistry, Faculty of Chemistry and BioTechMed Center, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
- Dpto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Samuel Blazquez
- Dpto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Dpto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Dpto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Felipe J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| |
Collapse
|
16
|
Phan A, Striolo A. Chemical Promoter Performance for CO 2 Hydrate Growth: A Molecular Perspective. ENERGY & FUELS : AN AMERICAN CHEMICAL SOCIETY JOURNAL 2023; 37:6002-6011. [PMID: 37114945 PMCID: PMC10123660 DOI: 10.1021/acs.energyfuels.3c00472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Revised: 03/25/2023] [Indexed: 06/19/2023]
Abstract
Carbon dioxide (CO2) hydrates, which contain a relatively large amount of captured CO2 (almost 30 wt % of CO2 with the balance being water), represent a promising CO2 sequestration option for climate change mitigation. To facilitate CO2 storage via hydrates, using chemical additives during hydrate formation might help to expedite formation/growth rates, provided the additives do not reduce the storage capacity. Implementing atomistic molecular dynamics, we study the impact of aziridine, pyrrolidine, and tetrahydrofuran (THF) on the kinetics of CO2 hydrate growth/dissociation. Our simulations are validated via reproducing experimental data for CO2 and CO2 + THF hydrates at selected operating conditions. The simulated results show that both aziridine and pyrrolidine could perform as competent thermodynamic and kinetic promoters. Furthermore, aziridine seems to exceed pyrrolidine and THF in expediting the CO2 hydrate growth rates under the same conditions. Our analysis unveils direct correlations between the kinetics of CO2 hydrate growth and a combination of the free energy barrier for desorption of CO2 from the hydrate surface and the binding free energy of chemical additives adsorbed at the growing hydrate substrate. The detailed thermodynamic analysis conducted in both hydrate and aqueous phases reveals molecular-level mechanisms by which CO2 hydrate promoters are active, which could help to enable the implementation of CO2 sequestration in hydrate-bearing reservoirs.
Collapse
Affiliation(s)
- Anh Phan
- School
of Chemistry and Chemical Engineering, Faculty of Engineering and
Physical Sciences, University of Surrey, Guildford, Surrey GU2
7XH, U.K.
| | - Alberto Striolo
- Department
of Chemical Engineering, University College
London, London WC1E 7JE, U.K.
- School
of Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, Oklahoma 73019, United States
| |
Collapse
|
17
|
Collado P, Piñeiro MM, Pérez-Rodríguez M. Molecular Simulation of CO 2 and H 2 Encapsulation in a Nanoscale Porous Liquid. NANOMATERIALS (BASEL, SWITZERLAND) 2023; 13:409. [PMID: 36770368 PMCID: PMC9920367 DOI: 10.3390/nano13030409] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Revised: 01/10/2023] [Accepted: 01/14/2023] [Indexed: 06/18/2023]
Abstract
In this study we analyse from a theoretical perspective the encapsulation of both gaseous H2 and CO2 at different conditions of pressure and temperature in a Type II porous liquid, composed by nanometric scale cryptophane-111 molecules dispersed in dichloromethane, using atomistic molecular dynamics. Gaseous H2 tends to occupy cryptophane-111's cavities in the early stages of the simulation; however, a remarkably greater selectivity of CO2 adsorption can be seen in the course of the simulation. Calculations were performed at ambient conditions first, and then varying temperature and pressure, obtaining some insight about the different adsorption found in each case. An evaluation of the host molecule cavities accessible volume was also performed, based on the guest that occupies the pore. Finally, a discussion between the different intermolecular host-guest interactions is presented, justifying the different selectivity obtained in the molecular simulation calculations. From the results obtained, the feasibility of a renewable separation and storage method for CO2 using these nanometric scale porous liquids is pointed out.
Collapse
|
18
|
Insights into the synergistic effects of metal particles (Ag, Cu, and Fe) and urea on CO2 clathrate hydrate growth using molecular dynamics simulations. Chem Eng Sci 2022. [DOI: 10.1016/j.ces.2022.118194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
19
|
Michalis VK, Economou IG, Stubos AK, Tsimpanogiannis IN. Phase equilibria molecular simulations of hydrogen hydrates via the direct phase coexistence approach. J Chem Phys 2022; 157:154501. [PMID: 36272800 DOI: 10.1063/5.0108738] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We report the three-phase (hydrate-liquid water-vapor) equilibrium conditions of the hydrogen-water binary system calculated with molecular dynamics simulations via the direct phase coexistence approach. A significant improvement of ∼10.5 K is obtained in the current study, over earlier simulation attempts, by using a combination of modifications related to the hydrogen model that include (i) hydrogen Lennard-Jones parameters that are a function of temperature and (ii) the water-guest energy interaction parameters optimized further by using the Lorentz-Berthelot combining rules, based on an improved description of the solubility of hydrogen in water.
Collapse
Affiliation(s)
| | - Ioannis G Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Athanasios K Stubos
- Environmental Research Laboratory, National Center for Scientific Research "Demokritos," 15310 Aghia Paraskevi Attikis, Greece
| | - Ioannis N Tsimpanogiannis
- Chemical Process and Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH), 57001 Thermi-Thessaloniki, Greece
| |
Collapse
|
20
|
Grabowska J, Blazquez S, Sanz E, Zerón IM, Algaba J, Míguez JM, Blas FJ, Vega C. Solubility of Methane in Water: Some Useful Results for Hydrate Nucleation. J Phys Chem B 2022; 126:8553-8570. [PMID: 36222501 DOI: 10.1021/acs.jpcb.2c04867] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this paper, the solubility of methane in water along the 400 bar isobar is determined by computer simulations using the TIP4P/Ice force field for water and a simple LJ model for methane. In particular, the solubility of methane in water when in contact with the gas phase and the solubility of methane in water when in contact with the hydrate has been determined. The solubility of methane in a gas-liquid system decreases as temperature increases. The solubility of methane in a hydrate-liquid system increases with temperature. The two curves intersect at a certain temperature that determines the triple point T3 at a certain pressure. We also determined T3 by the three-phase direct coexistence method. The results of both methods agree, and we suggest 295(2) K as the value of T3 for this system. We also analyzed the impact of curvature on the solubility of methane in water. We found that the presence of curvature increases the solubility in both the gas-liquid and hydrate-liquid systems. The change in chemical potential for the formation of hydrate is evaluated along the isobar using two different thermodynamic routes, obtaining good agreement between them. It is shown that the driving force for hydrate nucleation under experimental conditions is higher than that for the formation of pure ice when compared at the same supercooling. We also show that supersaturation (i.e., concentrations above those of the planar interface) increases the driving force for nucleation dramatically. The effect of bubbles can be equivalent to that of an additional supercooling of about 20 K. Having highly supersaturated homogeneous solutions makes possible the spontaneous formation of the hydrate at temperatures as high as 285 K (i.e., 10K below T3). The crucial role of the concentration of methane for hydrate formation is clearly revealed. Nucleation of the hydrate can be either impossible or easy and fast depending on the concentration of methane which seems to play the leading role in the understanding of the kinetics of hydrate formation.
Collapse
Affiliation(s)
- Joanna Grabowska
- Departamento Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain.,Department of Physical Chemistry, Faculty of Chemistry and BioTechMed Center, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Samuel Blazquez
- Departamento Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Iván M Zerón
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Jesús Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - José Manuel Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Felipe J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - Carlos Vega
- Departamento Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
21
|
Cabrera-Ramírez A, Prosmiti R. Modeling of Structure H Carbon Dioxide Clathrate Hydrates: Guest-Lattice Energies, Crystal Structure, and Pressure Dependencies. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:14832-14842. [PMID: 36110497 PMCID: PMC9465682 DOI: 10.1021/acs.jpcc.2c04140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Revised: 07/26/2022] [Indexed: 06/15/2023]
Abstract
We performed first-principles computations to investigate the complex interplay of molecular interaction energies in determining the lattice structure and stability of CO2@sH clathrate hydrates. Density functional theory computations using periodic boundary conditions were employed to characterize energetics and the key structural properties of the sH clathrate crystal under pressure, such as equilibrium lattice volume and bulk modulus. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. Structural relaxations of the fully CO2-filled and empty sH unit cells yield crystal structure and lattice energies, while their compressibility parameters were derived by including the pressure dependencies. The present quantum chemistry computations suggest anisotropy in the compressibility of the sH clathrate hydrates, with the crystal being less compressible along the a-axis direction than along the c-axis one, in distinction from nearly isotropic sI and sII structures. The detailed results presented here give insight into the complex nature of the underlying guest-host interactions, checking earlier assumptions, providing critical tests, and improving estimates. Such entries may eventually lead to better predictions of thermodynamic properties and formation conditions, with a direct impact on emerging hydrate-based technologies.
Collapse
Affiliation(s)
- Adriana Cabrera-Ramírez
- Institute
of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006, Madrid, Spain
- Doctoral
Programme in Theoretical Chemistry and Computational Modelling, Doctoral
School, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | - Rita Prosmiti
- Institute
of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006, Madrid, Spain
| |
Collapse
|
22
|
Fernández-Fernández ÁM, Conde MM, Pérez-Sánchez G, Pérez-Rodríguez M, Piñeiro MM. Molecular simulation of methane hydrate growth confined into a silica pore. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
23
|
Algaba J, Acuña E, Míguez JM, Mendiboure B, Zerón IM, Blas FJ. Simulation of the carbon dioxide hydrate-water interfacial energy. J Colloid Interface Sci 2022; 623:354-367. [PMID: 35594594 DOI: 10.1016/j.jcis.2022.05.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Revised: 05/04/2022] [Accepted: 05/05/2022] [Indexed: 10/18/2022]
Abstract
HYPOTHESIS Carbon dioxide hydrates are ice-like nonstoichiometric inclusion solid compounds with importance to global climate change, and gas transportation and storage. The thermodynamic and kinetic mechanisms that control carbon dioxide nucleation critically depend on hydrate-water interfacial free energy. Only two independent indirect experiments are available in the literature. Interfacial energies show large uncertainties due to the conditions at which experiments are performed. Under these circumstances, we hypothesize that accurate molecular models for water and carbon dioxide combined with computer simulation tools can offer an alternative but complementary way to estimate interfacial energies at coexistence conditions from a molecular perspective. CALCULATIONS We have evaluated the interfacial free energy of carbon dioxide hydrates at coexistence conditions (three-phase equilibrium or dissociation line) implementing advanced computational methodologies, including the novel Mold Integration methodology. Our calculations are based on the definition of the interfacial free energy, standard statistical thermodynamic techniques, and the use of the most reliable and used molecular models for water (TIP4P/Ice) and carbon dioxide (TraPPE) available in the literature. FINDINGS We find that simulations provide an interfacial energy value, at coexistence conditions, consistent with the experiments from its thermodynamic definition. Our calculations are reliable since are based on the use of two molecular models that accurately predict: (1) The ice-water interfacial free energy; and (2) the dissociation line of carbon dioxide hydrates. Computer simulation predictions provide alternative but reliable estimates of the carbon dioxide interfacial energy. Our pioneering work demonstrates that is possible to predict interfacial energies of hydrates from a truly computational molecular perspective and opens a new door to the determination of free energies of hydrates.
Collapse
Affiliation(s)
- Jesús Algaba
- Department of Chemical Engineering, South Kensington Campus, Imperial College London, SW7 2AZ London, United Kingdom
| | - Esteban Acuña
- Laboratorio de Simulacion Molecular y Quimica Computacional, CIQSO-Centro de Investigacion en Quimica Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - José Manuel Míguez
- Laboratorio de Simulacion Molecular y Quimica Computacional, CIQSO-Centro de Investigacion en Quimica Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - Bruno Mendiboure
- Laboratoire des Fluides Complexes et Leurs Reservoirs, UMR5150, Universite de Pau et des Pays de l'Adour, B. P. 1155, Pau Cedex 64014, France
| | - Iván M Zerón
- Laboratorio de Simulacion Molecular y Quimica Computacional, CIQSO-Centro de Investigacion en Quimica Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain
| | - Felipe J Blas
- Laboratorio de Simulacion Molecular y Quimica Computacional, CIQSO-Centro de Investigacion en Quimica Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21007 Huelva, Spain.
| |
Collapse
|
24
|
Kang DW, Lee W, Ahn YH, Lee JW. Exploring tuning phenomena of THF-H2 hydrates via molecular dynamics simulations. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118490] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
25
|
Study the effect of Ag nanoparticles on the kinetics of CO2 hydrate growth by molecular dynamics simulation. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117668] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
26
|
Fernández-Fernández ÁM, Pérez-Rodríguez M, Piñeiro MM. Molecular dynamics of fluoromethane type I hydrates. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116720] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
27
|
Jiao L, Wang Z, Li J, Zhao P, Wan R. Stability and dissociation studies of CO2 hydrate under different systems using molecular dynamic simulations. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116788] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
28
|
Honda K, Fujikawa R, Ma X, Yamamoto N, Fujiwara K, Kaneko A, Abe Y. The formation and growth model of a
CO
2
hydrate layer based on molecular dynamics. AIChE J 2021. [DOI: 10.1002/aic.17406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Kota Honda
- Department of Engineering Mechanics and Energy University of Tsukuba Ibaraki Japan
| | - Rintaro Fujikawa
- Department of Engineering Mechanics and Energy University of Tsukuba Ibaraki Japan
| | - Xiao Ma
- Department of Engineering Mechanics and Energy University of Tsukuba Ibaraki Japan
| | - Norifumi Yamamoto
- Department of Applied Chemistry, Faculty of Engineering Chiba Institute of Technology Chiba Japan
| | - Kota Fujiwara
- Department of Engineering Mechanics and Energy University of Tsukuba Ibaraki Japan
| | - Akiko Kaneko
- Department of Engineering Mechanics and Energy University of Tsukuba Ibaraki Japan
| | - Yutaka Abe
- Department of Engineering Mechanics and Energy University of Tsukuba Ibaraki Japan
| |
Collapse
|
29
|
Tsimpanogiannis IN. A novel hybrid method for the calculation of methane hydrate-water interfacial tension along the three-phase (hydrate-liquid water-vapor) equilibrium line. J Chem Phys 2021; 155:024702. [PMID: 34266278 DOI: 10.1063/5.0051383] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
We use a novel hybrid method to explore the temperature dependence of the solid-liquid interfacial tension of a system that consists of solid methane hydrate and liquid water. The calculated values along the three-phase (hydrate-liquid water-vapor) equilibrium line are obtained through the combination of available experimental measurements and computational results that are based on approaches at the atomistic scale, including molecular dynamics and Monte Carlo. An extensive comparison with available experimental and computational studies is performed, and a critical assessment and re-evaluation of previously reported data is presented.
Collapse
Affiliation(s)
- Ioannis N Tsimpanogiannis
- Chemical Process & Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH), 57001 Thermi-Thessaloniki, Greece
| |
Collapse
|
30
|
Cabrera-Ramírez A, Arismendi-Arrieta DJ, Valdés Á, Prosmiti R. Exploring CO 2 @sI Clathrate Hydrates as CO 2 Storage Agents by Computational Density Functional Approaches. Chemphyschem 2021; 22:359-369. [PMID: 33368985 DOI: 10.1002/cphc.202001035] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2020] [Indexed: 12/21/2022]
Abstract
The formation of specific clathrate hydrates and their transformation at given thermodynamic conditions depends on the interactions between the guest molecule/s and the host water lattice. Understanding their structural stability is essential to control structure-property relations involved in different technological applications. Thus, the energetic aspects relative to CO2 @sI clathrate hydrate are investigated through the computation of the underlying interactions, dominated by hydrogen bonds and van der Waals forces, from first-principles electronic structure approaches. The stability of the CO2 @sI clathrate is evaluated by combining bottom-up and top-down approaches. Guest-free and CO2 guest-filled aperiodic cages, up to the gradually CO2 occupation of the entire sI periodic unit cells were considered. Saturation, cohesive and binding energies for the systems are determined by employing a variety of density functionals and their performance is assessed. The dispersion corrections on the non-covalent interactions are found to be important in the stabilization of the CO2 @sI energies, with the encapsulation of the CO2 into guest-free/empty cage/lattice being always an energetically favorable process for most of the functionals studied. The PW86PBE functional with XDM or D3(BJ) dispersion corrections predicts a lattice constant in accord to the experimental values available, and simultaneously provides a reliable description for the guest-host interactions in the periodic CO2 @sI crystal, as well as the energetics of its progressive single cage occupancy process. It has been found that the preferential orientation of the single CO2 in the large sI crystal cages has a stabilizing effect on the hydrate, concluding that the CO2 @sI structure is favored either by considering the individual building block cages or the complete sI unit cell crystal. Such benchmark and methodology cross-check studies benefit new data-driven model research by providing high-quality training information, with new insights that indicate the underlying factors governing their structure-driven stability, and triggering further investigations for controlling the stabilization of these promising long-term CO2 storage materials.
Collapse
Affiliation(s)
| | - Daniel J Arismendi-Arrieta
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018, Donostia-San Sebastián, Spain
| | - Álvaro Valdés
- Escuela de Física, Universidad Nacional de Colombia, Sede Medellín, A. A., 3840, Medellíın, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006, Madrid, Spain
| |
Collapse
|
31
|
Cabrera-Ramírez A, Yanes-Rodríguez R, Prosmiti R. Computational density-functional approaches on finite-size and guest-lattice effects in CO 2@sII clathrate hydrate. J Chem Phys 2021; 154:044301. [PMID: 33514100 DOI: 10.1063/5.0039323] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
We performed first-principles computations to investigate guest-host/host-host effects on the encapsulation of the CO2 molecule in sII clathrate hydrates from finite-size clusters up to periodic 3D crystal lattice systems. Structural and energetic properties were first computed for the individual and first-neighbors clathrate-like sII cages, where highly accurate ab initio quantum chemical methods are available nowadays, allowing in this way the assessment of the density functional (DFT) theoretical approaches employed. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. On this basis, structural relaxations of the CO2-filled and empty sII unit cells yield lattice and compressibility parameters comparable to experimental and previous theoretical values available for sII hydrates. According to these data, the CO2 enclathration in the sII clathrate cages is a stabilizing process, either by considering both guest-host and host-host interactions in the complete unit cell or only the guest-water energies for the individual clathrate-like sII cages. CO2@sII clathrates are predicted to be stable whatever the dispersion correction applied and in the case of single cage occupancy are found to be more stable than the CO2@sI structures. Our results reveal that DFT approaches could provide a good reasonable description of the underlying interactions, enabling the investigation of formation and transformation processes as a function of temperature and pressure.
Collapse
Affiliation(s)
| | | | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
32
|
Cea-Klapp E, Míguez JM, Gómez-Álvarez P, Blas FJ, Quinteros-Lama H, Garrido JM. Molecular modelling techniques for predicting liquid-liquid interfacial properties of methanol plus alkane ( n-hexane, n-heptane, n-octane) mixtures. Phys Chem Chem Phys 2020; 22:27121-27133. [PMID: 33225339 DOI: 10.1039/d0cp04823b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this work, the liquid-liquid interfacial properties of methanol plus n-alkane (n-hexane, n-heptane, n-octane) mixtures are investigated at atmospheric pressure by two complementary molecular modelling techniques; namely, molecular dynamic simulations (MD) and density gradient theory (DGT) coupled with the PC-SAFT (perturbed-chain statistical associating fluid theory) equation of state. Furthermore, two molecular models of methanol are used, which are based on a non-polarisable three site approach. On the one hand, is the original (flexible) TraPPE-UA model force field. On the other hand, is the rigid approximation denoted as OPLS/2016. In both cases, n-alkanes are modelled using the TraPPE-UA model. Simulations are performed using the direct coexistence technique in the ensemble. Special attention is paid to the comparison between the estimations obtained from different methanol models, the available experimental data and theoretical calculations. In all cases, the rigid model is capable of predicting the experimental phase equilibrium and interfacial properties accurately. Unsurprisingly, the methanol-rich density and interfacial tension are overestimated using the TraPPE model combined with Lorentz-Berthelot mixing rules for predicting the mixture behaviour. Accurate comparison between MD and DGT plus PC-SAFT requires consideration of the cross-interactions between individual influence parameters and fitting the βij values. This latter aspect is particularly important because it allows the exploitation of the link between the EOS model and the direct molecular simulation of the corresponding fluid. At the same time, it was demonstrated that the key property defining the interfacial tension value is the absolute concentration of methanol in the methanol-rich phase. This behaviour indicates that there are more hydrogens bonded with each other, and they interact favourably with an increasing number of carbon atoms in the alkane.
Collapse
Affiliation(s)
- Esteban Cea-Klapp
- Department of Chemical Engineering, Universidad de Concepción, Concepción 4070386, Chile.
| | | | | | | | | | | |
Collapse
|
33
|
Cabrera-Ramírez A, Arismendi-Arrieta DJ, Valdés Á, Prosmiti R. Structural Stability of the CO 2 @sI Hydrate: a Bottom-Up Quantum Chemistry Approach on the Guest-Cage and Inter-Cage Interactions. Chemphyschem 2020; 21:2618-2628. [PMID: 33001534 DOI: 10.1002/cphc.202000753] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Revised: 09/30/2020] [Indexed: 12/27/2022]
Abstract
Through reliable first-principles computations, we have demonstrated the impact of CO2 molecules enclathration on the stability of sI clathrate hydrates. Given the delicate balance between the interaction energy components (van der Waals, hydrogen bonds) present on such systems, we follow a systematic bottom-up approach starting from the individual 512 and 512 62 sI cages, up to all existing combinations of two-adjacent sI crystal cages to evaluate how such clathrate-like models perform on the evaluation of the guest-host and first-neighbors inter-cage effects, respectively. Interaction and binding energies of the CO2 occupation of the sI cages were computed using DF-MP2 and different DFT/DFT-D electronic structure methodologies. The performance of selected DFT functionals, together with various semi-classical dispersion corrections schemes, were validated by comparison with reference ab initio DF-MP2 data, as well as experimental data from x-ray and neutron diffraction studies available. Our investigation confirms that the inclusion of the CO2 in the cage/s is an energetically favorable process, with the CO2 molecule preferring to occupy the large 512 62 sI cages compared to the 512 ones. Further, the present results conclude on the rigidity of the water cages arrangements, showing the importance of the inter-cage couplings in the cluster models under study. In particular, the guest-cage interaction is the key factor for the preferential orientation of the captured CO2 molecules in the sI cages, while the inter-cage interactions seems to cause minor distortions with the CO2 guest neighbors interactions do not extending beyond the large 512 62 sI cages. Such findings on these clathrate-like model systems are in accord with experimental observations, drawing a direct relevance to the structural stability of CO2 @sI clathrates.
Collapse
Affiliation(s)
| | - Daniel J Arismendi-Arrieta
- Donostia International Physics Center (DIPC), Paseo, Manuel de Lardizabal 4, 20018, Donostia-San Sebastián, Spain
| | - Álvaro Valdés
- Escuela de Física, Universidad Nacional de Colombia, Sede, A. A., 3840, Medellín, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), Serrano 123, Madrid, Spain
| |
Collapse
|
34
|
Tsimpanogiannis IN, Michalis VK, Economou IG. Novel methodology for the calculation of the enthalpy of enclathration of methane hydrates using molecular dynamics simulations. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1711976] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Chemical Process & Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH) Thermi-Thessaloniki, Greece
- Environmental Research Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Vasileios K. Michalis
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Ioannis G. Economou
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| |
Collapse
|
35
|
Sizova AA, Sizov VV, Brodskaya EN. Molecular Dynamics Simulation of the Stability of Spherical Nanoclusters of Methane and Carbon Dioxide Hydrates. COLLOID JOURNAL 2020. [DOI: 10.1134/s1061933x2002012x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
36
|
Castillo-Borja F, Bravo-Sánchez UI, Vázquez-Román R, Díaz-Ovalle CO. Biogas purification via sII hydrates in the presence of THF and DMSO solutions using MD simulations. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2019.111904] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
37
|
Yi L, Zhou X, He Y, Cai Z, Zhao L, Zhang W, Shao Y. Molecular Dynamics Simulation Study on the Growth of Structure II Nitrogen Hydrate. J Phys Chem B 2019; 123:9180-9186. [PMID: 31609605 DOI: 10.1021/acs.jpcb.9b06386] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Crystal growth of N2 hydrate in a three-phase system consisting of N2 hydrate, liquid water, and gaseous N2 was performed by molecular dynamics simulation at 260 K. Pressure influence on hydrate growth was evaluated. The kinetic properties including the growth rates and cage occupancies of the newly formed hydrate and the diffusion coefficient and concentration of N2 molecules in liquid phase were measured. The results showed that the growth of N2 hydrate could be divided into two stages where N2 molecules in gas phase had to dissolve in liquid phase and then form hydrate cages at the liquid-hydrate interface. The diffusion coefficient and concentration of N2 in liquid phase increased linearly with increasing pressure. As the pressure rose from 50 to 100 MPa, the hydrate growth rate kept increasing from 0.11 to 0.62 cages·ns-1·Å-2 and then dropped down to around 0.40 cages·ns-1·Å-2 once the pressure surpassed 100 MPa. During the hydrate formation, the initial sII N2 hydrate phase set in the system served as a template for the subsequent growth of N2 hydrate so that no new crystal structure was found. Analysis on the cage occupancies revealed that the amount of cages occupied by two N2 molecules increased evidently when the pressure was above 100 MPa, which slowed down the growth rate of hydrate cages. Additionally, a small fraction of defective cages including two N2 molecules trapped in 51265 cages and three N2 molecules trapped 51268 cages was observed during the hydrate growth.
Collapse
Affiliation(s)
| | - Xuebing Zhou
- Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences , Guangzhou 510640 , China
| | | | | | | | | | | |
Collapse
|
38
|
Michl J, Sega M, Dellago C. Phase stability of the ice XVII-based CO 2 chiral hydrate from molecular dynamics simulations. J Chem Phys 2019; 151:104502. [PMID: 31521081 DOI: 10.1063/1.5116540] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We computed the phase diagram of CO2 hydrates at high pressure (HP), from 0.3 to 20 kbar, by means of molecular dynamics simulations. The two CO2 hydrates known to occur in this pressure range are the cubic structure I (sI) clathrate and the HP hydrate, whose water framework is the recently discovered ice XVII. We investigated the stability of both hydrates upon heating (melting) as well as the phase changes upon compression. The CO2-filled ice XVII is found to be more stable than the sI clathrate and than the mixture of ice VI and dry ice at pressure values ranging from 6 to 18 kbar and in a wide temperature range, although a phenomenological correction suggests that the stability should more realistically range from 6.5 to 13.5 kbar. Our simulation results support the current hypothesis that the HP hydrate is stable at temperatures above the melting curve of ice VI.
Collapse
Affiliation(s)
- Jakob Michl
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090 Vienna, Austria
| | - Marcello Sega
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090 Vienna, Austria
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090 Vienna, Austria
| |
Collapse
|
39
|
Fernández-Fernández A, Pérez-Rodríguez M, Comesaña A, Piñeiro M. Three-phase equilibrium curve shift for methane hydrate in oceanic conditions calculated from Molecular Dynamics simulations. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2018.10.146] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
40
|
Pérez-Rodríguez M, Otero-Fernández J, Comesaña A, Fernández-Fernández ÁM, Piñeiro MM. Simulation of Capture and Release Processes of Hydrogen by β-Hydroquinone Clathrate. ACS OMEGA 2018; 3:18771-18782. [PMID: 31458440 PMCID: PMC6644111 DOI: 10.1021/acsomega.8b01798] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2018] [Accepted: 12/14/2018] [Indexed: 05/27/2023]
Abstract
Using molecular simulation techniques, we investigate the storage capabilities of H2 gas by the clathrate of hydroquinone (HQ). Quantum mechanics calculations have been used to assess structure and interactions at the atomic scale and molecular dynamics to model the HQ clathrate at successive equilibriums during the processes of capture and release of H2, as well as the diffusion of H2 inside the clathrate structure. The thermodynamic conditions of the simulations performed try to reproduce closely the corresponding experimental procedures, with results that are in good agreement with literature observed trends. The results obtained contribute to depict a more complete and better substantiated image of the mechanisms involved in stability and in the processes of capture and release of H2 by the HQ clathrate.
Collapse
|
41
|
Qiu N, Bai X, Sun N, Yu X, Yang L, Li Y, Yang M, Huang Q, Du S. Grand Canonical Monte Carlo Simulations on Phase Equilibria of Methane, Carbon Dioxide, and Their Mixture Hydrates. J Phys Chem B 2018; 122:9724-9737. [PMID: 30278135 DOI: 10.1021/acs.jpcb.8b04551] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The cage occupancy plays a crucial role in the thermodynamic stability of clathrate hydrates and is an important quantity for understanding the CO2-CH4 replacement phenomenon. In this work, the occupancy isotherms of pure CH4, pure CO2, and their mixture in sI and sII hydrates are studied by GCMC + MD simulations. The adsorption of CH4 and CO2 + CH4 in the sI and sII hydrates can be categorized as the one-site Langmuir type. The calculated occupancy ratio θL/θS and the abundance ratio of CO2 to CH4 vary with the temperature and pressure, which provide the prerequisite information for the prediction of CH4 recovery yield at different conditions in the CO2-CH4 gas exchange process. The phase equilibria of clathrate hydrates of pure gases and mixtures are explored and the corresponding heat of dissociation and hydration numbers are determined. The current investigation provides new perspectives to understand the mechanism behind the gas adsorption behavior of clathrate hydrates.
Collapse
Affiliation(s)
| | | | | | - Xiaohui Yu
- National Laboratory for Condensed Matter Physics, Institute of Physics , Chinese Academy of Sciences , Beijing 100190 , China
| | - Longbin Yang
- College of Power and Energy Engineering , Harbin Engineering University , Harbin , Heilongjiang 150001 , China
| | - Yanjun Li
- College of Power and Energy Engineering , Harbin Engineering University , Harbin , Heilongjiang 150001 , China
| | | | | | | |
Collapse
|
42
|
Arismendi-Arrieta DJ, Valdés Á, Prosmiti R. A Systematic Protocol for Benchmarking Guest-Host Interactions by First-Principles Computations: Capturing CO 2 in Clathrate Hydrates. Chemistry 2018; 24:9353-9363. [PMID: 29600599 DOI: 10.1002/chem.201800497] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2018] [Indexed: 01/19/2023]
Abstract
Clathrate hydrates of CO2 have been proposed as potential molecular materials in tackling important environmental problems related to greenhouse gases capture and storage. Despite the increasing interest in such hydrates and their technological applications, a molecular-level understanding of their formation and properties is still far from complete. Modeling interactions is a challenging and computationally demanding task, essential to reliably determine molecular properties. First-principles calculations for the CO2 guest in all sI, sII, and sH clathrate cages were performed, and the nature of the guest-host interactions, dominated by both hydrogen-bond and van der Waals forces, was systematically investigated. Different families of density functionals, as well as pairwise CO2 @H2 O model potentials versus wavefunction-based quantum approaches were studied for CO2 clathrate-like systems. Benchmark energies for new distance-dependent datasets, consisting of potential energy curves sampling representative configurations of the systems at the repulsive, near-equilibrium, and asymptotic/long-range regions of the full-dimensional surface, were generated, and a general protocol was proposed to assess the accuracy of such conventional and modern approaches at minimum and non-minimum orientations. Our results show that dispersion interactions are important in the guest-host stabilization energies of such clathrate cages, and the encapsulation of the CO2 into guest-free clathrate cages is always energetically favorable. In addition, the orientation of CO2 inside each cage was explored, and the ability of current promising approaches to accurately describe non-covalent CO2 @H2 O guest-host interactions in sI, sII, and sH clathrates was discussed, providing information for their applicability to future multiscale computer simulations.
Collapse
Affiliation(s)
| | - Álvaro Valdés
- Departamento de Física, Universidad Nacional de Colombia, Calle 26, Cra 39, Edificio, 404, Bogotá, Colombia
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006, Madrid, Spain
| |
Collapse
|
43
|
Conde M, Rovere M, Gallo P. Molecular dynamics simulations of freezing-point depression of TIP4P/2005 water in solution with NaCl. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.03.126] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
44
|
Comesaña A, Pérez-Rodríguez M, Fernández-Fernández AM, Piñeiro MM. A description of hydroquinone clathrates using molecular dynamics: Molecular model and crystalline structures for CH4and CO2guests. J Chem Phys 2018; 148:244502. [PMID: 29960310 DOI: 10.1063/1.5027807] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- A. Comesaña
- Universidade de Vigo, Departamento de Física Aplicada, Campus Lagoas-Marcosende, E36310 Vigo, Spain
| | - M. Pérez-Rodríguez
- Universidade de Vigo, Departamento de Física Aplicada, Campus Lagoas-Marcosende, E36310 Vigo, Spain
| | | | - M. M. Piñeiro
- Universidade de Vigo, Departamento de Física Aplicada, Campus Lagoas-Marcosende, E36310 Vigo, Spain
| |
Collapse
|
45
|
Waage MH, Trinh TT, van Erp TS. Diffusion of gas mixtures in the sI hydrate structure. J Chem Phys 2018; 148:214701. [PMID: 29884064 DOI: 10.1063/1.5026385] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Replacing methane with carbon dioxide in gas hydrates has been suggested as a way of harvesting methane, while at the same time storing carbon dioxide. Experimental evidence suggests that this process is facilitated if gas mixtures are used instead of pure carbon dioxide. We studied the free energy barriers for diffusion of methane, carbon dioxide, nitrogen, and hydrogen in the sI hydrate structure using molecular simulation techniques. Cage hops between neighboring cages were considered with and without a water vacancy and with a potential inclusion of an additional gas molecule in either the initial or final cage. Our results give little evidence for enhanced methane and carbon dioxide diffusion if nitrogen is present as well. However, the inclusion of hydrogen seems to have a substantial effect as it diffuses rapidly and can easily enter occupied cages, which reduces the barriers of diffusion for the gas molecules that co-occupy a cage with hydrogen.
Collapse
Affiliation(s)
- Magnus H Waage
- Department of Chemistry, Norwegian University of Science and Technology, Høgskoleringen 5, 7491 Trondheim, Norway
| | - Thuat T Trinh
- Department of Civil and Environmental Engineering, Norwegian University of Science and Technology, Høgskoleringen 7A, 7491 Trondheim, Norway
| | - Titus S van Erp
- Department of Chemistry, Norwegian University of Science and Technology, Høgskoleringen 5, 7491 Trondheim, Norway
| |
Collapse
|
46
|
Tsimpanogiannis IN, Costandy J, Kastanidis P, El Meragawi S, Michalis VK, Papadimitriou NI, Karozis SN, Diamantonis NI, Moultos OA, Romanos GE, Stubos AK, Economou IG. Using clathrate hydrates for gas storage and gas-mixture separations: experimental and computational studies at multiple length scales. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1471224] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Joseph Costandy
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| | - Panagiotis Kastanidis
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Sally El Meragawi
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| | - Vasileios K. Michalis
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Nikolaos I. Papadimitriou
- Environmental Research Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Stylianos N. Karozis
- Environmental Research Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | | | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - George E. Romanos
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Athanassios K. Stubos
- Environmental Research Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | | |
Collapse
|
47
|
Yagasaki T, Matsumoto M, Tanaka H. Adsorption of Kinetic Hydrate Inhibitors on Growing Surfaces: A Molecular Dynamics Study. J Phys Chem B 2018; 122:3396-3406. [PMID: 29278335 DOI: 10.1021/acs.jpcb.7b10356] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We investigate the mechanism of a typical kinetic hydrate inhibitor (KHI), polyvinylcaprolactam (PVCap), which has been applied to prevent hydrate plugs from forming in gas pipe lines, using molecular dynamics simulations of crystal growth of ethylene oxide hydrate. Water-soluble ethylene oxide is chosen as a guest species to avoid problems associated with the presence of the gas phase in the simulation cell such as slow crystal growth. A PVCap dodecamer adsorbs irreversibly on the hydrate surface which grows at supercooling of 3 K when the hydrophobic part of two pendent groups are trapped in open cages at the surface. The amide hydrogen bonds make no contribution to the adsorption. PVCap can adsorb on various crystallographic planes of sI hydrate. This is in contrast to antifreeze proteins, each of which prefers a specific plane of ice. The trapped PVCap gives rise to necessarily the concave surface of the hydrate. The crystal growth rate decreases with increasing surface curvature, indicating that the inhibition by PVCap is explained by the Gibbs-Thomson effect.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Research Institute for Interdisciplinary Science , Okayama University , Okayama 700-8530 , Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science , Okayama University , Okayama 700-8530 , Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science , Okayama University , Okayama 700-8530 , Japan
| |
Collapse
|
48
|
Tsimpanogiannis IN, Economou IG. Monte Carlo simulation studies of clathrate hydrates: A review. J Supercrit Fluids 2018. [DOI: 10.1016/j.supflu.2017.12.017] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
49
|
Conde MM, Rovere M, Gallo P. High precision determination of the melting points of water TIP4P/2005 and water TIP4P/Ice models by the direct coexistence technique. J Chem Phys 2017; 147:244506. [DOI: 10.1063/1.5008478] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- M. M. Conde
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - M. Rovere
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - P. Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| |
Collapse
|
50
|
Different Mechanism Effect between Gas-Solid and Liquid-Solid Interface on the Three-Phase Coexistence Hydrate System Dissociation in Seawater: A Molecular Dynamics Simulation Study. ENERGIES 2017. [DOI: 10.3390/en11010006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|