Next Article in Journal
Pedot:PSS/Graphene Oxide (GO) Ternary Nanocomposites for Electrochemical Applications
Previous Article in Journal
Metabolomic Analysis of Phytochemical Compounds from Ethanolic Extract of Lime (Citrus aurantifolia) Peel and Its Anti-Cancer Effects against Human Hepatocellular Carcinoma Cells
Previous Article in Special Issue
Undirected C-H Bond Activation in Aluminium Hydrido Enaminonates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Dynamics Study of Clathrate-like Ordering of Water in Supersaturated Methane Solution at Low Pressure

by
Rodion V. Belosludov
1,
Kirill V. Gets
2,3,
Ravil K. Zhdanov
2,3,
Yulia Y. Bozhko
2,3,
Vladimir R. Belosludov
2,3,*,
Li-Jen Chen
4 and
Yoshiyuki Kawazoe
5,6,7
1
Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan
2
Nikolaev Institute of Inorganic Chemistry, Siberian Branch, Russian Academy of Sciences, 630090 Novosibirsk, Russia
3
Department of Physics, Novosibirsk State University, 630090 Novosibirsk, Russia
4
Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
5
New Industry Creation Hatchery Center, Tohoku University, Sendai 980-8579, Japan
6
Department of Physics and Nanotechnology, SRM Institute of Science and Technology, Kattankurathur 603203, India
7
School of Physics, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(7), 2960; https://doi.org/10.3390/molecules28072960
Submission received: 23 February 2023 / Revised: 23 March 2023 / Accepted: 24 March 2023 / Published: 26 March 2023
(This article belongs to the Special Issue New Insights into Intramolecular Hydrogen Bonds)

Abstract

:
Using molecular dynamics, the evolution of a metastable solution for “methane + water” was studied for concentrations of 3.36, 6.5, 9.45, 12.2, and 14.8 mol% methane at 270 K and 1 bar during 100 ns. We have found the intriguing behavior of the system containing over 10,000 water molecules: the formation of hydrate-like structures is observed at 6.5 and 9.45 mol% concentrations throughout the entire solution volume. This formation of “blobs” and the following amorphous hydrate were studied. The creation of a metastable methane solution through supersaturation is the key to triggering the collective process of hydrate formation under low pressure. Even the first stage (0–1 ns), before the first fluctuating cavities appear, is a collective process of H-bond network reorganization. The formation of fluctuation cavities appears before steady hydrate growth begins and is associated with a preceding uniform increase in the water molecule’s tetrahedrality. Later, the constantly presented hydrate cavities become the foundation for a few independent hydrate nucleation centers, this evolution is consistent with the labile cluster and local structure hypotheses. This new mechanism of hydrogen-bond network reorganization depends on the entropy of the cavity arrangement of the guest molecules in the hydrate lattice and leads to hydrate growth.

1. Introduction

Gas hydrates are crystalline compounds with hydrogen-bonded water molecules acting as the host network and guest molecules placed in the host cavities. Typically, for gas hydrate formation, high pressure and low temperature are required. The most common hydrate structures in nature are the cubic structure I (sI), the cubic structure II (sII), and the hexagonal structure H (sH). The formation of the different hydrate structures depends on the size of the guest molecule [1].
Gas hydrate formation is an important process for a number of different applications, such as gas separation, energy storage, energy transport, cold energy storage, sequestration, and desalination. Hydrate formation from the liquid phase involves the crucial step of hydrate nucleation, which is very important for developing efficient hydrates for practical applications [2].
Hydrophobic hydration is important for understanding the dissolution of hydrophobic molecules and the formation of clathrate hydrate. The solubility of hydrophobic molecules in water under normal conditions is very low and does not affect the structure of the solution as a whole but changes it locally. For the formation of a hydrate, it is necessary to dissolve several orders of magnitude of additional molecules. In this case, it is not a local effect; as a result of hydrophobic hydration, the entire structure of the solution is transformed.
The typical time evolution of gas consumption during the hydrate formation process can be divided into three main phases: the dissolution phase, the metastable supersaturated phase, and the growth phase [2]. The dissolution phase starts with increasing gas uptake due to the dissolution into the liquid phase.
Experimental molecular-level description of the hydrate nucleation process in water is challenging not only due to the extremely small length/time scale involved but also due to the stochastic nature of the nucleation process, which is only possible with state-of-the-art measurements. It should also be noted that the anomalous behavior of water is also important. Density fluctuations of water with dimensions of ~0.6 nm have been observed [3] and later confirmed by theoretical investigations [4]. The properties of water are determined by its dynamic hydrogen-bond network [5], which influences the stochastic nature of the process of hydrate nucleation.
Over the past few decades numerous theoretical molecular dynamics (MD) simulations have been performed in which the authors unraveled crucial aspects of crystal nucleation in liquids (see reviews [2,6,7]). These MD simulations have been used to obtain a detailed understanding of the nucleation process in different systems, where the time required for the evolution of a liquid into a crystal is more or less accurately reproduced.
One hypothesis for hydrate formation [1] involves hydrate precursor structures, in which the water molecules form individual clathrate cages (hydrate precursors) around the dissolved guest molecules prior to the hydrate formation process. Then, if favorable conditions are met, these cages will agglomerate and form a solid clathrate hydrate. Using MD simulations of the hydration of hydrophobic molecules, volumetric numerical data describing the details of the process was obtained. The creation of structured hydration shells around the dissolved hydrophobic molecules was shown [8,9,10,11,12,13,14,15,16,17,18,19]. These clusters (often called “postcursors”) are preserved for an extended period of time in the liquid phase [20]. The degree of supercooling influences these shells, and more amorphous structures are formed at lower temperatures [21]. Conversely, in other studies [22,23] the authors did not mention any enhanced structure around the hydrophobic solutes and did not find clathrate cages in solutions with long lifetimes. Experimental studies were carried out to reveal the water structure around non-polar solutes [24,25,26,27,28]. These results are consistent with previous computer simulation data, where there are loosely organized hydration cages around the hydrophobic solutes but without evidence of the regular clathrate hydrate cages in solution that are present in water. For example, an NMR investigation of tetrahydrofuran (THF) hydrate formation observed THF hydration shells with no evidence of hydrate precursor formation [29].
There have been many research studies of the mechanism of gas hydrate nucleation [30,31,32,33,34,35,36]. Currently, there are two main molecular mechanisms for hydrate nucleation. The first is the labile cluster hypothesis (LCH) [30], and the second is the local structure hypothesis (LSH) [31]. According to the LCH, prior to hydrate formation, water molecules form labile, relatively long-lived, unstable clusters resembling the solid hydrate cavities. The forming structure is an agglomerate of cavities of various hydrate types [34,37]. Then, under favorable conditions, this agglomerate recombined into a solid clathrate hydrate phase [30], whose cage type ratio is closer to that of a hydrate crystal. This recombination, i.e., the formation of stable crystalline sI hydrate, passes through the polycrystal structure that consists of sI and sII nanocrystals connected by 51263 [34]. Conversely, LSH suggests that guest molecules are arranged in a configuration similar to that of the crystal due to concentration fluctuations, and then the water molecules rearrange around the guest molecules to form the hydrate structure proposed by Radhakrishnan and Trout [31]. The concept of LSH, in which the hydrate can appear first with the ordering of guest molecules followed by water molecules, differs from LCH. Moreover, the authors of LSH argue against LCH by showing that it is thermodynamically favorable for the labile clusters to disintegrate rather than agglomerate to form a larger cluster.
In the hydrate phase, the concentration of hydrophobic gas molecules such as methane is more than two orders of magnitude higher than the solubility of methane in liquid water [36]. In order to increase methane solubility, it is possible to lower the temperature and/or increase the pressure of the mixture. In one study [38], the authors examined metastable methane solubility using MD simulations. It was shown that to increase the methane concentration in the water phase by two orders of magnitude at a fixed temperature, it is also necessary to increase the pressure by more than two orders of magnitude. It follows from the results of another study [39], at T = 258.5 K and P = 1.2 bar, the equilibrium concentration of methane in water is xM = 0.0005 molar fraction, and at P = 945 bar, xM = 0.05. In that study, the critical solubility, above which metastable solutions spontaneously form a hydrate phase, was found to be about 0.05 molar fraction, comparable to the concentration of methane in the hydrate (xM = 0.07) with half-filling of all the hydrate cavities [39].
The very first stage of clathrate hydrate formation is characterized by gas consumption into the water phase with the formation of metastable supersaturated water. This stage is often called the induction time. Increasing gas concentration is observed until the critical nuclei of a new hydrate phase are formed. A critical nucleus can be defined as the minimum amount of a new phase that is capable of existing independently [39]. In water + gas systems, relaxation processes involving the structural redistribution of dissolved guest molecules occur during the induction time, leading to the initial stages of growth phase formation. The metastable supersaturation phase remains a key step toward a complete understanding of hydrate nucleation [40].
The direct production of a supersaturated aqueous hydrate gas solution is still a challenge. However, indirect methods of obtaining supersaturated solutions do exist. Thus, for example, solid amorphous solutions of N2, O2, CO, and Ar have been obtained [41,42], from which gas hydrates have been produced subsequently at unexpectedly low external pressures.
A study of the initial steps of hydrate formation [43] showed that a dodecahedral water cage (DWC) immersed in the bulk water phase can efficiently adsorb dissolved methane molecules, which causes the DWC lifetime to increase exponentially with the number of methane molecules adsorbed. At the same time, adsorbed methane molecules reduce the number of hydrogen bonds that tend to break the DWC and change the hydrogen-bond network surrounding the DWC from a liquid water structure toward a cage-like hydrate structure. Such behavior argues in favor of considering hydrate nucleation as a cooperative process, where the guest molecules can stabilize the hydrate cages as well as the hydrate structure itself [15,43]. Another approach to describe the initial evolution of water + gas systems toward hydrate formation was reported by Guo et al. [44], who proposed the cage adsorption hypothesis (CAH). This theory predicts a critical concentration of ~0.04 molar fraction for methane aqueous solution, above which methane hydrate begins to nucleate [39].
Several MD techniques and methods have been developed to obtain deeper insight into the mechanisms behind the initial stages of the metastable supersaturation phase and the growth of the hydrate phase under a high driving force (high excess pressure and low temperature supercooling). The key parameters found are order parameters [45,46]. Nucleation pathways associated with both host and guest ordering and the mutual influence of guest and host ordering are being explored. The resulting parameters are used to describe the dynamics and mechanisms of hydrate formation and other key properties.
A number of different simulation studies [47,48,49,50,51] have revealed that nano-bubbles in the gas phase trapped inside the water phase could act as centers of hydrate nucleation at the liquid-gas interface that function as heterogeneous nucleation centers. Additionally, studies from the Uchida group [52,53] have shown that the induction time of C2H6 and Xe hydrates is significantly lower if the water solution is made using a hydrate with nano-/micro- bubbles that act as nucleation centers.
To produce hydrate structure formation in the presence of a gas-water interface, MD simulation studies are carried out at high pressures (from 500 to 2000 bar) and/or low temperatures [14,15,16,17,18,19,23,32,33,34,38,39,54,55,56,57,58]. In the referenced works, the nucleation time varied from several hundreds of nanoseconds up to several microseconds depending on the exact thermodynamic (P-T) conditions, selection of guest type, and water molecule model.
A different approach was applied by Sarupria and Debeneditti [59] to study a supersaturated homogeneous aqueous methane solution. The authors simulated a methane solution with a 0.07 mol fraction of methane at T = 240 K and P = 200 bar and reported a several hundred nanosecond nucleation time for a 0.07 molar concentration of methane and found structures similar to those reported in the case of a gas-water interface at similar temperature [14,15,16,17,18,19,32,33,34,38,39]. The nucleation time is drastically reduced at higher methane concentrations (0.08) [60]. Moreover, increasing the degree of subcooling leads to the formation of more ordered hydrate structures. Thus, homogeneous supersaturated methane aqueous solutions make it possible to model hydrate growth under low overpressure or undercooling [21,60,61] in contrast to models with separated phases.
Later simulations showed four phases of nucleation [62]. There is an increase in solvated methane concentration in the aqueous domain via diffusion through the methane−water interface, the formation of unstructured clusters of methane molecules solvated in water, decreasing the water content of the nucleus to a value compatible with the type II methane clathrate hydrate composition, and the final reordering of solvated methane and water as in the ‘blob’ hypothesis [34].
The structures of methane, propane, and methane + propane hydrates formed from water mixtures with high structural order and low computational time have been obtained using MD simulations [60,63]. In these works, a moderate subcooling level has been applied, and the overall thermodynamic conditions considered are within the P and T range of practical interest. Other important steps to revealing the hydrate formation mechanism were made by other studies [35,55,56,59,64], where the initial stages of the metastable supersaturation phase and growth of the hydrate phase were studied at high excess pressure and low temperature supercooling with an aqueous phase prepared by melting the hydrate phase. The study of hydrate nucleation in a “methane + water + THF” system [65,66] showed the dependence of induction time, growth rate, and cage structure on the P-T conditions and THF concentration.
We later studied a dynamic hydrogen-bond network in the liquid water phase [5,67]. Based on this, we developed efficient methods to analyze hydrogen-bond networks and a method for determining the hydrate nucleus. In previous work [68], we showed the possibility of methane hydrate growth from a bulk water + methane solution or in the presence of sea salt and hydrate seed at excess pressure without notable induction time. However, the nucleation mechanism was not studied.
In the present work, we have focused on an examination of the initial stage of the methane hydrate phase growth mechanism from a supersaturated methane solution without excess pressure and supercooling applied using MD simulation methods. The initial solutions are based on a uniform distribution of gas molecules inside the bulk water phase, which results in an ideally stirred solution. The aim is to study the dynamics of hydrogen-bond network reconfiguration and the influence of fluctuating hydrate nuclei on the hydrate growth process in supersaturated methane solution, with an emphasis on studying the behavior of water molecules with a high degree of local order at the earliest steps.

