Sulfur hexafluoride's liquid-vapor coexistence curve, interfacial properties, and diffusion coefficients as predicted by a simple rigid model

We present here molecular-dynamics simulation results of the vapor-liquid coexistence curve, surface tension, and self-diffusion coefficients of sulfur hexafluoride. Sulfur hexafluoride is modeled as a rigid molecule, following the model proposed by Pawley [Mol. Phys. 43, 1321 (1981)]. Vapor-liquid coexistence curve and surface tension are obtained through direct molecular-dynamic simulations in the NVT ensemble. Simulation results are able to reproduce the qualitative shape of the vapor-liquid envelope. However, lower densities, a higher critical temperature, and an overestimated surface tension are obtained here. Those deviations are explained on the basis of the rigidity of the molecular model used. Self-diffusion coefficients are calculated from simulations in the NVE ensemble for different gas states at atmospheric pressure. The rigid model performs better for dynamical properties since simulation results provide very good agreement with available experimental data in this case.     

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.     

Surface tension of water and acid gases from Monte Carlo simulations

We report direct Monte Carlo (MC) simulations on the liquid-vapor interfaces of pure water, carbon dioxide, and hydrogen sulfide. In the case of water, the recent TIP4P/2005 potential model used with the MC method is shown to reproduce the experimental surface tension and to accurately describe the coexistence curves. The agreement with experiments is also excellent for CO(2) and H(2)S with standard nonpolarizable models. The surface tensions are calculated by using the mechanical and the thermodynamic definitions via profiles along the direction normal to the surface. We also discuss the different contributions to the surface tension due to the repulsion-dispersion and electrostatic interactions. The different profiles of these contributions are proposed in the case of water.     

Surface tension of associating fluids by Monte Carlo simulations

Canonical Monte Carlo (NVT-MC) simulations were performed to obtain surface tension and coexistence densities at the liquid-vapor interface of one-site associating Lennard-Jones and hard-core Yukawa fluids, as functions of association strength and temperature. The method to obtain the components of the pressure tensor from NVT-MC simulations was validated by comparing the equation of state of the associative hard sphere system with that coming from isothermal-isobaric Monte Carlo simulations. Surface tension of the associative Lennard-Jones fluid determined from NVT-MC is compared with previously reported results obtained by molecular dynamics simulations of a pseudomixture model of monomers and dimers. A good agreement was found between both methods. Values of surface tension of associative hard-core Yukawa fluids are presented here for the first time.     

Calculation of phase coexistence properties and surface tensions of n-alkanes with grand-canonical transition-matrix monte carlo simulation and finite-size scaling

Grand-canonical transition-matrix Monte Carlo is combined with configurational-bias and expanded ensemble Monte Carlo techniques to obtain saturated densities and vapor pressures of select n-alkanes. Surface tension values for butane, hexane, and octane are also computed via the finite-size scaling method of Binder. The exponential-6 model of Errington and Panagiotopoulos is used to describe the molecular interactions. The effect of the number of configurational-bias trial conformations on the efficiency of phase equilibra calculations is studied. We find that a broad range of trial conformation numbers give reasonable performance, with the optimal value increasing with decreasing temperature for a fixed chain length. Phase coexistence properties are in good agreement with literature values and are obtained with very reasonable computing resources. Similar to other recently developed n-alkane force fields, the exponential-6 model overestimates the surface tension relative to experimental values. Statistical uncertainties for coexistence properties obtained with the current approach are relatively small compared to existing methods.     

Toward a density functional description of liquid pH2

A finite-temperature density functional approach to describe the properties of parahydrogen in the liquid-vapor coexistence region is presented. The first proposed functional is zero-range, where the density-gradient term is adjusted so as to reproduce the surface tension of the liquid-vapor interface at low temperature. The second functional is finite-range and, while it is fitted to reproduce bulk pH(2) properties only, it is shown to yield surface properties in good agreement with experiments. These functionals are used to study the surface thickness of the liquid-vapor interface, the wetting transition of parahydrogen on a planar Rb model surface, and homogeneous cavitation in bulk liquid pH(2).     

Liquid-vapor interface of square-well fluids of variable interaction range

Properties of the liquid-vapor interface of square-well fluids with ranges of interaction lambda=1.5, 2.0, and 3.0 are obtained by Monte Carlo simulations and from square-gradient theories that combine the Carnahan-Starling equation of state for hard spheres with the second and third virial coefficients. The predicted surface tensions show good agreement with the simulation results for lambda=2 and for lambda=3 in a temperature range reasonably close to the critical point, 0.8相似文献   

Surface tension curvature correction for fluid–vapour interface by thermodynamics of curved boundary layers

Thermodynamics of curved boundary layers is used to study surface properties of Lennard–Jones fluids and an equation is obtained for the surface tension of the spherical interfaces. The work of cavity formation inside a fluid is utilised to calculate the surface tension. The scaled particle theory (SPT) equation of Pierotti and Stillinger are properly modified. The calculations by the modified Stillinger equation predicts the surface tension in any radius of curvature and the results are in good agreement with the experimental data. But for Pierotti equation, due to low sensitivity to temperature, good agreement obtained only at low temperatures.     

