Monte Carlo study of surface phase transitions in the three-dimensional Ising model

We present the results of extensive Monte Carlo simulations of phase transitions and critical behavior at the surface of a simple cubic Ising model. Profiles of the magnetization and internal energy are determined as a function of the distance from the surface, and we extract surface and bulk properties as a function of temperature and surface coupling Jₛ. The surface-bulk multicritical point is located with improved precision, Jₛ/J=1.52±{}0.02, and crossover behavior is studied. New estimates for critical exponents are extracted, γ₁=0.78±{}0.06, γ₁ᵐ=1.41±{}0.14, γ₁₁ᵐ=0.96±{}0.09, and along with critical amplitude ratios these are compared with theoretically predicted values.