Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation. II. Effect of multiple occupancy.

In this work, we determine the dissociation line of the nitrogen (N2) hydrate by computer simulation using the TIP4P/Ice model for water and the TraPPE force field for N2. This work is the natural extension of Paper I, in which the dissociation temperature of the N2 hydrate has been obtained at 500, 1000, and 1500 bar [Algaba et al., J. Chem. Phys. 159, 224707 (2023)] using the solubility method and assuming single occupancy. We extend our previous study and determine the dissociation temperature of the N2 hydrate at different pressures, from 500 to 4500 bar, taking into account the single and double occupancy of the N2 molecules in the hydrate structure. We calculate the solubility of N2 in the aqueous solution as a function of temperature when it is in contact with a N2-rich liquid phase and when in contact with the hydrate phase with single and double occupancy via planar interfaces. Both curves intersect at a certain temperature that determines the dissociation temperature at a given pressure. We observe a negligible effect of occupancy on the dissociation temperature. Our findings are in very good agreement with the experimental data taken from the literature. We have also obtained the driving force for the nucleation of the hydrate as a function of temperature and occupancy at several pressures. As in the case of the dissociation line, the effect of occupancy on the driving force for nucleation is negligible. To the best of our knowledge, this is the first time that the effect of the occupancy on the driving force for nucleation of a hydrate that exhibits sII crystallographic structure is studied from computer simulation.

Simulation of the THF hydrate-water interfacial free energy from computer simulation.

In this work, the tetrahydrofuran (THF) hydrate-water interfacial free energy is determined at 500 bar, at one point of the univariant two-phase coexistence line of the THF hydrate, by molecular dynamics simulation. The mold integration-host methodology, an extension of the original mold integration technique to deal with hydrate-fluid interfaces, is used to calculate the interfacial energy. Water is described using the well-known TIP4P/Ice model, and THF is described using a rigid version of the TraPPE model. We have recently used the combination of these two models to accurately describe the univariant two-phase dissociation line of the THF hydrate in a wide range of pressures from computer simulation [Algaba et al., J. Chem. Phys. 160, 164718 (2024)]. The THF hydrate-water interfacial free energy predicted in this work is compared with the only experimental data available in the literature. The value obtained, 27(2) mJ/m2, is in excellent agreement with the experimental data taken from the literature, 24(8) mJ/m2. To the best of our knowledge, this is the first time that the THF hydrate-water interfacial free energy is predicted from computer simulation. This work confirms that the mold integration technique can be used with confidence to predict the solid-fluid interfaces of complex structures, including hydrates that exhibit sI and sII crystallographic structures.

Prediction of the univariant two-phase coexistence line of the tetrahydrofuran hydrate from computer simulation.

In this work, the univariant two-phase coexistence line of the tetrahydrofuran (THF) hydrate is determined from 100 to 1000 bar by molecular dynamics simulations. This study is carried out by putting in contact a THF hydrate phase with a stoichiometric aqueous solution phase. Following the direct coexistence technique, the pressure is fixed, and the coexistence line is determined by analyzing if the hydrate phase grows or melts at different values of temperature. Water is described using the well-known TIP4P/Ice model. We have used two different models of THF based on the transferable parameters for phase equilibria-united atom approach (TraPPE-UA), the original (flexible) TraPPe-UA model and a rigid and planar version of it. Overall, at high pressures, small differences are observed in the results obtained by both models. However, large differences are observed in the computational efforts required by the simulations performed using both models, being the rigid and planar version much faster than the original one. The effect of the unlike dispersive interactions between the water and THF molecules is also analyzed at 250 bar using the rigid and planar THF model. In particular, we modify the Berthelot combining rule via a parameter ξO-THF that controls the unlike water-THF dispersive interactions. We analyze the effect on the dissociation temperature of the hydrate when ξO-THF is modified from 1.0 (original Berthelot combining rule) to 1.4 (modified Berthelot combining rule). We use the optimized value ξO-THF = 1.4 and the rigid THF model in a transferable way to predict the dissociation temperatures at other pressures. We find excellent agreement between computer simulation predictions and experimental data taken from the literature.

Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation.

