  We examine the performance of several molecular simulation techniques aimedat evaluation of the surface tension through its thermodynamic definition. Forall methods explored, the surface tension is calculated by approximating thechange in Helmholtz free energy associated with a change in interfacial areathrough simulation of a liquid slab at constant particle number, volume, andtemperature. The methods explored fall within three general classes:free-energy perturbation, the Bennett acceptance-ratio scheme, and the expandedensemble technique. Calculations are performed for both the truncatedLennard-Jones and square-well fluids at select temperatures spaced along theirrespective liquid-vapor saturation lines. Overall, we find that Bennett andexpanded ensemble approaches provide the best combination of accuracy andprecision. All of the methods, when applied using sufficiently small areaperturbation, generate equivalent results for the Lennard-Jones fluid. However,single-stage free-energy-perturbation methods and the closely related test-areatechnique recently introduced by Gloor et al. [J. Chem. Phys. 123, 134703(2005)], generate surface tension values for the square-well fluid that are notconsistent with those obtained from the more robust expanded ensemble andBennett approaches, regardless of the size of the area perturbation.Single-stage perturbation methods fail also for the Lennard-Jones system whenapplied using large area perturbations. Here an analysis of phase-space overlapproduces a quantitative explanation of the observed inaccuracy, and shows thatthe satisfactory results obtained in these cases from the test-area methodarise from a cancellation of errors that cannot be expected in general.
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.