1
|
Torrejón MJ, Algaba J, Blas FJ. Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation. II. Effect of multiple occupancy. J Chem Phys 2024; 161:054712. [PMID: 39092957 DOI: 10.1063/5.0220098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2024] [Accepted: 07/17/2024] [Indexed: 08/04/2024] Open
Abstract
In this work, we determine the dissociation line of the nitrogen (N2) hydrate by computer simulation using the TIP4P/Ice model for water and the TraPPE force field for N2. This work is the natural extension of Paper I, in which the dissociation temperature of the N2 hydrate has been obtained at 500, 1000, and 1500 bar [Algaba et al., J. Chem. Phys. 159, 224707 (2023)] using the solubility method and assuming single occupancy. We extend our previous study and determine the dissociation temperature of the N2 hydrate at different pressures, from 500 to 4500 bar, taking into account the single and double occupancy of the N2 molecules in the hydrate structure. We calculate the solubility of N2 in the aqueous solution as a function of temperature when it is in contact with a N2-rich liquid phase and when in contact with the hydrate phase with single and double occupancy via planar interfaces. Both curves intersect at a certain temperature that determines the dissociation temperature at a given pressure. We observe a negligible effect of occupancy on the dissociation temperature. Our findings are in very good agreement with the experimental data taken from the literature. We have also obtained the driving force for the nucleation of the hydrate as a function of temperature and occupancy at several pressures. As in the case of the dissociation line, the effect of occupancy on the driving force for nucleation is negligible. To the best of our knowledge, this is the first time that the effect of the occupancy on the driving force for nucleation of a hydrate that exhibits sII crystallographic structure is studied from computer simulation.
Collapse
Affiliation(s)
- Miguel J Torrejó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
| | - 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
|
2
|
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
|
3
|
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
|
4
|
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
|
5
|
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
|
6
|
Walsh MR. Comparing brute force to transition path sampling for gas hydrate nucleation with a flat interface: comments on time reversal symmetry. Phys Chem Chem Phys 2024; 26:5762-5772. [PMID: 38214888 DOI: 10.1039/d3cp05059a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
Fluid to solid nucleation is often investigated with the rare event method transition path sampling (TPS). I claim that the inherent irreversibility of solid nucleation, even at stationary conditions, calls into question TPS's applicability for determining solid nucleation mechanisms, especially for pre-critical behavior. Even when applied to a phenomenon which displays time reversal asymmetry like solid nucleation, TPS is a good means of exploring phase space and giving trends in post-critical structure, and its ability to facilitate nucleation rate and free energy calculations remains outstanding. Forward-only splitting and ratcheting methods such as forward flux sampling are more attractive for understanding nucleation mechanisms as they do not require time reversal symmetry, but at low driving forces may suffer from the same limitations as brute force: they may never make it to the first ratchet. Here I briefly summarize the TPS method and gas hydrate nucleation simulation literature, focusing on topics within both to facilitate a comparison of brute force hydrate nucleation to transition path sampling of hydrate nucleation. Perhaps anecdotally, the brute force technique results in more crystalline trajectories despite having higher driving forces than TPS. I maintain this difference is because of the inherent irreversibility of hydrate nucleation, meaning its pre-critical behavior cannot accurately be determined by the melting trajectories that comprise approximately half of the configurations in TPS's path ensemble.
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] [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
|
8
|
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 PMCID: PMC9623592 DOI: 10.1021/acs.jpcb.2c04867] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2022] [Revised: 09/16/2022] [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
|
9
|
Experimental characterization of CH4 and CO2 hydrates formation in presence of porous quartz and Cu gas-atomized particles: Thermodynamic analyses and evidences about the feasibility of CH4/CO2 reverse replacement. Chem Eng Res Des 2022. [DOI: 10.1016/j.cherd.2022.08.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
10
|
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]
|
11
|
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
|
12
|
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]
|
13
|
Gambelli AM, Stornelli G, Di Schino A, Rossi F. Methane and Carbon Dioxide Hydrate Formation and Dissociation in Presence of a Pure Quartz Porous Framework Impregnated with CuSn12 Metallic Powder: An Experimental Report. MATERIALS 2021; 14:ma14175115. [PMID: 34501204 PMCID: PMC8433810 DOI: 10.3390/ma14175115] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Revised: 08/20/2021] [Accepted: 09/03/2021] [Indexed: 11/16/2022]
Abstract
Hydrate formation and dissociation processes were carried out in the presence of a pure quartz porous medium impregnated with a metallic powder made with a CuSn12 alloy. Experiments were firstly made in the absence of that powder; then, different concentrations were added to the porous medium: 4.23 wt.%, 18.01 wt.%, and 30.66 wt.%. Then, the hydrate dissociation values were compared with those present in the literature. The porous medium was found to act as an inhibitor in the presence of carbon dioxide, while it did not alter methane hydrate, whose formation proceeded similarly to the ideal trend. The addition of CuSn12 promoted the process significantly. In particular, in concentrations of up to 18.01 wt.%, CO2 hydrate formed at milder conditions until it moved below the ideal equilibrium curve. For methane, the addition of 30.66 wt.% of powder significantly reduced the pressure required to form hydrate, but in every case, dissociation values remained below the ideal equilibrium curve.
Collapse
Affiliation(s)
- Alberto Maria Gambelli
- Engineering Department, University of Perugia, Via G. Duranti 93, 06125 Perugia, Italy;
- Correspondence: (A.M.G.); (A.D.S.)
| | - Giulia Stornelli
- Industrial Engineering Department, University of Rome “Tor Vergata”, Via del Politecnico 1, 00133 Rome, Italy;
| | - Andrea Di Schino
- Engineering Department, University of Perugia, Via G. Duranti 93, 06125 Perugia, Italy;
- Correspondence: (A.M.G.); (A.D.S.)
| | - Federico Rossi
- Engineering Department, University of Perugia, Via G. Duranti 93, 06125 Perugia, Italy;
| |
Collapse
|
14
|
Experimental investigation on the possibility of defining the feasibility of CO2/CH4 exchange into a natural gas hydrate marine reservoir via fast analysis of sediment properties. Chem Eng Res Des 2021. [DOI: 10.1016/j.cherd.2021.05.016] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
15
|
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
|