2. Results and Discussion

2.1. Fluctuation Character of the H-Bond Network Structure at the Hydrate-like Structure Pre-Nucleation Interval

The significant effect of guest molecules on the properties of the H-bond network, leading to the formation of small 512 fluctuation cavities at the initial stage of the hydrate-like structure (amorphous hydrate) nucleation process in a methane solution at low pressure, has been shown.
The characteristic time dependence of the number of large (51262) and small (512) sI hydrate cavities in solution with C = 0.6 and P = 1 bar is presented in Figure 1 for a single trajectory at which the formation and dissipation of sI cavities were observed over the time period 0–10 ns. The simulation could be divided into several intervals according to the features of small and large cavity formation. In the first interval (up to 1 ns) no cavities were formed. The second interval, 1–4.5 ns, is characterized by the formation and dissipation of single hydrate cavities (fluctuation cavities) in water solution (see Figure 1b), which have lifetimes of about one picosecond and could be the predecessors of hydrate nuclei. In the next interval (4.5–6 ns), at least one cavity exists in the volume, but because of fluctuations, it is possible to observe up to three different separate small cavities (see Figure 1c). Then after 6 ns, transformation of the hydrogen-bond network is observed, which leads to the formation of 2-3-4 small cavities as shown in Figure 1d, and only after a simulation time of 7 ns, does the fluctuation of large sI cavities become detectable and then persistent after 7.5 ns, although their number fluctuates (see Figure 1e). After a 10 ns simulation time, steady growth of the amorphous hydrate is observed. The number of cavities grows, and fluctuations do not lead to the disappearance of a significant fraction of the cavities. During the modeling of system transformations, there is no particular spatial ordering of the cavities, and the entire hydrogen-bond network of water molecules is involved. In this simulation, we focused only on the small (512) and large (51262) cavities that form the sI hydrate, but as noted earlier, other cavity types are present in this system and show similar behavior.
Figure 2a,b shows the time dependence of the order parameters F3 and F4 in the first 10 ns of simulation and comparisons to the bulk hydrate, bulk water, and cavity-forming water molecules in solution. The difference between F3 and F4 for bulk water and F3_liquid and F4_liquid is due to the fact that the TIP4P model series does not have an H-O-H angle equal to the right tetrahedral angle of 109.47 degrees, which is used in some other potentials. The value of the order parameters for bulk water does not change over time because the water is in an equilibrium state. The sharp decrease in the F3 parameter (an increase in tetrahedrality compared to pure water) occurs during the first 1 ns, which shows steady H-bond network reorganization (the difference is ~0.01 and gradually increases). Qualitatively similar results are presented for the F4 local torsion angle order parameter: whole solution order is close to pure water; however, the order parameters of cavity forming molecules in solution for both F3 and F4 are very close to those of a crystal in the same temperature and pressure region. Thus, the scattered cavities formed in the H-bond network are topologically close to hydrates. In general, the presence of dissolved methane molecules has an ordering effect on the structure of water. This water structure ordering in the presence of high concentrations of guest molecules is similar to the behavior of water molecules described in [69], where clathrate-like structures were shown to form under ambient conditions in the presence of confinement and a hydrophobic environment. The authors noted that as hydrophobic compounds, many different molecules can be chosen, such as alkanes such as methane, ethane, propane, etc. In this case, one may see the earliest stages of hydrate formation [34,62], which can be caused by a high local concentration of methane and an ordered arrangement of gas molecules [70], which is initially ensured by the homogeneity of the solution. Unlike many other works, here the distributed volume cavity formation can be observed. The reason for this is connected with the collective rearrangement of solution water molecules bonded by the hydrogen-bond network, which is discussed next.
Figure 2c shows the time dependence of the portion of highly ordered water molecules for which the parameter F3 is lower than 0.025 (highly ordered water molecules, HOWMs). This value was chosen because the water molecules forming the cavities have an average value of F3 < 0.025. As expected, the order parameters of these molecules are much more “crystalline” than for other molecules. The increase in the number of HOWMs compared to bulk water confirms the ordering of the structure prior to the formation of regular hydrate cavities. The fluctuations in the number of HOWMs indicate the dynamic nature of the local environment of all molecules united by the H-bond network.
Figure 3 shows the visualization of HOWMs during the first nanoseconds. At t = 0 ns (Figure 3a), the number of HOWMs is insignificant; however, at t = 0.1 ns (Figure 3b), this number is about ~20–25% of all molecules. It is important to note that HOWMs are delocalized in this time period (Figure 3a–f) in the entire system volume and do not form a single local cluster. This allows the cavity to fluctuate and form with any group of H-bond network HOWMs.

2.2. Radial Distribution Function of HOWMs

The radial distribution functions (RDF) were calculated to characterize the reorganization of the H-bond network structure. Figure 4a–d shows the oxygen–oxygen RDF for HOWMs in methane solution, bulk water, and bulk hydrate systems. The solution model peaks transform from liquid type with three peaks to hydrate type with four peaks over the time period. The position of the third peak shifts closer to the hydrate type. The high value of the first peak for the solution could indicate that HOWMs form groups that are uniformly distributed and form hydrate-like structures/fragments. At a larger time scale, further convergence of the RDF lines for the HOWMs and the general structure of water can be expected. The gas–oxygen RDF is discussed in the Supplementary Materials and is presented in Figure S2.

