Computer simulation of a gas–liquid surface. Part 1

Abstract

The gas–liquid surface of a system of Lennard-Jones (12, 6) molecules has been simulated by Monte Carlo and by Molecular Dynamic methods at temperatures which span most of the liquid range. For systems of 255 molecules the two methods lead to similar results and this agreement confirms that the density profile, as a function of height, falls monotonically from the density of the bulk liquid to that of the gas. The thickness of the surface layer is sensitive to the surface area, and appears to approach its thermodynamic limit for surface areas of 400σ² for a system of 4080 molecules. The density profile can be represented by a hyperbolic tangent of an appropriately scaled height. The thickness of the surface is of the order of two molecular diameters at temperatures near the triple point and increases rapidly as the critical point is approached. The computed surfacetens ions agree well with those calculated by statistical perturbation theory.

Monte Carlo and Molecular Dynamic simulation of a binary mixture shows clearly the adsorption of the component of higher vapour pressure; the amount absorbed agrees with that calculated from Gibbs's isotherm.