1
|
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
|
2
|
Li K, Chen B, Yang M, Song Y, Sum AK. Methane hydrate phase equilibrium considering dissolved methane concentrations and interfacial geometries from molecular simulations. J Chem Phys 2023; 159:244505. [PMID: 38153154 DOI: 10.1063/5.0174705] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 12/05/2023] [Indexed: 12/29/2023] Open
Abstract
Natural gas hydrates, mainly existing in permafrost and on the seabed, are expected to be a new energy source with great potential. The exploitation technology of natural gas hydrates is one of the main focuses of hydrate-related studies. In this study, a large-size liquid aqueous solution wrapping a methane hydrate system was established and molecular dynamics simulations were used to investigate the phase equilibrium conditions of methane hydrate at different methane concentrations and interfacial geometries. It is found that the methane concentration of a solution significantly affects the phase equilibrium of methane hydrates. Different methane concentrations at the same temperature and pressure can lead to hydrate formation or decomposition. At the same temperature and pressure, in a system reaching equilibrium, the size of spherical hydrate clusters is coupled to the solution concentration, which is proportional to the Laplace pressure at the solid-liquid interface. Lower solution concentrations reduce the phase equilibrium temperature of methane hydrates at the same pressure; as the concentration increases, the phase equilibrium temperature gradually approaches the actual phase equilibrium temperature. In addition, the interfacial geometry of hydrates affects the thermodynamic stability of hydrates. The spherical hydrate particles have the highest stability for the same volume. Through this study, we provide a stronger foundation to understand the principles driving hydrate formation/dissociation relevant to the exploitation of methane hydrates.
Collapse
Affiliation(s)
- Kehan Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, China
- Phases to Flow Laboratory, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| | - Bingbing Chen
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, China
| | - Mingjun Yang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, China
| | - Yongchen Song
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian, China
| | - Amadeu K Sum
- Phases to Flow Laboratory, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| |
Collapse
|
3
|
Li K, Chen B, Li M, Jiang L, Song Y, Yang M. Facilitation of Hydrate Dissociation and Structural Evolution by Major Marine Anions under Static Electric Fields. J Phys Chem B 2023; 127:10447-10457. [PMID: 37991934 DOI: 10.1021/acs.jpcb.3c06012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2023]
Abstract
Electric fields have been proven to be capable of significantly affecting the equilibrium state of hydrates. In this study, the thermodynamic properties and structural changes of methane hydrate (MH) in various anion solutions in an electric field at 0.7 V/nm were investigated by molecular dynamics simulations. The presence of anions significantly enhances the instability of methane hydrates under electric fields, leading to a staged dissociation process. First, the anions coexist with MH to form a temporary metastable structure under the action of an electric field. Then, the migration of anions causes the dissociation of nearby hydrates and the formation of flow channels in the hydrate layer, which leads to the complete dissociation of MH after a period. The promotive effects of F-, Br-, I-, and Cl- ions were close, while SO42- was relatively weak. The anions are still in hydration shells in the MH phase, but the structure of the hydration shells differs slightly from that in solution (the coordination numbers of I- and SO42- ions increased). The migration resistances of multiple anions to cross the surface of the hydrate layer are similar. However, inside the hydrate phase, the anions with a larger radius have a higher migration resistance. It is difficult for SO42- ions to migrate inside the hydrate phase, and they tend to form a metastable structure on the hydrate surface. Combining our previous studies, SrCl2 solution has the best hydrate promotion under an electric field environment.
Collapse
Affiliation(s)
- Kehan Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Bingbing Chen
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Mingjun Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Lanlan Jiang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Yongchen Song
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| | - Mingjun Yang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian University of Technology, Dalian 116024, China
| |
Collapse
|
4
|
Sun N, Li Y, Qiu N, Francisco JS, Du S. Effects of X (X = H2S, SO2 and N2O) mole fractions on adsorption behavior and phase equilibrium properties of CO2 + X mixed gas hydrate. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121661] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
|
5
|
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
|
6
|
Arjun A, Bolhuis PG. Homogeneous nucleation of crystalline methane hydrate in molecular dynamics transition paths sampled under realistic conditions. J Chem Phys 2023; 158:044504. [PMID: 36725504 DOI: 10.1063/5.0124852] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Methane hydrates are important from a scientific and industrial perspective, and form by nucleation and growth from a supersaturated aqueous solution of methane. Molecular simulation is able to shed light on the process of homogeneous nucleation of hydrates, using straightforward molecular dynamics or rare event enhanced sampling techniques with atomistic and coarse grained force fields. In our previous work [Arjun, T. A. Berendsen, and P. G. Bolhuis, Proc. Natl. Acad. Sci. U. S. A. 116, 19305 (2019)], we performed transition path sampling (TPS) simulations using all atom force fields under moderate driving forces at high pressure, which enabled unbiased atomistic insight into the formation of methane hydrates. The supersaturation in these simulations was influenced by the Laplace pressure induced by the spherical gas reservoir. Here, we investigate the effect of removing this influence. Focusing on the supercooled, supersaturated regime to keep the system size tractable, our TPS simulations indicate that nuclei form amorphous structures below roughly 260 K and crystalline sI structures above 260 K. For these temperatures, the average transition path lengths are significantly longer than in our previous study, pushing the boundaries of what can be achieved with TPS. The temperature to observe a critical nucleus of certain size was roughly 20 K lower compared to a spherical reservoir due to the lower concentration of methane in the solution, yielding a reduced driving force. We analyze the TPS results using a model based on classical nucleation theory. The corresponding free energy barriers are estimated and found to be consistent with previous predictions, thus adding to the overall picture of the hydrate formation process.
Collapse
Affiliation(s)
- A Arjun
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| | - Peter G Bolhuis
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
Yanes-Rodríguez R, Prosmiti R. Assessment of DFT approaches in noble gas clathrate-like clusters: stability and thermodynamics. Phys Chem Chem Phys 2021; 24:1475-1485. [PMID: 34935011 DOI: 10.1039/d1cp04935f] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have assessed the performance and accuracy of different wavefunction-based electronic structure methods, such as DFMP2 and domain-based local pair-natural orbital (DLPNO-CCSD(T)), as well as a variety of density functional theory (DFT) approaches on He@(H2O)N cage systems. We have selected representative clathrate-like structures corresponding to the building blocks present in each of the sI, sII and sH natural gas clathrate hydrates, and we have carefully studied the interaction between a He atom with each of their individual cages. We reported well-converged DFMP2 and DLPNO-CCSD(T) reference data, together with interaction and cohesive energies of four different density functionals (two GGA, revPBE and PW86PBE, and two hybrids, B3LYP and PBE0), including diverse dispersion correction schemes (D3(0), D3(BJ), D4 and XDM) for both He-filled and empty clathrate-like cages. After the analysis of the results, we came to the conclusion that the PW86PBE functional, with both XDM and D4 corrections, and the PBE0-D4 functional present reasonably adequate approaches to describe the guest-host noncovalent interactions that take place in such He clathrate hydrates. Taking into account that the He@sII is the only helium clathrate that scientists have been able to synthesize recently, we have performed a thermodynamic study on the individual 512 and 51264 cages present in the sII crystal. We determined the change in enthalpy, ΔH, and in Gibbs free energy, ΔG, at various temperatures and pressures, and we found out that in the range of experimental conditions the reactions associated with the encapsulation of the He atom inside the cages are exothermic and spontaneous. Finally, we highlighted the importance of an accurate description of the interaction in He@water mixtures, as a crucial component in construction of reliable data-driven models.
Collapse
Affiliation(s)
- Raquel Yanes-Rodríguez
- 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, Spain
| | - Rita Prosmiti
- Institute of Fundamental Physics (IFF-CSIC), CSIC, Serrano 123, 28006 Madrid, Spain.
| |
Collapse
|
10
|
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]
|
11
|
Sun Y, Li K, Su Y, Zhao J. Phase Diagrams for sII Clathrate Hydrates of CO 2 from First-Principles Thermodynamics. J Phys Chem A 2021; 125:5956-5962. [PMID: 34229440 DOI: 10.1021/acs.jpca.1c04673] [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/28/2022]
Abstract
Clathrate hydrates are crystalline solid compounds consisting of a water caged framework and guest molecules such as CH4, C2H6, and CO2. Understanding the phase equilibrium conditions of hydrates is significantly important for the industrial exploitation and experimental synthesis of hydrates. Based on the correct description of the intermolecular noncovalent interactions of clathrate hydrates with vdW-DF2, we studied the crystal structures and the chemical potential phase diagrams of sII hydrates encapsulated with CO2 molecules to provide a deep understanding of the stability mechanism of hydrates. Under the given p-T conditions, the partially occupied hydrates (136H2O·1CO2 and 136H2O·16CO2) and fully occupied hydrates (136H2O·24CO2) are thermodynamically stable, and the equilibrium temperature decreases as the relative CO2 chemical potential increases at the same pressure. We expect that the present study may provide vital information on the stability conditions of CO2 hydrates and trigger new experiments to establish an effective replacement strategy for CO2/CH4.
Collapse
Affiliation(s)
- Yuanze Sun
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
| | - Keyao Li
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
| | - Yan Su
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
| | - Jijun Zhao
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams, Dalian University of Technology, Ministry of Education, Dalian 116024, China
| |
Collapse
|
12
|
Arjun A, Bolhuis PG. Homogenous nucleation rate of CO 2 hydrates using transition interface sampling. J Chem Phys 2021; 154:164507. [PMID: 33940852 DOI: 10.1063/5.0044883] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Carbon dioxide and water can form solid clathrate structures in which water cages encapsulate the gas molecules. Such hydrates have sparked much interest due to their possible application in CO2 sequestration. How the solid structure forms exactly from the liquid phase via a homogenous nucleation process is still poorly understood. This nucleation event is rare on the molecular timescale even under moderate undercooling or supersaturation conditions because of the large free energy barrier toward crystallization, rendering a brute force simulation of hydrate nucleation unfeasible for moderate undercooling or supersaturation. Here, we perform transition interface sampling simulations to quantify the homogenous nucleation rate for CO2 hydrate formation using accurate atomistic force fields at 500 bars for three different temperatures between 260 and 273 K. Collecting more than 100 000 pathways comprising roughly two milliseconds of simulation time, we computed a nucleation rate in the amorphous phase of ∼1021 nuclei s-1 cm-3 for a temperature of 260 K and a rate of ∼1012 nuclei s-1 cm-3 for a temperature of 265 K. For a temperature of 273 K, we find that the hydrate forms an sI crystalline phase with a rate of order of ∼101 nuclei s-1 cm-3. We compare these rates to classical nucleation theory estimates as well as experiments, and to nucleation rate estimates for methane hydrates and discuss possible causes of the observed differences. Our findings shed light on the kinetics of this important clathrate and should assist in future hydrate formation investigation.
Collapse
Affiliation(s)
- A Arjun
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| | - Peter G Bolhuis
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| |
Collapse
|
13
|
Rahbari A, Hens R, Ramdin M, Moultos OA, Dubbeldam D, Vlugt TJH. Recent advances in the continuous fractional component Monte Carlo methodology. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1828585] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- A. Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - R. Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - M. Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - O. A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - D. Dubbeldam
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, Netherlands
| | - T. J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| |
Collapse
|
14
|
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]
|
15
|
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]
|
16
|
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
|
17
|
Veesam SK, Ravipati S, Punnathanam SN. Recent advances in thermodynamics and nucleation of gas hydrates using molecular modeling. Curr Opin Chem Eng 2019. [DOI: 10.1016/j.coche.2019.01.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
18
|
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]
|
19
|
Choudhary N, Chakrabarty S, Roy S, Kumar R. A comparison of different water models for melting point calculation of methane hydrate using molecular dynamics simulations. Chem Phys 2019. [DOI: 10.1016/j.chemphys.2018.08.036] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
20
|
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
|
21
|
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
|
22
|
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
|
23
|
Yuhara D, Brumby PE, Wu DT, Sum AK, Yasuoka K. Analysis of three-phase equilibrium conditions for methane hydrate by isometric-isothermal molecular dynamics simulations. J Chem Phys 2018; 148:184501. [DOI: 10.1063/1.5016609] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Daisuke Yuhara
- Department of Mechanical Engineering, Keio University, 3-4-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - Paul E. Brumby
- Department of Mechanical Engineering, Keio University, 3-4-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| | - David T. Wu
- Chemistry Department, Colorado School of Mines, Golden, Colorado 80401, USA
| | - Amadeu K. Sum
- Hydrates Energy Innovation Lab, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| | - Kenji Yasuoka
- Department of Mechanical Engineering, Keio University, 3-4-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
| |
Collapse
|
24
|
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]
|
25
|
Papadimitriou NI, Tsimpanogiannis IN, Economou IG, Stubos AK. Monte Carlo simulations of the separation of a binary gas mixture (CH4 + CO2) using hydrates. Phys Chem Chem Phys 2018; 20:28026-28038. [DOI: 10.1039/c8cp02050g] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The current study employs Grand Canonical Monte Carlo simulations in order to calculate the process efficiency of separating CH4 + CO2 gas mixtures by utilizing structure sI clathrate hydrates.
Collapse
Affiliation(s)
- Nikolaos I. Papadimitriou
- National Center for Scientific Research “Demokritos”
- Environmental Research Laboratory
- 15310 Aghia Paraskevi Attikis
- Greece
- Public Power Corporation
| | - Ioannis N. Tsimpanogiannis
- National Center for Scientific Research “Demokritos”
- Environmental Research Laboratory
- 15310 Aghia Paraskevi Attikis
- Greece
- National Center for Scientific Research “Demokritos”
| | - Ioannis G. Economou
- Texas A&M University at Qatar
- Chemical Engineering Program
- Education City
- Doha
- Qatar
| | - Athanassios K. Stubos
- National Center for Scientific Research “Demokritos”
- Environmental Research Laboratory
- 15310 Aghia Paraskevi Attikis
- Greece
| |
Collapse
|
26
|
Tanaka H, Yagasaki T, Matsumoto M. On the Thermodynamic Stability of Clathrate Hydrates VI: Complete Phase Diagram. J Phys Chem B 2017; 122:297-308. [DOI: 10.1021/acs.jpcb.7b10581] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hideki Tanaka
- Research Institute for Interdisciplinary
Science, Okayama University, Okayama 700-8530, Japan
| | - 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
| |
Collapse
|
27
|
Jin D, Coasne B. Molecular Simulation of the Phase Diagram of Methane Hydrate: Free Energy Calculations, Direct Coexistence Method, and Hyperparallel Tempering. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:11217-11230. [PMID: 28793774 DOI: 10.1021/acs.langmuir.7b02238] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Different molecular simulation strategies are used to assess the stability of methane hydrate under various temperature and pressure conditions. First, using two water molecular models, free energy calculations consisting of the Einstein molecule approach in combination with semigrand Monte Carlo simulations are used to determine the pressure-temperature phase diagram of methane hydrate. With these calculations, we also estimate the chemical potentials of water and methane and methane occupancy at coexistence. Second, we also consider two other advanced molecular simulation techniques that allow probing the phase diagram of methane hydrate: the direct coexistence method in the Grand Canonical ensemble and the hyperparallel tempering Monte Carlo method. These two direct techniques are found to provide stability conditions that are consistent with the pressure-temperature phase diagram obtained using rigorous free energy calculations. The phase diagram obtained in this work, which is found to be consistent with previous simulation studies, is close to its experimental counterpart provided the TIP4P/Ice model is used to describe the water molecule.
Collapse
Affiliation(s)
- Dongliang Jin
- Laboratoire Interdisciplinaire de Physique (LIPhy), CNRS and Université Grenoble Alpes , F-38000 Grenoble, France
| | - Benoit Coasne
- Laboratoire Interdisciplinaire de Physique (LIPhy), CNRS and Université Grenoble Alpes , F-38000 Grenoble, France
| |
Collapse
|