2.3. Growth Outside the Stability Region

Using the MD method, we examined the evolution of two sets of independently generated water + methane systems with different methane concentrations in the system at P = 1 bar and T = 270 K. The total number of water molecules was 10,001 and 1200 (for C = 0.6-small). The chosen concentrations are much higher than the equilibrium methane concentration in water, and thus the supersaturated homogeneous solutions obtained here are metastable. The pressure and temperature fluctuations are presented in Figures S3a,b and S4a,b in the Supplementary Materials. The simulation for C = 0.2 over 100 ns (Figure 5a–d) shows neither hydrate-like structure formation nor phase separation but a metastable supersaturated solution with homogeneously distributed gas molecules. At methane concentrations C = 0.4 (Figure 5e–h) and C = 0.6 (Figure 5i–l), the simulation results in the appearance of a volume region with a different methane concentration. The region with higher methane concentrations corresponds to a hydrate-like structure phase. The other region is a water solution with a greatly decreased methane concentration in comparison with the initial supersaturated phase. Increasing concentration (C = 0.8, Figure 5m–p) leads to the appearance of the third region of high methane concentration—gas nanobubbles. The formation of a gas phase at C = 0.8 is a competing process with the amorphous hydrate formation process: unlike the other systems where all runs provide similar results, for C = 0.8 in the first run, behavior similar to C = 1.0 is observed, whereas in the second run we observed a hydrate-like structure growth process but with a nano-bubble formation. At C = 1.0 (Figure 5q–t), no hydrate-like structure formation was observed.
In the simulations of the small system for 700 ns, a similar mechanism of hydrate formation was observed. Longer simulation times showed a significant decrease in the gas concentration in solution and the transition of the gas to the amorphous hydrate phase (Figure 5u–x). In this case, a fluctuation rearrangement of the amorphous hydrate structure is observed, which, however, qualitatively does not undergo significant changes with further modeling. Figure shows the separated amorphous hydrate layer and liquid solution layer. Thus, the same system evolution path at different sizes can be seen and will be discussed later.
Complete 100 ns evolutions of metastable supersaturated solutions with C = 0.2, 0.4, 0.6, 0.8, and 1.0 are presented in Supplementary Videos S1–S5, respectively. The structures bonded by stable hydrogen bonds and incorporated into these structures can be seen at early stages that indicates the presence of “blobs” that is one step in the way towards hydrate formation [34]. According to a previous study [38], nanobubbles formed in the water phase at C = 0.8 and C = 1.0 could be considered as hydrate nucleation centers at the microsecond timescale [71].
Experimentally obtaining supersaturated liquid solutions at atmospheric pressure is a task that does not have a direct solution, but the formation of a solid solution seems to be possible through the deposition of water vapor and methane on a cooled substrate [41,42].
The chemical potential of water molecules in the hydrate (or hydrate-like structure) phase depends not only on the interaction between water molecules in the hydrate lattice but also on the number of guest molecules in the hydrate cavities and exists as an entropy contribution, as shown in our previous work [72]. The value of this contribution at a certain guest molecule concentration exceeds the difference between the chemical potential of water molecules in the liquid phase and the empty hydrate lattice framework. In this case, due to a sufficiently large number of dissolved methane molecules, the filling of the hydrate-like structures is large enough to make it energetically preferable with respect to the liquid phase even at atmospheric pressure. We observed this phenomenon at concentrations of 0.4 (6.5 mol%) and 0.6 (9.45 mol%) during MD simulations. The formation of cavities may be preceded by the distribution and concentration of methane molecules that provoke the growth of hydrate, as stated in [31,70]. In this case, the formation of an amorphous hydrate is observed, similar to the works [34,37], but with a difference, which lies in the delocalization and the formation of cavities in the entire system volume. In the case presented here, the systems go through the stages of ‘blob’ and amorphous hydrate formation, which could be followed by the formation of crystalline hydrate.
The molecular dynamics method is not designed to study thermodynamic stability; the focus is on studying the dynamics of system change [73], such as viscosity, diffusion, and other transport properties. Therefore, the study of thermodynamic stability requires the use of other methods or techniques, especially in the region close to phase equilibrium, where the influence of fluctuations cannot be neglected.
In molecular dynamics simulations, the characteristic cell size rarely exceeds 100 Å, that is quite sufficient to describe the phenomena of interest. Nevertheless, the formation of hydrate nanocrystals with a characteristic size of ~200 Å, which can be in an aqueous gas solution without causing crystallization of the solution, has been experimentally shown [74]. This once again confirms the unpredictability of the behavior of small systems, which is somewhat similar to manifestations of the confinement effect. The authors of [75] showed that nanoporous media make it possible to increase gas solubility/uptake under milder conditions. So, it was shown that the confinement effect makes it possible to obtain methane hydrate at a lower pressure [76] and even at 1 bar [77].
When the concentration increases to more than 0.8 (12.2 mol%), phase separation occurs, which is in agreement with theoretical work [60]. This leads to a significant decrease in methane concentration in the water phase and disruption of this solid structure formation mechanism. With a decrease to C = 0.2 (3.36 mol%) or lower methane concentration, the amount of dissolved methane molecules in the water phase is insufficient for an adequately large entropy contribution to the water chemical potential, which was confirmed by the presence of a minimum methane concentration in the solution to initiate any solid structure formation [39].
Molecular dynamics simulations using the TIP4P/Ice water model [55] showed the growth of methane hydrate at T = 245 K and P = 0.2 bar, which is outside the methane hydrate stability area [78]. The authors hypothesized that unexpected hydrate growth is due to the spherical shape of the initially separated gas phase, which can create additional pressure and cause stable hydrate growth [55].
Previous results [21,60,61] assume that using a supersaturated methane solution as the initial structure allows one to observe hydrate formation at moderate conditions such as lower pressure and/or lower sub-cooling. Our results show that the effect of supersaturation of an aqueous solution of methane has greater importance than was previously suggested and allows hydrate-like structure formation at 1 bar.
Figure 6a shows the smoothed dependence of the potential energy of systems with different contents of guest molecules depending on the simulation time. It can be seen that for a system with a concentration C = 0.2, the energy value does not change with time, which indicates the absence of any evolution of the system. However, at higher concentrations (C = 0.4, 0.6), a monotonic decrease in potential energy is observed, which is accompanied by a rearrangement of the hydrogen-bond network of the system into an amorphous hydrate. At a concentration of C = 0.8, the behavior of the potential energy indicates two parallel processes: the process of rearranging the network of hydrogen bonds into a hydrate-like structure at the initial stage of simulation (up to 20 ns) and the process of isolating guest molecules into a separate phase, which remains the only one after 20 ns of simulation. A further increase in the concentration of guest molecules (C = 1.0) leads to the immediate release of the gas phase of guest molecules, which is expressed in an increase in potential energy. These graphs show the advantage of the formation of hydrate-like structures at certain concentrations of methane from the point of view of potential energy.
The construction of a solution model with a concentration of guest molecules significantly exceeding the thermodynamically stable one creates the conditions for the formation of amorphous hydrate. It can be assumed that while the concentration of guest molecules in the liquid phase of water is high, the process of hydrate formation will continue. However, as the guest molecules in the liquid phase are exhausted, the driving force of hydrate formation caused by their presence will gradually decrease, which will subsequently lead, first, to a slowdown in growth and then to the beginning of its dissociation. In our modeling, we observe only the initial stage of this process, where only amorphous hydrate growth is observed. With continued modeling, one may expect “amorphous hydrate—methane solution” phase separation (as it was shown in Figure 5u–x) and the formation of hydrate, whose structure will be closer to crystalline, taking into account amorphous cavities on its surface. Although the decomposition of hydrate-like structures is eventually expected due to thermodynamic conditions, it can be seen that the resulting structures are quite viable, as shown by the modeling of a small system (Figure 5u–x). The stability of the detected phases and their evolution with time are not fully studied here and can be the subject of a separate study.
Figure 6b shows the dependencies of the total number of 512, 51262, 51263, and 51264 hydrate, and 4151062, 4151063, and 4151064 hydrate-like cavities, forming an amorphous hydrate agglomeration, in the systems. A continuing increase in the number of cavities is seen for systems C = 0.4 and 0.6. For the C = 0.4 system, about ~13% of methane molecules are captured by cavities of various types during the formation of amorphous hydrate, and further growth can be expected. In the case of the system C = 0.6, by the end of the first 100 ns of simulation, more than 25% of methane molecules are captured. More molecules could be captured by unfinished cavities (i.e., cavity caps). Further, the rate of cavity formation slows down significantly, which is consistent with the fact that the formation of cavities at C = 0.2 is not observed in the considered time interval, and the composition of the liquid phase of the solution is still far from equilibrium. Due to phase separation, no growth is observed in the system C = 0.8 after 50 ns of simulation; however, the decay (dissociation) of the formed structure is also not observed. The increase in potential energy shown in Figure 6a for this system is due to the continued release of methane from solution.
The volume distribution of hydrate cavities for the cavity forming systems C = 0.4, 0.6, 0.8 is presented at Figure 7. Initial steps of delocalized hydrate formation in the entire system volume are clearly observed. Fluctuation formations and the decay of cavities were observed throughout the simulation time. Due to the relatively low methane concentration at C = 0.4, some hydrate nuclei are dissipating or merging into the larger one. Higher methane concentrations (C = 0.6 and 0.8) lead to more sustainable growth in the number of nucleation centers. As a result, the formation of an amorphous hydrate can be seen. The decrease in cavity number for C = 0.8 that may look like localization is associated with the release of gas from the solution shown at Figure 5m–p. The figures clearly show that the cavities in the resulting structures are connected to each other to form an amorphous structure. So, it is not possible to distinguish a clear crystal structure of sI, sII, or other types on the given time scale. For the sake of clarity, the 4151062, 4151063, and 4151064 cavities were not shown; however, these cavities are directly connected with the presented cavities, and due to their topology, they can easily transform to the hydrate type.
The dependence of the number of cavities on time for systems with different methane concentrations is presented in Figure 8. Methane forms an sI hydrate that consists only of 512 (Figure 8a) and 51262 (Figure 8b) cavities. Nevertheless, during the structure formation process, other cavity types could be formed because of fluctuations, surface effects, and other effects. These cavity types were also drawn in the plots (Figure 8c–g). For systems with concentrations of C = 0.4, 0.6, and 0.8, the number of observed cavities increases with time. It should be noted that the formation of large (51262) cavities occurs after the onset of the formation of small (512) ones. An increase in the number of different types of cavities is strongly correlated with methane concentrations in the water phase. The system with C = 0.6 shows the fastest speed, and, therefore, this concentration is the most favorable for hydrate-like structure formation under these conditions. We believe that at different pressure/temperature conditions, the most favorable methane concentration should be nearly the same. In other words, amorphous hydrate formation is faster if the gas/water molar relation is within a certain range. Similar dependencies for series A and several smaller models are given in the Supplementary Materials (Figures S5a–g and S6a–g).
The sI hydrate unit cell consists of six large and two small cavities, but in our calculations, this relation is not met due to a large number of other cavity types (51263, 51264, 4151062, 4151063, and 4151064) that reconfigure in later formation phases. Thus, the 4151062 cavity is topologically very close to the large 51262 cavity present in the sI hydrate and can easily reconfigure into a proper large 51262 cavity. Using such an assumption, the resulting relationship is not far from the target 3/1 value. Disintegration (change of structure/size) of the cavities is observed only as fluctuations in time. In the time interval of 0–100 ns, growth of the amorphous hydrate (C = 0.4 and 0.6) was observed. The case of C = 0.8 is exceptional, where slow amorphous hydrate decomposition is caused by gas phase separation. However, the resulting phases (0.4 and 0.6) are not thermodynamically stable, but under the given p-T conditions, they are metastable. Additional driving forces (such as supercooling, high pressure, promoters, or a helper gas) could be used to stabilize the structure.
The dependence of the normalized H-bond number on time and concentration is shown in Figure 9a and confirms previous results. The growth of NH-bond/Nmol indicates an increasing number of molecules involved in the formation of four H-bonds that are connected with crystal-like ordering, where each molecule forms exactly two donor and two acceptor H-bonds (NH-bond/Nmol = 2). At C = 0.2, the local structure of the solution does not change in terms of mean values. The decrease in NH-bond/Nmol at C = 1.0 clearly indicates that a high proportion of molecules form the phase boundary. The NH-bond/Nmol steadily increases for solutions with initial C = 0.4 and 0.6, and in the case of C = 0.8, during the first 25 ns (amorphous hydrate growth). In this case, the growth rate of the hydrate-like structure from solution with C = 0.8 is higher than at 0.6 and much higher than at C = 0.4.
The distribution of average NH-bond/Nmol in dependence on methane concentration makes (Figure 9b) it possible to estimate the methane concentration boundaries for successful hydrate-like structure formation. At C > 0.2 and C < 0.8 a stable increase in the NH-bond/Nmol number is observed. Thus, this range could be used for hydrate-like structure formation at low pressure.
The F3 and F4 order parameters were calculated for a 0–100 ns time interval (Figure 9c,d, respectively) for the considered solutions. Solutions with C = 0.4 and 0.6 show an increase in the gradual ordering of the H-bond network during the entire simulation time. The difference from the ideal value is due to the fact that a significant fraction of water remained in the liquid phase. The solution with C = 0.8 revealed ordering until t = 23 ns; however, the number of hydrate cavities increased. The discrepancy between these results is explained by an increase in methane nano-bubbles and the influence of the bubble surface. The disordering of an H-bond network solution with C = 1.0 caused by methane bubble formation is expressed as an increase in the F3 parameter and a decrease in the F4 parameter. The solution with C = 0.2 shows no significant change over time. The value F4 ≈ –0.03 for C = 0.2 and 1.0 corresponds to a liquid (bulk water), and the small difference from F4_liquid is not significant and can be caused by the features of the geometry of the chosen water model.
The obtained results indicate the cooperative nature of the hydrate-like structure nucleation process, which is determined by the rearrangement of the hydrogen-bond network, which is associated with an individual change in the tetrahedral order parameter of each water molecule in solution.
A change in the distribution of gas molecules also indicates a partial transition of the solution to the amorphous hydrate. Due to the initial uniform mixing of gas molecules, the peaks in the regions of 3–5 Å (first peak) and 6–8 Å (second peak) are present (Figure 10a–e) at t = 1 ns. The comparison with bulk hydrate and the increase in the first and second peaks of g(r) with time, as well as the appearance of the third peak, clearly confirm amorphous hydrate growth in solutions with C = 0.4 (Figure 10b), 0.6 (Figure 10c), and 0.8 (Figure 10d). The value of g(r) for the second peak exceeds the height of the first peak (0.4 and 0.6) or reaches a close value (0.8); however, in the case of 0.4 and 0.6, the distribution becomes constant with distance, indicating the presence of a significant amount of methane mixed with water even after t = 100 ns. However, in a system with a methane concentration of 0.8, due to the formation of a bubble (phase separation), the distribution value does not reach a constant value with increasing gas-gas distance. Moreover, at t = 20 ns, the growth of the second peak is insignificant. This is in agreement with theoretical works on carbon dioxide [79] and methane [62] hydrate formation.
The second peak of g(r) for the C = 1.0 solution (Figure 10e) smoothens over 20 ns, which is faster than the disordering of order parameters taking place during 30–40 ns and indicates the separation of phases is nearly complete, but residual methane segregates slowly. A decrease in g(r) with r clearly confirms the phase separation.
For the solution with C = 0.2 (Figure 10a), where no hydrate growth or phase separation is observed over 100 ns, no qualitative differences appear over time. The constant value of g(r) at high r values indicates that the methane concentration is not high enough to break the H-bond network.

