Curvature dependence of surface free energy of liquid drops and bubbles: A simulation study.

We study the excess free energy due to phase coexistence of fluids by Monte Carlo simulations using successive umbrella sampling in finite L×L×L boxes with periodic boundary conditions. Both the vapor-liquid phase coexistence of a simple Lennard-Jones fluid and the coexistence between A-rich and B-rich phases of a symmetric binary (AB) Lennard-Jones mixture are studied, varying the density ρ in the simple fluid or the relative concentration x(A) of A in the binary mixture, respectively. The character of phase coexistence changes from a spherical droplet (or bubble) of the minority phase (near the coexistence curve) to a cylindrical droplet (or bubble) and finally (in the center of the miscibility gap) to a slablike configuration of two parallel flat interfaces. Extending the analysis of Schrader et al., [Phys. Rev. E 79, 061104 (2009)], we extract the surface free energy γ(R) of both spherical and cylindrical droplets and bubbles in the vapor-liquid case and present evidence that for R→∞ the leading order (Tolman) correction for droplets has sign opposite to the case of bubbles, consistent with the Tolman length being independent on the sign of curvature. For the symmetric binary mixture, the expected nonexistence of the Tolman length is confirmed. In all cases and for a range of radii R relevant for nucleation theory, γ(R) deviates strongly from γ(∞) which can be accounted for by a term of order γ(∞)/γ(R)-1∝R(-2). Our results for the simple Lennard-Jones fluid are also compared to results from density functional theory, and we find qualitative agreement in the behavior of γ(R) as well as in the sign and magnitude of the Tolman length.

Anisotropic interfacial tension, contact angles, and line tensions: a graphics-processing-unit-based Monte Carlo study of the Ising model.

As a generic example for crystals where the crystal-fluid interface tension depends on the orientation of the interface relative to the crystal lattice axes, the nearest-neighbor Ising model on the simple cubic lattice is studied over a wide temperature range, both above and below the roughening transition temperature. Using a thin-film geometry L(x)×L(y)×L(z) with periodic boundary conditions along the z axis and two free L(x)×L(y) surfaces at which opposing surface fields ±H(1) act, under conditions of partial wetting, a single planar interface inclined under a contact angle θ<π/2 relative to the yz plane is stabilized. In the y direction, a generalization of the antiperiodic boundary condition is used that maintains the translational invariance in the y direction despite the inhomogeneity of the magnetization distribution in this system. This geometry allows a simultaneous study of the angle-dependent interface tension, the contact angle, and the line tension (which depends on the contact angle, and on temperature). All these quantities are extracted from suitable thermodynamic integration procedures. In order to keep finite-size effects as well as statistical errors small enough, rather large lattice sizes (of the order of 46 million sites) were found to be necessary, and the availability of very efficient code implementation of graphics processing units was crucial for the feasibility of this study.