Search

We investigate the effect of pressure on the carbon dioxide (CO2) hydrate–water interfacial free energy along its dissociation line using advanced computer simulation techniques. In previous works, we have determined the interfacial energy of the hydrate at 400 bars using the TIP4P/Ice and TraPPE molecular models for water and CO2, respectively, in combination with two different extensions of the Mold Integration technique [J. Colloid Interface Sci. 623, 354 (2022) and J. Chem. Phys. 157, 134709 (2022)]. Results obtained from computer simulation, 29(2) and 30(2) mJ/m2, are found to be in excellent agreement with the only two measurements that exist in the literature, 28(6) mJ/m2 determined by Uchida et al. [J. Phys. Chem. B 106, 8202 (2002)] and 30(3) mJ/m2 determined by Anderson et al. [J. Phys. Chem. B 107, 3507 (2002)]. Since the experiments do not allow to obtain the variation of the interfacial energy along the dissociation line of the hydrate, we extend our previous studies to quantify the effect of pressure on the interfacial energy at different pressures. Our results suggest that there exists a correlation between the interfacial free energy values and the pressure, i.e., it decreases with the pressure between 100 and 1000 bars. We expect that the combination of reliable molecular models and advanced simulation techniques could help to improve our knowledge of the thermodynamic parameters that control the interfacial free energy of hydrates from a molecular perspective. [ABSTRACT FROM AUTHOR]

The growth pattern and nucleation rate of carbon dioxide hydrate critically depend on the precise value of the hydrate–water interfacial free energy. There exist in the literature only two independent experimental measurements of this thermodynamic magnitude: one obtained by Uchida et al. [J. Phys. Chem. B 106, 8202 (2002)], 28(6) mJ/m2, and the other by Anderson and co-workers [J. Phys. Chem. B 107, 3507 (2003)], 30(3) mJ/m2. Recently, Algaba et al. [J. Colloid Interface Sci. 623, 354 (2022)] have extended the mold integration method proposed by Espinosa and co-workers [J. Chem. Phys. 141, 134709 (2014)] to deal with the CO2 hydrate–water interfacial free energy (mold integration–guest or MI-H). Computer simulations predict a value of 29(2) mJ/m2, in excellent agreement with experimental data. The method is based on the use of a mold of attractive wells located at the crystallographic positions of the oxygen atoms of water molecules in equilibrium hydrate structures to induce the formation of a thin hydrate slab in the liquid phase at coexistence conditions. We propose here a new implementation of the mold integration technique using a mold of attractive wells located now at the crystallographic positions of the carbon atoms of the CO2 molecules in the equilibrium hydrate structure. We find that the new mold integration–guest methodology, which does not introduce positional or orientational information of the water molecules in the hydrate phase, is able to induce the formation of CO2 hydrates in an efficient way. More importantly, this new version of the method predicts a CO2 hydrate–water interfacial energy value of 30(2) mJ/m2, in excellent agreement with experimental data, which is also fully consistent with the results obtained using the previous methodology. [ABSTRACT FROM AUTHOR]

As a first step of an ongoing study of thermodynamic properties and adsorption of complex fluids in confined media, we present a new theoretical description for spherical monomers using the Statistical Associating Fluid Theory for potential of Variable Range (SAFT-VR) and a Non-Local Density Functional Theory (NLDFT) with Weighted Density Approximations (WDA). The well-known Modified Fundamental Measure Theory is used to describe the inhomogeneous hard-sphere contribution as a reference for the monomer and two WDA approaches are developed for the dispersive terms from the high-temperature Barker and Henderson perturbation expansion. The first approach extends the dispersive contributions using the scalar and vector weighted densities introduced in the Fundamental Measure Theory (FMT) and the second one uses a coarse-grained (CG) approach with a unique weighted density. To test the accuracy of this new NLDFT/SAFT-VR coupling, the two versions of the theoretical model are compared with Grand Canonical Monte Carlo (GCMC) molecular simulations using the same molecular model. Only the version with the "CG" approach for the dispersive terms provides results in excellent agreement with GCMC calculations in a wide range of conditions while the "FMT" extension version gives a good representation solely at low pressures. Hence, the "CG" version of the theoretical model is used to reproduce methane adsorption isotherms in a Carbon Molecular Sieve and compared with experimental data after a characterization of the material. The whole results show an excellent agreement between modeling and experiments. Thus, through a complete and consistent comparison both with molecular simulations and with experimental data, the NLDFT/SAFT-VR theory has been validated for the description of monomers. [ABSTRACT FROM AUTHOR]

In a first part, interfacial properties of a pure monoatomic fluid interacting through the Mie n-6 potential (n=8, 10, 12, and 20) have been studied using extensive molecular simulations. Monte Carlo and molecular dynamics simulations have been employed, using, respectively, the test area approach and the mechanic route. In order to yield reference values, simulations have been performed with a cutoff radius equal to 10σ, which is shown to be sufficient to avoid long range corrections. It is shown that both approaches provide results consistent with each other. Using the molecular simulations results, it is demonstrated that a unique scaling law is able to provide an accurate estimation of the surface tension whatever the repulsive exponent n, even far from the critical point. Furthermore, it is shown that the surface tension of the Mie n-6 fluid is as well accurately described by a unique Parachor’s law. Density profiles are shown to be well represented by the tanh mean field profile, with slight deviations for the lowest temperatures and the smallest n. In addition, the interfacial width is shown to increase when n decreases (for a given reduced temperature) and to follow the usual scaling behavior for not too low temperature. In a second part, interfacial properties of the Mie n-6 fluid computed by the gradient theory, coupled with an equation of state based on the Barker–Henderson perturbation theory, have been compared with those obtained by molecular simulations. It is demonstrated that, even far from the critical point, the gradient theory is efficient to compute surface tensions and density profiles of this model fluid, provided the equation of state accurately model the phase behavior of the fluid involved (which is not the case for n=8 in this study). [ABSTRACT FROM AUTHOR]

The growth pattern and nucleation rate of carbon dioxide hydrate critically depend on the precise value of the hydrate-water interfacial free energy. There exist in the literature only two independent experimental measurements of this thermodynamic magnitude: one obtained by Uchida et al. [J. Phys. Chem. B 106, 8202 (2002)], 28(6) mJ/m 2 , and the other by Anderson and co-workers [J. Phys. Chem. B 107, 3507 (2003)], 30(3) mJ/m 2 . Recently, Algaba et al. [J. Colloid Interface Sci. 623, 354 (2022)] have extended the mold integration method proposed by Espinosa and co-workers [J. Chem. Phys. 141, 134709 (2014)] to deal with the CO 2 hydrate-water interfacial free energy (mold integration-guest or MI-H). Computer simulations predict a value of 29(2) mJ/m 2 , in excellent agreement with experimental data. The method is based on the use of a mold of attractive wells located at the crystallographic positions of the oxygen atoms of water molecules in equilibrium hydrate structures to induce the formation of a thin hydrate slab in the liquid phase at coexistence conditions. We propose here a new implementation of the mold integration technique using a mold of attractive wells located now at the crystallographic positions of the carbon atoms of the CO 2 molecules in the equilibrium hydrate structure. We find that the new mold integration-guest methodology, which does not introduce positional or orientational information of the water molecules in the hydrate phase, is able to induce the formation of CO 2 hydrates in an efficient way. More importantly, this new version of the method predicts a CO 2 hydrate-water interfacial energy value of 30(2) mJ/m 2 , in excellent agreement with experimental data, which is also fully consistent with the results obtained using the previous methodology.