Joanna Grabowska

Understanding the mechanisms by which gas hydrates nucleate from aqueous solutions is essential for predicting and controlling their formation. Experimental observations have suggested that hydrate nucleation occurs preferentially at interfaces between the aqueous phase and hydrate-forming reservoirs. However, recent simulation studies have raised questions regarding the general validity of this assumption. Here, we investigate the spontaneous nucleation of H2S hydrates at 500 bar using molecular dynamics simulations employing the TIP4P/Ice water model and a four-site H2S potential, with the aim of clarifying the competition between bulk and interfacial nucleation pathways. To ensure rigorous control of the thermodynamic conditions, we first determine the three-phase coexistence temperature (T3) from the intersection of the aqueous solubility and hydrate–solution equilibrium curves; our calculated value of T3 = 314(3) K shows excellent agreement with an experimentally inferred value of approximately 316 K. Then, under deeply supercooled conditions, we directly compare nucleation kinetics across two distinct configurations: an interfacial system with an aqueous solution in contact with an H2S reservoir and a standalone bulk solution at identical saturation. Our results show that nucleation rates in bulk systems are comparable to those obtained in the presence of an interface, indicating that the interface does not promote hydrate nucleation under the conditions explored. This behavior is consistent with our previous findings for CO2 hydrates and suggests that bulk nucleation can be the prevailing pathway under deep supercooling conditions.

Gas hydrates are of great relevance to both the oil industry and the environment. Understanding how these solid structures nucleate from aqueous solutions is essential for controlling their formation. Experimental studies have often suggested that hydrate nucleation originates at the interface between the aqueous phase and the guest-molecule reservoir. To assess this hypothesis, we perform molecular dynamics simulations of CO₂ hydrate nucleation. First, we place hydrate seeds at different positions relative to the interface and monitor their evolution, finding that seeds embedded in the bulk are more likely to grow than those located near or at the interface. Second, we analyze spontaneous nucleation simulations with and without an interface. Our previous work showed that nucleation rates are indistinguishable in both systems, strongly indicating that the interface does not play a role. Here, trajectory analysis reveals that hydrates nucleate in regions of locally high CO₂ concentration, which arise spontaneously in the bulk and are not associated with the interface. Our results indicate that hydrate nucleation does not preferentially occur at the interface, at least under the deep supercooling conditions explored in this work. Further work at higher temperatures, and considering alternative nucleation locations, is needed to reconcile experiments and simulations, and thereby reach a deep understanding of the mechanism of hydrate formation.

We investigate the nucleation of carbon dioxide (CO2) hydrates from carbon dioxide aqueous solutions by means of molecular dynamics simulations using the TIP4P/Ice and the TraPPE models for water and CO2, respectively. We work at 400 bar and different temperatures and CO2 concentrations. We use brute force molecular dynamics when the supersaturation or the supercooling is so high so that nucleation occurs spontaneously and Seeding otherwise. We use both methods for a particular state and found an excellent agreement when using a linear combination of q̄3 and q̄12 order parameters to identify critical clusters. With such order parameter, we get a rate of 1025 m−3 s−1 for nucleation in a CO2 saturated solution at 255 K (35 K of supercooling). By comparison with our previous work on methane hydrates, we conclude that nucleation of CO2 hydrates is several orders of magnitude faster due to a lower interfacial free energy between the crystal and the solution. By combining our nucleation studies with a recent calculation of the hydrate–solution interfacial free energy at coexistence [Algaba et al., J. Colloid Interface Sci. 623, 354–367 (2022)], we obtain a prediction of the nucleation rate temperature dependence for CO2-saturated solutions (the experimentally relevant concentration). On the one hand, we open the window for comparison with experiments for supercooling larger than 25 K. On the other hand, we conclude that homogeneous nucleation is impossible for supercooling lower than 20 K. Therefore, nucleation must be heterogeneous in typical experiments where hydrate formation is observed at low supercooling. To assess the hypothesis that nucleation occurs at the solution-CO2 interface, we run spontaneous nucleation simulations in two-phase systems and find, by comparison with single-phase simulations, that the interface does not affect hydrate nucleation, at least at the deep supercooling at which this study was carried out (40 and 45 K). Overall, our work sheds light on molecular and thermodynamic aspects of hydrate nucleation.

Using molecular dynamics simulations, we show that a molecule of moderately active antifreeze protein (type III AFP, QAE HPLC-12 isoform) is able to interact with ice in an indirect manner. This interaction occurs between the ice binding site (IBS) of the AFP III molecule and the surface of ice, and it is mediated by liquid water, which separates these surfaces. As a result, the AFP III molecule positions itself at a specific orientation and distance relative to the surface of ice, which enables the effective binding (via hydrogen bonds) of the molecule with the nascent ice surface. Our results show that the final adsorption of the AFP III molecule on the surface of ice is not achieved by chaotic diffusion movements, but it is preceded by a remote, water-mediated interaction between the IBS and the surface of ice. The key factor that determines the existence of this interaction is the ability of water molecules to spontaneously form large, high-volume aggregates that can be anchored to both the IBS of the AFP molecule and the surface of ice. The results presented in this work for AFP III are in full agreement with the ones obtained by us previously for hyperactive CfAFP, which indicates that the mechanism of the remote interaction of these molecules with ice remains unchanged despite significant differences in the molecular structure of their ice binding sites. For that reason, we can expect that also other types of AFPs interact with the ice surface according to an analogous mechanism.

