Wetting transitions near the bulk critical point: Monte Carlo simulations for the Ising model

Abstract
Critical, tricritical, and first-order wetting transitions are studied near the bulk critical point of a simple cubic nearest-neighbor Ising model by extensive Monte Carlo simulations. The model applies an exchange J in the bulk and exchange Js in the surface planes, where surface fields H1 also act in addition to a possible bulk field H. Lattices in a thin-film geometry L×L×D are used, with two free L×L surfaces (with L up to 256) and film thickness D up to 160, applying a very fast fully vectorizing multispin coding program. Our results present the first quantitative evidence for the scaling theory due to Nakanishi and Fisher, which links wetting behavior and surface critical behavior. In particular, we show that the tricritical wetting line (Jst/J) merges into the surface-bulk multicritical point with associated critical field H1t∼(1-T/Tc )Δ1m, while the critical field for critical wetting H1c vanishes as &, where Δ1(Δ1m) are the ‘‘gap’’ exponents for the surface-layer magnetization at the ordinary (or surface-bulk multicritical) transition. The mean-field character of critical wetting in this model is again confirmed.