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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
3
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
4
|
Debenedetti PG, Kim YY, Meldrum FC, Tanaka H. Special Topic Preface: Nucleation-Current understanding approaching 150 years after Gibbs. J Chem Phys 2024; 160:100401. [PMID: 38465675 DOI: 10.1063/5.0203119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Accepted: 02/14/2024] [Indexed: 03/12/2024] Open
Affiliation(s)
- Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Yi-Yeoun Kim
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Fiona C Meldrum
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Hajime Tanaka
- Department of Fundamental Engineering, The University of Tokyo, Tokyo 153-8505, Japan
- Research Center for Advanced Science and Technology, The University of Tokyo, Tokyo 153-8904, Japan
| |
Collapse
|
5
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
6
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
7
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|