Each multi‐electron state is characterized by the symmetry of the electron distribution as expressed by the symbol of the corresponding irreducible representation Γ and by the spin multiplicity M = 2S + 1 with the sum S of unpaired electrons. The coefficients z1, z2, and z3 depend on the number of d electrons and have been calculated and tabulated for each state (Tanabe & Sugano, 1954a). Figure 3.4 shows how the 3d electron configurations d5 (Fe3+) and d6 (Fe2+) give rise to multi‐electron states and how their energies vary as a function of the ratio Δ/B (Tanabe and Sugano, 1954b). We will see below that the ratio Δ/B increases with compression.
To incorporate the changes in state energies that result from compression, I assume the crystal‐field parameters to scale with the volume ratio V0/V as:
An electrostatic point charge model for an octahedrally coordinated transition metal cation suggests that δ = 5 (Burns, 1993). High‐pressure spectroscopic measurements, however, have shown that real materials may deviate from this prediction (Burns, 1985; Drickamer & Frank, 1973). Enhanced covalent bonding might also lead to a reduction of the Racah B parameter with increasing compression (Abu‐Eid & Burns, 1976; Keppler et al., 2007; Stephens & Drickamer, 1961a, 1961b) and would imply b < 0.
Most multi‐electron states are degenerate and allow a number m > 1 of different electron configurations. The total degeneracy of a given state MΓ is then given by the product mM and contributes a configurational entropy equal to kln(mM), where k is the Boltzmann constant. The Helmholtz free energy of a given state can then be expressed as (Badro et al., 2005; Sturhahn et al., 2005):
Figure 3.4 Energy diagrams for multi‐electron states that arise from the 3d electron configurations d5 (Fe3+) (a) and d6 (Fe2+) (b) in the limit of strong octahedral crystal fields. Note the change in electronic ground state for each electron configuration when the ratio Δ/B reaches the value marked by a vertical line. The distributions of d electrons over the crystal‐field orbitals t2 and e are illustrated for the three states that have the lowest energies at the change in electronic ground state.
At T > 0, and in particular at temperatures of Earth’s mantle, the d electron configuration can be thermally excited into states of higher energies. The thermal population of higher energy states becomes most important when the energies of other states approach the energy of the ground state, for example, as a result of compression. A spin transition results from a compression‐induced change of the electronic ground state from a high‐spin to a low‐spin state (Figure 3.4). At finite temperatures, the broadening of the transition reflects the incipient population of the low‐spin state before and the decaying population of the high‐spin state after their energies cross over as long as the energies of both states remain close enough to allow for thermally induced transitions between both states (Holmström & Stixrude, 2015; Sturhahn et al., 2005; Tsuchiya et al., 2006). However, states other than the high‐spin or low‐spin state may get close enough in energy to the ground state to become thermally populated by a significant fraction of d electrons. The thermal distribution of electrons over accessible multi‐electron states contributes to the overall configurational entropy. The excess contribution to the Helmholtz free energy per transition metal cation that results from compression‐ and temperature‐induced changes in the occupation of electronic states can then be identified as:
where Δϕ = ϕ(V, T) − ϕ0(V0, T0) are the changes in the fractions ϕ of d electrons that occupy each electronic state at a given volume and temperature with respect to the fractions ϕ0 at ambient conditions. Note that the configurational entropy term is multiplied by the number d of d electrons per transition metal cation. The fraction of d electrons that occupy a given state can be found by applying the equilibrium condition (∂FEL/∂ni)V,T,N = 0 to a micro‐canonical ensemble with the absolute numbers ni = ϕiN of electrons in each state and the total number N of electrons (Sturhahn et al., 2005):
The denominator sums over all multi‐electron states being considered.
The excess contributions to pressure and bulk modulus then follow from the definitions:
and
All excess quantities need to be multiplied by the number and fraction of crystallographic sites that are occupied by the transition metal cation.
The parametrization of spin transitions presented above builds on the formulation proposed by Sturhahn et al. (2005) and resembles earlier attempts to predict the effect of spin transitions in Fe2+ on thermodynamic properties of mantle minerals (Badro et al., 2005; Gaffney, 1972; Gaffney & Anderson, 1973; Ohnishi, 1978). Unlike these forward modeling approaches, I will use the formulation here in a semi‐empirical way with adjustable parameters that will be constrained by experimental observations. In contrast to formulations based on the Gibbs free energy with pressure as variable and elastic compliances as parameters that have proven useful in describing the results of DFT computations (Shukla et al., 2016; Wentzcovitch et al., 2009; Wu et al., 2013), a formulation of excess contributions in terms of volume and temperature is more consistent with finite‐strain theory based on the assumption of homogeneous and isotropic finite strain and can, in principle, be generalized to evaluate the effect on individual components of the elastic stiffness tensor and hence on the shear modulus. The impact of spin transitions on the shear modulus, however, appears to be minor (Fu et al., 2018; Marquardt et al., 2009b; Shukla et al., 2016; Wu et al., 2013). Note that by accounting for changes in the electronic configurations in terms of excess properties it is no longer necessary to derive sets of finite‐strain parameters for each individual electronic configuration as required by earlier formulations based on high‐spin and low‐spin states only (Chen et al., 2012; Speziale et al., 2007; Wu et al., 2013, 2009), since the excess contributions are simply added to the (cold) elastic and thermal contributions. Although I assume here that transition metal cations with different electronic configurations mix ideally among each other and with other cations, it is in principle possible to add further terms to account for nonideal mixing behavior and interactions between transition metal cations (Holmström and Stixrude, 2015; Ohnishi and Sugano, 1981; Sturhahn et al., 2005). To some extent and for