Articulo
Heterogeneous nucleation of a droplet pinned at a chemically inhomogeneous substrate: a simulation study of the two-dimensional Ising case
Autor
Trobo, Marta Liliana
Albano, Ezequiel Vicente
Binder, Kurt
Institución
Resumen
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<sub>b</sub> 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<sub>b</sub> < H<sub>b</sub><sup>crit</sup>, the critical droplet radius is so large that a critical droplet having the contact angle θ<sub>c</sub> 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<sub>b</sub><sup>crit</sup> < H<sub>b</sub> < H<sub>b</sub>*, 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 θ<sub>eff</sub> in the regime θ<sub>c</sub> < θ<sub>eff</sub> < π/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 θ<sub>eff</sub> and the excess density of the droplets, scaled by b, are functions of the product bH<sub>b</sub> but do not depend on both variables separately. Since the free energy barrier for the "depinning" of the droplet (i.e., growth of θ<sub>eff</sub> to π - θ<sub>c</sub>) vanishes when θ<sub>eff</sub> approaches π/2, in practice only angles θ<sub>eff</sub> up to about θ<sub>eff</sub><sup>max</sup> ≃ 70 were observed. For larger fields (H<sub>b</sub> > H<sub>b</sub>*), the droplets nucleated at the chemical inhomogeneity grow to the full system size. While the relaxation time for the growth scales as τ<sub>G</sub> ∝ H<sub>b</sub><sup>-1</sup>, the nucleation time τ<sub>N</sub> scales as ln <sub>N</sub> ∝ H<sub>b</sub><sup>-1</sup>. 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 θ<sub>c</sub> is used. Facultad de Ciencias Exactas Facultad de Ingeniería Instituto de Física de Líquidos y Sistemas Biológicos