Determination of surface tension in binary mixtures using transition-matrix Monte Carlo

We present a methodology based on grand-canonical transition-matrix Monte Carlo and finite-size scaling analysis to calculate surface tensions in binary mixtures. In particular, mixture transition-matrix Monte Carlo is first used to calculate apparent, system-size-dependent free-energy barriers separating coexisting fluid phases. Finite-size scaling is then used to extrapolate these values to the infinitely large system limit to determine the true thermodynamic surface tension. A key distinction of the methodology is that it yields the entire isothermal surface-tension curve for a binary mixture in a relatively small number of simulations. We demonstrate the utility of the method by calculating surface-tension curves for three binary Lennard-Jones mixtures. While we have only examined the surface tension of simple fluids in this work, the method is general and can be extended to molecular fluids as well as to determine interfacial tensions of liquid-liquid interfaces.     

Computation of surface tensions using expanded ensemble simulations

A method for the direct simulation of the surface tension is examined. The technique is based on the thermodynamic route to the interfacial tension and makes use of the expanded ensemble simulation method for the calculation of the free energy difference between two inhomogeneous systems with the same number of particles, temperature, and volume, but different interfacial area. The method is completely general and suitable for systems with either continuous or discontinuous interactions. The adequacy of the expanded ensemble method is assessed by computing the interfacial tension of the planar vapor-liquid interface of Lennard-Jones, Lennard-Jones dimers, Gay-Berne, and square-well model fluids; in the latter, the interactions are discontinuous and the present method does not exhibit the asymmetry of other related methods, such as the test area. The expanded ensemble simulation results are compared with simulation data obtained from other techniques (mechanical and test area) with overall good agreement.     

Vapor-liquid interfacial properties of rigid-linear Lennard-Jones chains

We have obtained the interfacial properties of short rigid-linear chains formed from tangentially bonded Lennard-Jones monomeric units from direct simulation of the vapour-liquid interface. The full long-range tails of the potential are accounted for by means of an improved version of the inhomogeneous long-range corrections of Janec?ek [J. Phys. Chem. B 110, 6264-6269 (2006)] proposed recently by MacDowell and Blas [J. Chem. Phys. 131, 074705 (2009)] valid for spherical as well as for rigid and flexible molecular systems. Three different model systems comprising of 3, 4, and 5 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 obtain 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. The surface tension has been scaled by critical properties and represented as a function of the difference between coexistence densities relative to the critical density.     

Calculation of surface tension via area sampling

We examine the performance of several molecular simulation techniques aimed at evaluation of the surface tension through its thermodynamic definition. For all methods explored, the surface tension is calculated by approximating the change in Helmholtz free energy associated with a change in interfacial area through simulation of a liquid slab at constant particle number, volume, and temperature. The methods explored fall within three general classes: free-energy perturbation, the Bennett acceptance-ratio scheme, and the expanded ensemble technique. Calculations are performed for both the truncated Lennard-Jones and square-well fluids at select temperatures spaced along their respective liquid-vapor saturation lines. Overall, we find that Bennett and expanded ensemble approaches provide the best combination of accuracy and precision. All of the methods, when applied using sufficiently small area perturbation, generate equivalent results for the Lennard-Jones fluid. However, single-stage free-energy-perturbation methods and the closely related test-area technique recently introduced by Gloor et al. [J. Chem. Phys. 123, 134703 (2005)] generate surface tension values for the square-well fluid that are not consistent with those obtained from the more robust expanded ensemble and Bennett approaches, regardless of the size of the area perturbation. Single-stage perturbation methods fail also for the Lennard-Jones system when applied using large area perturbations. Here an analysis of phase-space overlap produces a quantitative explanation of the observed inaccuracy and shows that the satisfactory results obtained in these cases from the test-area method arise from a cancellation of errors that cannot be expected in general. We also briefly analyze the variation in method performance with respect to the adjustable parameters inherent to the techniques.     