In this work, we determine the dissociation line of the nitrogen (N2) hydrate by computer simulation using the TIP4P/Ice model for water and the TraPPE force field for N2. We use the solubility method proposed recently by some of us to evaluate the dissociation temperature of the hydrate at different pressures, from 500 to 1500 bar. Particularly, we calculate the solubility of N2 in the aqueous solution when it is in contact with a N2-rich liquid phase and when in contact with the hydrate phase via planar interfaces as functions of temperature. Since the solubility of N2 decreases with temperature in the first case and increases with temperature in the second case, both curves intersect at a certain temperature that determines the dissociation temperature at a given pressure. We find a good agreement between the predictions obtained in this work and the experimental data taken from the literature in the range of pressures considered in this work. From our knowledge of the solubility curves of N2 in the aqueous solution, we also determine the driving force for nucleation of the hydrate, as a function of temperature, at different pressures. In particular, we use two different thermodynamic routes to evaluate the change in chemical potential for hydrate formation. Although the driving force for nucleation slightly decreases (in absolute value) when the pressure is increased, our results indicate that the effect of pressure can be considered negligible in the range of pressures studied in this work. To the best of our knowledge, this is the first time the driving force for nucleation of a hydrate that exhibits crystallographic structure sII, along its dissociation line, is studied from computer simulation.

Effect of pressure on the carbon dioxide hydrate-water interfacial free energy along its dissociation line.

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.

Solubility of carbon dioxide in water: Some useful results for hydrate nucleation.

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.

Molecular dynamics of liquid-liquid equilibrium and interfacial properties of aqueous solutions of methyl esters.

In this work, the liquid-liquid phase equilibria and interfacial properties of methyl ester + water binary mixtures are determined at atmospheric pressure and from 278 to 358 K combining the direct coexistence technique and molecular dynamics simulations. Methyl esters are modelled using new parametrizations based on the united atom TraPPE model force field proposed recently by us [E. Feria, J. Algaba, J. M. Míguez, A. Mejía, P. Gómez-Álvarez and F. J. Blas, Phys. Chem. Chem. Phys., 2019, 22, 4974-4983] that are able to predict the vapour-liquid interfacial properties of pure methyl esters with high accuracy. In the case of water, we consider the well-known TIP4P/2005 model, the most popular rigid and non-polarizable model to describe the interfacial properties of pure water. The simulations are performed using the direct coexistence technique in the isothermal-isobaric or NPzT ensemble in combination with molecular dynamics. We obtain density profiles, temperature-densities and temperature-composition projections of the phase diagrams, and interfacial tensions. The liquid-liquid interfacial tension is calculated from the normal and tangential components of the pressure tensor according to the mechanical virial route. We pay attention particularly to the ability of the molecular models in predicting the experimental behavior of the systems. Simulation results are able to account for the liquid-liquid phase equilibria of these binary mixtures, in good agreement with the experimental data taken from the literature. Unfortunately, experimental values for interfacial tensions are substantially overestimated by predictions from computer simulations in all cases. To our knowledge, this is the first time that the liquid-liquid phase equilibrium and interfacial properties of methyl ester + water mixtures have been predicted from computer simulations.

Simulation of the CO2 hydrate-water interfacial energy: The mold integration-guest methodology.

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.

Vapour-liquid phase equilibria and interfacial properties of fatty acid methyl esters from molecular dynamics simulations.

