Quick Search:    advanced search

GSW Home   GeoRef Home   My GSW Alerts   Contact GSW   About GSW   Journals List   Help

This Article Figures Only Full Text Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Palin, E. J. Articles by Harrison, R. J. Search for Related Content GeoRef GeoRef Citation

A Monte Carlo investigation of the thermodynamics of cation ordering in 2-3 spinels

Department of Earth Sciences, University of Cambridge, Downing Street, Cambridge CB2 3EQ, U.K.

Correspondence: ^* E-mail: ejp24{at}cam.ac.uk

The Monte Carlo (MC) simulation technique is a powerful tool for the investigation of thermodynamic and kinetic phenomena in minerals, and is especially well suited to the study of cation ordering. We have performed MC simulations of eight end-member 2-3 spinels (X²⁺ = Mg, Fe, Zn, Ni; X³⁺ = Al, Fe) using pair interaction parameters, J_i, and chemical potentials, µ, derived from atomistic simulations. The J_i values for all but one of these spinels are remarkably similar, despite their different character (normal vs. inverse). The sign of µ, and hence the tendency to form a normal or inverse spinel, was correctly predicted in all cases. Agreement between the simulated and observed cation distributions as a function of temperature is good for the normal spinels and poor for the inverse spinels. Agreement could be greatly improved for the inverse spinels through relatively modest adjustments to the simulation parameters (usually increasing the strength of the tetrahedral-octahedral, T-O, interactions, and decreasing the magnitude of µ).

We have developed an atomistic random-mixing model for cation ordering in spinels and compared it with the macroscopic O’Neill-Navrotsky model. In so doing, we have determined the relative contributions of µ, tetrahedral-tetrahedral (T-T), octahedral-octahedral (O-O), and T-O interactions to the O’Neill-Navrotsky coefficients and ß. We found that the value of ß depends on the relative enthalpy contributions of (T-T + O-O) vs. T-O interactions, a useful insight considering the large spread of values found experimentally to be taken by ß.

We used the thermodynamic integration technique to quantify the entropy, and hence the amount of short-range order, present in the spinels studied. We found that there is virtually no short-range order in the normal spinels. There is significant short-range order in the inverse spinels, though in the experimentally accessible temperature range, the contribution of this short-range order to the entropy is comparatively small. At very low temperatures, we find that the octahedral cations in the inverse spinels become ordered, reducing the symmetry to P4₁22, in agreement with other simulated findings for 2-3 spinels and experimental findings for 4-2 spinels.

Key Words: Spinel • thermodynamics • Monte Carlo simulation • O’Neill-Navrotsky model • short-range order

This article has been cited by other articles:

				E. J. Palin, A. M. Walker, and R. J. Harrison A computational study of order-disorder phenomena in Mg2TiO4 spinel (qandilite) American Mineralogist, August 1, 2008; 93(8-9): 1363 - 1372. [Abstract] [Full Text] [PDF]

				E. J. Palin and R. J. Harrison A computational investigation of cation ordering phenomena in the binary spinel system MgAl2O4-FeAl2O4 Mineralogical Magazine, December 1, 2007; 71(6): 611 - 624. [Abstract] [Full Text] [PDF]