Calculation of the surface tension and pressure components from a non-exponential perturbation method of the thermodynamic route

Surface tension is probably the most important interfacial property and a large number of techniques have been devoted to its calculation. Usually, this calculation is carried out using mechanical or thermodynamic definitions. The mechanical route uses an arbitrary choice to affect the contribution of the pairwise force. To overcome this arbitrariness, a thermodynamic route based on the area perturbation (test-area (TA) method) has been developed for the calculation of surface tension. The volume perturbation (VP) method provides an original route to compute the components of the pressure tensor. These two routes are developed from the perturbation theory leading to working expressions using exponential averages of energy. The use of exponential averages makes the calculation strongly dependent on the occurrence of low values of ΔU. Additionally, the decomposition of the energy to obtain local surface tension is nontrivial. From the explicit derivation of the partition function the exponential average is avoided providing an interesting alternative to TA, VP, and mechanical methods. To make a consistent comparison, we study the profiles of the surface tension along the direction normal to the surface for the different definitions and techniques in the cases of liquid-vapor interfaces of acids gases, binary, and apolar systems.     

Vapor-liquid surface tension of strong short-range Yukawa fluid

The thermodynamic properties of strong short-range attractive Yukawa fluids, κ = 10, 9, 8, and 7, are determined by combining the slab technique with the standard and the replica exchange Monte Carlo (REMC) methods. A good agreement was found among the coexistence curves of these systems calculated by REMC and those previously reported in the literature. However, REMC allows exploring the coexistence at lower temperatures, where dynamics turns glassy. To obtain the surface tension we employed, for both methods, a procedure that yields the pressure tensor components for discontinuous potentials. The surface tension results obtained by the standard MC and REMC techniques are in good agreement.     

Surface tension and vapor-liquid phase coexistence of confined square-well fluid

Phase equilibria of a square-well fluid in planar slit pores with varying slit width are investigated by applying the grand-canonical transition-matrix Monte Carlo (GC-TMMC) with the histogram-reweighting method. The wall-fluid interaction strength was varied from repulsive to attractive such that it is greater than the fluid-fluid interaction strength. The nature of the phase coexistence envelope is in agreement with that given in literature. The surface tension of the vapor-liquid interface is calculated via molecular dynamics simulations. GC-TMMC with finite size scaling is also used to calculate the surface tension. The results from molecular dynamics and GC-TMMC methods are in very good mutual agreement. The vapor-liquid surface tension, under confinement, was found to be lower than the bulk surface tension. However, with the increase of the slit width the surface tension increases. For the case of a square-well fluid in an attractive planar slit pore, the vapor-liquid surface tension exhibits a maximum with respect to wall-fluid interaction energy. We also report estimates of critical properties of confined fluids via the rectilinear diameter approach.     

Molecular simulation study of the effect of pressure on the vapor-liquid interface of the square-well fluid

We examine a model system to study the effect of pressure on the surface tension of a vapor-liquid interface. The system is a two-component mixture of spheres interacting with the square-well (A-A) and hard-sphere (B-B) potentials and with unlike (A-B) interactions ranging (for different cases) from hard sphere to strongly attractive square well. The bulk-phase and interfacial properties are measured by molecular dynamics simulation for coexisting vapor-liquid phases for various mixture compositions, pressures, and temperatures. The variation of the surface tension with pressure compares well to values given by surface-excess formulas derived from thermodynamic considerations. We find that surface tension increases with pressure only for the case of an inert solute (hard-sphere A-B interactions) and that the presence of A-B attractions strongly promotes a decrease of surface tension with pressure. An examination of density and composition profiles is made to explain these effects in terms of surface-adsorption arguments.     

Interfacial properties of cyclic hydrocarbons: a Monte Carlo study