In this paper, the solubility of carbon dioxide (CO2) in water along the isobar of 400 bar is determined by computer simulations using the well-known TIP4P/Ice force field for water and the TraPPE model for CO2. In particular, the solubility of CO2 in water when in contact with the CO2 liquid phase and the solubility of CO2 in water when in contact with the hydrate have been determined. The solubility of CO2 in a liquid–liquid system decreases as the temperature increases. The solubility of CO2 in a hydrate–liquid system increases with temperature. The two curves intersect at a certain temperature that determines the dissociation temperature of the hydrate at 400 bar (T3). We compare the predictions with T3 obtained using the direct coexistence technique in a previous work. The results of both methods agree, and we suggest 290(2) K as the value of T3 for this system using the same cutoff distance for dispersive interactions. We also propose a novel and alternative route to evaluate the change in chemical potential for the formation of hydrates along the isobar. The new approach is based on the use of the solubility curve of CO2 when the aqueous solution is in contact with the hydrate phase. It considers rigorously the non-ideality of the aqueous solution of CO2, providing reliable values for the driving force for nucleation of hydrates in good agreement with other thermodynamic routes used. It is shown that the driving force for hydrate nucleation at 400 bar is larger for the methane hydrate than for the carbon dioxide hydrate when compared at the same supercooling. We have also analyzed and discussed the effect of the cutoff distance of dispersive interactions and the occupancy of CO2 on the driving force for nucleation of the hydrate.

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas–liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log₁₀ (J/(m³ s)) = −7(5).

In this paper, the solubility of methane in water along the 400 bar isobar is determined by computer simulations using the TIP4P/Ice force field for water and a simple LJ model for methane. In particular, the solubility of methane in water when in contact with the gas phase and the solubility of methane in water when in contact with the hydrate has been determined. The solubility of methane in a gas-liquid system decreases as temperature increases. The solubility of methane in a hydrate-liquid system increases with temperature. The two curves intersect at a certain temperature that determines the triple point T₃ at a certain pressure. We also determined T₃ by the three-phase direct coexistence method. The results of both methods agree, and we suggest 295(2) K as the value of T₃ for this system. We also analyzed the impact of curvature on the solubility of methane in water. We found that the presence of curvature increases the solubility in both the gas-liquid and hydrate-liquid systems. The change in chemical potential for the formation of hydrate is evaluated along the isobar using two different thermodynamic routes, obtaining good agreement between them. It is shown that the driving force for hydrate nucleation under experimental conditions is higher than that for the formation of pure ice when compared at the same supercooling. We also show that supersaturation (i.e., concentrations above those of the planar interface) increases the driving force for nucleation dramatically. The effect of bubbles can be equivalent to that of an additional supercooling of about 20 K. Having highly supersaturated homogeneous solutions makes possible the spontaneous formation of the hydrate at temperatures as high as 285 K (i.e., 10K below T₃). The crucial role of the concentration of methane for hydrate formation is clearly revealed. Nucleation of the hydrate can be either impossible or easy and fast depending on the concentration of methane which seems to play the leading role in the understanding of the kinetics of hydrate formation.

Using computer simulations, the structural properties of solvation water of three model hydrophobic molecules, methane and two fullerenes (C60 and C80), were studied. Systems were simulated at temperatures in the range of 250–298 K. By analyzing both the local ordering of the molecules of water in the solvation layers and the structure of hydrogen bond network, it is shown that in the solvation layer of hydrophobic molecules, ordered aggregates consisting of water molecules are formed. Even though it is difficult to define the exact structure of these aggregates, their existence alone is clearly noticeable. Moreover, these aggregates become more pronounced with the decrease of temperature. The existence of the ordered aggregates around the hydrophobic solutes complies with the concept of “icebergs” proposed by Frank and Evans.

Most protein molecules do not adsorb onto ice, one of the exceptions being so-called antifreeze proteins. In this paper, we describe that there is a force pushing an antifreeze protein molecule away from the ice surface when it is not oriented with its ice-binding plane toward the ice and that this pushing force may be also present even when the protein is oriented with its ice-binding plane toward the ice. This force is absent only when certain specific distance criteria are met, regarding the surface of ice and the protein. It acts at early stages of adsorption, prior to the solidification of water between the ice and the protein molecule nearby. We propose the water-originating mechanism of the generation of this force and also the mechanism of remote attachment of an antifreeze molecule to the ice surface. In liquid water, there exist locally favored structures, ordered and of high specific volume. The presence of a protein molecule usually shifts the equilibrium that exists in liquid water toward increasing the number of high-density, disordered structures and diminishing the number of low-density structures. Creation of the locally favored structures may be hampered not only near the non-ice-binding surfaces but also between the ice surface and the protein surface, if the distance between these surfaces does not allow these structures to develop because the available space is not sufficient for their proper formation. This conclusion is supported by the analysis of the mean geometry of a single hydrogen bond, as well as of the hydrogen bond network in the solvation layer and a structural order parameter that characterizes the separation between the first and second solvation shells of a water molecule.