Monte Carlo Simulation of a Random-Field Ising Antiferromagnet
  Phase transitions in the three-dimensional diluted Ising antiferromagnet in an applied magnetic field are analyzed numerically. It is found that random magnetic field in a system with spin concentration below a certain threshold induces a crossover from second-order phase transition to first-order transition to a new phase characterized by a spin-glass ground state and metastable energy states at finite temperatures. 
V. V. Prudnikov, V. N. Borodikhin