3. Calculation Methods

Five homogeneous methane + water mixture models with different compositions (see Table 1) were constructed. The selected methane concentrations (C) correspond to 20%, 40%, 60%, 80%, and 100% rates of single filling of small (512) and large (51262) cavities of sI hydrate composed of 10,001 and 1200 (C = 0.6-small) water molecules. In contrast with previous works [35,55,56,59,64], where the aqueous phase was (homogeneously) supersaturated with gas and prepared by melting the hydrate phase, in the current study the initial supersaturated solution was similar to a different study [60]. The water molecules were placed at the nodes of a virtual cubic lattice with a step of 3 Å and a random orientation in space. The number of lattice sites slightly exceeded the number of water molecules, so some sites remained empty, but the model remained cubic. The gas molecules were also randomly arranged at the nodes of another independent virtual lattice with the same step size. For verification purposes, two series of simulations (A and B) were carried out with different distributions of water and gas molecules. By default, we present the results corresponding to series B, and the qualitative differences between series A and B are described. Before initiating the simulation, molecular non-overlap was confirmed.
All simulations were conducted at T = 270 K and P = 1 bar in the NPT-ensemble with a 1 fs timestep using LAMMPS software [80]. The time constants of the Nose–Hoover thermostat and barostat [81,82] implemented in LAMMPS were 300 fs and 3000 fs, respectively. Details of selected parameters are presented in the Supplementary Materials. Long-range interaction was calculated using the PPPM method [83]. The cut-off distance of the Coulomb and van der Waals forces was set to 13 Å. The SHAKE algorithm [84] was used in order to maintain the geometry of the water molecules. Periodic boundary conditions were applied for every model.
Snapshots for analysis purposes were written every 100 ps for all long-term simulations (more than 10 ns). In order to study the pre-nucleation processes and early steps of hydrate-like structure nucleation itself for the large 0.6 model (10,001 H2O + 1043 CH4) during the first 10 ns, snapshots were taken every 100 fs. In order to reveal structure changes due to nucleation, the pure water model containing 10,001 H2O molecules and the bulk methane hydrate model containing 6 × 6 × 6 sI unit cells (9936 H2O + 1728 CH4) were simulated at the same conditions over 10 ns and 1 ns, respectively.
The TIP4P/Ice potential [85] (σ = 3.1668 Å, ε = 0.21095 kcal/mol, qO = −1.1794|e|, qH = +0.5897|e|, O-M distance = 0.1577 Å, O-H distance = 0.9572 Å, and H-O-H angle = 104.52°) was used to describe the water molecules due to its high melting point accuracy in comparison to other simple 3-point and 4-point potentials. Water molecules were considered rigid. Methane molecules were considered single-point particles, with potential parameters described elsewhere [86].
The structures of the homogeneous mixture (series A and B, series of smaller systems) and the simulation run parameters with force-field data are presented in Supplementary Materials.
The analysis of the hydrogen-bond (H-bond) network was carried out as described in our previous works [5,67]. The geometric criterion for H-bonding was chosen as follows: ROO ≤ 3.2 Å and ∠HOH ≤ 30°. The local (short-range order) tetrahedrality order parameter F3 [45] and torsion angle order parameter F4 [46] were used in this study. Reference points of the F3 parameter for the ideal water crystal and bulk liquid are F3_crystal = 0 and F3_liquid ≈ 0.1. Reference points of the F4 parameter for the hydrate, bulk liquid, and ice are F4_hydrate ≈ 0.7, F4_liquid ≈ 0, and F4_ice ≈ −0.3.
To find the cavities in a water solution, a graph of the connectivity of water molecules by H-bonds was created based on the geometric criterion of H-bonding. The closed chains of 4, 5, and 6 H2O molecules were found using this graph. These chains were considered the faces of a polygon. Subsequently, a connection graph of the faces was created through their common (adjacent) sides. Using this graph, the closed polyhedral structures were searched, among which sI hydrate cavities (512 and 51262) and topologically close hydrate ones 4151062, 4151063, and 4151064 cavities were selected and studied. The obtained sets of coordinates of water molecules were visualized.

