Molecular dynamics study of the surface tension of a binary immiscible fluid

J Chem Phys. 2004 Feb 8;120(6):2893-900. doi: 10.1063/1.1637340.

Abstract

The planar interface between two liquids having two degrees of affinity to mix has been studied by molecular dynamics simulations. The surface tension is calculated from the normal, PN, and transverse, PT, components of the pressure tensor P for a wide range of temperatures. An unusual increase in surface tension with increasing temperature is attributed to a pressure induced void transfer mechanism that is justified by basic thermodynamic arguments. This effect is diminished on the addition of a modest attractive potential between the two species, and there is a turnover point at higher temperatures beyond which the surface tension decreases with increasing temperature. An order parameter is identified as the gradient of the mole fraction distribution through the interfacial region. An additional effect is the dramatic inversion of the kinetic and potential contributions to the PN profile as the temperature is varied. It is found that a commonly used approximation for P, the Irving-Kirkwood 1 or IK1 method, results in a relatively modest unphysical variability in PN that weakly violates the condition of local mechanical stability. However, this artifact does not prevent the IK1 method from producing an interfacial tension which is nearly identical to that derived from the complete IK formula with no additional approximations.