1
|
Molecular Insights into the Effect of Nitrogen Bubbles on the Formation of Tetrahydrofuran Hydrates. Molecules 2022; 27:molecules27154945. [PMID: 35956899 PMCID: PMC9370114 DOI: 10.3390/molecules27154945] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Revised: 07/26/2022] [Accepted: 07/30/2022] [Indexed: 11/17/2022] Open
Abstract
In this work, a molecular dynamics simulation was conducted to study the microscopic mechanism of how nitrogen bubbles affect the formation of THF hydrates at the molecular level. The results obtained reveal that the nitrogen bubble can promote the formation of THF hydrates. In the system with a nitrogen bubble, more THF-filled cages were generated, and the crystal structure was more orderly. The promotion of nitrogen bubbles on hydrate crystallization comes from the dissolution of nitrogen molecules. Some of dissolved nitrogen molecules can be enclosed in small hydrate cages near the nitrogen bubble, which can serve as stable sites for hydrate crystal growth, resulting in the fact that THF-filled cages connected with N2-filled cages are much more stable and have a long lifetime. The results in this work can help to understand the promotion effect of micro- and nano-air bubbles on the crystallization of THF hydrates.
Collapse
|
2
|
Ismail NA, Koh CA. Growth Rate and Morphology Study of Tetrahydrofuran Hydrate Single Crystals and the Effect of Salt. CrystEngComm 2022. [DOI: 10.1039/d2ce00176d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
For an advanced gas hydrate management strategy in flow assurance as well as the application of a gas hydrate technology (e.g., desalination), it is important to understand the kinetics and...
Collapse
|
3
|
Atakoohi SE, Naeiji P, Peyvandi K, Mollashahi Sanatgar S. The experimental study and molecular dynamic simulation of THF hydrate growth kinetics in the presence of Arabic and Guar gum: New approaches in promotion of THF hydrate formation. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.115249] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
4
|
Nada H. Melt crystallization mechanism analyzed with dimensional reduction of high-dimensional data representing distribution function geometries. Sci Rep 2020; 10:15465. [PMID: 32963268 PMCID: PMC7508891 DOI: 10.1038/s41598-020-72455-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2020] [Accepted: 09/02/2020] [Indexed: 11/21/2022] Open
Abstract
Melt crystallization is essential to many industrial processes, including semiconductor, ice, and food manufacturing. Nevertheless, our understanding of the melt crystallization mechanism remains poor. This is because the molecular-scale structures of melts are difficult to clarify experimentally. Computer simulations, such as molecular dynamics (MD), are often used to investigate melt structures. However, the time evolution of the structural order in a melt during crystallization must be analyzed properly. In this study, dimensional reduction (DR), which is an unsupervised machine learning technique, is used to evaluate the time evolution of structural order. The DR is performed for high-dimensional data representing an atom–atom pair distribution function and the distribution function of the angle formed by three nearest neighboring atoms at each period during crystallization, which are obtained by an MD simulation of a supercooled Lennard–Jones melt. The results indicate that crystallization occurs via the following activation processes: nucleation of a crystal with a distorted structure and reconstruction of the crystal to a more stable structure. The time evolution of the local structures during crystallization is also evaluated with this method. The present method can be applied to studies of the mechanism of crystallization from a disordered system for real materials, even for complicated multicomponent materials.
Collapse
|
5
|
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: 11] [Impact Index Per Article: 2.2] [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
|
6
|
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
|
7
|
Nada H. Anisotropy in geometrically rough structure of ice prismatic plane interface during growth: Development of a modified six-site model of H 2O and a molecular dynamics simulation. J Chem Phys 2017; 145:244706. [PMID: 28049310 DOI: 10.1063/1.4973000] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
This paper presents a modified version of the six-site model of H2O [H. Nada and J. P. J. M. van der Eerden, J. Chem. Phys. 118, 7401 (2003)]. Although the original six-site model was optimized by assuming the cut-off of the Coulomb interaction at an intermolecular distance of 10 Å, the modified model is optimized by using the Ewald method for estimating the Coulomb interaction. Molecular dynamics (MD) simulations of an ice-water interface suggest that the melting point of ice at 1 atm in the modified model is approximately 274.5 K, in good agreement with the real melting point of 273.15 K. MD simulations of bulk ice and water suggest that the modified model reproduces not only the structures and density curves of ice and water, but also the diffusion coefficient of water molecules in water near the melting point at 1 atm. Using the modified model, a large-scale MD simulation of the growth at an ice-water interface of the prismatic plane is performed to elucidate the anisotropy in the interface structure during growth. Simulation results indicate that the geometrical roughness of the ice growth front at the interface is greater in the c-axis direction than in the direction normal to the c-axis when it is analyzed along the axes parallel to the prismatic plane. In addition, during the growth at the interface, the transient appearance of specific crystallographic planes, such as a {202¯1} pyramidal plane, occurs preferentially at the ice growth front. The effect of different ensembles with different simulation systems on the anisotropy in the interface structure is also investigated.
Collapse
Affiliation(s)
- Hiroki Nada
- National Institute of Advanced Industrial Science and Technology (AIST), 16-1 Onogawa, Tsukuba, Ibaraki 305-8569, Japan
| |
Collapse
|
8
|
Sujith KS, Ramachandran CN. Natural Gas Evolution in a Gas Hydrate Melt: Effect of Thermodynamic Hydrate Inhibitors. J Phys Chem B 2016; 121:153-163. [DOI: 10.1021/acs.jpcb.6b07782] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- K. S. Sujith
- Department of Chemistry, Indian Institute of Technology Roorkee, Roorkee 247667, India
| | - C. N. Ramachandran
- Department of Chemistry, Indian Institute of Technology Roorkee, Roorkee 247667, India
| |
Collapse
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
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
|
13
|
Yagasaki T, Matsumoto M, Andoh Y, Okazaki S, Tanaka H. Dissociation of Methane Hydrate in Aqueous NaCl Solutions. J Phys Chem B 2014; 118:11797-804. [DOI: 10.1021/jp507978u] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [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
| | - Yoshimichi Andoh
- Department
of Applied Chemistry, Nagoya University, Nagoya 464-8603, Japan
| | - Susumu Okazaki
- Department
of Applied Chemistry, Nagoya University, Nagoya 464-8603, 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
|
14
|
Seo M, Jang E, Kim K, Choi S, Kim JS. Understanding anisotropic growth behavior of hexagonal ice on a molecular scale: A molecular dynamics simulation study. J Chem Phys 2012; 137:154503. [DOI: 10.1063/1.4759113] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
|
15
|
Sarupria S, Debenedetti PG. Molecular Dynamics Study of Carbon Dioxide Hydrate Dissociation. J Phys Chem A 2011; 115:6102-11. [DOI: 10.1021/jp110868t] [Citation(s) in RCA: 90] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Sapna Sarupria
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo G. Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
16
|
SUN C, LI W, YANG X, LI F, YUAN Q, MU L, CHEN J, LIU B, CHEN G. Progress in Research of Gas Hydrate. Chin J Chem Eng 2011. [DOI: 10.1016/s1004-9541(09)60192-0] [Citation(s) in RCA: 69] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
17
|
Bauer M, Többens DM, Mayer E, Loerting T. Pressure-amorphized cubic structure II clathrate hydrate: crystallization in slow motion. Phys Chem Chem Phys 2011; 13:2167-71. [DOI: 10.1039/c0cp01351j] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
18
|
Finney AR, Rodger PM. Applying the Z method to estimate temperatures of melting in structure II clathrate hydrates. Phys Chem Chem Phys 2011; 13:19979-87. [DOI: 10.1039/c1cp21919g] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
19
|
Nada H, Furukawa Y. Growth inhibition at the ice prismatic plane induced by a spruce budworm antifreeze protein: a molecular dynamics simulation study. Phys Chem Chem Phys 2011; 13:19936-42. [DOI: 10.1039/c1cp21929d] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
|