Extension of the Test-Area methodology for calculating solid-fluid interfacial tensions in cylindrical geometry

J Chem Phys. 2013 Apr 7;138(13):134701. doi: 10.1063/1.4795836.

Abstract

We extend the well-known Test-Area methodology of Gloor et al. [J. Chem. Phys. 123, 134703 (2005)], originally proposed to evaluate the surface tension of planar fluid-fluid interfaces along a computer simulation in the canonical ensemble, to deal with the solid-fluid interfacial tension of systems adsorbed on cylindrical pores. The common method used to evaluate the solid-fluid interfacial tension invokes the mechanical relation in terms of the tangential and normal components of the pressure tensor relative to the interface. Unfortunately, this procedure is difficult to implement in the case of cylindrical geometry, and particularly complex in case of nonspherical molecules. Following the original work of Gloor et al., we perform free-energy perturbations due to virtual changes in the solid-fluid surface. In this particular case, the radius and length of the cylindrical pore are varied to ensure constant-volume virtual changes of the solid-fluid surface area along the simulation. We apply the modified methodology for determining the interfacial tension of a system of spherical Lennard-Jones molecules adsorbed inside cylindrical pores that interact with fluid molecules through the generalized 10-4-3 Steele potential recently proposed by Siderius and Gelb [J. Chem. Phys. 135, 084703 (2011)]. We analyze the effect of pore diameter, density of adsorbed molecules, and fluid-fluid cutoff distance of the Lennard-Jones intermolecular potential on the solid-fluid interfacial tension. This extension, as the original Test-Area formulation, offers clear advantages over the classical mechanical route of computational efficiency, easy of implementation, and generality.