Stress-induced changes in magnetite: Insights from a numerical analysis of the verwey transition

Magnetic susceptibility behaviour around the Verwey transition of magnetite (≈ 125 K) is known to be sensitive to stress, composition and oxidation. From the isotropic point (≈ 130 K) to room temperature, decreasing magnetic susceptibility indicates an increase in magnetocrystalline anisotropy. In this study, we present a model which numerically analyses low-temperature magnetic susceptibility curves (80 to 280 K) of an experimentally shocked (up to 30 GPa) and later heated (973 K) magnetite ore. To quantify variations of the transition shape caused by both shock and heating, the model statistically describes local variations in the Verwey transition temperature within bulk magnetite. For the description, Voigt profiles are used, which indicate variations between a Gaussian and a Lorentzian character. These changes are generally interpreted as variations in the degree of correlation between observed events, i.e. between local transition temperatures in the model. Shock pressures exceeding the Hugoniot elastic limit of magnetite ($ \ge $ 5 GPa) cause an increase in transition width and Verwey transition temperature, which is partially recovered by heat treatment. Above the Verwey transition temperature, susceptibility variations related to the magnetocrystalline anisotropy are described with an exponential approach. The room temperature magnetic susceptibility relative to the maximum near the isotropic point is reduced after shock, which is related to grain size reduction. Since significant oxidation and cation substitution can be excluded for the studied samples, variations are only attributed to changes in elastic strain associated with shock-induced deformation and annealing due to heat treatment. The shocked magnetite shows a high correlation between local transition temperatures which is reduced by heat treatment. The model allows a quantitative description of low-temperature magnetic susceptibility curves of experimentally shocked and subsequently heat-treated polycrystalline magnetite around the Verwey transition temperature. The curves are accurately reproduced within the experimental uncertainties. Further applications for analysing magnetite-bearing rocks seem possible if model parameters, such as for oxidation are included into the model.