4. Conclusions

The collective nature of the methane hydrate nucleation process in the dynamic H-bond network of TIP4P/Ice water outside the region of thermodynamic stability was considered. At T = 270 K and P = 1 bar, the process of amorphous hydrate formation depended on the methane concentration in water was observed. At methane concentrations of 0.4 (6.5 mol%) and 0.6 (9.45 mol%), growth of the amorphous hydrate phase was observed, which showed the intriguing behavior of methane-containing hydrate-like structure formation outside the region of thermodynamic stability of methane hydrate. At a concentration of 0.8 (12.2 mol%), the competing process of the formation of methane bubbles in water was found. Early stages of hydrate formation are observed: ordering of the hydrogen-bond network, formation of ‘blobs’ throughout the entire solution volume, and amorphous hydrate. It can be assumed that the nature of this process is connected with the sufficiently high entropy contribution of the guest molecules to the chemical potential of the water molecules in the hydrate phase. This contribution allows the transformation of the H-bond network of a metastable liquid solution into the H-bond network of a more energetically preferable amorphous hydrate structure. The further behavior of obtained hydrate-like structures needs additional research at different thermodynamic conditions, however, the small system simulation showed the dynamic amorphous hydrate with a certain number of hydrate cavities exists at the microsecond time scale.
The collective character of H-bond network reorganization in methane solution was revealed as an increase in water short-range order tetrahedrality at the initial stage of simulation (0–1 ns) before the appearance of the first fluctuation of sI hydrate cavities. Molecules with a high degree of tetrahedrality (low F3 order parameter values) are distributed uniformly. A significant change in F3 preceded the formation of cavities, and thus the entire network of hydrogen bonds was involved in the amorphous hydrate formation process. Fluctuating cavities appeared in the interval 1.5–6 ns. Later, the fluctuating and stable cavities were constantly present in the system, and an increase in their number indicated the initial stage of the formation of nanohydrate phases.
The formation of hydrate cavities and hydrate nanostructures under the described conditions is possible using real-time Raman spectrometry as in Ref. [87].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28072960/s1, Barostat and thermostat time parameter choice: discussion; Figure S1: Spatial distribution of small and large cavities; Figure S2: Unformed by the number of G-O pairs in the methane-oxygen RDF of 0.6 methane solution in comparison to bulk hydrate data under the same conditions; Figure S3: Time dependence of system pressure at different methane concentrations.; Figure S4: Time dependence of system temperature at different methane concentrations; Table S1: Designation of methane + water solution compositions for large series A systems and composition of smaller systems; Video S1: Evolution of methane + water solution at a concentration of 0.2 (3.36 mol%) during 100 ns of modeling; Video S2: Evolution of a methane + water solution at a concentration of 0.4 (6.5 mol%) during 100 ns of modeling; Video S3: Evolution of a methane + water solution at concentration of 0.6 (9.45 mol%) during 100 ns of modeling; Video S4: Evolution of methane + water solution at concentration of 0.8 (12.2 mol%) during 100 ns of modeling; Video S5: Evolution of a methane + water solution at a concentration of 1.0 (14.8 mol%) during 100 ns of modeling; Figure S5: Time dependence of the number of 512, 51262, 4151062, 51263, 4151063, 51264, 4151064 type cavities in series A simulations; Figure S6: Time dependence of the number of 512, 51262, 4151062, 51263, 4151063, 51264, 4151064 type cavities in smaller model; Simulation data: starting homogeneous solution structures for series A, B, and smaller systems, and LAMMPS run scripts. References [88,89,90,91,92] are cited in the supplementary materials.

Author Contributions

Conceptualization, V.R.B. and R.V.B.; methodology, V.R.B.; software, K.V.G. and R.K.Z.; validation, K.V.G., R.K.Z. and Y.Y.B.; formal analysis, L.-J.C. and Y.K.; investigation, R.V.B.; data curation, V.R.B.; writing—original draft preparation, K.V.G. and R.K.Z.; writing—review and editing, R.V.B., V.R.B., Y.Y.B., L.-J.C. and Y.K.; visualization, K.V.G. and R.K.Z.; supervision, V.R.B.; project administration, V.R.B. All authors have read and agreed to the published version of the manuscript.

Funding

V.R.B., K.V.G., Z.R.K., Y.Y.B. and L.-J.C. are grateful to the Russian Foundation for Basic Research (project number 21-52-52001) jointly with the Ministry of Science and Technology of Taiwan (grant no. 110-2923-E-002-008-MY3) for funding this research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated and analyzed during the current study are available from the corresponding author by request.

Acknowledgments