We have determined the phase equilibria and interfacial properties of a methyl ester homologous series (from methyl acetate to methyl heptanoate) using direct simulations of the vapour-liquid interfaces. The methyl esters are modelled using the united atom approach in combination with transferable parameters for phase equilibria (TraPPE) force fields for alkanes, alkenes, carbon dioxide, ethers, and carboxylic acids in a transferable way. This allows us to take into account explicitly both dispersive and coulombic interactions, as well as the repulsive Pauli-exclusion interactions. Simulations are performed in the NVT or canonical ensemble using molecular dynamics. Vapour-liquid surface tension is determined using the virial route, i.e., evaluating the normal and tangential components of the pressure tensor along the simulation box. We have also calculated density profiles, coexistence densities, vapour pressures, surface entropies and enthalpies, and interfacial thickness as functions of temperature, as well as the normal boiling temperatures and the critical temperatures, densities, and pressures for each member of the series. Special attention is paid to the comparison between experimental data taken from the literature and our results obtained using molecular dynamics simulations. We also analyze the effect of increasing the molecular weight of the methyl esters (at fixed temperature) on all the properties considered, with special emphasis on phase equilibria envelopes and surface tension. The TraPPE force fields transferred from other molecules and chemical families are able to predict very accurately the experimental vapour-liquid phase envelopes of methyl esters. We also compare the results obtained from simulations of the surface tension, with experimental data taken from the literature. To our knowledge, this is the first time that vapour-liquid phase equilibria and interfacial properties, and particularly surface tension, of this methyl ester homologous series are obtained using computer simulation.

Molecular modelling techniques for predicting liquid-liquid interfacial properties of methanol plus alkane (n-hexane, n-heptane, n-octane) mixtures.

In this work, the liquid-liquid interfacial properties of methanol plus n-alkane (n-hexane, n-heptane, n-octane) mixtures are investigated at atmospheric pressure by two complementary molecular modelling techniques; namely, molecular dynamic simulations (MD) and density gradient theory (DGT) coupled with the PC-SAFT (perturbed-chain statistical associating fluid theory) equation of state. Furthermore, two molecular models of methanol are used, which are based on a non-polarisable three site approach. On the one hand, is the original (flexible) TraPPE-UA model force field. On the other hand, is the rigid approximation denoted as OPLS/2016. In both cases, n-alkanes are modelled using the TraPPE-UA model. Simulations are performed using the direct coexistence technique in the ensemble. Special attention is paid to the comparison between the estimations obtained from different methanol models, the available experimental data and theoretical calculations. In all cases, the rigid model is capable of predicting the experimental phase equilibrium and interfacial properties accurately. Unsurprisingly, the methanol-rich density and interfacial tension are overestimated using the TraPPE model combined with Lorentz-Berthelot mixing rules for predicting the mixture behaviour. Accurate comparison between MD and DGT plus PC-SAFT requires consideration of the cross-interactions between individual influence parameters and fitting the ßij values. This latter aspect is particularly important because it allows the exploitation of the link between the EOS model and the direct molecular simulation of the corresponding fluid. At the same time, it was demonstrated that the key property defining the interfacial tension value is the absolute concentration of methanol in the methanol-rich phase. This behaviour indicates that there are more hydrogens bonded with each other, and they interact favourably with an increasing number of carbon atoms in the alkane.

An accurate density functional theory for the vapor-liquid interface of chain molecules based on the statistical associating fluid theory for potentials of variable range for Mie chainlike fluids.

A new Helmholtz free energy density functional is presented to predict the vapor-liquid interface of chainlike molecules. The functional is based on the last version of the statistical associating fluid theory for potentials of variable range for homogeneous Mie chainlike fluids (SAFT-VR Mie). Following the standard formalism, the density functional theory (SAFT-VR Mie DFT) is constructed using a perturbative approach in which the free energy density contains a reference term to describe all the short-range interactions treated at the local level, and a perturbative contribution to account for the attractive perturbation which incorporates the long-range dispersive interactions. In this first work, we use a mean-field version of the theory in which the pair correlations are neglected in the attractive term. The SAFT-VR Mie DFT formalism is used to examine the effect of molecular chain length and the repulsive exponent of the intermolecular potential on density profiles and surface tension of linear chains made up of up to six Mie (λr - 6) segments with different values of the repulsive exponent of the intermolecular potential. Theoretical predictions from the theory are compared directly with molecular simulation data for density profiles and surface tension of Mie chainlike molecules taken from the literature. Agreement between theory and simulation data is good for short-chain molecules under all thermodynamic conditions of coexistence considered. Once the theory has proven that it is able to predict the interfacial properties, and particularly interfacial tension, the SAFT-VR Mie DFT formalism is used to predict the interfacial behavior of two new coarse-grained models for carbon dioxide and water recently proposed in the literature. In particular, the theoretical formalism, in combination with the coarse-grained models for carbon dioxide and water, is able to predict the interfacial properties of these important substances in a reasonable way.

