On interfacial properties of tetrahydrofuran: Atomistic and coarse-grained models from molecular dynamics simulation.

We have determined the interfacial properties of tetrahydrofuran (THF) from direct simulation of the vapor-liquid interface. The molecules are modeled using six different molecular models, three of them based on the united-atom approach and the other three based on a coarse-grained (CG) approach. In the first case, THF is modeled using the transferable parameters potential functions approach proposed by Chandrasekhar and Jorgensen [J. Chem. Phys. 77, 5073 (1982)] and a new parametrization of the TraPPE force fields for cyclic alkanes and ethers [S. J. Keasler et al., J. Phys. Chem. B 115, 11234 (2012)]. In both cases, dispersive and coulombic intermolecular interactions are explicitly taken into account. In the second case, THF is modeled as a single sphere, a diatomic molecule, and a ring formed from three Mie monomers according to the SAFT-γ Mie top-down approach [V. Papaioannou et al., J. Chem. Phys. 140, 054107 (2014)]. Simulations were performed in the molecular dynamics canonical ensemble and the vapor-liquid surface tension is evaluated from the normal and tangential components of the pressure tensor along the simulation box. In addition to the surface tension, we have also obtained density profiles, coexistence densities, critical temperature, density, and pressure, and interfacial thickness as functions of temperature, paying special attention to the comparison between the estimations obtained from different models and literature experimental data. The simulation results obtained from the three CG models as described by the SAFT-γ Mie approach are able to predict accurately the vapor-liquid phase envelope of THF, in excellent agreement with estimations obtained from TraPPE model and experimental data in the whole range of coexistence. However, Chandrasekhar and Jorgensen model presents significant deviations from experimental results. We also compare the predictions for surface tension as obtained from simulation results for all the models with experimental data. The three CG models predict reasonably well (but only qualitatively) the surface tension of THF, as a function of temperature, from the triple point to the critical temperature. On the other hand, only the TraPPE united-atoms models are able to predict accurately the experimental surface tension of the system in the whole temperature range.

Liquid-liquid interfacial properties of a symmetrical Lennard-Jones binary mixture.

We determine the interfacial properties of a symmetrical binary mixture of equal-sized spherical Lennard-Jones molecules, σ11 = σ22, with the same dispersive energy between like species, ϵ11 = ϵ22, but different dispersive energies between unlike species low enough to induce phase separation. We use the extensions of the improved version of the inhomogeneous long-range corrections of Janec̆ek [J. Phys. Chem. B 110, 6264 (2006)], presented recently by MacDowell and Blas [J. Chem. Phys. 131, 074705 (2009)] and Martínez-Ruiz et al. [J. Chem. Phys. 141, 184701 (2014)], to deal with the interaction energy and microscopic components of the pressure tensor. We perform Monte Carlo simulations in the canonical ensemble to obtain the interfacial properties of the symmetrical mixture with different cut-off distances rc and in combination with the inhomogeneous long-range corrections. The pressure tensor is obtained using the mechanical (virial) and thermodynamic route. The liquid-liquid interfacial tension is also evaluated using three different procedures, the Irving-Kirkwood method, the difference between the macroscopic components of the pressure tensor, and the test-area methodology. This allows to check the validity of the recent extensions presented to deal with the contributions due to long-range corrections for intermolecular energy and pressure tensor in the case of binary mixtures that exhibit liquid-liquid immiscibility. In addition to the pressure tensor and the surface tension, we also obtain density profiles and coexistence densities and compositions as functions of pressure, at a given temperature. According to our results, the main effect of increasing the cut-off distance rc is to sharpen the liquid-liquid interface and to increase the width of the biphasic coexistence region. Particularly interesting is the presence of a relative minimum in the total density profiles of the symmetrical mixture. This minimum is related with a desorption of the molecules at the interface, a direct consequence of a combination of the weak dispersive interactions between unlike species of the symmetrical binary mixture, and the presence of an interfacial region separating the two immiscible liquid phases in coexistence.

Effect of dispersive long-range corrections to the pressure tensor: the vapour-liquid interfacial properties of the Lennard-Jones system revisited.

