where δ δ j , k is 1 if j and k are the same site and zero otherwise. J is a prefactor with dimensions of energy, and g is another coupling coefficient that determines the relative strength of the external field compared to the nearest neighbour interaction.