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.

Comparison of the vapour–liquid surface tensions for substances studied in this work as obtained from SAFT-VR Mie DFT and experiments.

Determination of liquid-vapor equilibrium and critical properties of fatty acids for biodiesel production through molecular dynamics

In recent years, biodiesel production has emerged as an option for renewable and green fuel generation due to the constant reduction of fossil fuel reservoirs. Biofuels as biodiesel also show valuable attributes, environmentally speaking, due to their low environmental impact, contributing to the achievement of sustainability. However, costs are not allowable for large-scale production. Thereby, several novel processes have been proposed (e.g., reactive distillation) to solve this issue. An inconvenience for the development of these processes is the little information in the literature about the critical properties of fatty acids, which are precursors of biodiesel. Determination of critical properties for fatty acids through experimentation is difficult. The reason is that fatty acids tend to self-associate (to dimerize) due to carboxylic groups presence through hydrogen bonds, and consequently, have higher boiling points than other compounds of similar molecular mass (e.g., hydrocarbons, esters). Therefore, alternative methods for this determination are required. One choice is the group-contribution method, which is based on the structure of the molecule; however, results can significantly vary among different group-contribution approaches. Another alternative (and the focus of this research) for the determination of these properties is molecular simulation techniques. In this work, the liquid-vapor equilibrium as a function of temperature and the surface tension of three pure fatty acids of long chain (linoleic, oleic, and palmitic acid) have been calculated. Simulations have been performed by molecular dynamics using the method of direct determination of phase coexistence with the software GROMACS; in which the transferable potentials for phase equilibria united atom forcefield (TraPPE-UA) have been implemented for these specific molecules. Orthobaric densities and surface tension values have been reported at temperatures near the critical point (from 650 K to 800 K). Critical properties (temperature, pressure, density) have been extrapolated from trajectories obtained in these simulations using scaling law relations. Critical properties for these compounds are not available experimentally, therefore, group contribution calculations from the literature were used as a reference. In this comparison, the palmitic acid properties calculated in this work, show the best agreement among the three substances investigated.

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 NP_z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.

Preferential Orientations and Anomalous Interfacial Tensions in Aqueous Solutions of Alcohols

Literature studies on interfacial tension versus temperature between normal alcohols and water show that it increases with temperature and exhibits a maximum value at a given temperature depending on the molecular weight of the alcohol. This very unusual behavior is supposedly accompanied by the formation of monolayers of alcohol molecules oriented preferentially at the interface, a structural issue not confirmed until now. We use molecular-based models for water and alcohols in combination with molecular dynamics simulations to provide physical insights, from a molecular perspective, into the structural and thermodynamic behavior at the liquid-liquid interfaces of aqueous solutions of alcohols.