Energy minimization methods were used to simulate the migration of Zr, Si,
and O vacancies in zircon (ZrSiO4). Two sets of interatomic potentials were
employed for comparison: one with O-Si-O three-body terms for the SiO4, an
d one without. Results for Si were inconclusive. but consistent with mainta
ining the integrity of the SiO4 molecular units. Both Zr and O vacancies ca
n migrate on three-dimensional sublattice networks, thus supporting the exp
erimentally observed diffusional isotropy. The predicted Zr vacancy migrati
on energy (1.16-1.38 eV) was in good agreement with experiment if supplemen
ted by Zr vacancy formation via Schottky or Frenkel defects (6.21-12.28 eV/
defect). Oxygen vacancy migration energies were predicted to be 0.99-1.16 e
V, somewhat lower than the experimental Value of 4.64 eV measured in natura
l zircons, which thus may include significant contributions from vacancy fo
rmation mechanisms at 3.31-6.52 eV/defect. (C) 1999 Elsevier Science B.V. A
ll rights reserved.