The Monte Carlo technique is used to study the vapor-liquid interface of cyclopentane, cyclohexane, and benzene. The OPLS and TraPPE potential fields are compared in the temperature range from 298.15 to 348.15 K (273.15-298.15 K for C5H10). A new method for the treatment of the long-range interactions in inhomogeneous simulations is used. When this new method is employed, the obtained values of saturated liquid density and of enthalpy of vaporization are equal to those obtained using the bulk isothermal-isobaric Monte Carlo technique. The values of surface tension become independent of the cutoff distance and they are significantly larger than those when only simple spherical truncation of intermolecular interactions is used.     

Monte Carlo simulation of equilibrium reactions at vapor-liquid interfaces

Chemical reactions are known to behave differently, depending upon their local environment. While the interactions with neighboring molecules may alter both the kinetics of chemical reactions and the overall equilibrium conversion, we have performed simulations of the latter. The particular environment that we address is the vapor-liquid interface, since only a few, limited studies have explored the influence of an interface on equilibrium reaction behavior. Simple dimerization reactions are modeled, as well as more complex multicomponent reactions, using the reactive Monte Carlo (RxMC) simulation technique. We find that the conversion of a reaction can be markedly different at an interface as compared to the bulk vapor and liquid phases, and these trends are analyzed with respect to specific intermolecular interactions. In conjunction, we calculate the surface tension of the reacting fluids at the interface, which is found to have unusual scaling behavior, with respect to the system temperature.     

On the long-range corrections to computer simulation results for the Lennard-Jones vapor-liquid interface

The long-range corrections (LRCs) to the configurational energy have been taken into consideration in the Monte Carlo simulation of the vapor-liquid interface for a pure Lennard-Jones (LJ) fluid. The simulated bulk densities agree satisfactorily with those obtained from the Gibbs ensemble method, and the simulated surface tension values agree reasonably well with those reported in the literature for a larger number of molecules and a larger cut-off distance. To compare the influence of the potential forms on the simulation results, a truncated LJ potential, and a shifted and truncated LJ potential have been examined. Although the bulk densities and surface tensions calculated for different model fluids are strongly affected by the LRC, the different potentials essentially lead to similar density values and similar surface tension values when the respective calculated values are compared on the basis of a reduced temperature scale.     

The nature of the calculation of the pressure in molecular simulations of continuous models from volume perturbations

We consider some fundamental aspects of the calculation of the pressure from simulations by performing volume perturbations. The method, initially proposed for hard-core potentials by Eppenga and Frenkel [Mol. Phys.52, 1303 (1984)] and then extended to continuous potentials by Harismiadis et al. [J. Chem. Phys. 105, 8469 (1996)], is based on the numerical estimate of the change in Helmholtz free energy associated with the perturbation which, in turn, can be expressed as an ensemble average of the corresponding Boltzmann factor. The approach can be easily generalized to the calculation of components of the pressure tensor and also to ensembles other than the canonical ensemble. The accuracy of the method is assessed by comparing simulation results obtained from the volume-perturbation route with those obtained from the usual virial expression for several prototype fluid models. Monte Carlo simulation data are reported for bulk fluids and for inhomogeneous systems containing a vapor-liquid interface.     

A general perturbation approach for equation of state development: applications to simple fluids, ab initio potentials, and fullerenes

A new perturbation scheme based on the Barker-Henderson perturbation theory [J. Chem. Phys. 47, 4714 (1967)] is proposed to predict the thermodynamic properties of spherical molecules. Accurate predictions of second virial coefficients and vapor-liquid coexistence properties are obtained for a large variety of potential functions (square well, Yukawa, Sutherland, Lennard-Jones, Buckingham, Girifalco). New Gibbs ensemble Monte Carlo simulations of the generalized exp-m Buckingham potential are reported. An extension of the perturbation approach to mixtures is proposed, and excellent predictions of vapor-liquid equilibria are obtained for Lennard-Jones mixtures. The perturbation scheme can be applied to complex potential functions fitted to ab initio data to predict the properties of real molecules such as neon. The new approach can also be used as an auxiliary tool in molecular simulation studies, to efficiently optimize an intermolecular potential on macroscopic properties or match force fields based on different potential functions.     