Vapour-liquid interfacial properties of square-well chains from density functional theory and Monte Carlo simulation.

The statistical associating fluid theory for attractive potentials of variable range (SAFT-VR) density functional theory (DFT) developed by [Gloor et al., J. Chem. Phys., 2004, 121, 12740-12759] is used to predict the interfacial behaviour of molecules modelled as fully-flexible square-well chains formed from tangentially-bonded monomers of diameter σ and potential range λ = 1.5σ. Four different model systems, comprising 4, 8, 12, and 16 monomers per molecule, are considered. In addition to that, we also compute a number of interfacial properties of molecular chains from direct simulation of the vapour-liquid interface. The simulations are performed in the canonical ensemble, and the vapour-liquid interfacial tension is evaluated using the wandering interface (WIM) method, a technique based on the thermodynamic definition of surface tension. Apart from surface tension, we also obtain density profiles, coexistence densities, vapour pressures, and critical temperature and density, paying particular attention to the effect of the chain length on these properties. According to our results, the main effect of increasing the chain length (at fixed temperature) is to sharpen the vapour-liquid interface and to increase the width of the biphasic coexistence region. As a result, the interfacial thickness decreases and the surface tension increases as the molecular chains get longer. The interfacial thickness and surface tension appear to exhibit an asymptotic limiting behaviour for long chains. A similar behaviour is also observed for the coexistence densities and critical properties. Agreement between theory and simulation results indicates that SAFT-VR DFT is only able to predict qualitatively the interfacial properties of the model. Our results are also compared with simulation data taken from the literature, including the vapour-liquid coexistence densities, vapour pressures, and surface tension.

Understanding the interfacial behavior in isopycnic Lennard-Jones mixtures by computer simulations.

The physical characterization of the singular interfacial behavior of heterogeneous fluid systems is a very important step in preliminary stages of the design process, and also in the subsequent procedures for the determination of the optimal operating conditions. Molar isopycnicity or molar density inversion is a special case of phase equilibrium behavior that directly affects the relative position of phases in heterogeneous mixtures, without being affected by gravitational fields. This work is dedicated to characterize the impact of molar density inversion on the interfacial properties of Lennard-Jones binary mixtures. The results and specific trends of the molar density inversion phenomena on the peculiar calculated composition profiles across the interface and interfacial tensions are explored by using canonical molecular dynamics simulations of the Lennard-Jones binary mixtures. Our results show that the density inversion causes drastic changes in the density profiles of the mixtures. In particular, symmetrical and equal-sized Lennard-Jones mixtures always exhibit desorption along the interfacial zone, i.e. the interfacial concentration profiles show a relative minimum at the interface of the total density profiles that increases when the dispersive energy parameter (ε(ij)) between unlike species decreases. However, as the asymmetry of the Lennard-Jones mixtures increases (σ(i) ≠ σ(j)), the concentration profiles display a relative maximum at the interface, which implies the adsorption of the total density profiles along the interfacial zone.

Density functional theory for the description of spherical non-associating monomers in confined media using the SAFT-VR equation of state and weighted density approximations.

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.

Extension of the Test-Area methodology for calculating solid-fluid interfacial tensions in cylindrical geometry.

