Multicomponent diffusion in ionic crystals: theoretical model and application to combined tracer- and interdiffusion in alkali feldspar

We present a model for multicomponent diffusion in ionic crystals. The model accounts for vacancy-mediated diffusion on a sub-lattice and for diffusion due to binary exchange of different ionic species without involvement of vacancies on the same sub-lattice. The diffusive flux of a specific ionic species depends on the self-diffusion coefficients, on the diffusion coefficients related to the binary exchanges, and on the site fractions of all ionic species. The model delivers explicit expressions for these dependencies, which lead to a set of coupled non-linear diffusion equations. We applied the model to diffusion of [Formula: see text] Na, [Formula: see text] K, and [Formula: see text] K in alkali feldspar. To this end, gem-quality crystals of alkali feldspar were used together with [Formula: see text] K doped KCl salt as diffusion couples, which were annealed at temperatures between 800[Formula: see text] and 950[Formula: see text] C. Concentration-distance data for [Formula: see text] Na, [Formula: see text] K, and [Formula: see text] K were obtained by Time of Flight Secondary Ion Mass Spectrometry. Over the entire investigated temperature range the Na self-diffusion coefficient is by a factor of [Formula: see text] higher than the K self-diffusion coefficient. Diffusion mediated by binary [Formula: see text] K–[Formula: see text] K exchange is required for obtaining satisfactory fits of the model curves to the experimental data, and the respective kinetic coefficient is well constrained.