Rodion V. Belosludov is grateful to the E-IMR Center at the Institute for Materials Research, Tohoku University. The authors are also thankful to Tina M. Briere for carefully reading this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Sloan, E.D., Jr.; Koh, C.A. Clathrate Hydrates of Natural Gases, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2007; p. 752. [Google Scholar] [CrossRef]
  2. Khurana, M.; Yin, Z.; Linga, P. A Review of Clathrate Hydrate Nucleation. ACS Sustain. Chem. Eng. 2017, 5, 11176–11203. [Google Scholar] [CrossRef]
  3. Bosio, L.; Teixeira, J.; Stanley, H.E. Enhanced density fuctuations in supercooled H2O, D2O, and ethanol-water solutions: Evidence from small-angle x-ray scattering. Phys. Rev. Lett. 1981, 46, 597–600. [Google Scholar] [CrossRef]
  4. Geiger, A.; Mausbach, P.; Schnitker, J.; Blumberg, R.L.; Stanley, H.E. Structure and dynamics of the hydrogen bond network in water by computer simulations. J. Phys. Colloq. 1984, 45, C7-13–C7-30. [Google Scholar] [CrossRef] [Green Version]
  5. Belosludov, V.; Gets, K.; Zhdanov, R.; Malinovsky, V.; Bozhko, Y.; Belosludov, R.; Surovtsev, N.; Subbotin, O.; Kawazoe, Y. The nano-structural inhomogeneity of dynamic hydrogen bond network of TIP4P/2005 water. Sci. Rep. 2020, 10, 7323. [Google Scholar] [CrossRef]
  6. Sosso, G.C.; Chen, J.; Cox, S.J.; Fitzner, M.; Pedevilla, P.; Zen, A.; Michaelides, A. Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations. Chem. Rev. 2016, 116, 7078–7116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. English, N.J.; MacElroy, J.M.D. Perspectives on molecular simulation of clathrate hydrates: Progress, prospects and challenges. Chem. Eng. Sci. 2015, 121, 133–156. [Google Scholar] [CrossRef]
  8. Owicki, J.C.; Scheraga, H.A. Monte Carlo calculation in isothermal-isobaric ensemble.2. Dilute aqueous solution of methane. J. Am. Chem. Soc. 1977, 99, 7413–7418. [Google Scholar] [CrossRef]
  9. Swaminathan, S.; Harrison, S.W.; Beveridge, D.L. Monte Carlo studies on the structure of dilute aqueous solution of methane. J. Am. Chem. Soc. 1978, 100, 5705–5712. [Google Scholar] [CrossRef]
  10. Geiger, A.; Rahman, A.; Stillinger, F.H. Molecular dynamics study of the hydration of Lennard-Jones solutes. J. Chem. Phys. 1979, 70, 263–276. [Google Scholar] [CrossRef] [Green Version]
  11. Guillot, B.; Guissani, Y. A computer simulation study of the temperature dependence of the hydrophobic hydration. J. Chem. Phys. 1993, 99, 8075–8094. [Google Scholar] [CrossRef]
  12. Chau, P.L.; Mancera, R.L. Computer simulation of structural effect of pressure on the hydrophobic hydration of methane. Mol. Phys. 1999, 96, 109–122. [Google Scholar] [CrossRef]
  13. Long, J.; Sloan, E.D. Quantized water cluster around apolar molecules. Mol. Simul. 1993, 11, 145–161. [Google Scholar] [CrossRef]
  14. Skipper, N.T. Computer simulation of methane-water solutions. Evidence for a temperature-dependent hydrophobic attraction. Chem. Phys. Lett. 1993, 207, 424–429. [Google Scholar] [CrossRef]
  15. Mancera, R.L.; Buckingham, A.D.; Skipper, N.T. The aggregation of methane in aqueous solution. J. Chem. Soc. Faraday Trans. 1997, 93, 2263–2267. [Google Scholar] [CrossRef]
  16. Silverstein, K.A.T.; Haymet, A.D.J.; Dill, K.A. Molecular model of hydrophobic solvation. J. Chem. Phys. 1999, 111, 8000–8009. [Google Scholar] [CrossRef]
  17. Silverstein, K.A.T.; Haymet, A.D.J.; Dill, K.A. A simple model of water and the hydrophobic effect. J. Am. Chem. Soc. 1998, 120, 3166–3175. [Google Scholar] [CrossRef]
  18. Mancera, R.L. Computer simulation of effect of salt on the hydrophobic effect. J. Chem. Soc. Faraday Trans. 1998, 94, 3549–3559. [Google Scholar] [CrossRef]
  19. Lazaridis, T. Solvent reorganization energy and entropy in hydrophobic hydration. J. Phys. Chem. B 2000, 104, 4964–4979. [Google Scholar] [CrossRef]
  20. Lederhos, J.P.; Long, J.P.; Sum, A.; Christiansen, R.L.; Sloan, E.D. Effective kinetic inhibitors for natural gas hydrates. Chem. Eng. Sci. 1996, 51, 1221–1229. [Google Scholar] [CrossRef]
  21. Arjun Berendsen, T.A.; Bolhuis, P.G. Unbiased atomistic insight in the competing nucleation mechanisms of methane hydrates. Proc. Natl. Acad. Sci. USA 2019, 116, 19305–19310. [Google Scholar] [CrossRef] [Green Version]
  22. Meng, E.C.; Kollman, P.A. Molecular dynamics studies of the properties of water around simple organic solutes. J. Phys. Chem. 1996, 100, 11460–11470. [Google Scholar] [CrossRef]
  23. Hernandez-Cobos, J.; Mackie, A.D.; Vega, L.F. The hydrophobic hydration of methane as a function of temperature from histogram reweighting Monte Carlo simulations. J. Chem. Phys. 2001, 114, 7527–7535. [Google Scholar] [CrossRef]
  24. Broadbent, R.D.; Neilson, G.W. The interatomic structure of argon in water. J. Chem. Phys. 1994, 100, 7543–7547. [Google Scholar] [CrossRef]
  25. De Jong, P.H.K.; Wilson, J.E.; Neilson, G.W.; Buckingham, A.D. Hydrophobic hydration of methane. Mol. Phys. 1997, 91, 99–104. [Google Scholar] [CrossRef]
  26. Koh, C.A.; Wisbey, R.P.; Wu, X.; Westacott, R.E.; Scoper, A.K. Water ordering around methane during hydrate formation. J. Chem. Phys. 2000, 113, 6390–6397. [Google Scholar] [CrossRef]
  27. Filipponi, A.; Bowron, D.T.; Lobban, C.; Finney, J.L. Structural determination of hydrophobic hydration shell of Kr. Phys. Rev. Lett. 1997, 79, 1293–1296. [Google Scholar] [CrossRef]
  28. Bowron, D.T.; Filipponi, A.; Roberts, M.A.; Finney, J.L. Hydrophobic hydration and the formation of a clathrate hydrate. Phys. Rev. Lett. 1998, 81, 4164–4167. [Google Scholar] [CrossRef]
  29. Gao, S.; Chapman, W.G.; House, W. NMR and Viscosity Investigation of Clathrate Hydrate Formation and Dissociation. Ind. Eng. Chem. Res. 2005, 44, 7373–7379. [Google Scholar] [CrossRef]
  30. Sloan, E.D.; Fleyfel, F.A. A molecular mechanism for gas hydrate nucleation from ice. AIChE J. 1991, 37, 1281–1292. [Google Scholar] [CrossRef]
  31. Radhakrishnan, R.; Trout, B.L. A new approach for studying nucleation phenomena using molecular simulations: Application to CO2 hydrate clathrates. J. Chem. Phys. 2002, 117, 1786–1796. [Google Scholar] [CrossRef] [Green Version]
  32. Moon, C.; Taylor, P.C.; Rodger, P.M. Molecular Dynamics Study of Gas Hydrate Formation. J. Am. Chem. Soc. 2003, 125, 4706–4707. [Google Scholar] [CrossRef] [PubMed]
  33. Hawtin, R.W.; Quigley, D.; Rodger, P.M. Gas hydrate nucleation and cage formation at a water/methane interface. Phys. Chem. Chem. Phys. 2008, 10, 4853–4864. [Google Scholar] [CrossRef]
  34. Jacobson, L.C.; Hujo, W.; Molinero, V. Amorphous Precursors in the Nucleation of Clathrate Hydrates. J. Am. Chem. Soc. 2010, 132, 11806–11811. [Google Scholar] [CrossRef]
  35. Vatamanu, J.; Kusalik, P.G. Observation of two-step nucleation in methane hydrates. Phys. Chem. Chem. Phys. 2010, 12, 15065–15072. [Google Scholar] [CrossRef] [PubMed]
  36. Guo, G.-J.; Zhang, Y.-G.; Liu, C.-J.; Li, K.-H. Using the face-saturated incomplete cage analysis to quantify the cage compositions and cage linking structures of amorphous phase hydrates. Phys. Chem. Chem. Phys. 2011, 13, 12048–12057. [Google Scholar] [CrossRef] [PubMed]
  37. Jacobson, L.C.; Hujo, W.; Molinero, V. Nucleation pathways of clathrate hydrates: Effect of guest size and solubility. J. Phys. Chem. B 2010, 114, 13796–13807. [Google Scholar] [CrossRef]
  38. Walsh, M.R.; Koh, C.A.; Sloan, E.D.; Sum, A.K.; Wu, D.T. Microsecond simulations of spontaneous methane hydrate nucleation and growth. Science 2009, 326, 1095–1098. [Google Scholar] [CrossRef] [PubMed]
  39. Guo, G.-J.; Rodger, P.M. Solubility of Aqueous Methane under Metastable Conditions: Implications for Gas Hydrate Nucleation. J. Phys. Chem. B 2013, 117, 6498–6504. [Google Scholar] [CrossRef]
  40. Rodger, P. Methane hydrate: Melting and memory. Ann. N. Y. Acad. Sci. 2000, 912, 474–482. [Google Scholar] [CrossRef]
  41. Hallbrucker, A.; Mayer, E. Unexpectedly stable nitrogen, oxygen, carbon monoxide and argon clathrate hydrates from vapour-deposited amorphous solid water: An X-ray and two-step differential scanning calorimetry study. J. Chem. Soc. Faraday Trans. 1990, 86, 3785–3792. [Google Scholar] [CrossRef]
  42. Hallbrucker, A.; Mayer, E. Unexpectedly stable clathrate hydrates formed from microporous vapor-deposited amorphous solid water at low “external” guest pressures and their astrophysical implications. Icarus 1991, 90, 176–180. [Google Scholar] [CrossRef]
  43. Guo, G.-J.; Zhang, Y.-G.; Liu, H. Effect of Methane Adsorption on the Lifetime of a Dodecahedral Water Cluster Immersed in Liquid Water: A Molecular Dynamics Study on the Hydrate Nucleation Mechanisms. J. Phys. Chem. C 2007, 111, 2595–2606. [Google Scholar] [CrossRef]
  44. Guo, G.-J.; Li, M.; Zhang, Y.-G.; Wu, C.-H. Why can water cages adsorb aqueous methane? A potential of mean force calculation on hydrate nucleation mechanisms. Phys. Chem. Chem. Phys. 2009, 11, 10427–10437. [Google Scholar] [CrossRef] [PubMed]
  45. Baez, L.A.; Clancy, P. Computer Simulation of the Crystal Growth and Dissolution of Natural Gas Hydrates. Ann. N. Y. Acad. Sci. 1994, 715, 177–186. [Google Scholar] [CrossRef]
  46. Rodger, P.M.; Forester, T.R.; Smith, W. Simulations of the methane hydrate/methane gas interface near hydrate forming conditions conditions. Fluid Ph. Equilib. 1996, 116, 326–332. [Google Scholar] [CrossRef]
  47. Yagasaki, T.; Matsumoto, M.; Andoh, Y.; Okazaki, S.; Tanaka, H. Effect of bubble formation on the dissociation of methane hydrate in water: A molecular dynamics study. J. Phys. Chem. B 2014, 118, 1900–1906. [Google Scholar] [CrossRef]
  48. Bagherzadeh, S.A.; Alavi, S.; Ripmeester, J.A.; Englezos, P. Formation of methane nano-bubbles during hydrate decomposition and their effect on hydrate growth. J. Chem. Phys. 2015, 142, 214701. [Google Scholar] [CrossRef] [PubMed]
  49. Knott, B.C.; LaRue, J.L.; Wodtke, A.M.; Doherty, M.F.; Peters, B. Communication: Bubbles, crystals, and laser-induced nucleation. J. Chem. Phys. 2011, 134, 171102. [Google Scholar] [CrossRef] [PubMed]
  50. Matsumoto, M.; Wada, Y.; Oonaka, A.; Onoe, K. Polymorph control of glycine by antisolvent crystallization using nitrogen minute-bubbles. J. Cryst. Growth 2013, 373, 73–77. [Google Scholar] [CrossRef]
  51. Zhang, Y.; Zhao, L.; Deng, S.; Zhao, R.; Nie, X.; Liu, Y. Effect of Nanobubble Evolution on Hydrate Process: A Review. J. Therm. Sci. 2019, 28, 948–961. [Google Scholar] [CrossRef]
  52. Uchida, T.; Yamazaki, K.; Gohara, K. Generation of micro-and nano-bubbles in water by dissociation of gas hydrates. Korean J. Chem. Eng. 2016, 33, 1749–1755. [Google Scholar] [CrossRef] [Green Version]
  53. Uchida, T.; Yamazaki, K.; Gohara, K. Gas nanobubbles as nucleation acceleration in the gas-hydrate memory effect. J. Phys. Chem. C 2016, 120, 26620–26629. [Google Scholar] [CrossRef]
  54. Ripmeester, J.A.; Alavi, S. Some current challenges in clathrate hydrate science: Nucleation, decomposition and the memory effect. Curr. Opin. Solid State Mater. Sci. 2016, 20, 344–351. [Google Scholar] [CrossRef] [Green Version]
  55. Walsh, M.R.; Beckham, G.T.; Koh, C.A.; Sloan, E.D.; Wu, D.T.; Sum, A.K. Methane Hydrate Nucleation Rates from Molecular Dynamics Simulations: Effects of Aqueous Methane Concentration, Interfacial Curvature, and System Size. J. Phys. Chem. C 2011, 115, 21241–21248. [Google Scholar] [CrossRef]
  56. Walsh, M.R.; Rainey, J.D.; Lafond, P.G.; Park, D.H.; Beckham, G.T.; Jones, M.D.; Lee, K.-H.; Koh, C.A.; Sloan, E.D.; Wu, D.T.; et al. The Cages, Dynamics, and Structuring of Incipient Methane Clathrate Hydrates. Phys. Chem. Chem. Phys. 2011, 13, 19951–19959. [Google Scholar] [CrossRef] [Green Version]
  57. Liang, S.; Kusalik, P.G. Crystal growth simulations of H2S hydrate. J. Phys. Chem. B 2010, 114, 9563–9571. [Google Scholar] [CrossRef]
  58. Liang, S.; Kusalik, P.G. Exploring nucleation of H2S hydrates. Chem. Sci. 2011, 2, 1286–1292. [Google Scholar] [CrossRef]
  59. Sarupria, S.; Debenedetti, P.G. Homogeneous Nucleation of Methane Hydrate in Microsecond Molecular Dynamics Simulations. J. Phys. Chem. Lett. 2012, 3, 2942–2947. [Google Scholar] [CrossRef]
  60. Jimenez-Angeles, F.; Firoozabadi, A. Nucleation of methane hydrates at moderate subcooling by molecular dynamics simulations. J. Phys. Chem. C 2014, 118, 11310–11318. [Google Scholar] [CrossRef]
  61. Choudhary, N.; Kushwaha, O.S.; Bhattacharjee, G.; Chakrabarty, S.; Kumar, R. Macro and Molecular Level Insights on Gas Hydrate Growth in the Presence of Hofmeister Salts. Ind. Eng. Chem. Res. 2020, 59, 20591–20600. [Google Scholar] [CrossRef]
  62. Lauricella, M.; Meloni, S.; English, N.J.; Peters, B.; Ciccotti, G. Methane clathrate hydrate nucleation mechanism by advanced molecular simulations. J. Phys. Chem. C 2014, 118, 22847–22857. [Google Scholar] [CrossRef]
  63. Jiménez-Ángeles, F.; Firoozabadi, A. Nucleation pathways of gas hydrates from molecular dynamics simulations. In Proceedings of the 8th International Conference on Gas Hydrates (ICGH8), Beijing, China, 28 July–1 August 2014. [Google Scholar] [CrossRef]
  64. Barnes, B.C.; Sum, A.K. Advances in molecular simulations of clathrate hydrates. Curr. Opin. Chem. Eng. 2013, 2, 184–190. [Google Scholar] [CrossRef]
  65. Wu, J.Y.; Chen, L.J.; Chen, Y.P.; Lin, S.T. Molecular dynamics study on the growth mechanism of methane plus tetrahydrofuran mixed hydrates. J. Phys. Chem. C 2015, 119, 19883–19890. [Google Scholar] [CrossRef]
  66. Wu, J.Y.; Chen, L.J.; Chen, Y.P.; Lin, S.T. Molecular dynamics study on the nucleation of methane+ tetrahydrofuran mixed guest hydrate. Phys. Chem. Chem. Phys. 2016, 18, 9935–9947. [Google Scholar] [CrossRef]
  67. Gets, K.; Belosludov, V.; Zhdanov, R.; Bozhko, Y.; Belosludov, R.; Subbotin, O.; Marasanov, N.; Kawazoe, Y. Transformation of hydrogen bond network during CO2 clathrate hydrate dissociation. Appl. Surf. Sci. 2020, 499, 143644. [Google Scholar] [CrossRef] [Green Version]
  68. Belosludov, V.R.; Gets, K.V.; Zhdanov, R.K.; Bozhko, Y.Y.; Belosludov, R.V.; Chen, L.J. Collective Effect of Transformation of a Hydrogen Bond Network at the Initial State of Growth of Methane Hydrate. JETP Lett. 2022, 115, 124–129. [Google Scholar] [CrossRef]
  69. 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. [Google Scholar] [CrossRef]
  70. Lauricella, M.; Ciccotti, G.; English, N.J.; Peters, B.; Meloni, S. Mechanisms and nucleation rate of methane hydrate by dynamical nonequilibrium molecular dynamics. J. Phys. Chem. C 2017, 121, 24223–24234. [Google Scholar] [CrossRef]
  71. Uchida, T.; Miyoshi, H.; Sugibuchi, R.; Suzuta, A.; Yamazaki, K.; Gohara, K. Contribution of ultra-fine bubbles to promoting effect on propane hydrate formation. Front. Chem. 2020, 8, 480. [Google Scholar] [CrossRef]
  72. Belosludov, V.R.; Subbotin, O.S.; Krupskii, D.S.; Belosludov, R.V.; Kawazoe, Y.; Kudoh, J.I. Physical and chemical properties of gas hydrates: Theoretical aspects of energy storage application. Mater. Trans. 2007, 48, 704–710. [Google Scholar] [CrossRef] [Green Version]
  73. Das, S.; Tadepalli, K.M.; Roy, S.; Kumar, R. A review of clathrate hydrate nucleation, growth and decomposition studied using molecular dynamics simulation. J. Mol. Liquid. 2022, 348, 118025. [Google Scholar] [CrossRef]
  74. Lehmkühler, F.; Paulus, M.; Sternemann, C.; Lietz, D.; Venturini, F.; Gutt, C.; Tolan, M. The carbon dioxide−water interface at conditions of gas hydrate formation. J. Am. Chem. Soc. 2009, 131, 585–589. [Google Scholar] [CrossRef] [PubMed]
  75. Both, A.K.; Gao, Y.; Zeng, X.C.; Cheung, C.L. Gas hydrates in confined space of nanoporous materials: New frontier in gas storage technology. Nanoscale 2021, 13, 7447–7470. [Google Scholar] [CrossRef] [PubMed]
  76. Casco, M.E.; Silvestre-Albero, J.; Ramírez-Cuesta, A.J.; Rey, F.; Jordá, J.L.; Bansode, A.; Urakawa, A.; Peral, I.; Martínez-Escandell, M.; Kaneko, K.; et al. Methane hydrate formation in confined nanospace can surpass nature. Nat. Commun. 2015, 6, 6432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Yu, K.B.; Yazaydin, A.O. Does confinement enable methane hydrate growth at low pressures? Insights from molecular dynamics simulations. J. Phys. Chem. C 2020, 124, 11015–11022. [Google Scholar] [CrossRef] [PubMed]
  78. Makogon, T.Y.; Sloan, E.D.J. Phase equilibrium for methane hydrate from 190 to 262K. J. Chem. Eng. Data 1994, 39, 351–353. [Google Scholar] [CrossRef]
  79. Hirai, S.; Okazaki, K.; Tabe, Y.; Kawamura, K. CO2 clathrate-hydrate formation and its mechanism by molecular dynamics simulation. Energy Convers. Manag. 1997, 38, S301–S306. [Google Scholar] [CrossRef]
  80. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  81. Nosé, S. A molecular dynamics method for simulations in the canonical ensemble. Mol. Phys. 1984, 52, 255–268. [Google Scholar] [CrossRef]
  82. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [Green Version]
  83. Hockney, R.W.; Eastwood, J.W. Computer Simulation Using Particles, 1st ed.; CRC Press: Boca Raton, FL, USA, 1988; p. 540. [Google Scholar]
  84. Ryckaert, J.P.; Ciccotti, G.; Berendsen, H.J. Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 1977, 23, 327–341. [Google Scholar] [CrossRef] [Green Version]
  85. Abascal, J.L.F.; Sanz, E.; Fernández, R.G.; Vega, C. A potential model for the study of ices and amorphous water: TIP4P/Ice. J. Chem. Phys. 2005, 122, 234511. [Google Scholar] [CrossRef] [Green Version]
  86. Martin, M.G.; Siepmann, J.I. Transferable Potentials for Phase Equilibria. 1. United-atom Description of n-Alkanes. J. Phys. Chem. B 1998, 102, 2569–2577. [Google Scholar] [CrossRef]
  87. Kang, K.C.; Hong, S.Y.; Cho, S.J.; Kim, D.H.; Lee, J.D. Evaluation of desalination by nanostructured hydrate formation and pellet production process. J. Nanosci. Nanotechnol. 2017, 17, 4059–4062. [Google Scholar] [CrossRef]
  88. Mi, F.; He, Z.; Fang, B.; Ning, F.; Jiang, G. Molecular insights into the effects of surface property and pore size of non-swelling clay on methane hydrate formation. Fuel 2022, 311, 122607. [Google Scholar] [CrossRef]
  89. Tung, Y.T.; Chen, L.J.; Chen, Y.P.; Lin, S.T. Molecular dynamics study on the growth of structure I methane hydrate in aqueous solution of sodium chloride. J. Phys. Chem. B 2012, 116, 14115–14125. [Google Scholar] [CrossRef]
  90. DeFever, R.S.; Sarupria, S. Nucleation mechanism of clathrate hydrates of water-soluble guest molecules. J. Chem. Phys. 2017, 147, 204503. [Google Scholar] [CrossRef]
  91. English, N.J.; John, S.T.; Carey, D.J. Mechanisms for thermal conduction in various polymorphs of methane hydrate. Phys. Rev. B 2009, 80, 134306. [Google Scholar] [CrossRef]
  92. Conde, M.M.; Rovere, M.; Gallo, P. High precision determination of the melting points of water TIP4P/2005 and water TIP4P/Ice models by the direct coexistence technique. J. Chem. Phys. 2017, 147, 244506. [Google Scholar] [CrossRef]
