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 ${\mathit{J}}_{\mathit{s}}$. The surface-bulk multicritical point is located with improved precision, ${\mathit{J}}_{\mathit{s}}$/J=1.52±0.02, and crossover behavior is studied. New estimates for critical exponents are extracted, ${\gamma }_1$=0.78±0.06, ${\gamma }_1^{\mathit{m}}$=1.41±0.14, ${\gamma }_{11}^{\mathit{m}}$=0.96±0.09, and along with critical amplitude ratios these are compared with theoretically predicted values.