Simulating the critical behaviour of complex fluids

ALTHOUGH the liquid–gas phase equilibria of simple fluids have been studied extensively since the seminal work of van der Waals, the properties of fluids with more complex molecular structures, such as polymers, present a less tractable problem both theoretically and experimentally. The phase behaviour of hydrocarbons is of particular importance for the petrochemical industry. But despite significant experimental and theoretical efforts, the phase diagrams of the straight-chain alkanes longer than decane (C10) are known only partially, and even qualitative aspects such as the chain-length dependence of the critical properties are poorly understood. Until recently it was considered impossible to estimate the critical properties of such complex fluids using computer simulations. Here we report Monte Carlo simulations of the phase diagrams of alkanes with unbranched carbon chains as long as C48, up to the vicinity of the liquid–vapour critical points. Our calculations show that, in contrast to the traditional view, the critical density of the long-chain alkanes decreases rather than increases with carbon number. This work indicates that simulations can be used as an 'engineering tool' to estimate properties that are not readily accessible experimentally.

[1]  C. Tsonopoulos Critical constants of normal alkanes from methane to polyethylene , 1987 .

[2]  D. Frenkel,et al.  Computer simulations in the Gibbs ensemble , 1989 .

[3]  Amyn S. Teja,et al.  The critical temperatures and densities of the n-alkanes from pentane to octadecane , 1990 .

[4]  A. Panagiotopoulos Direct determination of phase coexistence properties of fluids by Monte Carlo simulation in a new ensemble , 1987 .

[5]  Juan J. de Pablo,et al.  Simulation of phase equilibria for chain molecules , 1992 .

[6]  Herman J. C. Berendsen,et al.  MOLECULAR-DYNAMICS SIMULATION OF A BILAYER-MEMBRANE , 1982 .

[7]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[8]  Athanassios Z. Panagiotopoulos,et al.  Phase equilibria by simulation in the Gibbs ensemble , 1988 .

[9]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid hydrocarbons , 1984 .

[10]  I. R. Mcdonald,et al.  Monte Carlo study of the properties of self-assembled monolayers formed by adsorption of CH3(CH2)15SH on the (111) surface of gold , 1993 .

[11]  B. Smit,et al.  Vapour-liquid equilibria for quadrupolar Lennard-Jones fluids , 1990 .

[12]  Daan Frenkel,et al.  Configurational bias Monte Carlo: a new sampling scheme for flexible chains , 1992 .

[13]  Berend Smit,et al.  Novel scheme to study structural and thermal properties of continuously deformable molecules , 1992 .

[14]  C. Tsonopoulos,et al.  The Critical Constants of Normal Alkanes From Methane to Polyethylene: II. Application of the Flory Theory , 1993 .

[15]  Søren Toxvaerd,et al.  Molecular dynamics calculation of the equation of state of alkanes , 1990 .