We propose an extension of the improved version of the inhomogeneous long-range corrections of Janecek [J. Phys. Chem. B 110, 6264-6269 (2006)], presented recently by MacDowell and Blas [J. Chem. Phys. 131, 074705 (2009)] to account for the intermolecular potential energy of spherical, rigid, and flexible molecular systems, to deal with the contributions to the microscopic components of the pressure tensor due to the dispersive long-range corrections. We have performed Monte Carlo simulations in the canonical ensemble to obtain the interfacial properties of spherical Lennard-Jones molecules with different cutoff distances, r(c) = 2.5, 3, 4, and 5σ. In addition, we have also considered cutoff distances r(c) = 2.5 and 3σ in combination with the inhomogeneous long-range corrections proposed in this work. The normal and tangential microscopic components of the pressure tensor are obtained using the mechanical or virial route in combination with the recipe of Irving and Kirkwood, while the macroscopic components are calculated using the Volume Perturbation thermodynamic route proposed by de Miguel and Jackson [J. Chem. Phys. 125, 164109 (2006)]. The vapour-liquid interfacial tension is evaluated using three different procedures, the Irving-Kirkwood method, the difference between the macroscopic components of the pressure tensor, and the Test-Area methodology. In addition to the pressure tensor and the surface tension, we also obtain density profiles, coexistence densities, vapour pressure, critical temperature and density, and interfacial thickness as functions of temperature, paying particular attention to the effect of the cutoff distance and the long-range corrections on these properties. According to our results, the main effect of increasing the cutoff distance (at fixed temperature) is to sharpen the vapour-liquid interface, to decrease the vapour pressure, and to increase the width of the biphasic coexistence region. As a result, the interfacial thickness decreases, the width of the tangential microscopic component of the pressure tensor profile increases, and the surface tension increases as the cutoff distance is larger. We have also checked the effect of the impulsive contribution to the pressure due to the discontinuity of the intermolecular interaction potential when it is cut. If this contribution is not accounted for in the calculation of the microscopic components of the pressure tensor, incorrect values of both components as well as a wrong structure along the vapour-liquid interface are obtained.

Effect of molecular flexibility of Lennard-Jones chains on vapor-liquid interfacial properties.

We have determined the interfacial properties of short fully flexible chains formed from tangentially bonded Lennard-Jones monomeric units from direct simulation of the vapor-liquid interface. The results obtained are compared with those corresponding to rigid-linear chains formed from the same chain length, previously determined in the literature [F. J. Blas, A. I. M.-V. Bravo, J. M. Míguez, M. M. Piñeiro, and L. G. MacDowell, J. Chem. Phys. 137, 084706 (2012)]. The full long-range tails of the potential are accounted for by means of an improved version of the inhomogeneous long-range corrections of Janecek [J. Phys. Chem. B 129, 6264 (2006)] proposed recently by MacDowell and Blas [J. Chem. Phys. 131, 074705 (2008)] valid for spherical as well as for rigid and flexible molecular systems. Three different model systems comprising of 3, 5, and 6 monomers per molecule are considered. The simulations are performed in the canonical ensemble, and the vapor-liquid interfacial tension is evaluated using the test-area method. In addition to the surface tension, we also obtained density profiles, coexistence densities, critical temperature and density, and interfacial thickness as functions of temperature, paying particular attention to the effect of the chain length and rigidity on these properties. According to our results, the main effect of increasing the chain length (at fixed temperature) is to sharpen the vapor-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. Comparison between predictions for fully flexible and rigid-linear chains, formed by the same number of monomeric units, indicates that the main effects of increasing the flexibility, i.e., passing from a rigid-linear to a fully flexible chain, are: (a) to decrease the difference between the liquid and vapor densities; (b) to decrease the critical temperature and to increase the critical density; (c) to smooth the density profiles along the interfacial region; (d) to increase the interfacial thickness; and (e) to decrease the vapor-liquid surface tension.

Universal scaling behaviour of surface tension of molecular chains.

We use and extend the universal relationship recently proposed by Galliero [G. Galliero, J. Chem. Phys. 133, 074705 (2010)], based on a combination of the corresponding-states principle of Guggenheim [E. A. Guggenheim, J. Chem. Phys. 13, 253 (1945)] and the parachor approach of Macleod [J. Macleod, Trans. Faraday Soc. 19, 38 (1923)], to predict the vapour-liquid surface tension of fully flexible chainlike Lennard-Jones molecules. In the original study of Galliero, the reduced surface tension of short-chain molecules formed by up to five monomers is expressed as a unique function of the difference between the liquid and vapour coexistence densities. In this work, we extend the applicability of the recipe and demonstrate that it is also valid for predicting the surface tension of two different chainlike molecular models, namely, linear tangent chains that interact through the Lennard-Jones intermolecular potential and fully flexible chains formed by spherical segments interacting through the square-well potential. Computer simulation data for vapour-liquid surface tension of fully flexible and rigid linear Lennard-Jones, and fluid flexible square-well chains is taken from our previous works. Our results indicate that the universal scaling relationship is able to correlate short- and long-chain molecules with different degrees of flexibility and interacting through different intermolecular potentials.

On interfacial tension calculation from the test-area methodology in the grand canonical ensemble.

We propose the extension of the test-area methodology, 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 slitlike pores using the grand canonical ensemble. In order to check the adequacy of the proposed extension, we apply the method for determining the density profiles and interfacial tension of spherical molecules adsorbed in slitlike pore with different pore sizes and solid-fluid dispersive energy parameters along the same simulation. We also calculate the solid-fluid interfacial tension using the original test-area method in the canonical ensemble. Agreement between the results obtained from both methods indicate that both methods are fully equivalent. The advantage of the new methodology is that allows to calculate simultaneously the density profiles and the amount of molecules adsorbed onto a slitlike pore, as well as the solid-fluid interfacial tension. This ensures that the chemical potential at which all properties are evaluated during the simulation is exactly the same since simulations can be performed in the grand canonical ensemble, mimicking the conditions at which the adsorption experiments are most usually carried out in the laboratory.