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

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 Jₛ in the surface planes, where surface fields H₁ 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 (Jₛᵗ/J) merges into the surface-bulk multicritical point with associated critical field H₁ₜ{~}(1-T/Tc{)}^{{Δ}1m}$, while the critical field for critical wetting ${H}1c$ vanishes as , where ${{Δ}}₁$(${{Δ}}₁ₘ$) 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.