Dynamics and hysteresis in square lattice artificial spin ice

Dynamical effects under geometrical frustration are considered in a model for artificial spin ice on a square lattice in two dimensions. Each island of the spin ice has a three-component Heisenberg-like dipole moment subject to shape anisotropies that influence its direction. The model has real dynamics, including rotation of the magnetic degrees of freedom, going beyond the Ising-type models of spin ice. The dynamics is studied using a Langevin equation solved via a second-order Heun algorithm. Thermodynamic properties such as the specific heat are presented for different couplings. A peak in specific heat is related to a type of melting-like phase transition present in the model. Hysteresis in an applied magnetic field is calculated for model parameters where the system is able to reach thermodynamic equilibrium.