Magnetic hysteresis in an Ising-like dipole-dipole model

Using zero-temperature Monte Carlo simulations we have studied the magnetic hysteresis in a three-dimensional Ising model with nearest-neighbor exchange and dipolar interaction. The average magnetization of spins located inside a sphere on a cubic lattice is determined as a function of magnetic field varied periodically. The simulations have justified the appearance of hysteresis and allowed us to have a deeper insight into the series of metastable states developed during this process.