Isotherms of phospholipid monolayers measured by a pendant drop technique

The Axisymmetric Drop Shape Analysis (ADSA) has been used to study the surface pressure/area isotherms of insoluble surfactant monolayers. The continuous measurement of surface tension as a function of surface area by increasing and decreasing the drop volume allows to investigate the phase transitions in monolayers. The isotherms of two phospholipids, dipalmitoyl phosphatidyl choline (DPPC) and dimyristoyl phosphatidyl ethanolamine (DMPE), show good agreement with those measured by using a conventional Langmuir-Blodgett film balance, except in the coexistence region. The observed disagreements are discussed in terms of differences in compression rate, curvature of the surface and effect of impurities. Evidence of possible geometric effects on monolayer domain formation and growth is given on the basis of BAM images.Due to the small total surface area, the ADSA technique provides advantages as regards homogeneity of temperature, surface pressure, surface concentration and the symmetry of area changes.     

Thin-thick surface phase coexistence and boundary tension of the square-well fluid on a weak attractive surface

Prewetting transition is studied for the square-well fluid of attractive-well diameter lambda(ff)sigma(ff)=1.5 in the presence of a homogeneous surface modeled by the square-well potential of attractive well from 0.8sigma(ff) to 1.8sigma(ff). We investigate surface phase coexistence of thin-thick film transition using grand-canonical transition matrix Monte Carlo (GC-TMMC) and histogram reweighting techniques. Molecular dynamics (MD) and GC-TMMC are utilized to predict the properties of the fluid for various surface fluid affinities. Occurrences of prewetting transition with the variation of surface affinity are observed for a domain of reduced temperature from T(*)=0.62 to 0.75. We have used MD and GC-TMMC+finite size scaling (FSS) simulations to calculate the boundary tension as a function of temperature as well as surface affinity. Boundary tensions via MD and GC-TMMC+FSS methods are in good agreement. The boundary tension increases with the decrease of wall-fluid affinity. Prewetting critical properties are calculated using rectilinear diameter approach and scaling analysis. We found that critical temperature and density increase with the decrease of wall-fluid affinity.     

Density functional theory of size-dependent surface tension of Lennard-Jones fluid droplets using a double well type Helmholtz free energy functional

A double well type Helmholtz free energy density functional and a model density profile for a two phase vapor-liquid system are used to obtain the size-dependent interfacial properties of the vapor-liquid interface at coexistence condition along the lines of van der Waals and Cahn and Hilliard density functional formalism of the interface. The surface tension, temperature-density curve, density profile, and thickness of the interface of Lennard-Jones fluid droplet-vapor equilibrium, as predicted in this work are reported. The planar interfacial properties, obtained from consideration of large radius of the liquid drop, are in good agreement with the results of other earlier theories and experiments. The same free energy model has been tested by solving the equations numerically, and the results compare well with those from the use of model density profile.     

Comparative study of the effect of tail corrections on surface tension determined by molecular simulation

We report results from a comparative study of the influence of tail corrections on the surface tension of the Lennard-Jones fluid. We find that cutoff-independent surface tensions can be obtained by applying a set of tail corrections recently introduced by Janecek at each step of an interfacial Monte Carlo (MC) or molecular dynamics (MD) simulation. The effect of tail corrections on an alternative methodology for calculating surface tension, the combination of finite-size scaling and grand-canonical transition-matrix Monte Carlo (FSS/GC-TMMC), was also investigated. Using this indirect method, surface tensions were calculated with standard (bulk-fluid) tail corrections and lattice sums, the latter usually considered more accurate but computationally more intensive than the former. With standard tail corrections, we find that the surface tension decreases with increasing cutoff distance, reaching a limiting value corresponding to the maximum cutoff possible, namely half the simulation box length. In contrast, surface tension values obtained with the lattice summation were cutoff-independent. More importantly, these values were equivalent to those surface tension values obtained using standard tail corrections and a cutoff distance of half the box length. We also find that the surface tension values obtained here are in agreement with those found in the literature. Last, we find that surface tension values obtained by MD and FSS/GC-TMMC are in decent agreement so long as the appropriate tail correction schemes are used, and that the relative uncertainties in the surface tensions calculated by MD are generally an order of magnitude greater than those calculated by FSS/GC-TMMC. However, the time required by MD on a single central processing unit is less than that required by FSS/GC-TMMC.     

Thermodynamic properties of short-range square well fluid

The interfacial properties of short-range square well fluid with lambda=1.15, 1.25, and 1.375 were determined by using single canonical Monte Carlo simulations. Simulations were carried out in the vapor-liquid region. The coexistence curves of these models were calculated and compared to those previously reported in the literature and good agreement was found among them. We found that the surface tension curves for any potential model of short range form a single master curve when we plot gamma* vs TT(c). It is demonstrated that the critical reduced second virial coefficient B(2)* as a function of interaction range or T(c)* is not constant.     

