Direct Coexistence Simulations Research Articles

Note: A simple correlation to locate the three phase coexistence line in methane-hydrate simulations

In this Note, the three phases coexistence (hydrate-water-methane) temperature is determined by using direct coexistence simulations for a force field in which methane is described by a single Lennard-Jones site, and water is described by either SPC or SPC/E models. Results were obtained for two pressures: 100 and 400 bars. The results of this note along with those obtained in our previous work [M. M. Conde and C. Vega, J. Chem. Phys. 133, 064507 (2010)]10.1063/1.3466751 allow to present the three phases coexistence line for a number of water models: SPC, SPC/E, TIP4P, TIP4P/2005, and TIP4P/Ice. A correlation between the three phases coexistence temperature T3 of the methane hydrate and the melting point of ice at normal pressure Tm of the water model is found. Models with a good prediction of the melting temperature of ice Ih, have also good predictions of the three phases coexistence line. Models with a bad prediction of the melting temperature of ice Ih, provide poor estimates of T3. As a rule of thumb, the three phase coexistence line of the methane hydrate at 100 bars for a water-methane force field is located approximately at about 15(7) K above the melting temperature of ice Ih of the considered water model. In summary, when modeling the methane-hydrate, if you are looking for a three phases coexistence line in good agreement with experiment, you need first a water model that describes accurately the melting temperature of ice Ih.

Solubility of NaCl in water by molecular simulation revisited

In this paper, the solubility of NaCl in water is evaluated by using computer simulations for three different force fields. The condition of chemical equilibrium (i.e., equal chemical potential of the salt in the solid and in the solution) is obtained at room temperature and pressure to determine the solubility of the salt. We used the same methodology that was described in our previous work [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] although several modifications were introduced to improve the accuracy of the calculations. It is found that the predictions of the solubility are quite sensitive to the details of the force field used. Certain force fields underestimate the experimental solubility of NaCl in water by a factor of four, whereas the predictions of other force fields are within 20% of the experimental value. Direct coexistence molecular dynamic simulations were also performed to determine the solubility of the salt. Reasonable agreement was found between the solubility obtained from free energy calculations and that obtained from direct coexistence simulations. This work shows that the evaluation of the solubility of salts in water can now be performed in computer simulations. The solubility depends on the ion-ion, ion-water, and water-water interactions. For this reason, the prediction of the solubility can be quite useful in future work to develop force fields for ions in water.

Melting and superheating of sI methane hydrate: Molecular dynamics study

Melting and decay of the superheated sI methane structure are studied using molecular dynamics simulation. The melting curve is calculated by the direct coexistence simulations in a wide range of pressures up to 5000 bar for the SPC/E, TIP4P/2005 and TIP4P/Ice water models and the united-atom model for methane. We locate the kinetic stability boundary of the superheated metastable sI structure that is found to be surprisingly high comparing with the predictions based on the classical nucleation theory.

Applying the Z method to estimate temperatures of melting in structure II clathrate hydrates

Equilibrium melting temperatures for structure II THF hydrate and argon/xenon (Ar/Xe) binary hydrate have been calculated using molecular dynamics using two melting techniques, namely the Z method [Belonoshko et al., Phys. Rev. B, 2006, 73, 012201] (applied for the first time to complex molecular solids) and direct phase coexistence simulations. The two methods give results in moderate agreement: calculations with the Z method give T(fus) to be 250.7 K (0.77 katm) for THF and 244.3 K (1.86 katm) for Ar/Xe hydrate respectively; the corresponding direct phase coexistence calculations give T(fus) in the range 235-240 K (0.77 katm) for THF and 240-252.5 K (1.86 katm) for Ar/Xe hydrate. The Z method was found to define the key thermodynamic states with high precision, although required long simulation times with these multicomponent molecular systems to ensure the complete melting required by the method. In contrast, the direct phase coexistence method did bracket the equilibrium temperature with little difficulty, but small thermodynamic driving forces close to phase equilibrium generated long-lived fluctuations, that obscured the precise value of phase coexistence conditions within the bracketed range.

Determining the three-phase coexistence line in methane hydrates using computer simulations