Figure 1. (a) Time dependence of proper large and small sI cavity number. Spatial distribution of small and large cavities at different simulation time points: (b) t = 1.500 ns; (c) t = 4.651 ns; (d) t = 8.365 ns; (e) t = 10.0 ns. Intermediate time steps are presented in Figure S1 in Supplementary Materials.
Figure 1. (a) Time dependence of proper large and small sI cavity number. Spatial distribution of small and large cavities at different simulation time points: (b) t = 1.500 ns; (c) t = 4.651 ns; (d) t = 8.365 ns; (e) t = 10.0 ns. Intermediate time steps are presented in Figure S1 in Supplementary Materials.
Molecules 28 02960 g001
Figure 2. Time dependence of (a) F3 and (b) F4 order parameters of 0.6 H2O + CH4 solution (black line), pure water (orange line), and bulk hydrate (blue line) phases (time period corresponding to 0–1 ns and stretched to 0–10 ns for comparison). The green open circles represent the corresponding data for cavity-forming molecules. The data are smoothed; (c) Time dependence of the portion of all molecules having molecular parameter F3 < 0.025 in 0.6 methane solution, bulk hydrate, and pure water. The data are smoothed.
Figure 2. Time dependence of (a) F3 and (b) F4 order parameters of 0.6 H2O + CH4 solution (black line), pure water (orange line), and bulk hydrate (blue line) phases (time period corresponding to 0–1 ns and stretched to 0–10 ns for comparison). The green open circles represent the corresponding data for cavity-forming molecules. The data are smoothed; (c) Time dependence of the portion of all molecules having molecular parameter F3 < 0.025 in 0.6 methane solution, bulk hydrate, and pure water. The data are smoothed.
Molecules 28 02960 g002
Figure 3. Spatial distribution of highly ordered water molecules (HOWMs) in 0.6 solution at different simulation time points: (a) t = 0.0 ns; (b) t = 0.1 ns; (c) t = 1.0 ns; (d) t = 3.0 ns; (e) t = 6.0 ns; (f) t = 10.0 ns.
Figure 3. Spatial distribution of highly ordered water molecules (HOWMs) in 0.6 solution at different simulation time points: (a) t = 0.0 ns; (b) t = 0.1 ns; (c) t = 1.0 ns; (d) t = 3.0 ns; (e) t = 6.0 ns; (f) t = 10.0 ns.
Molecules 28 02960 g003
Figure 4. The highly ordered water molecule (HOWM) oxygen–oxygen (eigen F3 < 0.025) RDF g(r) for the 0.6 methane solution, pure water, and bulk phase methane hydrate at different simulation time points: (a) t = 0.1 ns; (b) t = 4.0 ns; (c) t = 6.0 ns; (d) t = 10.0 ns. Water and hydrate data are averaged over the simulation time. The rOO distance distribution step is 0.01 Å.
Figure 4. The highly ordered water molecule (HOWM) oxygen–oxygen (eigen F3 < 0.025) RDF g(r) for the 0.6 methane solution, pure water, and bulk phase methane hydrate at different simulation time points: (a) t = 0.1 ns; (b) t = 4.0 ns; (c) t = 6.0 ns; (d) t = 10.0 ns. Water and hydrate data are averaged over the simulation time. The rOO distance distribution step is 0.01 Å.
Molecules 28 02960 g004aMolecules 28 02960 g004b
Figure 5. Spatial distribution of H-bonds and gas molecules in the methane + water solution at concentrations of (ad) C = 0.2, (eh) C = 0.4, (il) C = 0.6, (mp) C = 0.8 and (qt) C = 1.0 at simulation time points: 1 ns, 20 ns, 40 ns, and 100 ns and (ux) for replicated 3 × 1 × 1 times C = 0.6-small at time points 1, 50, 100, and 700 ns. Coordinates of water molecules are averaged over 1 ns in order to show the stable hydrogen bonds (red dashed lines). Hydrate-like structure growth at C = 0.8 is not characteristic.
Figure 5. Spatial distribution of H-bonds and gas molecules in the methane + water solution at concentrations of (ad) C = 0.2, (eh) C = 0.4, (il) C = 0.6, (mp) C = 0.8 and (qt) C = 1.0 at simulation time points: 1 ns, 20 ns, 40 ns, and 100 ns and (ux) for replicated 3 × 1 × 1 times C = 0.6-small at time points 1, 50, 100, and 700 ns. Coordinates of water molecules are averaged over 1 ns in order to show the stable hydrogen bonds (red dashed lines). Hydrate-like structure growth at C = 0.8 is not characteristic.
Molecules 28 02960 g005aMolecules 28 02960 g005b
Figure 6. Time dependence of (a) system potential energy and (b) overall cavity number. Solid lines represent the smoothed data of the cavity number.
Figure 6. Time dependence of (a) system potential energy and (b) overall cavity number. Solid lines represent the smoothed data of the cavity number.
Molecules 28 02960 g006
Figure 7. Spatial distribution of 512 (orange), 51262 (blue), 51263 (green), and 51264 (magenta) cavities for (ac) C = 0.4, (df) C = 0.6, and (gi) C = 0.8 systems at simulation time points: 20 ns, 40 ns, and 100 ns. Hydrate-like structure growth at a concentration of 0.8 is not characteristic; the drawing angle is the same as in Figure 5.
Figure 7. Spatial distribution of 512 (orange), 51262 (blue), 51263 (green), and 51264 (magenta) cavities for (ac) C = 0.4, (df) C = 0.6, and (gi) C = 0.8 systems at simulation time points: 20 ns, 40 ns, and 100 ns. Hydrate-like structure growth at a concentration of 0.8 is not characteristic; the drawing angle is the same as in Figure 5.
Molecules 28 02960 g007
Figure 8. Time dependence of the number of (a) 512; (b) 51262; (c) 4151062; (d) 51263; (e) 4151063; (f) 51264; (g) 4151064 type cavities of amorphous methane hydrate agglomeration. Solid lines represent the smoothed data.
Figure 8. Time dependence of the number of (a) 512; (b) 51262; (c) 4151062; (d) 51263; (e) 4151063; (f) 51264; (g) 4151064 type cavities of amorphous methane hydrate agglomeration. Solid lines represent the smoothed data.
Molecules 28 02960 g008
Figure 9. (a) Time dependence of H-bond number per water molecule and (b) normalized H-bond number as a function of gas concentration at different time points. Time dependence of (c) F3 and (d) F4 order parameters. When t = 0 ns, F3 = 0.1 for all methane concentrations. Solid lines at (a,c,d) represent the smoothed data.
Figure 9. (a) Time dependence of H-bond number per water molecule and (b) normalized H-bond number as a function of gas concentration at different time points. Time dependence of (c) F3 and (d) F4 order parameters. When t = 0 ns, F3 = 0.1 for all methane concentrations. Solid lines at (a,c,d) represent the smoothed data.
Molecules 28 02960 g009
Figure 10. Gas–gas RDFs g(r) in (a) C = 0.2; (b) C = 0.4; (c) C = 0.6; (d) C = 0.8; (e) C = 1.0 solutions at different time points in comparison with the bulk hydrate phase. The rGG distance distribution step is 0.01 Å.
Figure 10. Gas–gas RDFs g(r) in (a) C = 0.2; (b) C = 0.4; (c) C = 0.6; (d) C = 0.8; (e) C = 1.0 solutions at different time points in comparison with the bulk hydrate phase. The rGG distance distribution step is 0.01 Å.
Molecules 28 02960 g010
Table 1. Designation of concentration, C, for methane (NM) + water (NW) solution compositions and correspondence to molar concentration (xM).
Table 1. Designation of concentration, C, for methane (NM) + water (NW) solution compositions and correspondence to molar concentration (xM).
CNMNWxM
0.234810,0013.36 mol%
0.469610,0016.5 mol%
0.6104310,0019.45 mol%
0.8139110,00112.2 mol%
1.0173910,00114.8 mol%
0.6-small12512009.43 mol%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Belosludov, R.V.; Gets, K.V.; Zhdanov, R.K.; Bozhko, Y.Y.; Belosludov, V.R.; Chen, L.-J.; Kawazoe, Y. Molecular Dynamics Study of Clathrate-like Ordering of Water in Supersaturated Methane Solution at Low Pressure. Molecules 2023, 28, 2960. https://doi.org/10.3390/molecules28072960

AMA Style

Belosludov RV, Gets KV, Zhdanov RK, Bozhko YY, Belosludov VR, Chen L-J, Kawazoe Y. Molecular Dynamics Study of Clathrate-like Ordering of Water in Supersaturated Methane Solution at Low Pressure. Molecules. 2023; 28(7):2960. https://doi.org/10.3390/molecules28072960

Chicago/Turabian Style

Belosludov, Rodion V., Kirill V. Gets, Ravil K. Zhdanov, Yulia Y. Bozhko, Vladimir R. Belosludov, Li-Jen Chen, and Yoshiyuki Kawazoe. 2023. "Molecular Dynamics Study of Clathrate-like Ordering of Water in Supersaturated Methane Solution at Low Pressure" Molecules 28, no. 7: 2960. https://doi.org/10.3390/molecules28072960

APA Style

Belosludov, R. V., Gets, K. V., Zhdanov, R. K., Bozhko, Y. Y., Belosludov, V. R., Chen, L. -J., & Kawazoe, Y. (2023). Molecular Dynamics Study of Clathrate-like Ordering of Water in Supersaturated Methane Solution at Low Pressure. Molecules, 28(7), 2960. https://doi.org/10.3390/molecules28072960

Article Metrics

Back to TopTop