1
|
Tanaka H, Matsumoto M, Yagasaki T, Takeuchi M, Mori Y, Kono T. Stability mechanism of crystalline CO2 and Xe. J Chem Phys 2024; 161:084501. [PMID: 39177089 DOI: 10.1063/5.0223879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2024] [Accepted: 08/08/2024] [Indexed: 08/24/2024] Open
Abstract
We explore the phase behaviors of simple molecular crystals in order to investigate the molecular basis of the stability mechanism relative to their liquid counterparts. The free energies of the face centered cubic crystals of Xe and CO2 are calculated as a collection of oscillators, and those of the liquids are from an equation of state via molecular dynamics simulations. The vibrational free energy in the solid is separated into the harmonic and anharmonic terms. The harmonic free energies decrease harshly with the expansion of the volume manifested as the large positive Grüneisen parameters, but the anharmonic free energies are positive and increase with volume, both of which originate from the deviation of the potential surface from the parabolic curve. The anharmonic free energies, though less significant in magnitude and destabilize the solids thermodynamically, serve to enhance their mechanical stability. The solid-liquid phase boundaries cannot be settled correctly without the exquisite balance between the two opposing contributions. A sharp contrast regarding the solid free energy is found in low-pressure ice, where the harmonic free energy does not decrease monotonically with volume and its anharmonic free energy is negative.
Collapse
Affiliation(s)
- Hideki Tanaka
- Toyota Physical and Chemical Research Institute, Nagakute 480-1192, Japan
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Takuma Yagasaki
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Osaka 560-8531, Japan
| | - Munetaka Takeuchi
- Ochanomizu University, 2-1-1 Ohtsuka, Bunkyo-ku, Tokyo 112-8610, Japan
| | - Yoshihito Mori
- Ochanomizu University, 2-1-1 Ohtsuka, Bunkyo-ku, Tokyo 112-8610, Japan
| | - Takumi Kono
- Engineering Advancement Association of Japan, 1-11-9 Azabudai, Minato-ku, Tokyo 105-0001, Japan
| |
Collapse
|
2
|
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
|
3
|
Liu L, Guan D, Lu Y, Sun M, Liu Y, Zhao J, Yang L. A Molecular Dynamics Study on Xe/Kr Separation Mechanisms Using Crystal Growth Method. ACS OMEGA 2024; 9:25822-25831. [PMID: 38911791 PMCID: PMC11191100 DOI: 10.1021/acsomega.4c00108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 03/31/2024] [Accepted: 05/23/2024] [Indexed: 06/25/2024]
Abstract
The separation of xenon/krypton gas mixtures is a valuable but challenging endeavor in the gas industry due to their similar physical characteristics and closely sized molecules. To address this, we investigated the effectiveness of the hydrate-based gas separation method for mixed Xe-Kr gas via molecular dynamics (MD) simulations. The formation process of hydrates facilitates the encapsulation of guest molecules within hydrate cages, offering a potential strategy for gas separation. Higher temperatures and pressures are advantageous for accelerating the hydrate growth rate. The final occupancy of guest molecules and empty cages within 512, 51264, and all hydrate cages were thoroughly examined. An increase in the pressure and temperature enhanced the occupancy rates of Xe in both 512 and 51264 cages, whereas elevated pressure alone improved the occupancy of Kr in 51264 cages. However, the impact of temperature and pressure on Kr occupancy within 512 cages was found to be minimal. Elevated temperature and pressure resulted in a reduced occupancy of empty cages. Predominantly, 51264 cages were occupied by Xe, whereas Kr showed a propensity to occupy the 512 cages. With increasing simulated pressure, the final occupancy of Xe molecules in all cages rose from 0.37 to 0.41 for simulations at 260 K, while the final occupancy of empty cages decreased from 0.24 to 0.2.
Collapse
Affiliation(s)
- Liangliang Liu
- Shenyang
Aircraft Design Institute Shenyang 110042, China
| | - Dawei Guan
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Yi Lu
- Shenyang
Aircraft Design Institute Shenyang 110042, China
| | - Mingrui Sun
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Yu Liu
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Jiafei Zhao
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| | - Lei Yang
- Key
Laboratory of Ocean Energy Utilization and Energy Conservation of
Ministry of Education, Dalian University
of Technology, Dalian 116024, China
| |
Collapse
|
4
|
Blazquez S, Algaba J, Míguez JM, Vega C, Blas FJ, Conde MM. Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate. J Chem Phys 2024; 160:164721. [PMID: 38686998 DOI: 10.1063/5.0201295] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/01/2024] [Indexed: 05/02/2024] Open
Abstract
Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field.
Collapse
Affiliation(s)
- S Blazquez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - J M Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y del Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| |
Collapse
|
5
|
Algaba J, Blazquez S, Míguez JM, Conde MM, Blas FJ. Three-phase equilibria of hydrates from computer simulation. III. Effect of dispersive interactions in the methane and carbon dioxide hydrates. J Chem Phys 2024; 160:164723. [PMID: 38686999 DOI: 10.1063/5.0201309] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/02/2024] [Indexed: 05/02/2024] Open
Abstract
In this work, the effect of the range of dispersive interactions in determining the three-phase coexistence line of the CO2 and CH4 hydrates has been studied. In particular, the temperature (T3) at which solid hydrate, water, and liquid CO2/gas CH4 coexist has been determined through molecular dynamics simulations using different cutoff values (from 0.9 to 1.6 nm) for dispersive interactions. The T3 of both hydrates has been determined using the direct coexistence simulation technique. Following this method, the three phases in equilibrium are put together in the same simulation box, the pressure is fixed, and simulations are performed at different temperatures T. If the hydrate melts, then T > T3. Conversely, if the hydrate grows, then T < T3. The effect of the cutoff distance on the dissociation temperature has been analyzed at three different pressures for CO2 hydrate: 100, 400, and 1000 bar. Then, we have changed the guest and studied the effect of the cutoff distance on the dissociation temperature of the CH4 hydrate at 400 bar. Moreover, the effect of long-range corrections for dispersive interactions has been analyzed by running simulations with homo- and inhomogeneous corrections and a cutoff value of 0.9 nm. The results obtained in this work highlight that the cutoff distance for the dispersive interactions affects the stability conditions of these hydrates. This effect is enhanced when the pressure is decreased, displacing the T3 about 2-4 K depending on the system and the pressure.
Collapse
Affiliation(s)
- J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J M Míguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y del Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| |
Collapse
|
6
|
Fernández-Fernández ÁM, Bárcena Á, Conde MM, Pérez-Sánchez G, Pérez-Rodríguez M, Piñeiro MM. Modeling oceanic sedimentary methane hydrate growth through molecular dynamics simulation. J Chem Phys 2024; 160:144107. [PMID: 38591679 DOI: 10.1063/5.0203116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Accepted: 03/23/2024] [Indexed: 04/10/2024] Open
Abstract
The crystallization process of methane hydrates in a confined geometry resembling seabed porous silica sedimentary conditions has been studied using molecular dynamics simulations. With this objective in mind, a fully atomistic quartz silica slit pore has been designed, and the temperature stability of a methane hydrate crystalline seed in the presence of water and guest molecule methane has been analyzed. NaCl ion pairs have been added in different concentrations, simulating salinity conditions up to values higher than average oceanic conditions. The structure obtained when the hydrate crystallizes inside the pore is discussed, paying special attention to the presence of ionic doping inside the hydrate and the subsequent induced structural distortion. The shift in the hydrate stability conditions due to the increasing water salinity is discussed and compared with the case of unconfined hydrate, concluding that the influence of the confinement geometry and pore hydrophilicity produces a larger deviation in the confined hydrate phase equilibria.
Collapse
Affiliation(s)
| | - Álvaro Bárcena
- Dpto. de Física Aplicada, Univ. de Vigo, Vigo 36310, Spain
| | - María M Conde
- Dpto. de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, Madrid 28006, Spain
| | | | | | | |
Collapse
|
7
|
Cao P, Wu J, Ning F. Mechanical properties of amorphous CO 2 hydrates: insights from molecular simulations. Phys Chem Chem Phys 2024; 26:9388-9398. [PMID: 38444360 DOI: 10.1039/d4cp00203b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Understanding physicochemical properties of amorphous gas hydrate systems is of great significance to reveal structural stabilities of polycrystalline gas hydrate systems. Furthermore, amorphous gas hydrates can occur ordinarily in the nucleation events of gas hydrate systems. Herein, the mechanical properties of amorphous carbon dioxide hydrates are examined by means of all-atom classical molecular dynamic simulations. Our molecular simulation results reveal that mechanical strengths of amorphous carbon dioxide hydrates are evidently governed by temperatures, confining pressures, and ratios of water to carbon dioxide molecules. Notably, under compressive loads, amorphous carbon dioxide hydrates firstly exhibit monotonic strain hardening, followed by an interesting distinct phenomenon characterized by a steady flow stress at further large deformation strains. Furthermore, structural evolutions of amorphous carbon dioxide hydrates are analyzed on the basis of the N-Hbond DOP order parameter. These important findings can not only contribute to our understanding of the structural stabilities of amorphous gas hydrate systems, but also help to develop fundamental understandings about grain boundaries of gas hydrate systems.
Collapse
Affiliation(s)
- Pinqiang Cao
- School of Resource and Environmental Engineering, Wuhan University of Science and Technology, Wuhan, Hubei 430081, China.
| | - Jianyang Wu
- Department of Physics, Research Institute for Biomimetics and Soft Matter, Fujian Provincial Key Laboratory for Soft Functional Materials Research, Xiamen University, Xiamen 361005, China.
| | - Fulong Ning
- Faculty of Engineering, China University of Geosciences, Wuhan, Hubei 430074, China.
| |
Collapse
|
8
|
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
|
9
|
Hao X, Li C, Meng Q, Sun J, Huang L, Bu Q, Li C. Molecular Dynamics Simulation of the Three-Phase Equilibrium Line of CO 2 Hydrate with OPC Water Model. ACS OMEGA 2023; 8:39847-39854. [PMID: 37901483 PMCID: PMC10601413 DOI: 10.1021/acsomega.3c05673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Accepted: 09/26/2023] [Indexed: 10/31/2023]
Abstract
The three-phase coexistence line of the CO2 hydrate was determined using molecular dynamics (MD) simulations. By using the classical and modified Lorentz-Berthelot (LB) parameters, the simulations were carried out at 10 different pressures from 3 to 500 MPa. For the OPC water model, simulations with the classic and the modified LB parameters both showed negative deviations from the experimental values. For the TIP4P/Ice water model, good agreement with experimental equilibrium data can be achieved when the LB parameter is adjusted based on the solubility of CO2 in water. Our results also show that the influence of the water model on the equilibrium prediction is much larger than the CO2 model. Current simulations indicated that the H2O-H2O and H2O-CO2 cross-interactions' parameters might contribute equally to the accurate prediction of T3. According to our simulations, the prediction of T3 values showed relatively higher accuracy while using the combination of TIP4P/Ice water and EPM2 CO2 with modified LB parameter. Furthermore, varied χ values are recommended for accurate T3 estimation over a wide pressure range. The knowledge obtained in this study will be helpful for further accurate MD simulation of the process of CO2/CH4 replacement.
Collapse
Affiliation(s)
- Xiluo Hao
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Chengfeng Li
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Qingguo Meng
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Jianye Sun
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Li Huang
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Qingtao Bu
- Key
Laboratory of Gas Hydrate, Ministry of Natural Resources, Qingdao Institute of Marine Geology, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| | - Congying Li
- Center
of Deep Sea Research, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
- Laboratory
for Marine Mineral Resources, Laoshan Laboratory, Qingdao 266071, China
| |
Collapse
|
10
|
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
|
11
|
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
|
12
|
Michalis VK, Economou IG, Stubos AK, Tsimpanogiannis IN. Phase equilibria molecular simulations of hydrogen hydrates via the direct phase coexistence approach. J Chem Phys 2022; 157:154501. [PMID: 36272800 DOI: 10.1063/5.0108738] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We report the three-phase (hydrate-liquid water-vapor) equilibrium conditions of the hydrogen-water binary system calculated with molecular dynamics simulations via the direct phase coexistence approach. A significant improvement of ∼10.5 K is obtained in the current study, over earlier simulation attempts, by using a combination of modifications related to the hydrogen model that include (i) hydrogen Lennard-Jones parameters that are a function of temperature and (ii) the water-guest energy interaction parameters optimized further by using the Lorentz-Berthelot combining rules, based on an improved description of the solubility of hydrogen in water.
Collapse
Affiliation(s)
| | - Ioannis G Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Athanasios K Stubos
- Environmental Research Laboratory, National Center for Scientific Research "Demokritos," 15310 Aghia Paraskevi Attikis, Greece
| | - Ioannis N Tsimpanogiannis
- Chemical Process and Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH), 57001 Thermi-Thessaloniki, Greece
| |
Collapse
|
13
|
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]
|
14
|
Mohr S, Pétuya R, Sarria J, Purkayastha N, Bodnar S, Wylde J, Tsimpanogiannis IN. Assessing the effect of a liquid water layer on the adsorption of hydrate anti-agglomerants using molecular simulations. J Chem Phys 2022; 157:094703. [DOI: 10.1063/5.0100260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We have performed Molecular Dynamics simulations to study the adsorption of ten hydrate anti-agglomerants onto a mixed methane-propane sII hydrate surface covered by layers of liquid water of various thickness. As a general trend, we found that the more liquid water is present on the hydrate surface the less favorable the adsorption becomes, even though there are considerable differences between the individual molecules, indicating that the presence and thickness of this liquid water layer is a crucial parameter for anti-agglomerant adsorption studies. Additionally, we found that there exists an optimal thickness of the liquid water layer favoring hydrate growth due to the presence of both liquid water and hydrate-forming guest molecules. For all other cases of liquid water layer thickness, hydrate growth is slower due to the limited availability of hydrate-forming guests close to the hydrate formation front. Finally, we investigated the connection between the thickness of the liquid water layer and the degree of subcooling, and found a very good agreement between our Molecular Dynamics simulations and theoretical predictions.
Collapse
Affiliation(s)
| | | | | | | | - Scot Bodnar
- Clariant Oil Services, United States of America
| | | | | |
Collapse
|
15
|
Hao X, Li C, Liu C, Meng Q, Sun J. The performance of OPC water model in prediction of the phase equilibria of methane hydrate. J Chem Phys 2022; 157:014504. [DOI: 10.1063/5.0093659] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular dynamics (MD) simulations were performed to determine the three-phase coexistence line of sI methane hydrates. The MD simulations were carried out at four different pressures (4, 10, 40 and 100 MPa) by using direct phase coexistence method. In current simulations, water was described by either TIP4P/Ice or OPC models and methane was described as a simple Lennard-Jones (LJ) interaction site. Lorentz-Berthelot combining rules were used to calculate the parameters of the cross interactions. For OPC model, positive deviations from the energetic Lorentz-Berthelot rule were also considered based on the solubility of methane in water. For TIP4P/Ice water model, the obtained three phase coexistence temperatures showed good agreement with experiment data at higher pressures, which is consistent with previous predictions. For OPC water model, simulations using the classic and the modified LB parameters both showed negative deviations to the experimental values. Our results also indicated that the deviation of the T3 prediction by OPC model not much correlated with the predicted melting point of ice. At 4 MPa, the modified OPC model showed outstanding prediction of hydrate equilibrium temperature, even better than the prediction by TIP4P/Ice. The relative higher accuracy in biomolecular MD of OPC model suggests that this model may have a better performance in hydrate MD simulations of biomolecule-based additives.
Collapse
Affiliation(s)
- Xiluo Hao
- Qingdao Institute of Marine Geology, China
| | | | | | | | - Jianye Sun
- Qingdao Institute of Marine Geology, China
| |
Collapse
|
16
|
Yu Z, Li J, Zhang X. A new hypothesis for cavitation nucleation in gas saturated solutions: clustering of gas molecules lowers significantly the surface tension. Chin J Chem Eng 2022. [DOI: 10.1016/j.cjche.2022.06.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
17
|
Molecular insights into the heterogeneous crystal growth of tetrahydrofuran hydrate: Kinetic and interfacial properties. J Mol Graph Model 2022; 115:108205. [DOI: 10.1016/j.jmgm.2022.108205] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2022] [Revised: 04/17/2022] [Accepted: 04/23/2022] [Indexed: 11/22/2022]
|
18
|
Mileo PGM, Rogge SMJ, Houlleberghs M, Breynaert E, Martens JA, Van Speybroeck V. Interfacial study of clathrates confined in reversed silica pores. JOURNAL OF MATERIALS CHEMISTRY. A 2021; 9:21835-21844. [PMID: 34707871 PMCID: PMC8491980 DOI: 10.1039/d1ta03105h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Accepted: 08/22/2021] [Indexed: 05/08/2023]
Abstract
Storing methane in clathrates is one of the most promising alternatives for transporting natural gas (NG) as it offers similar gas densities to liquefied and compressed NG while offering lower safety risks. However, the practical use of clathrates is limited given the extremely low temperatures and high pressures necessary to form these structures. Therefore, it has been suggested to confine clathrates in nanoporous materials, as this can facilitate clathrate's formation conditions while preserving its CH4 volumetric storage. Yet, the choice of nanoporous materials to be employed as the clathrate growing platform is still rather arbitrary. Herein, we tackle this challenge in a systematic way by computationally exploring the stability of clathrates confined in alkyl-grafted silica materials with different pore sizes, ligand densities and ligand types. Based on our findings, we are able to propose key design criteria for nanoporous materials favoring the stability of a neighbouring clathrate phase, namely large pore sizes, high ligand densities, and smooth pore walls. We hope that the atomistic insight provided in this work will guide and facilitate the development of new nanomaterials designed to promote the formation of clathrates.
Collapse
Affiliation(s)
- Paulo G M Mileo
- Center for Molecular Modeling (CMM), Ghent University Technologiepark 46 B-9052 Zwijnaarde Belgium
| | - Sven M J Rogge
- Center for Molecular Modeling (CMM), Ghent University Technologiepark 46 B-9052 Zwijnaarde Belgium
| | - Maarten Houlleberghs
- Center for Surface Chemistry and Catalysis, Katholieke Universiteit Leuven Celestijnenlaan 200F 3001 Heverlee Belgium
| | - Eric Breynaert
- Center for Surface Chemistry and Catalysis, Katholieke Universiteit Leuven Celestijnenlaan 200F 3001 Heverlee Belgium
| | - Johan A Martens
- Center for Surface Chemistry and Catalysis, Katholieke Universiteit Leuven Celestijnenlaan 200F 3001 Heverlee Belgium
| | - Veronique Van Speybroeck
- Center for Molecular Modeling (CMM), Ghent University Technologiepark 46 B-9052 Zwijnaarde Belgium
| |
Collapse
|
19
|
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]
|
20
|
Abstract
Water and methane can stay together under low temperature and high pressure in the forms of liquid solutions and crystalline solids. From liquid and gaseous states to crystalline solids or the contrary processes, amorphous methane hydrates can occur in these evolution scenarios. Herein, mechanical properties of amorphous methane hydrates are explored for the first time to bridge the gap between mechanical responses of monocrystalline and polycrystalline methane hydrates. Our results demonstrate that mechanical properties of amorphous methane hydrates are strongly governed by our original proposed order parameter, namely, normalized hydrogen-bond directional order parameter. Followed by this important achievement, a multistep deformation mechanism core is proposed to explain mechanical properties of amorphous methane hydrates. Through an extensive detailed analysis of amorphous methane hydrates, our simulation results not only greatly enlarge our fundamental understanding for mechanical responses of amorphous methane hydrates in geological systems but also offer a fresh perspective in structure-property topics of solid materials in future science and technology.
Collapse
Affiliation(s)
- Pinqiang Cao
- School of Resource and Environmental Engineering, Wuhan University of Science and Technology, Wuhan, Hubei 430081, China
| |
Collapse
|
21
|
Tsimpanogiannis IN. A novel hybrid method for the calculation of methane hydrate-water interfacial tension along the three-phase (hydrate-liquid water-vapor) equilibrium line. J Chem Phys 2021; 155:024702. [PMID: 34266278 DOI: 10.1063/5.0051383] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
We use a novel hybrid method to explore the temperature dependence of the solid-liquid interfacial tension of a system that consists of solid methane hydrate and liquid water. The calculated values along the three-phase (hydrate-liquid water-vapor) equilibrium line are obtained through the combination of available experimental measurements and computational results that are based on approaches at the atomistic scale, including molecular dynamics and Monte Carlo. An extensive comparison with available experimental and computational studies is performed, and a critical assessment and re-evaluation of previously reported data is presented.
Collapse
Affiliation(s)
- Ioannis N Tsimpanogiannis
- Chemical Process & Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH), 57001 Thermi-Thessaloniki, Greece
| |
Collapse
|
22
|
Guo Q, Hu W, Zhang Y, Zhang K, Dong B, Qin Y, Li W. Molecular dynamics simulation of the interfacial properties of methane-water and methane-brine systems. MOLECULAR SIMULATION 2021. [DOI: 10.1080/08927022.2021.1929969] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Qiuyi Guo
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian, People’s Republic of China
| | - Wenfeng Hu
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian, People’s Republic of China
| | - Yue Zhang
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian, People’s Republic of China
| | - Kun Zhang
- School of Ocean and Civil Engineering, Dalian Ocean University, Dalian, People’s Republic of China
| | - Bo Dong
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian, People’s Republic of China
- School of Energy and Power Engineering, Dalian University of Technology, Dalian, People’s Republic of China
| | - Yan Qin
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian, People’s Republic of China
- School of Energy and Power Engineering, Dalian University of Technology, Dalian, People’s Republic of China
| | - Weizhong Li
- Key Laboratory of Ocean Energy Utilization and Energy Conservation of Ministry of Education, Dalian, People’s Republic of China
- School of Energy and Power Engineering, Dalian University of Technology, Dalian, People’s Republic of China
| |
Collapse
|
23
|
Cruz FJAL, Mota JPB. Structure and thermodynamics of empty clathrate hydrates below the freezing point of water. Phys Chem Chem Phys 2021; 23:16033-16043. [PMID: 34286770 DOI: 10.1039/d1cp00893e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recently prepared as a new H2O phase, ice XVI was obtained by degassing a Ne-sII clathrate hydrate under vacuum, however very little is known of that crystalline solid under temperatures (T ≤ 220 K) and pressures (p ≤ 5000 bar) relevant for the Earth's environment and geochemistry. In this work, atomically detailed calculations using long time-scale molecular simulations, seldom paralelled before, are employed to probe empty sII clathrate hydrates. It is found that the volumetric response to an applied pressure-temperature gradient is accurately described by the Parsafar and Mason equation of state with an accuracy of at least 99.7%. Structural deformation induced upon the crystals is interpreted by monitoring the unit cell length and isobaric thermal expansivity, whilst benchmarked against previous neutron diffraction measurements of ice XVI and hexagonal ice under room pressure conditions; a critical comparison is established with other sII guest occupied lattices (CH4, CO2 and CnH2n+2 with n = 2, 3, 4), often found in permafrost regions and in the margins of continental shelves. Such an analysis reveals that empty sII frameworks are slightly more stable to thermal deformation than their sI analogues and that hexagonal ice is the structurally most stable of the condensed H2O phases addressed here. Of paramount importance for the oil and natural gas industries, heat capacities obtained from enthalpy profiles are identical for the sI and sII empty clathrates up to 2000 bar and diverge by only ∼7.3% at 5000 bar. The canonical tetrahedral symmetry of water-bonded networks is analysed in terms of an angular and a distance order parameters, which are observed to decrease (increase) as pressure (temperature) increases (decreases). The results now being reported constitute a landmark for future studies dealing with high-pressure and very low-temperature conditions, characteristic of the Earth's permafrost environment and other planetary interiors, whilst contributing to expand our knowledge regarding the recently discovered ice XVI phase.
Collapse
Affiliation(s)
- Fernando J A L Cruz
- LAQV-REQUIMTE, Department of Chemistry, Faculdade de Ciências e Tecnologia, NOVA University Lisbon, 2829-516 Caparica, Portugal.
| | | |
Collapse
|
24
|
Mohr S, Pétuya R, Wylde J, Sarria J, Purkayastha N, Ward Z, Bodnar S, Tsimpanogiannis IN. Size dependence of the dissociation process of spherical hydrate particles via microsecond molecular dynamics simulations. Phys Chem Chem Phys 2021; 23:11180-11185. [DOI: 10.1039/d1cp01223a] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The dissociation process of spherical sII mixed methane–propane hydrate particles in liquid hydrocarbon was investigated via microsecond-long molecular dynamics simulations.
Collapse
Affiliation(s)
- Stephan Mohr
- Nextmol (Bytelab Solutions SL)
- Barcelona
- Spain
- Barcelona Supercomputing Center (BSC)
- Barcelona
| | - Rémi Pétuya
- Nextmol (Bytelab Solutions SL)
- Barcelona
- Spain
| | - Jonathan Wylde
- Clariant Oil Services, Clariant Corporation
- Houston
- USA
- Heriot Watt University
- Edinburgh
| | - Juan Sarria
- Clariant Produkte (Deutschland) GmbH
- Frankfurt
- Germany
| | | | - Zachary Ward
- Clariant Oil Services, Clariant Corporation
- Houston
- USA
| | - Scot Bodnar
- Clariant Oil Services, Clariant Corporation
- Houston
- USA
| | - Ioannis N. Tsimpanogiannis
- Chemical Process & Energy Resources Institute (CPERI)
- Centre for Research & Technology Hellas (CERTH)
- Thermi-Thessaloniki
- Greece
| |
Collapse
|
25
|
Tsimpanogiannis IN, Michalis VK, Economou IG. Novel methodology for the calculation of the enthalpy of enclathration of methane hydrates using molecular dynamics simulations. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1711976] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Chemical Process & Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH) Thermi-Thessaloniki, Greece
- Environmental Research Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Vasileios K. Michalis
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
| | - Ioannis G. Economou
- Institute of Nanoscience and Nanotechnology, Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research ‘Demokritos’, Aghia Paraskevi Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| |
Collapse
|
26
|
Yu KB, Yazaydin AO. Does Confinement Enable Methane Hydrate Growth at Low Pressures? Insights from Molecular Dynamics Simulations. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2020; 124:11015-11022. [PMID: 32582402 PMCID: PMC7304911 DOI: 10.1021/acs.jpcc.0c02246] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2020] [Revised: 04/30/2020] [Indexed: 05/23/2023]
Abstract
Natural methane hydrates are estimated to be the largest source of unexploited hydrocarbon fuel. The ideal conditions for methane hydrate formation are low temperatures and high pressures. On the other hand, recent experimental studies suggest that porous materials, thanks to their confinement effects, can enable methane hydrate formation at milder conditions, although there has not been a consensus on this. A number of studies have investigated methane hydrate growth in confinement by employing molecular simulations; however, these were carried out at either very high pressures or very low temperatures. Therefore, the effects of confinement on methane hydrate growth at milder conditions have not yet been elaborated by molecular simulations. In order to address this, we carried out a systematic study by performing molecular dynamics (MD) simulations of methane water systems. Using a direct phase coexistence approach, microsecond-scale MD simulations in the isobaric-isothermal (NPT) ensemble were performed in order to study the behavior of methane hydrates in the bulk and in confined nanospaces of hydroxylated silica pores at external pressures ranging from 1 to 100 bar and a simulation temperature corresponding to a 2 °C experimental temperature. We validated the combination of the TIP4P/ice water and TraPPE-UA methane models in order to correctly predict the behavior of methane hydrates in accordance to their phase equilibria. We also demonstrated that the dispersion corrections applied to short-range interactions lead to artificially induced hydrate growth. We observed that in the confinement of a hydroxylated silica pore, a convex-shaped methane nanobuble forms, and methane hydrate growth primarily takes place in the center of the pore rather than the surfaces where a thin water layer exists. Most importantly, our study showed that in the nanopores methane hydrate growth can indeed take place at pressures which would be too low for the growth of methane hydrates in the bulk.
Collapse
|
27
|
Tsimpanogiannis IN, Lichtner PC. Methane dissolution inside bulk or porous-medium-confined water at near-hydrate equilibrium conditions. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1699186] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Chemical Process & Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH), Thermi-Thessaloniki, Greece
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
| | - Peter C. Lichtner
- OFM Research Southwest, Santa Fe, NM, USA
- Center for Water & the Environment, University of New Mexico, Albuquerque, NM, USA
| |
Collapse
|
28
|
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]
|
29
|
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]
|
30
|
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
|
31
|
Sicard F, Bui T, Monteiro D, Lan Q, Ceglio M, Burress C, Striolo A. Emergent Properties of Antiagglomerant Films Control Methane Transport: Implications for Hydrate Management. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:9701-9710. [PMID: 30058809 DOI: 10.1021/acs.langmuir.8b01366] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
The relationship between collective properties and performance of antiagglomerants (AAs) used in hydrate management is handled using molecular dynamics simulations and enhanced sampling techniques. A thin film of AAs adsorbed at the interface between one flat sII methane hydrate substrate and a fluid hydrocarbon mixture containing methane and n-dodecane is studied. The AA considered is a surface-active compound with a complex hydrophilic head that contains both amide and tertiary ammonium cation groups and hydrophobic tails. At a sufficiently high AA density, the interplay between the surfactant layer and the liquid hydrocarbon excludes methane from the interfacial region. In this scenario, we combine metadynamics and umbrella sampling frameworks to study accurately the free-energy landscape and the equilibrium rates associated with the transport of one methane molecule across the AA film. We observe that the local configurational changes of the liquid hydrocarbon packed within the AA film are associated with high free-energy barriers for methane transport. The time scales estimated for the transport of methane across the AA film can be, in some cases, comparable to those reported in the literature for the growth of hydrates, suggesting that one possible mechanism by which AAs delay the formation of hydrate plugs could be providing a barrier to methane transport. Considering the interplay between the structural design and collective properties of AAs might be of relevance to improve their performance in flow assurance.
Collapse
Affiliation(s)
- François Sicard
- Department of Chemical Engineering , University College London , WC1E 7JE London , U.K
| | - Tai Bui
- Department of Chemical Engineering , University College London , WC1E 7JE London , U.K
| | | | - Qiang Lan
- Halliburton , Houston , 77032 Texas , United States
| | - Mark Ceglio
- Halliburton , Houston , 77032 Texas , United States
| | | | - Alberto Striolo
- Department of Chemical Engineering , University College London , WC1E 7JE London , U.K
| |
Collapse
|
32
|
Bui T, Sicard F, Monteiro D, Lan Q, Ceglio M, Burress C, Striolo A. Antiagglomerants Affect Gas Hydrate Growth. J Phys Chem Lett 2018; 9:3491-3496. [PMID: 29870264 DOI: 10.1021/acs.jpclett.8b01180] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
In gas clathrate hydrates, inclusion gas molecules stabilize crystalline water structures. In addition to being fundamentally interesting, gas hydrates attract significant practical attention because of their possible application in various high-tech technologies. However, gas hydrates pose health, safety, and environmental risks when they form within oil and gas pipelines, as well as within hydrocarbon-producing and treatment facilities. Among available strategies to control and sometimes prevent hydrate plug formation is the use of surface-active low-molecular-weight compounds, known as antiagglomerants (AAs). AAs prevent the agglomeration of small hydrate particles into large plugs. It is not clear whether AAs promote or frustrate hydrate growth. We present two molecular mechanisms by which AAs promote and frustrate, respectively, hydrate growth. Our results could lead to innovative methodologies for managing hydrates in high-tech applications, as well as for securing the safety of oil and gas operations.
Collapse
Affiliation(s)
- Tai Bui
- Department of Chemical Engineering , University College London , WC1 E7JE London , United Kingdom
| | - Francois Sicard
- Department of Chemical Engineering , University College London , WC1 E7JE London , United Kingdom
| | | | - Qiang Lan
- Halliburton , Houston , Texas 77032 , United States
| | - Mark Ceglio
- Halliburton , Houston , Texas 77032 , United States
| | | | - Alberto Striolo
- Department of Chemical Engineering , University College London , WC1 E7JE London , United Kingdom
| |
Collapse
|
33
|
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
|
34
|
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
|
35
|
Yagasaki T, Matsumoto M, Tanaka H. Adsorption of Kinetic Hydrate Inhibitors on Growing Surfaces: A Molecular Dynamics Study. J Phys Chem B 2018; 122:3396-3406. [PMID: 29278335 DOI: 10.1021/acs.jpcb.7b10356] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We investigate the mechanism of a typical kinetic hydrate inhibitor (KHI), polyvinylcaprolactam (PVCap), which has been applied to prevent hydrate plugs from forming in gas pipe lines, using molecular dynamics simulations of crystal growth of ethylene oxide hydrate. Water-soluble ethylene oxide is chosen as a guest species to avoid problems associated with the presence of the gas phase in the simulation cell such as slow crystal growth. A PVCap dodecamer adsorbs irreversibly on the hydrate surface which grows at supercooling of 3 K when the hydrophobic part of two pendent groups are trapped in open cages at the surface. The amide hydrogen bonds make no contribution to the adsorption. PVCap can adsorb on various crystallographic planes of sI hydrate. This is in contrast to antifreeze proteins, each of which prefers a specific plane of ice. The trapped PVCap gives rise to necessarily the concave surface of the hydrate. The crystal growth rate decreases with increasing surface curvature, indicating that the inhibition by PVCap is explained by the Gibbs-Thomson effect.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Research Institute for Interdisciplinary Science , Okayama University , Okayama 700-8530 , Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science , Okayama University , Okayama 700-8530 , Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science , Okayama University , Okayama 700-8530 , Japan
| |
Collapse
|
36
|
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]
|
37
|
Luis DP, Romero-Ramirez IE, González-Calderón A, López-Lemus J. The coexistence temperature of hydrogen clathrates: A molecular dynamics study. J Chem Phys 2018; 148:114503. [DOI: 10.1063/1.5017854] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- D. P. Luis
- CONACYT Research Fellow-Centro de Ingeniería y Desarrollo Industrial (Sede Campeche), Ave. Playa Pie de la Cuesta No. 702, Desarrollo San Pablo, Querétaro, Querétaro 76125, Mexico
| | - I. E. Romero-Ramirez
- CONACYT Research Fellow-Centro de Ingeniería y Desarrollo Industrial (Sede Campeche), Ave. Playa Pie de la Cuesta No. 702, Desarrollo San Pablo, Querétaro, Querétaro 76125, Mexico
| | - A. González-Calderón
- CONACYT Research Fellow-Centro de Ingeniería y Desarrollo Industrial (Sede Campeche), Ave. Playa Pie de la Cuesta No. 702, Desarrollo San Pablo, Querétaro, Querétaro 76125, Mexico
| | - J. López-Lemus
- Facultad de Ciencias, Universidad Autónoma del Estado de México, Toluca CP 50295, Mexico
| |
Collapse
|
38
|
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
|
39
|
Different Mechanism Effect between Gas-Solid and Liquid-Solid Interface on the Three-Phase Coexistence Hydrate System Dissociation in Seawater: A Molecular Dynamics Simulation Study. ENERGIES 2017. [DOI: 10.3390/en11010006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
40
|
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
|
41
|
13 The Role of Molecular Thermodynamics in Developing Industrial Processes and Novel Products That Meet the Needs for a Sustainable Future. ACTA ACUST UNITED AC 2017. [DOI: 10.1201/9781315153209-14] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
42
|
Waage MH, Vlugt TJH, Kjelstrup S. Phase Diagram of Methane and Carbon Dioxide Hydrates Computed by Monte Carlo Simulations. J Phys Chem B 2017; 121:7336-7350. [DOI: 10.1021/acs.jpcb.7b03071] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | | |
Collapse
|
43
|
Papadimitriou NI, Tsimpanogiannis IN, Economou IG, Stubos AK. Storage of H2 in Clathrate Hydrates: Evaluation of Different Force-Fields used in Monte Carlo Simulations. Mol Phys 2017. [DOI: 10.1080/00268976.2016.1277590] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Nikolaos I. Papadimitriou
- National Center for Scientific Research ‘Demokritos’, Environmental Research Laboratory, Aghia Paraskevi Attikis, Greece
| | - Ioannis N. Tsimpanogiannis
- National Center for Scientific Research ‘Demokritos’, Environmental Research Laboratory, Aghia Paraskevi Attikis, Greece
- Texas A&M University at Qatar, Chemical Engineering Program, Doha, Qatar
| | | | - Athanassios K. Stubos
- National Center for Scientific Research ‘Demokritos’, Environmental Research Laboratory, Aghia Paraskevi Attikis, Greece
| |
Collapse
|
44
|
Zhang Z, Guo GJ. The effects of ice on methane hydrate nucleation: a microcanonical molecular dynamics study. Phys Chem Chem Phys 2017; 19:19496-19505. [DOI: 10.1039/c7cp03649c] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The NVE simulations realize the ice shrinking when methane hydrate nucleates both heterogeneously and homogeneously.
Collapse
Affiliation(s)
- Zhengcai Zhang
- Key Laboratory of Earth and Planetary Physics
- Institute of Geology and Geophysics
- Chinese Academy of Sciences
- Beijing 100029
- China
| | - Guang-Jun Guo
- Key Laboratory of Earth and Planetary Physics
- Institute of Geology and Geophysics
- Chinese Academy of Sciences
- Beijing 100029
- China
| |
Collapse
|
45
|
Costandy J, Michalis VK, Tsimpanogiannis IN, Stubos AK, Economou IG. Molecular dynamics simulations of pure methane and carbon dioxide hydrates: lattice constants and derivative properties. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1241442] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Joseph Costandy
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| | | | - Ioannis N. Tsimpanogiannis
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
- Environmental Research Laboratory, National Center for Scientific Research (NCSR) “Demokritos”, Athens, Greece
| | - Athanassios K. Stubos
- Environmental Research Laboratory, National Center for Scientific Research (NCSR) “Demokritos”, Athens, Greece
| | | |
Collapse
|
46
|
Espinosa JR, Young JM, Jiang H, Gupta D, Vega C, Sanz E, Debenedetti PG, Panagiotopoulos AZ. On the calculation of solubilities via direct coexistence simulations: Investigation of NaCl aqueous solutions and Lennard-Jones binary mixtures. J Chem Phys 2016; 145:154111. [DOI: 10.1063/1.4964725] [Citation(s) in RCA: 61] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- J. R. Espinosa
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. M. Young
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - H. Jiang
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - D. Gupta
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - C. Vega
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E. Sanz
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - P. G. Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - A. Z. Panagiotopoulos
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
47
|
Affiliation(s)
- Alberto Striolo
- Department of Chemical Engineering, University College London, London, UK
| |
Collapse
|
48
|
Benavides AL, Aragones JL, Vega C. Consensus on the solubility of NaCl in water from computer simulations using the chemical potential route. J Chem Phys 2016; 144:124504. [PMID: 27036458 DOI: 10.1063/1.4943780] [Citation(s) in RCA: 83] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The solubility of NaCl in water is evaluated by using three force field models: Joung-Cheatham for NaCl dissolved in two different water models (SPC/E and TIP4P/2005) and Smith Dang NaCl model in SPC/E water. The methodology based on free-energy calculations [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] and [J. L. Aragones et al., J. Chem. Phys. 136, 244508 (2012)] has been used, except, that all calculations for the NaCl in solution were obtained by using molecular dynamics simulations with the GROMACS package instead of homemade MC programs. We have explored new lower molalities and made longer runs to improve the accuracy of the calculations. Exploring the low molality region allowed us to obtain an analytical expression for the chemical potential of the ions in solution as a function of molality valid for a wider range of molalities, including the infinite dilute case. These new results are in better agreement with recent estimations of the solubility obtained with other methodologies. Besides, two empirical simple rules have been obtained to have a rough estimate of the solubility of a certain model, by analyzing the ionic pairs formation as a function of molality and/or by calculating the difference between the NaCl solid chemical potential and the standard chemical potential of the salt in solution.
Collapse
Affiliation(s)
- A L Benavides
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J L Aragones
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
49
|
Papadimitriou NI, Tsimpanogiannis IN, Economou IG, Stubos AK. The effect of lattice constant on the storage capacity of hydrogen hydrates: a Monte Carlo study. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1202456] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Nikolaos I. Papadimitriou
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos”, Athens, Greece
| | - Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos”, Athens, Greece
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| | | | - Athanassios K. Stubos
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos”, Athens, Greece
| |
Collapse
|
50
|
Zhang Z, Liu CJ, Walsh MR, Guo GJ. Effects of ensembles on methane hydrate nucleation kinetics. Phys Chem Chem Phys 2016; 18:15602-8. [PMID: 27222203 DOI: 10.1039/c6cp02171a] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
By performing molecular dynamics simulations to form a hydrate with a methane nano-bubble in liquid water at 250 K and 50 MPa, we report how different ensembles, such as the NPT, NVT, and NVE ensembles, affect the nucleation kinetics of the methane hydrate. The nucleation trajectories are monitored using the face-saturated incomplete cage analysis (FSICA) and the mutually coordinated guest (MCG) order parameter (OP). The nucleation rate and the critical nucleus are obtained using the mean first-passage time (MFPT) method based on the FS cages and the MCG-1 OPs, respectively. The fitting results of MFPT show that hydrate nucleation and growth are coupled together, consistent with the cage adsorption hypothesis which emphasizes that the cage adsorption of methane is a mechanism for both hydrate nucleation and growth. For the three different ensembles, the hydrate nucleation rate is quantitatively ordered as follows: NPT > NVT > NVE, while the sequence of hydrate crystallinity is exactly reversed. However, the largest size of the critical nucleus appears in the NVT ensemble, rather than in the NVE ensemble. These results are helpful for choosing a suitable ensemble when to study hydrate formation via computer simulations, and emphasize the importance of the order degree of the critical nucleus.
Collapse
Affiliation(s)
- Zhengcai Zhang
- Key Laboratory of Earth and Planetary Physics, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China.
| | | | | | | |
Collapse
|