Molecular dynamics simulations have been performed to estimate the three-phase (solid hydrate-liquid water-gaseous methane) coexistence line for the water-methane binary mixture. The temperature at which the three phases are in equilibrium was determined for three different pressures, namely, 40, 100, and 400 bar by using direct coexistence simulations. In the simulations water was described by using either TIP4P, TIP4P/2005, or TIP4P/Ice models and methane was described as simple Lennard-Jones interaction site. Lorentz-Berthelot combining rules were used to obtain the parameters of the cross interactions. For the TIP4P/2005 model positive deviations from the energetic Lorentz-Berthelot rule were also considered to indirectly account for the polarization of methane when introduced in liquid water. To locate the three-phase coexistence point, two different global compositions were used, which yielded (to within statistical uncertainty) the same predictions for the three-phase coexistence temperatures, although with a somewhat different time evolution. The three-phase coexistence temperatures obtained at different pressures when using the TIP4P/Ice model of water were in agreement with the experimental results. The main reason for this is that the TIP4P/Ice model reproduces the melting point of ice I(h).

The stability of a crystal with diamond structure for patchy particles with tetrahedral symmetry

The phase diagram of model anisotropic particles with four attractive patches in a tetrahedral arrangement has been computed at two different values of the range of the potential, with the aim of investigating the conditions under which a diamond crystal can be formed. We find that the diamond phase is never stable for our longer-ranged potential. At low temperatures and pressures, the fluid freezes into a body-centered-cubic solid that can be viewed as two interpenetrating diamond lattices with a weak interaction between the two sublattices. Upon compression, an orientationally ordered face-centered-cubic crystal becomes more stable than the body-centered-cubic crystal, and at higher temperatures, a plastic face-centered-cubic phase is stabilized by the increased entropy due to orientational disorder. A similar phase diagram is found for the shorter-ranged potential, but at low temperatures and pressures, we also find a region over which the diamond phase is thermodynamically favored over the body-centered-cubic phase. The higher vibrational entropy of the diamond structure with respect to the body-centered-cubic solid explains why it is stable even though the enthalpy of the latter phase is lower. Some preliminary studies on the growth of the diamond structure starting from a crystal seed were performed. Even though the diamond phase is never thermodynamically stable for the longer-ranged model, direct coexistence simulations of the interface between the fluid and the body-centered-cubic crystal and between the fluid and the diamond crystal show that at sufficiently low pressures, it is quite probable that in both cases the solid grows into a diamond crystal, albeit involving some defects. These results highlight the importance of kinetic effects in the formation of diamond crystals in systems of patchy particles.

Melting point and phase diagram of methanol as obtained from computer simulations of the OPLS model

In this work, the melting point and the phase diagram of methanol is determined via computer simulations using the OPLS model. The three different solid structures that are found experimentally were considered. By computing the free energies of both the fluid phase and the three different solid structures (alpha, beta, gamma), the initial solid-solid and fluid-solid coexistence points were determined. By performing Gibbs-Duhem integration, the complete coexistence lines were evaluated. In this way, it was possible to compute, for the first time, the complete phase diagram for a potential model of methanol. It is found that the optimized potential model for liquid simulations (OPLS) provides reasonable predictions for the densities of the three solid polymorphs, although they tend to be somewhat low when compared with the experiment. Overall the model provides a qualitatively correct description of the phase diagram of methanol. The beta solid, which is thermodynamically stable in the experimental phase diagram of methanol, is found to be metastable in the phase diagram of the model. The alpha phase is stable at low pressures and the gamma phase is stable at high pressures, in agreement with experiment. Thus, the model is able to predict the existence of the gamma solid at high pressure. From free energy calculations we found that the melting point of the model at room pressure is 215 K. That was further confirmed by direct coexistence simulations. Thus, the model presents a melting point about 40 K above the experimental value of 175 K. Thus the OPLS model provides a reasonable description of the phase diagram of methanol, but it could probably be modified to improve the phase diagram predictions.

The melting point of ice Ih for common water models calculated from direct coexistence of the solid-liquid interface

In this work we present an implementation for the calculation of the melting point of ice I(h) from direct coexistence of the solid-liquid interface. We use molecular dynamics simulations of boxes containing liquid water and ice in contact. The implementation is based on the analysis of the evolution of the total energy along NpT simulations at different temperatures. We report the calculation of the melting point of ice I(h) at 1 bar for seven water models: SPC/E, TIP4P, TIP4P-Ew, TIP4P/ice, TIP4P/2005, TIP5P, and TIP5P-E. The results for the melting temperature from the direct coexistence simulations of this work are in agreement (within the statistical uncertainty) with those obtained previously by us from free energy calculations. By taking into account the results of this work and those of our free energy calculations, recommended values of the melting point of ice I(h) at 1 bar for the above mentioned water models are provided.