We report on a procedure developed to create stochastic models of hyperfine interactions for complex diffusion mechanisms and demonstrate its application to simulate perturbed angular correlation spectra for the divacancy and 6-jump cycle diffusion mechanisms in L1 2 -structured compounds.Diffusion plays an important role in materials processing and in determining suitable applications for materials. Therefore, it is of interest to refine methodologies capable of discerning individual atomic jumps that take place in the diffusion process, particularly in intermetallic compounds for which complex diffusion mechanisms can dominate the simple vacancy mechanism. Here, complex diffusion refers to cases where diffusing species change sublattices or multiple defect species cooperate in diffusion. In principle, hyperfine methods are ideal for discerning atomic-scale diffusion processes because of the short-range nature of the quadrupole interaction. More specifically, perturbed angular correlation (PAC) spectroscopy is attractive because it is sensitive to jump rates on the order of 0.1 to 1 GHz, which can be obtained in typical materials in the experimentally convenient temperature range between room temperature and 1200 • C. This work is funded in part by NSF grant DMR 06-06006 (Metals Program) and computational resources were provided in part by KY EPSCoR grant RSF 012-03.