We extend the well-known Test-Area methodology of Gloor et al. [J. Chem. Phys. 123, 134703 (2005)], originally proposed to evaluate the surface tension of planar fluid-fluid interfaces along a computer simulation in the canonical ensemble, to deal with the solid-fluid interfacial tension of systems adsorbed on cylindrical pores. The common method used to evaluate the solid-fluid interfacial tension invokes the mechanical relation in terms of the tangential and normal components of the pressure tensor relative to the interface. Unfortunately, this procedure is difficult to implement in the case of cylindrical geometry, and particularly complex in case of nonspherical molecules. Following the original work of Gloor et al., we perform free-energy perturbations due to virtual changes in the solid-fluid surface. In this particular case, the radius and length of the cylindrical pore are varied to ensure constant-volume virtual changes of the solid-fluid surface area along the simulation. We apply the modified methodology for determining the interfacial tension of a system of spherical Lennard-Jones molecules adsorbed inside cylindrical pores that interact with fluid molecules through the generalized 10-4-3 Steele potential recently proposed by Siderius and Gelb [J. Chem. Phys. 135, 084703 (2011)]. We analyze the effect of pore diameter, density of adsorbed molecules, and fluid-fluid cutoff distance of the Lennard-Jones intermolecular potential on the solid-fluid interfacial tension. This extension, as the original Test-Area formulation, offers clear advantages over the classical mechanical route of computational efficiency, easy of implementation, and generality.

Influence of the long-range corrections on the interfacial properties of molecular models using Monte Carlo simulation.

We analyze the influence of the long-range corrections, due to the dispersive term of the intermolecular potential energy, on the surface tension using direct simulation of the vapour-liquid interface of different molecular models. Although several calculation methods have been proposed recently to compute the fluid-fluid interfacial properties, the truncation of the intermolecular potential or the use of the tail corrections represents a contribution relevant from a quantitative perspective. In this work, a simplified model for methane, namely a spherical Lennard-Jones intermolecular potential, has been considered first, and afterwards other models including rigid non polarizable structures with both Lennard-Jones sites and point electric charges, representing some of the most popular models to describe water (namely the original TIP4P model, and the TIP4P/Ew and TIP4P/2005 versions), and carbon dioxide (MSM, EPM2, TraPPE, and ZD models) have been studied. Our results show that for all cases tested, including those in which the electrostatic interactions may be predominant, an incomplete account of the long-range corrections produces a systematic underestimation of the computed interfacial tension.

Monte Carlo simulation of flexible trimers: from square well chains to amphiphilic primitive models.

In this work, we present Monte Carlo computer simulation results of a primitive model of self-assembling system based on a flexible 3-mer chain interacting via square-well interactions. The effect of switching off the attractive interaction in an extreme sphere is analyzed, since the anisotropy in the molecular potential promotes self-organization. Before addressing studies on self-organization it is necessary to know the vapor liquid equilibrium of the system to avoid to confuse self-organization with phase separation. The range of the attractive potential of the model, λ, is kept constant and equal to 1.5σ, where σ is the diameter of a monomer sphere, while the attractive interaction in one of the monomers was gradually turned off until a pure hard body interaction was obtained. We present the vapor-liquid coexistence curves for the different models studied, their critical properties, and the comparison with the SAFT-VR theory prediction [A. Gil-Villegas, A. Galindo, P. J. Whitehead, S. J. Mills, G. Jackson, and A. N. Burgess, J. Chem. Phys. 106, 4168 (1997)]. Evidence of self-assembly for this system is discussed.

Semi-infinite boundary conditions for the simulation of interfaces: the Ar/CO2(s) model revisited.

We propose a method to account for the long tail corrections of dispersive forces in inhomogeneous systems. This method deals separately with the two interfaces that are usually present in a simulation setup, effectively establishing semi-infinite boundary conditions that are appropriate for the study of the interface between two infinite bulk phases. Using the wandering interface method, we calculate surface free energies of vapor-liquid, wall-liquid, and wall-vapor interfaces for a model of Lennard-Jones argon adsorbed on solid carbon dioxide. The results are employed as input to Young's equation, and the wetting temperature located. This estimate is compared with predictions from the method of effective interface potentials and good agreement is found. Our results show that truncating Ar-Ar interactions at two and a half molecular diameters results in a dramatic decrease of the wetting temperature of about 40%.

