Abstract
A numerical method involving the time integration of equations of motion is developed for calculating response functions of antiferromagnets. The method is particularly useful in calculating dynamic response functions S(k, omega ) and can also be used to obtain densities of states. The technique is applied to the two-dimensional random antiferromagnets Rb2Mn0.5Ni0.5F4 and Rb2MncZn1-cF4. In the former case, good agreement is found with the recently measured S(k, omega ) and in the latter some interesting clustering effects that should be accessible experimentally are found.