Droplets pinned at chemically inhomogenous substrates: A simulation study of the two-dimensional Ising case

Height of a liquid drop on a wetting stripe

Adsorption of liquid on a planar wall decorated by a hydrophilic stripe of width L is considered. Under the condition that the wall is only partially wet (or dry) while the stripe tends to be wet completely, a liquid drop is formed above the stripe. The maximum height ℓ_{m}(δμ) of the drop depends on the stripe width L and the chemical potential departure from saturation δμ where it adopts the value ℓ_{0}=ℓ_{m}(0). Assuming a long-range potential of van der Waals type exerted by the stripe, the interfacial Hamiltonian model is used to show that ℓ_{0} is approached linearly with δμ with a slope which scales as L^{2} over the region satisfying L≲ξ_{∥}, where ξ_{∥} is the parallel correlation function pertinent to the stripe. This suggests that near the saturation there exists a universal curve ℓ_{m}(δμ) to which the adsorption isotherms corresponding to different values of L all collapse when appropriately rescaled. Although the series expansion based on the interfacial Hamiltonian model can be formed by considering higher order terms, a more appropriate approximation in the form of a rational function based on scaling arguments is proposed. The approximation is based on exact asymptotic results, namely, that ℓ_{m}∼δμ^{-1/3} for L→∞ and that ℓ_{m} obeys the correct δμ→0 behavior in line with the results of the interfacial Hamiltonian model. All the predictions are verified by the comparison with a microscopic density functional theory (DFT) and, in particular, the rational function approximation-even in its simplest form-is shown to be in a very reasonable agreement with DFT for a broad range of both δμ and L.

Heterogeneous nucleation of a droplet pinned at a chemically inhomogeneous substrate: A simulation study of the two-dimensional Ising case

Heterogeneous nucleation is studied by Monte Carlo simulations and phenomenological theory, using the two-dimensional lattice gas model with suitable boundary fields. A chemical inhomogeneity of length b at one boundary favors the liquid phase, while elsewhere the vapor is favored. Switching on the bulk field H_b favoring the liquid, nucleation and growth of the liquid phase starting from the region of the chemical inhomogeneity are analyzed. Three regimes occur: for small fields, H_b<H_b^crit, the critical droplet radius is so large that a critical droplet having the contact angle θ_c required by Young's equation in the region of the chemical inhomogeneity does not yet "fit" there since the baseline length of the circle-cut sphere droplet would exceed b. For H_b^crit<H_b<H_b^*, such droplets fit inside the inhomogeneity and are indeed found in simulations with large enough observation times, but these droplets remain pinned to the chemical inhomogeneity when their baseline has grown to the length b. Assuming that these pinned droplets have a circle cut shape and effective contact angles θ_eff in the regime θ_c < θ_eff < π/2, the density excess due to these droplets can be predicted and is found to be in reasonable agreement with the simulation results. On general grounds, one can predict that the effective contact angle θ_eff and the excess density of the droplets, scaled by b, are functions of the product bH_b but do not depend on both variables separately. Since the free energy barrier for the "depinning" of the droplet (i.e., growth of θ_eff to π - θ_c) vanishes when θ_eff approaches π/2, in practice only angles θ_eff up to about θ_eff^max≃70° were observed. For larger fields (H_b>H_b^*), the droplets nucleated at the chemical inhomogeneity grow to the full system size. While the relaxation time for the growth scales as τ_G∝H_b^-1, the nucleation time τ_N scales as lnτ_N∝H_b^-1. However, the prefactor in the latter relation, as evaluated for our simulations results, is not in accord with an extension of the Volmer-Turnbull theory to two-dimensions, when the theoretical contact angle θ_c is used.