We investigate the quasistatic magnetic hysteresis of ferromagnetic thin films grown on a vicinal substrate, using Monte Carlo simulations within a two-dimensional XY model. Intrinsic in-plane anisotropy is assigned to surface sites according to their local symmetry. The simulated hysteresis loops show a strong anisotropy: the coercive field is the largest when the external field is along the step direction and vanishes when the external field is perpendicular to the step direction. In general, the coercivity increases with increasing step density, but displays a more complex dependence on film thickness. The simulations also suggest that the mechanism for the magnetization reversal is coherent rotation. These results are in good agreement with experiments. © 2002 American Institute of Physics. ͓DOI: 10.1063/1.1433179͔The magnetic properties of ferromagnetic thin films and multilayers have been extensively studied because of their potential impact on magnetic recording devices. As the thickness of a film is reduced, its properties are expected to be strongly influenced by surfaces and interfaces, which are inevitably rough at atomic scales. The ultimate goal of studies of the influences of surface/interface roughness on magnetic properties of thin films 1-17 is to engineer desirable magnetic properties by artificially creating and controlling the surface/ interface structure and morphology.Surface/interface roughness strongly influences both the static ͑e.g., magnetization͒ 17 and dynamic ͑e.g., magnetic hysteresis͒ 16 magnetic properties of thin films. Experiments have demonstrated that surface steps induce an in-plane uniaxial magnetic anisotropy, with the easy axis parallel to the step direction, in a variety of magnetic thin films 1-11 grown on metal as well as on semiconductor stepped surfaces. The measured hysteresis loops show that the coercive field decreases while the saturation field increases when the external field is turned away from the direction parallel to the steps to the direction perpendicular to the steps, and such uniaxial anisotropy increases with increasing step density. 10 Theoretical modeling and simulations have been carried out to investigate the effect of steps on magnetization reversal. 14 -16 In particular, Hyman et al. 16 have recently derived, within a two-dimensional XY model, a phase diagram of hysteresis loops in the parameter space of anisotropy strength and step density, for magnetic reversal on vicinal surfaces.In this article, we carry out Monte Carlo simulations of hysteresis loops of ferromagnetic thin films to extend these earlier theoretical studies. 14 -16 Our focus is to investigate systematically the effect of step density and film thickness on magnetic hysteresis anisotropy in ultrathin films grown on vicinal substrates. Our simulations show that the coercive ͑saturation͒ field decreases ͑increases͒ monotonically as the angle between the external field and the step direction is increases from 0°to 90°, exhibiting a strong hysteresis anisotropy. The coercivity, for ...