The Wolf method applied to the liquid-vapor interface of water

The Wolf method for the calculation of electrostatic interactions is applied in a liquid phase and at the liquid-vapor interface of water and its results are compared with those from the Ewald sums method. Molecular dynamics simulations are performed to calculate the radial distribution functions at room temperature. The interface simulations are used to obtain the coexisting densities and surface tension along the coexistence curve. The water model is a flexible version of the extended simple point charge model. The Wolf method gives good structural results, fair coexistence densities, and poor surface tensions as compared with those obtained using the Ewald sums method.     

The Prediction of Surface Tension and Thermodynamic Analysis of the Surface in Mixtures of Cryogenic Liquids

The surface tensions of 42 binary cryogenic mixtures at low temperature are correlated using the Shereshefsky model and excellent results are obtained. The average percent deviation is about ～ 1.08%. The Gibbs energy change in the surface region is calculated and is used to obtain the excess number of molecular layers in the surface region. Furthermore, the model is used to derive an equation for the standard Gibbs energy of adsorption. The experimental standard Gibbs energy of adsorption is obtained from surface tension data and compared with calculated data. The agreement between experimental and calculated data is found to be very good. The magnitude of the Gibbs energy change in the surface region and the standard Gibbs energy of adsorption are discussed in terms of nature and type of intermolecular interactions in binary mixtures.     

Surface tension of quantum fluids from molecular simulations

We present the first molecular simulations of the vapor-liquid surface tension of quantum liquids. The path integral formalism of Feynman was used to account for the quantum mechanical behavior of both the liquid and the vapor. A replica-data parallel algorithm was implemented to achieve good parallel performance of the simulation code on at least 32 processors. We have computed the surface tension and the vapor-liquid phase diagram of pure hydrogen over the temperature range 18-30 K and pure deuterium from 19 to 34 K. The simulation results for surface tension and vapor-liquid orthobaric densities are in very good agreement with experimental data. We have computed the interfacial properties of hydrogen-deuterium mixtures over the entire concentration range at 20.4 and 24 K. The calculated equilibrium compositions of the mixtures are in excellent agreement with experimental data. The computed mixture surface tension shows negative deviations from ideal solution behavior, in agreement with experimental data and predictions from Prigogine's theory. The magnitude of the deviations at 20.4 K are substantially larger from simulations and from theory than from experiments. We conclude that the experimentally measured mixture surface tension values are systematically too high. Analysis of the concentration profiles in the interfacial region shows that the nonideal behavior can be described entirely by segregation of H(2) to the interface, indicating that H(2) acts as a surfactant in H(2)-D(2) mixtures.     

Surface tension of the Widom-Rowlinson model

We consider the computation of the surface tension of the fluid-fluid interface for the Widom-Rowlinson [J. Chem. Phys. 52, 1670 (1970)] binary mixture from direct simulation of the inhomogeneous system. We make use of the standard mechanical route, in which the surface tension follows from the computation of the normal and tangential components of the pressure tensor of the system. In addition to the usual approach, which involves simulations of the inhomogeneous system in the canonical ensemble, we also consider the computation of the surface tension in an ensemble where the pressure perpendicular (normal) to the planar interface is kept fixed. Both approaches are seen to provide consistent values of the interfacial tension. The issue of the system-size dependence of the surface tension is addressed. In addition, simulations of the fluid-fluid coexistence properties of the mixture are performed in the semigrand canonical ensemble. Our results are compared with existing data of the Widom-Rowlinson mixture and are also examined in the light of the vapor-liquid equilibrium of the thermodynamically equivalent one-component penetrable sphere model.     

Effect of temperature and pressure on surface tension of polystyrene in supercritical carbon dioxide

The surface tension of polymers in a supercritical fluid is one of the most important physicochemical parameters in many engineering processes, such as microcellular foaming where the surface tension between a polymer melt and a fluid is a principal factor in determining cell nucleation and growth. This paper presents experimental results of the surface tension of polystyrene in supercritical carbon dioxide, together with theoretical calculations for a corresponding system. The surface tension is determined by Axisymmetric Drop Shape Analysis-Profile (ADSA-P), where a high pressure and temperature cell is designed and constructed to facilitate the formation of a pendent drop of polystyrene melt. Self-consistent field theory (SCFT) calculations are applied to simulate the surface tension of a corresponding system, and good qualitative agreement with experiment is obtained. The physical mechanisms for three main experimental trends are explained by using SCFT, and none of the explanations quantitatively depend on the configurational entropy of the polymer constituents. These calculations therefore rationalize the use of simple liquid models for the quantitative prediction of surface tensions of polymers. As pressure and temperature increase, the surface tension of polystyrene decreases. A linear relationship is found between surface tension and temperature, and between surface tension and pressure; the slope of surface tension change with temperature is dependent on pressure.