A linear/quadratic order parameter coupling description of the Verwey transition in magnetite, FeO.

The change in symmetry Fd3m → Cc at the Verwey transition in magnetite puts it in a class of phase transitions with linear/quadratic coupling between two separate order parameters. Direct coupling between an order parameter Q, to represent an electronic instability, and an order parameter Q, to represent cation charge ordering, has the form λQQ, with T(Q) < T(Q), but there must also be indirect coupling through the common strain, e, due to strain coupling terms λeQ and λeQ. Q has the symmetry of irrep Γ while the pattern of cation charge ordering of Fe and Fe on octahedral sites depends on some combination of irreps Δ, X, X, W and W. The software package ISOVIZ has been used to show how reported patterns of order for the simplified structure in space group P2/c can be understood in terms of a linearly dependent mix of patterns with symmetry Δ and X, so that Q can be treated in the first instance as though it has the symmetry of Δ. Spontaneous strains calculated from published lattice parameters and symmetry-adapted atomic displacements from previous structural refinements in space group Cc have been used to confirm that the two order parameters have different temperature dependences, consistent with this phenomenological treatment. The effect of chemical doping can be understood in terms of the development of local strain heterogeneity which acts to suppress the macroscopic strains and which appears to have a greater influence on charge ordering than on the electronic structure.