Efficient prediction of thermodynamic properties of quadrupolar fluids from simulation of a coarse-grained model: the case of carbon dioxide

Monte Carlo simulations are presented for a coarse-grained model of real quadrupolar fluids. Molecules are represented by particles interacting with Lennard-Jones forces plus the thermally averaged quadrupole-quadrupole interaction. The properties discussed include the vapor-liquid coexistence curve, the vapor pressure along coexistence, and the surface tension. The full isotherms are also accessible over a wide range of temperatures and densities. It is shown that the critical parameters (critical temperature, density, and pressure) depend almost linearly on a quadrupolar parameter q=Q(*4)T*, where Q* is the reduced quadrupole moment of the molecule and T* the reduced temperature. The model can be applied to a variety of small quadrupolar molecules. We focus on carbon dioxide as a test case, but consider nitrogen and benzene, too. Experimental critical temperature, density, and quadrupolar moment are sufficient to fix the parameters of the model. The resulting agreement with experiments is excellent and marks a significant improvement over approaches which neglect quadrupolar effects. The same coarse-grained model was also applied in the framework of perturbation theory in the mean spherical approximation. As expected, the latter deviates from the Monte Carlo results in the critical region, but is reasonably accurate at lower temperatures.     

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.     

Adsorption of n-alkane vapours at the water surface

Monte Carlo simulations are reported here to predict the surface tension of the liquid-vapour interface of water upon adsorption of alkane vapours (methane to hexane). A decrease of the surface tension has been established from n-pentane. A correlation has been evidenced between the decrease of the surface tension and the absence of specific arrangement at the water surface for n-pentane and n-hexane. The thermodynamic stability of the adsorption layer and the absence of film for longer alkanes have been checked through the calculation of a potential of mean force. This complements the work recently published [Ghoufi et al., Phys. Chem. Chem. Phys., 2010, 12, 5203] concerning the adsorption of methane at the water surface. The decrease of the surface tension has been interpreted in terms of the degree of hydrogen bonding of water molecules at the liquid-vapour interface upon adsorption.     

Computation of interfacial properties via grand canonical transition matrix Monte Carlo simulation

We examine two free-energy-based methods for studying the wetting properties of a fluid in contact with a solid substrate. Application of the first approach involves examination of the adsorption behavior of a fluid at a single substrate, while the second technique requires investigation of the properties of a system confined between two parallel substrates. Both of the techniques rely upon computation and analysis of the density dependence of a system's surface free energy and provide the contact angle and solid-vapor and solid-liquid interfacial tensions for substrate-fluid combinations within the partial wetting regime. Grand canonical transition matrix Monte Carlo simulation is used to obtain the required free-energy curves. The methods examined within this work are general and are applicable to a wide range of molecular systems. We probe the performance of the methods by computing the interfacial properties for two systems in which an atomistic fluid interacts with a fcc crystal. For both of the systems studied we find good agreement between our results and those obtained via the mechanical definition of the interfacial tension.     

Vapor-liquid interfacial properties of mutually saturated water/1-butanol solutions

Adsorption and ordering at the vapor-liquid interfaces of mutually saturated water/1-butanol solutions at a temperature of 298.15 K were investigated using configurational-bias Monte Carlo simulations in the Gibbs ensemble and compared to the surface properties of neat water and 1-butanol liquids. A dense 1-butanol monolayer is observed at the surface of the water-rich phase, which results in a substantial decrease of its surface tension. In contrast, there is no enrichment of water molecules at the surface of the butanol-rich phase, and its surface tension is not significantly changed. Analysis of the interfacial structures reveals that these systems exhibit orientational ordering and composition heterogeneity. Analysis of the hydrogen-bonding distributions suggests that the formation of the 1-butanol monolayer is driven by an excellent match between water and the primary alcohol; that is, additional hydrogen bonds are formed between the excess free hydrogens of surface water and the excess hydrogen-bond acceptor sites of 1-butanol.