Interaction of rhyolite melts with monazite,xenotime, and zircon surfaces

The interfacial contact region between a rhyolite melt and the accessory minerals monazite, xenotime, and zircon is investigated using molecular dynamics simulations. On all surfaces, major structural rearrangement extends about 1 nm into the melt from the interface. As evidenced by the structural perturbations in the ion distribution profiles, the affinity of the melt for the surface increases in going from monazite to xenotime to zircon. Alkali ions are enriched in the melt in contact with an inert wall, as well as at the mineral surfaces. Melt in contact with zircon has a particularly strong level of aluminum enrichment. In xenotime, the enrichment of aluminum is less than that in zircon, but still notable. In monazite, the aluminum enrichment in the contact layer is much less. It is expected that the relative surface energies of these accessory minerals will be a strong function of the aluminum content of the melt and that nucleation of zircon, in particular, would be easier for melts with higher aluminum concentration. The crystal growth rate for zircon is expected to be slower at a higher aluminum concentration because of the effectiveness of aluminum in solvating the zircon surface. The variable interfacial concentration profiles across the series of accessory minerals will likely affect the kinetics of trace element incorporation, as the trace elements must compete with the major elements for surface sites on the growing accessory minerals.