Experimental determination of Fe–Mg interdiffusion in orthopyroxene as a function of Fe content

We have measured the dependence of the Fe–Mg interdiffusion coefficient, D Fe-Mg , on the ferrosilite component in orthopyroxene, which so far has not been experimentally calibrated. Diffusion couples, consisting of approximately 1 µm thin-films were deposited by pulsed laser ablation on orthopyroxene crystals of En 91 Fs 9 composition. Diffusion experiments were carried out at atmospheric pressure in vertical gas mixing furnaces (CO-CO 2 ) at temperatures between 950 and 1100 °C at constant f O 2 = 10 –7 Pa. Using a focused ion beam-scanning electron microscope (FIB-SEM), FIB-foils were cut from diffusion couples before and after annealing. Diffusion profiles were extracted by using combined backscattered electron (BSE) imaging and energy dispersive X-ray (EDXS) mapping on FIB-foils which allowed to resolve concentration gradients within 20 nm. The microstructure of the diffusion experiments was investigated using transmission electron microscopy (TEM). Using this method, we obtained the first experimentally determined data on the dependence of D Fe-Mg on the ferrosilite content in orthopyroxene at different temperatures, appearing to increase with increasing a temperature. For the temperature range (950 – 1100 °C), log f O 2 = -7 Pa, along [001] in the composition Fs9, the log D Fe-Mg yields the following Arrhenius equation: DFe - Mg [m² /s] = 3.8 · 10- 9 exp [ - 261.07 ± 24[kJ/mol]/(R/T[K])] D F e - M g [ m 2 / s ] = 3.8 · 10 - 9 exp [ - 261.07 ± 24 [ k J / m o l ] / ( R / T [ K ] ) ] The effect of X Fe on D Fe-Mg, given by D(XFe ) = D(XFe = 0.09) · 10^m (XFe - 0.09), D ( X Fe ) = D ( X Fe = 0.09 ) · 10 m ( X Fe - 0.09 ) , can be calculated by the following parameterization where m follows a linear regression of m on temperature: m = - 2.711 · 10⁴ /T(K) + 23.5408 m = - 2.711 · 10 4 / T ( K ) + 23.5408 By considering the D Fe-Mg dependence on X Fe , the timescales of natural processes obtained from modelling the compositional zoning of natural crystal may considerably differ from previously estimated timescales.