Magnetic hysteresis in 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. 
Gyorgy Szabo and Gyorgy Kadar