1
|
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
|
2
|
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
|
3
|
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
|
4
|
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
|
5
|
Sharma M, Singh S. Carbon dioxide sequestration in natural gas hydrates - effect of flue and noble gases. Phys Chem Chem Phys 2023; 25:30211-30222. [PMID: 37830431 DOI: 10.1039/d3cp03777k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2023]
Abstract
Clean energy is one of the immediate requirements all over the world to tackle the global energy demands. Natural gas hydrates (NGHs) are one of the proposed alternatives that could be used to extract methane as clean energy and simultaneously sequestrate carbon dioxide. However, the formation of CH4-CO2 mixed hydrates and the first hydrate layer besides the interface reduces the rate of CO2 sequestration and methane extraction in NGHs, and thus, multistep extraction of methane is one of the proposed solutions. We report the atomic level factors that could enhance CO2 sequestration in the newly formed first hydrate layer besides the interface in the presence of flue and noble gases using DFT calculations and molecular dynamics simulations at 250 K and 0.15 kbar. The simulations show the formation of stable dual cages (large-large or small-large) that lead to the formation of a four-caged, Y-shaped cluster (growth synthon) which leads to the formation of a hydrate unit cell in heterogeneous medium. Among the flue and noble gases, only argon forms energetically favorable dual cages with itself and CO2 due to which enhanced CO2 sequestration is observed at different concentrations of Ar and CO2 where the CO2 : Ar (2.5 : 1.5) system shows the best CO2 sequestration in the first layer besides the interface. The results also provide understanding into the previously reported concentration dependent CO2 selectivity in sI hydrates in the presence of third gases (N2 and H2S).
Collapse
Affiliation(s)
- Manju Sharma
- School of Chemistry, University of Hyderabad, Hyderabad, 500 046, India.
| | - Satyam Singh
- School of Chemistry, University of Hyderabad, Hyderabad, 500 046, India.
| |
Collapse
|
6
|
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
|
7
|
Pineda M, Phan A, Koh CA, Striolo A, Stamatakis M. Stochastic Cellular Automata Modeling of CO 2 Hydrate Growth and Morphology. CRYSTAL GROWTH & DESIGN 2023; 23:4222-4239. [PMID: 37304394 PMCID: PMC10251419 DOI: 10.1021/acs.cgd.3c00045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/12/2023] [Revised: 05/04/2023] [Indexed: 06/13/2023]
Abstract
Carbon dioxide (CO2) hydrates are important in a diverse range of applications and technologies in the environmental and energy fields. The development of such technologies relies on fundamental understanding, which necessitates not only experimental but also computational studies of the growth behavior of CO2 hydrates and the factors affecting their crystal morphology. As experimental observations show that the morphology of CO2 hydrate particles differs depending on growth conditions, a detailed understanding of the relation between the hydrate structure and growth conditions would be helpful. To this end, this work adopts a modeling approach based on hybrid probabilistic cellular automata to investigate variations in CO2 hydrate crystal morphology during hydrate growth from stagnant liquid water presaturated with CO2. The model, which uses free energy density profiles as inputs, correlates the variations in growth morphology to the system subcooling ΔT, i.e., the temperature deficiency from the triple CO2-hydrate-water equilibrium temperature under a given pressure, and properties of the growing hydrate-water interface, such as surface tension and curvature. The model predicts that when ΔT is large, parabolic needle-like or dendrite crystals emerge from planar fronts that deform and lose stability. In agreement with chemical diffusion-limited growth, the position of such planar fronts versus time follows a power law. In contrast, the tips of the emerging parabolic crystals steadily grow in proportion to time. The modeling framework is computationally fast and produces complex growth morphology phenomena under diffusion-controlled growth from simple, easy-to-implement rules, opening the way for employing it in multiscale modeling of gas hydrates.
Collapse
Affiliation(s)
- Miguel Pineda
- Thomas
Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London WC1E 7JE, United Kingdom
- Institute
for Materials Discovery, University College
London, WC1H 0AJ, London, United Kingdom
| | - Anh Phan
- Department
of Chemical and Process Engineering, University
of Surrey, Guildford, GU2 7XH, United Kingdom
| | - Carolyn Ann Koh
- Center
for Hydrate Research, Chemical & Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, United States
| | - Alberto Striolo
- Thomas
Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London WC1E 7JE, United Kingdom
- School
of Chemical, Biological and Materials Engineering, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Michail Stamatakis
- Thomas
Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London WC1E 7JE, United Kingdom
| |
Collapse
|
8
|
Li J, Wang Z, Liang Z. Analysis of Influencing Factors and Kinetic Characteristics of Spherical Methane Hydrate Decomposition. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2023; 39:7122-7131. [PMID: 37167340 DOI: 10.1021/acs.langmuir.3c00421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Herein, several molecular systems are simulated by molecular dynamics to study the decomposition process and fluctuation-dissipation characteristics of spherical methane hydrates under different conditions. The spherical radius and the movement of the hydrate-liquid water interface during decomposition are measured. Different fitted formulas of the variation of methane numbers are obtained from the decomposition of spherical and bulk methane hydrates. Fluctuation-dissipation characteristics for spherical methane hydrates with different radii are analyzed, which show that increasing the scale of hydrates can increase the relaxation time and slow down the fluctuation process. The variations of the hydrogen bond and hydrogen-bond lifetime are calculated. For hydrate phase water, the peak of the hydrogen-bond lifetime lies between 8 and 10 ps. After complete decomposition, the hydrogen-bond lifetime mainly distributes in 0 and 2 ps and the peak disappears. The effects of temperature, cage occupancy, liquid phase environment, and spherical hydrate scale are explored. The decomposition activation energy for the spherical hydrate with a radius of 20 Å is calculated to be 52.23 kJ/mol. It can speed up the decomposition rate as well as the diffusion of methane and water molecules with a lower cage occupancy. For the effect of the liquid phase environment, it is found that the number of liquid water rarely affects the decomposition. However, when the Na+ and Cl- concentrations change from 0 to 10%, the decomposition time reduces from ∼511 to ∼369 ps, which indicates that there is an obviously positive impact on decomposition.
Collapse
Affiliation(s)
- Jia Li
- Energy and Power Department, Shandong University of Science and Technology, Qingdao 266590, Shandong, China
| | - Zhaoliang Wang
- Energy and Power Department, China University of Petroleum, Qingdao 266580, Shandong, China
| | - Zhenju Liang
- Energy and Power Department, China University of Petroleum, Qingdao 266580, Shandong, China
| |
Collapse
|
9
|
Li T, Liu N, Huang J. Effects of carbon nanotube on methane hydrate formation by molecular dynamics simulation. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120621] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
10
|
|
11
|
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
|
12
|
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]
|
13
|
A review of clathrate hydrate nucleation, growth and decomposition studied using molecular dynamics simulation. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118025] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
14
|
Study the effect of Ag nanoparticles on the kinetics of CO2 hydrate growth by molecular dynamics simulation. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117668] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
15
|
Wang PW, Wu DT, Lin ST. Promotion mechanism for the growth of CO 2 hydrate with urea using molecular dynamics simulations. Chem Commun (Camb) 2021; 57:5330-5333. [PMID: 33928959 DOI: 10.1039/d0cc06165d] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The role of urea as a kinetic promoter for the growth of CO2 hydrates is revealed for the first time using molecular dynamics simulations. Analysis of simulation trajectories shows that urea plays two important roles in the growth process: increasing mass transport of CO2 and catalyzing cage formation at the solid-liquid interface.
Collapse
Affiliation(s)
- Po-Wei Wang
- Department of Chemical Engineering, National Taiwan University, Taiwan.
| | - David T Wu
- Department of Chemistry and Department of Chemical and Biological Engineering, Colorado School of Mines, USA
| | - Shiang-Tai Lin
- Department of Chemical Engineering, National Taiwan University, Taiwan.
| |
Collapse
|
16
|
Liu N, Liu T. Different pathways for methane hydrate nucleation and crystallization at high supercooling: Insights from molecular dynamics simulations. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115466] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
17
|
Grabowska J, Kuffel A, Zielkiewicz J. Revealing the Frank-Evans "Iceberg" Structures within the Solvation Layer around Hydrophobic Solutes. J Phys Chem B 2021; 125:1611-1617. [PMID: 33539702 PMCID: PMC7898264 DOI: 10.1021/acs.jpcb.0c09489] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Using computer simulations, the structural properties of solvation water of three model hydrophobic molecules, methane and two fullerenes (C60 and C80), were studied. Systems were simulated at temperatures in the range of 250-298 K. By analyzing both the local ordering of the molecules of water in the solvation layers and the structure of hydrogen bond network, it is shown that in the solvation layer of hydrophobic molecules, ordered aggregates consisting of water molecules are formed. Even though it is difficult to define the exact structure of these aggregates, their existence alone is clearly noticeable. Moreover, these aggregates become more pronounced with the decrease of temperature. The existence of the ordered aggregates around the hydrophobic solutes complies with the concept of "icebergs" proposed by Frank and Evans.
Collapse
Affiliation(s)
- Joanna Grabowska
- Faculty of Chemistry, Department of Physical Chemistry, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Anna Kuffel
- Faculty of Chemistry, Department of Physical Chemistry, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Jan Zielkiewicz
- Faculty of Chemistry, Department of Physical Chemistry, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland
| |
Collapse
|
18
|
Chen C, Hu W, Yang L, Zhao J, Song Y. Gas supersaturation and diffusion joint controlled CH4 nanobubble evolution during hydrate dissociation. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114614] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
19
|
Liu T, Liu N, Chen L. Microscopic Insights into the Formation of Methane Hydrate in the Absence/Presence of Electrolyte Ions. ChemistrySelect 2020. [DOI: 10.1002/slct.201904693] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Tingsong Liu
- University of Shanghai for Science and Technology School of Energy and Power EngineeringUniversity of Shanghai for Science and Technology Shanghai China
| | - Ni Liu
- University of Shanghai for Science and Technology School of Energy and Power EngineeringUniversity of Shanghai for Science and Technology Shanghai China
| | - Litao Chen
- University of Shanghai for Science and Technology School of Energy and Power EngineeringUniversity of Shanghai for Science and Technology Shanghai China
| |
Collapse
|
20
|
Castillo-Borja F, Bravo-Sánchez UI, Vázquez-Román R, Díaz-Ovalle CO. Biogas purification via sII hydrates in the presence of THF and DMSO solutions using MD simulations. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2019.111904] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
21
|
Mirzaeifard S, Servio P, Rey AD. Characterization of nucleation of methane hydrate crystals: Interfacial theory and molecular simulation. J Colloid Interface Sci 2019; 557:556-567. [DOI: 10.1016/j.jcis.2019.09.056] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Revised: 09/06/2019] [Accepted: 09/17/2019] [Indexed: 01/18/2023]
|
22
|
Yi L, Zhou X, He Y, Cai Z, Zhao L, Zhang W, Shao Y. Molecular Dynamics Simulation Study on the Growth of Structure II Nitrogen Hydrate. J Phys Chem B 2019; 123:9180-9186. [PMID: 31609605 DOI: 10.1021/acs.jpcb.9b06386] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Crystal growth of N2 hydrate in a three-phase system consisting of N2 hydrate, liquid water, and gaseous N2 was performed by molecular dynamics simulation at 260 K. Pressure influence on hydrate growth was evaluated. The kinetic properties including the growth rates and cage occupancies of the newly formed hydrate and the diffusion coefficient and concentration of N2 molecules in liquid phase were measured. The results showed that the growth of N2 hydrate could be divided into two stages where N2 molecules in gas phase had to dissolve in liquid phase and then form hydrate cages at the liquid-hydrate interface. The diffusion coefficient and concentration of N2 in liquid phase increased linearly with increasing pressure. As the pressure rose from 50 to 100 MPa, the hydrate growth rate kept increasing from 0.11 to 0.62 cages·ns-1·Å-2 and then dropped down to around 0.40 cages·ns-1·Å-2 once the pressure surpassed 100 MPa. During the hydrate formation, the initial sII N2 hydrate phase set in the system served as a template for the subsequent growth of N2 hydrate so that no new crystal structure was found. Analysis on the cage occupancies revealed that the amount of cages occupied by two N2 molecules increased evidently when the pressure was above 100 MPa, which slowed down the growth rate of hydrate cages. Additionally, a small fraction of defective cages including two N2 molecules trapped in 51265 cages and three N2 molecules trapped 51268 cages was observed during the hydrate growth.
Collapse
Affiliation(s)
| | - Xuebing Zhou
- Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences , Guangzhou 510640 , China
| | | | | | | | | | | |
Collapse
|
23
|
Nakate P, Ghosh B, Das S, Roy S, Kumar R. Molecular dynamics study on growth of carbon dioxide and methane hydrate from a seed crystal. Chin J Chem Eng 2019. [DOI: 10.1016/j.cjche.2019.02.006] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
24
|
Parui S, Jana B. Factors Promoting the Formation of Clathrate-Like Ordering of Water in Biomolecular Structure at Ambient Temperature and Pressure. J Phys Chem B 2019; 123:811-824. [PMID: 30605607 DOI: 10.1021/acs.jpcb.8b11172] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Clathrate hydrate forms when a hydrophobic molecule is entrapped inside a water cage or cavity. Although biomolecular structures also have hydrophobic patches, clathrate-like water is found in only a limited number of biomolecules. Also, while clathrate hydrates form at low temperature and moderately higher pressure, clathrate-like water is observed in biomolecular structure at ambient temperature and pressure. These indicate presence of other factors along with hydrophobic environment behind the formation of clathrate-like water in biomolecules. In the current study, we presented a systematic approach to explore the factors behind the formation of clathrate-like water in biomolecules by means of molecular dynamics simulation of a model protein, maxi, which is a naturally occurring nanopore and has clathrate-like water inside the pore. Removal of either confinement or hydrophobic environment results in the disappearance of clathrate-like water ordering, indicating a coupled role of these two factors. Apart from these two factors, clathrate-like water ordering also requires anchoring groups that can stabilize the clathrate-like water through hydrogen bonding. Our results uncover crucial factors for the stabilization of clathrate-like ordering in biomolecular structure which can be used for the development of new biomolecular structure promoting clathrate formation.
Collapse
Affiliation(s)
- Sridip Parui
- School of Chemical Sciences , Indian Association for the Cultivation of Science , Jadavpur, Kolkata 700032 , India
| | - Biman Jana
- School of Chemical Sciences , Indian Association for the Cultivation of Science , Jadavpur, Kolkata 700032 , India
| |
Collapse
|
25
|
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]
|
26
|
Molecular dynamics simulation of methane hydrate formation in presence and absence of amino acid inhibitors. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.08.108] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
27
|
Parui S, Jana B. Molecular Insights into the Unusual Structure of an Antifreeze Protein with a Hydrated Core. J Phys Chem B 2018; 122:9827-9839. [PMID: 30286600 DOI: 10.1021/acs.jpcb.8b05350] [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/19/2022]
Abstract
The primary driving force for protein folding is the formation of a well-packed, anhydrous core. However, recently, the crystal structure of an antifreeze protein, maxi, has been resolved where the core of the protein is filled with water, which apparently contradicts the existing notion of protein folding. Here, we have performed standard molecular dynamics (MD) simulation, replica exchange MD (REMD) simulation, and umbrella sampling using TIP4P water at various temperatures (300, 260, and 240 K) to explore the origin of this unusual structural feature. It is evident from standard MD and REMD simulations that the protein is found to be stable at 240 K in its unusual state. The core of protein has two layers of semi-clathrate water separating the methyl groups of alanine residues from different helical strands. However, with increasing temperature (260 and 300 K), the stability decreases as the core becomes dehydrated, and methyl groups of alanine are tightly packed driven by hydrophobic interactions. Calculation of the potential of mean force by an umbrella sampling technique between a pair of model hydrophobes resembling maxi protein at 240 K shows the stabilization of second solvent-separated minima (SSM), which provides a thermodynamic rationale of the unusual structural feature in terms of weakening of the hydrophobic interaction. Because the stabilization of SSMs is implicated for cold denaturation, it suggests that the maxi protein is so designed by nature where the cold denatured-like state becomes the biologically active form as it works near or below the freezing point of water.
Collapse
Affiliation(s)
- Sridip Parui
- Department of Physical Chemistry , Indian Association for the Cultivation of Science , Jadavpur, Kolkata 700032 , India
| | - Biman Jana
- Department of Physical Chemistry , Indian Association for the Cultivation of Science , Jadavpur, Kolkata 700032 , India
| |
Collapse
|
28
|
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
|
29
|
Warrier P, Khan MN, Srivastava V, Maupin CM, Koh CA. Overview: Nucleation of clathrate hydrates. J Chem Phys 2018; 145:211705. [PMID: 28799342 DOI: 10.1063/1.4968590] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular level knowledge of nucleation and growth of clathrate hydrates is of importance for advancing fundamental understanding on the nature of water and hydrophobic hydrate formers, and their interactions that result in the formation of ice-like solids at temperatures higher than the ice-point. The stochastic nature and the inability to probe the small length and time scales associated with the nucleation process make it very difficult to experimentally determine the molecular level changes that lead to the nucleation event. Conversely, for this reason, there have been increasing efforts to obtain this information using molecular simulations. Accurate knowledge of how and when hydrate structures nucleate will be tremendously beneficial for the development of sustainable hydrate management strategies in oil and gas flowlines, as well as for their application in energy storage and recovery, gas separation, carbon sequestration, seawater desalination, and refrigeration. This article reviews various aspects of hydrate nucleation. First, properties of supercooled water and ice nucleation are reviewed briefly due to their apparent similarity to hydrates. Hydrate nucleation is then reviewed starting from macroscopic observations as obtained from experiments in laboratories and operations in industries, followed by various hydrate nucleation hypotheses and hydrate nucleation driving force calculations based on the classical nucleation theory. Finally, molecular simulations on hydrate nucleation are discussed in detail followed by potential future research directions.
Collapse
Affiliation(s)
- Pramod Warrier
- Center for Hydrate Research, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| | - M Naveed Khan
- Center for Hydrate Research, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| | - Vishal Srivastava
- Center for Hydrate Research, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| | - C Mark Maupin
- Center for Hydrate Research, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| | - Carolyn A Koh
- Center for Hydrate Research, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| |
Collapse
|
30
|
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
|
31
|
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
|
32
|
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
|
33
|
Xu J, Li L, Liu J, Wang X, Yan Y, Zhang J. The molecular mechanism of the inhibition effects of PVCaps on the growth of sI hydrate: an unstable adsorption mechanism. Phys Chem Chem Phys 2018. [DOI: 10.1039/c8cp00010g] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Adsorption and non-binding-hydrate sites produce synergistic effects that lead to unstable adsorption and inhibition effects (recurring destruction of hydrate).
Collapse
Affiliation(s)
- Jiafang Xu
- National Engineering Laboratory for Subsea Equipment Testing and Detection Technology
- School of Petroleum Engineering
- China University of Petroleum (East China)
- Qingdao
- People's Republic of China
| | - Liwen Li
- National Engineering Laboratory for Subsea Equipment Testing and Detection Technology
- School of Petroleum Engineering
- China University of Petroleum (East China)
- Qingdao
- People's Republic of China
| | - Jinxiang Liu
- School of Physics and Technology
- University of Jinan
- Jinan
- People's Republic of China
| | - Xiaopu Wang
- National Engineering Laboratory for Subsea Equipment Testing and Detection Technology
- School of Petroleum Engineering
- China University of Petroleum (East China)
- Qingdao
- People's Republic of China
| | - Youguo Yan
- College of Science
- China University of Petroleum (East China)
- Qingdao
- People's Republic of China
| | - Jun Zhang
- College of Science
- China University of Petroleum (East China)
- Qingdao
- People's Republic of China
| |
Collapse
|
34
|
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]
|
35
|
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
|
36
|
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]
|
37
|
Liu Y, Zhao L, Deng S, Bai D. Evolution of bubbles in decomposition and replacement process of methane hydrate. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1359745] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Yinan Liu
- Key Laboratory of Efficient Utilization of Low and Medium Grade Energy, Tianjin University, Ministry of Education of China, Tianjin, China
| | - Li Zhao
- Key Laboratory of Efficient Utilization of Low and Medium Grade Energy, Tianjin University, Ministry of Education of China, Tianjin, China
| | - Shuai Deng
- Key Laboratory of Efficient Utilization of Low and Medium Grade Energy, Tianjin University, Ministry of Education of China, Tianjin, China
| | - Dongsheng Bai
- Department of Chemistry, School of Science, Beijing Technology and Business University, Beijing, China
| |
Collapse
|
38
|
Sun N, Li Z, Qiu N, Yu X, Zhang X, Li Y, Yang L, Luo K, Huang Q, Du S. Ab Initio Studies on the Clathrate Hydrates of Some Nitrogen- and Sulfur-Containing Gases. J Phys Chem A 2017; 121:2620-2626. [PMID: 28304172 DOI: 10.1021/acs.jpca.6b11850] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Ab initio calculations are performed to investigate the host-guest interactions and multiple occupancies of some sulfur- (H2S, CS2) and nitrogen-containing (N2, NO, and NH3) molecules in dodecahedral, tetrakaidecahedral, and hexakaidecahedral water cages in this work. Five functionals in the framework of density functional theory are compared, and the M06-2X method appears to be the best to predict the binding energies as well as the geometries. Results show that N2 and NO molecules are more stable in the 51264 cage, while NH3 and H2S prefer to stabilize in the 51262 cage. This suggests that the sI hydrates of NH3 and H2S exhibit higher stability than the sII structures and that sII NO hydrate is more stable than sI NO hydrate. N2 is found to be more stable in type II structure with single occupancy and to form type I hydrate with multiple occupancy, which is consistent with the experimental observations. As to the guest molecule CS2, it may undergo severe structural deformation in the 512 and 51262 cage. For multiple occupancies, the 512, 51262, and 51264 water cages can trap up to two N2 molecules, and the 51264 water cage can accommodate two H2S molecules. This work is expected to provide new insight into the formation mechanism of clathrate hydrates for atmospherically important molecules.
Collapse
Affiliation(s)
- Ningru Sun
- School of Chemistry and Materials Science, Heilongjiang University , Harbin, Heilongjiang 150080, P. R. China
| | - Zewen Li
- School of Chemistry and Materials Science, Heilongjiang University , Harbin, Heilongjiang 150080, P. R. China
| | - Nianxiang Qiu
- Engineering Laboratory of Specialty Fibers and Nuclear Energy Materials, Ningbo Institute of Industrial Technology, Chinese Academy of Sciences , Ningbo, Zhejiang 315201, P. R. China
| | - Xiaohui Yu
- National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Science , Beijing 100190, P. R. China
| | - Xuran Zhang
- Engineering Laboratory of Specialty Fibers and Nuclear Energy Materials, Ningbo Institute of Industrial Technology, Chinese Academy of Sciences , Ningbo, Zhejiang 315201, P. R. China
| | - Yanjun Li
- College of Power and Energy Engineering, Harbin Engineering University , Harbin, Heilongjiang 150001, P. R. China
| | - Longbin Yang
- College of Power and Energy Engineering, Harbin Engineering University , Harbin, Heilongjiang 150001, P. R. China
| | - Kan Luo
- Engineering Laboratory of Specialty Fibers and Nuclear Energy Materials, Ningbo Institute of Industrial Technology, Chinese Academy of Sciences , Ningbo, Zhejiang 315201, P. R. China
| | - Qing Huang
- Engineering Laboratory of Specialty Fibers and Nuclear Energy Materials, Ningbo Institute of Industrial Technology, Chinese Academy of Sciences , Ningbo, Zhejiang 315201, P. R. China
| | - Shiyu Du
- Engineering Laboratory of Specialty Fibers and Nuclear Energy Materials, Ningbo Institute of Industrial Technology, Chinese Academy of Sciences , Ningbo, Zhejiang 315201, P. R. China
| |
Collapse
|
39
|
Yan KF, Li XS, Chen ZY, Xia ZM, Xu CG, Zhang Z. Molecular Dynamics Simulation of the Crystal Nucleation and Growth Behavior of Methane Hydrate in the Presence of the Surface and Nanopores of Porous Sediment. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2016; 32:7975-7984. [PMID: 27398713 DOI: 10.1021/acs.langmuir.6b01601] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The behavior of hydrate formation in porous sediment has been widely studied because of its importance in the investigation of reservoirs and in the drilling of natural gas hydrate. However, it is difficult to understand the hydrate nucleation and growth mechanism on the surface and in the nanopores of porous media by experimental and numerical simulation methods. In this work, molecular dynamics simulations of the nucleation and growth of CH4 hydrate in the presence of the surface and nanopores of clay are carried out. The molecular configurations and microstructure properties are analyzed for systems containing one H2O hydrate layer (System A), three H2O hydrate layers (System B), and six H2O hydrate layers (System C) in both clay and the bulk solution. It is found that hydrate formation is more complex in porous media than in the pure bulk solution and that there is cooperativity between hydrate growth and molecular diffusion in clay nanopores. The hydroxylated edge sites of the clay surface could serve as a source of CH4 molecules to facilitate hydrate nucleation. The diffusion velocity of molecules is influenced by the growth of the hydrate that forms a block in the throats of the clay nanopore. Comparing hydrate growth in different clay pore sizes reveals that the pore size plays an important role in hydrate growth and molecular diffusion in clay. This simulation study provides the microscopic mechanism of hydrate nucleation and growth in porous media, which can be favorable for the investigation of the formation of natural gas hydrate in sediments.
Collapse
Affiliation(s)
- Ke-Feng Yan
- Key Laboratory of Gas Hydrate, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences , Guangzhou 510640, China
- Guangdong Key Laboratory of New and Renewable Energy Research and Development , Guangzhou 510640, China
| | - Xiao-Sen Li
- Key Laboratory of Gas Hydrate, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences , Guangzhou 510640, China
- Guangdong Key Laboratory of New and Renewable Energy Research and Development , Guangzhou 510640, China
| | - Zhao-Yang Chen
- Key Laboratory of Gas Hydrate, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences , Guangzhou 510640, China
- Guangdong Key Laboratory of New and Renewable Energy Research and Development , Guangzhou 510640, China
| | - Zhi-Ming Xia
- Key Laboratory of Gas Hydrate, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences , Guangzhou 510640, China
- Guangdong Key Laboratory of New and Renewable Energy Research and Development , Guangzhou 510640, China
| | - Chun-Gang Xu
- Key Laboratory of Gas Hydrate, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences , Guangzhou 510640, China
- Guangdong Key Laboratory of New and Renewable Energy Research and Development , Guangzhou 510640, China
| | - Zhiqiang Zhang
- College of Mining Engineering, Taiyuan University of Technology , Taiyuan 030024, China
| |
Collapse
|
40
|
A Theoretical Study of the Hydration of Methane, from the Aqueous Solution to the sI Hydrate-Liquid Water-Gas Coexistence. Int J Mol Sci 2016; 17:ijms17060378. [PMID: 27240339 PMCID: PMC4926321 DOI: 10.3390/ijms17060378] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2015] [Revised: 02/24/2016] [Accepted: 02/26/2016] [Indexed: 11/30/2022] Open
Abstract
Monte Carlo and molecular dynamics simulations were done with three recent water models TIP4P/2005 (Transferable Intermolecular Potential with 4 Points/2005), TIP4P/Ice (Transferable Intermolecular Potential with 4 Points/ Ice) and TIP4Q (Transferable Intermolecular Potential with 4 charges) combined with two models for methane: an all-atom one OPLS-AA (Optimal Parametrization for the Liquid State) and a united-atom one (UA); a correction for the C–O interaction was applied to the latter and used in a third set of simulations. The models were validated by comparison to experimental values of the free energy of hydration at 280, 300, 330 and 370 K, all under a pressure of 1 bar, and to the experimental radial distribution functions at 277, 283 and 291 K, under a pressure of 145 bar. Regardless of the combination rules used for σC,O, good agreement was found, except when the correction to the UA model was applied. Thus, further simulations of the sI hydrate were performed with the united-atom model to compare the thermal expansivity to the experiment. A final set of simulations was done with the UA methane model and the three water models, to study the sI hydrate-liquid water-gas coexistence at 80, 230 and 400 bar. The melting temperatures were compared to the experimental values. The results show the need to perform simulations with various different models to attain a reliable and robust molecular image of the systems of interest.
Collapse
|
41
|
Luis D, López-Lemus J, Maspoch ML, Franco-Urquiza E, Saint-Martin H. Methane hydrate: shifting the coexistence temperature to higher temperatures with an external electric field. MOLECULAR SIMULATION 2016. [DOI: 10.1080/08927022.2016.1139704] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- D.P. Luis
- CONACYT Research Fellow-Centro de Ingeniería y Desarrollo Industrial, Querétaro, México
| | - J. López-Lemus
- Facultad de ciencias, Universidad Autónoma del Estado de México, Toluca, México
| | - M. Ll. Maspoch
- Centre Català del Plàstic, Universitat Politècnica de Catalunya, Terrassa, Spain
| | - E.A. Franco-Urquiza
- CONACYT Research Fellow-Centro de Ingeniería y Desarrollo Industrial, Querétaro, México
| | - H. Saint-Martin
- Instituto de Ciencias Físicas, Universidad Nacional Autónoma de México,
| |
Collapse
|
42
|
Wu JY, Chen LJ, Chen YP, Lin ST. Molecular dynamics study on the nucleation of methane + tetrahydrofuran mixed guest hydrate. Phys Chem Chem Phys 2016; 18:9935-47. [PMID: 26750660 DOI: 10.1039/c5cp06419h] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The nucleation of methane (CH4), tetrahydrofuran (THF), and CH4 + THF hydrates are investigated by microsecond MD simulations. These three systems exhibit distinct structural developments in the aqueous phase quantified by the formation of cage structures of hydrogen bonded water molecules. The development of a cluster of cages in the CH4 system is limited by the scarce CH4 molecules in the solution, while in the THF system it is limited by the short lifetime of cages. In the CH4 + THF mixed guest system, a small cluster of caged CH4 molecules can be rapidly stabilized by abundant neighboring cages of THF molecules. Therefore, the induction time of the CH4 + THF mixed guest system is found to be significantly shorter than that of the pure CH4 and pure THF systems. Furthermore, the structure of cages found in the initially formed cage clusters are often different from the typical 5(12)6(n) (n = 0, 2, 3, 4) cages observed in clathrate hydrate systems. The cluster of cages may grow or transform into structure I or II clathrate hydrate in the later stages.
Collapse
Affiliation(s)
- Jyun-Yi Wu
- Department of Chemical Engineering, National Taiwan University, Taipei, 10617, Taiwan.
| | | | | | | |
Collapse
|
43
|
Yagasaki T, Matsumoto M, Tanaka H. Effects of thermodynamic inhibitors on the dissociation of methane hydrate: a molecular dynamics study. Phys Chem Chem Phys 2015; 17:32347-57. [PMID: 26587576 DOI: 10.1039/c5cp03008k] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate the effects of methanol and NaCl, which are known as thermodynamic hydrate inhibitors, on the dissociation kinetics of methane hydrate in aqueous solutions by using molecular dynamics simulations. It is shown that the dissociation rate is not constant but changes with time. The dissociation rate in the initial stage is increased by methanol whereas it is decreased by NaCl. This difference arises from the opposite effects of the two thermodynamic inhibitors on the hydration free energy of methane. The dissociation rate of methane hydrate is increased by the formation of methane bubbles in the aqueous phase because the bubbles absorb surrounding methane molecules. It is found that both methanol and NaCl facilitate the bubble formation. However, their mechanisms are completely different from each other. The presence of ions enhances the hydrophobic interactions between methane molecules. In addition, the ions in the solution cause a highly non-uniform distribution of dissolved methane molecules. These two effects result in the easy formation of bubbles in the NaCl solution. In contrast, methanol assists the bubble formation because of its amphiphilic character.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Department of Chemistry, Faculty of Science, Okayama University, Okayama, 700-8530, Japan
| | - Masakazu Matsumoto
- Department of Chemistry, Faculty of Science, Okayama University, Okayama, 700-8530, Japan
| | - Hideki Tanaka
- Department of Chemistry, Faculty of Science, Okayama University, Okayama, 700-8530, Japan and Research Center of New Functional Materials for Energy Production, Storage and Transport, Okayama, 700-8530, Japan.
| |
Collapse
|
44
|
Luis DP, Herrera-Hernández EC, Saint-Martin H. A theoretical study of the dissociation of the sI methane hydrate induced by an external electric field. J Chem Phys 2015; 143:204503. [PMID: 26627964 DOI: 10.1063/1.4936214] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular dynamics simulations in the equilibrium isobaric-isothermal (NPT) ensemble were used to examine the strength of an external electric field required to dissociate the methane hydrate sI structure. The water molecules were modeled using the four-site TIP4P/Ice analytical potential and methane was described as a simple Lennard-Jones interaction site. A series of simulations were performed at T = 260 K with P = 80 bars and at T = 285 K with P = 400 bars with an applied electric field ranging from 1.0 V nm(-1) to 5.0 V nm(-1). For both (T,P) conditions, applying a field greater than 1.5 V nm(-1) resulted in the orientation of the water molecules such that an ice Ih-type structure was formed, from which the methane was segregated. When the simulations were continued without the external field, the ice-like structures became disordered, resulting in two separate phases: gas methane and liquid water.
Collapse
Affiliation(s)
- D P Luis
- CONACYT Research Fellow-Centro de Ingeniería y Desarrollo Industrial, Av. Playa Pie de la Cuesta No. 702, Desarrollo San Pablo, Querétaro, Querétaro 76125, Mexico
| | - E C Herrera-Hernández
- CONACYT Research Fellow-Centro de Ingeniería y Desarrollo Industrial, Av. Playa Pie de la Cuesta No. 702, Desarrollo San Pablo, Querétaro, Querétaro 76125, Mexico
| | - H Saint-Martin
- Instituto de Ciencias Físicas, Universidad Autónoma de México, Apartado Postal 48-3, Cuernavaca, Morelos 62251, Mexico
| |
Collapse
|
45
|
Yagasaki T, Matsumoto M, Tanaka H. Adsorption Mechanism of Inhibitor and Guest Molecules on the Surface of Gas Hydrates. J Am Chem Soc 2015; 137:12079-85. [DOI: 10.1021/jacs.5b07417] [Citation(s) in RCA: 104] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Takuma Yagasaki
- Department
of Chemistry, Faculty of Science, Okayama University, Okayama, 700-8530, Japan
| | - Masakazu Matsumoto
- Department
of Chemistry, Faculty of Science, Okayama University, Okayama, 700-8530, Japan
| | - Hideki Tanaka
- Department
of Chemistry, Faculty of Science, Okayama University, Okayama, 700-8530, Japan
- Research Center of New Functional Materials for Energy Production, Storage and Transport, Okayama, 700-8530, Japan
| |
Collapse
|
46
|
Bagherzadeh SA, Alavi S, Ripmeester J, Englezos P. Formation of methane nano-bubbles during hydrate decomposition and their effect on hydrate growth. J Chem Phys 2015; 142:214701. [DOI: 10.1063/1.4920971] [Citation(s) in RCA: 62] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
47
|
Jimenez-Angeles F, Firoozabadi A. Enhanced Hydrate Nucleation Near the Limit of Stability. J Phys Chem A 2015. [DOI: 10.1021/acs.jpca.5b01869] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
48
|
Jimenez-Angeles F, Firoozabadi A. Enhanced Hydrate Nucleation Near the Limit of Stability. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2015; 119:8798-8804. [PMID: 25822221 DOI: 10.1021/acs.jpcc.5b01869] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Clathrate hydrates are crystalline structures composed of small guest molecules trapped into cages formed by hydrogen-bonded water molecules. In hydrate nucleation, water and the guest molecules may stay in a metastable fluid mixture for a long period. Metastability is broken if the concentration of the guest is above certain limit. We perform molecular dynamics (MD) simulations of supersaturated water-propane solutions close to the limit of stability. We show that hydrate nucleation can be very fast in a very narrow range of composition at moderate temperatures. Propane density fluctuations near the fluid-fluid demixing are coupled with crystallization producing en- hanced nucleation rates. This is the first report of propane-hydrate nucleation by MD simulations. We observe motifs of the crystalline structure II in line with experiments and new hydrate cages not reported in the literature. Our study relates nucleation to the fluid-fluid spinodal decomposition and demonstration that the enhanced nucleation phenomenon is more general than short range attractive interactions as suggested in nucleation of proteins.
Collapse
|
49
|
Michalis VK, Costandy J, Tsimpanogiannis IN, Stubos AK, Economou IG. Prediction of the phase equilibria of methane hydrates using the direct phase coexistence methodology. J Chem Phys 2015; 142:044501. [DOI: 10.1063/1.4905572] [Citation(s) in RCA: 90] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Vasileios K. Michalis
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23847, Doha, Qatar
| | - Joseph Costandy
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23847, Doha, Qatar
| | - Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research NCSR “Demokritos,” Aghia Paraskevi, Attiki GR-15310, Greece
| | - Athanassios K. Stubos
- Environmental Research Laboratory, National Center for Scientific Research NCSR “Demokritos,” Aghia Paraskevi, Attiki GR-15310, Greece
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23847, Doha, Qatar
| |
Collapse
|
50
|
Das S, Baghel VS, Roy S, Kumar R. A molecular dynamics study of model SI clathrate hydrates: the effect of guest size and guest–water interaction on decomposition kinetics. Phys Chem Chem Phys 2015; 17:9509-18. [PMID: 25767053 DOI: 10.1039/c5cp00678c] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
One of the options suggested for methane recovery from natural gas hydrates is molecular replacement of methane by suitable guests like CO2 and N2.
Collapse
Affiliation(s)
- Subhadip Das
- Physical Chemistry Division
- National Chemical Laboratory
- Pune-411008
- India
| | - Vikesh Singh Baghel
- Chemical Engineering and Process Development
- National Chemical Laboratory
- Pune-411008
- India
| | - Sudip Roy
- Physical Chemistry Division
- National Chemical Laboratory
- Pune-411008
- India
| | - Rajnish Kumar
- Chemical Engineering and Process Development
- National Chemical Laboratory
- Pune-411008
- India
| |
Collapse
|