Density functional theory for the prediction of interfacial properties of molecular fluids within the SAFT-γ coarse-grained approach.

Recently, we have proposed the SAFT-VR Mie MF DFT approach [Algaba et al., Phys. Chem. Chem. Phys., 2019, 21, 11937-11948] to investigate systems that exhibit fluid-fluid interfaces. This formalism is based on the combination of the Statistical Associating Fluid Theory for attractive potentials of variable range using Mie intermolecular potential (SAFT-VR Mie) and a Density Functional Theory (DFT) treatment of the free energy. A mean-field approach is used to evaluate the attractive term, neglecting the pair correlations associated to attractions. This theory has been combined with reported SAFT-γ Coarse-Grained (CG) Mie force fields to provide an excellent description of the vapor-liquid interface of carbon dioxide and water pure fluids. The present work is a natural and necessary extension of this previous study. We assess the adequacy of the proposed methodology for dealing with inhomogeneous fluid systems of large complex molecules, in particular carbon tetrafluoride and sulfur hexafluoride greenhouse gases, the refrigerant 2,3,3,3-tetrafluoro-1-propene, and the long-chain n-decane and n-eicosane hydrocarbons. The obvious diversity of these fluids, their chemical and industrial interest, and the fact of that SAFT-γ CG Mie force fields have been reported for them justify such choice. With the aim of testing the theory, we perform Molecular Dynamics simulations in the canonical ensemble using the direct coexistence technique for the same models. We focus both on bulk, such as coexistence diagrams and vapor pressure curves, as well as interfacial properties, including surface tension. The comparison of the theoretical predictions with the computational results as well as with experimental data taken from the literature demonstrates the reliability and generalization of this method for dealing simultaneously with vapor-liquid equilibrium and interfacial phenomena. Hence, it appears as a potential tool for the interface analysis, with the main advantage over molecular simulation of low computational cost, and solving the experimental difficulties in accurately measuring the surface tension of certain systems.

Simulation of the carbon dioxide hydrate-water interfacial energy.

HYPOTHESIS: Carbon dioxide hydrates are ice-like nonstoichiometric inclusion solid compounds with importance to global climate change, and gas transportation and storage. The thermodynamic and kinetic mechanisms that control carbon dioxide nucleation critically depend on hydrate-water interfacial free energy. Only two independent indirect experiments are available in the literature. Interfacial energies show large uncertainties due to the conditions at which experiments are performed. Under these circumstances, we hypothesize that accurate molecular models for water and carbon dioxide combined with computer simulation tools can offer an alternative but complementary way to estimate interfacial energies at coexistence conditions from a molecular perspective. CALCULATIONS: We have evaluated the interfacial free energy of carbon dioxide hydrates at coexistence conditions (three-phase equilibrium or dissociation line) implementing advanced computational methodologies, including the novel Mold Integration methodology. Our calculations are based on the definition of the interfacial free energy, standard statistical thermodynamic techniques, and the use of the most reliable and used molecular models for water (TIP4P/Ice) and carbon dioxide (TraPPE) available in the literature. FINDINGS: We find that simulations provide an interfacial energy value, at coexistence conditions, consistent with the experiments from its thermodynamic definition. Our calculations are reliable since are based on the use of two molecular models that accurately predict: (1) The ice-water interfacial free energy; and (2) the dissociation line of carbon dioxide hydrates. Computer simulation predictions provide alternative but reliable estimates of the carbon dioxide interfacial energy. Our pioneering work demonstrates that is possible to predict interfacial energies of hydrates from a truly computational molecular perspective and opens a new door to the determination of free energies of hydrates.