Surface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water

The test-area (TA) perturbation approach has been gaining popularity as a methodology for the direct computation of the interfacial tension in molecular simulation. Though originally implemented for planar interfaces, the TA approach has also been used to analyze the interfacial properties of curved liquid interfaces. Here, we provide an interpretation of the TA method taking the view that it corresponds to the change in free energy under a transformation of the spatial metric for an affine distortion. By expressing the change in configurational energy of a molecular configuration as a Taylor expansion in the distortion parameter, compact relations are derived for the interfacial tension and its energetic and entropic components for three different geometries: planar, cylindrical, and spherical fluid interfaces. While the tensions of the planar and cylindrical geometries are characterized by first-order changes in the energy, that of the spherical interface depends on second-order contributions. We show that a greater statistical uncertainty is to be expected when calculating the thermodynamic properties of a spherical interface than for the planar and cylindrical cases, and the evaluation of the separate entropic and energetic contributions poses a greater computational challenge than the tension itself. The methodology is employed to determine the vapour-liquid interfacial tension of TIP4P/2005 water at 293 K by molecular dynamics simulation for planar, cylindrical, and spherical geometries. A weak peak in the curvature dependence of the tension is observed in the case of cylindrical threads of condensed liquid at a radius of about 8 Å, below which the tension is found to decrease again. In the case of spherical drops, a marked decrease in the tension from the planar limit is found for radii below ∼ 15 Å; there is no indication of a maximum in the tension with increasing curvature. The vapour-liquid interfacial tension tends towards the planar limit for large system sizes for both the cylindrical and spherical cases. Estimates of the entropic and energetic contributions are also evaluated for the planar and cylindrical geometries and their magnitudes are in line with the expectations of our simple analysis.