Real-space determination of orbital states driving successive phase transitions in FeV2O4

"Direct experimental access to orbital states in strongly correlated materials remains a major challenge, despite their central role in driving coupled structural and magnetic phase transitions. In systems where electronic correlations, electron–lattice coupling, and relativistic spin–orbit interactions compete on comparable energy scales, even first-principles calculations often yield multiple metastable solutions, hindering the unambiguous identification of the ground state. Here, we demonstrate that the orbital states of the spinel oxide FeV2O4, which possesses active orbital degrees of freedom on both Fe and V ions, are uniquely resolved by combining valence electron density (VED) analysis based on state-of-the-art synchrotron x-ray diffraction with spin-polarized density-functional-theory calculations. Our results reveal that temperature-dependent rearrangements of orbital occupations drive successive structural transitions that accompany collinear and noncoplanar ferrimagnetic orders, establishing a direct correspondence between orbital anisotropy and spin structure. More broadly, this work shows that experimentally determined VED provides a decisive real-space constraint on competing theoretical solutions, offering a powerful and broadly applicable framework for elucidating the microscopic mechanisms of complex phase transitions in strongly correlated electron systems."

"The anisotropy of valence electrons is fundamentally intertwined with a variety of emergent quantum phenomena, including charge and orbital order(1–5), unconventional superconductivity(6, 7), quantum spin liquid(8), and multipole order(9). In such systems, orbital configuration directly governs the electronic states and physical properties. However, since external fields that couple directly to orbital degrees of freedom are limited, orbital states are considerably less accessible to direct experimental observation than lattice distortions or magnetic structures. Traditionally, orbital states have been inferred indirectly from Jahn–Teller distortions, which manifest as structural deformations of metal–ligand polyhedra through electron-phonon coupling. While such structural information generally provides valuable insight in strongly correlated electron systems, it is rarely sufficient on its own. Electron correlation, relativistic spin–orbit interaction, and finite-temperature effects act on comparable energy scales, obscuring the simple correspondence between lattice distortions and orbital occupations."

"A prototypical example in which orbital degrees of freedom give rise to exotic electronic states is the spinel-type iron vanadium oxide FeV2O4, which crystallizes in the space group Fd3̅ m at high temperatures [Fig. 1A]. In this compound, each Fe2+ ion (3d6) is coordinated by four O2- ions to form a FeO4 regular tetrahedron, while each V3+ ion (3d2) is surrounded by six O2- ions to form a VO6 octahedron with trigonal distortions. In contrast to AV2O4 (A = Mg, Zn, Cd, Mn)(10–15), both the Fe and V ions host orbital degeneracies in the high-temperature cubic phase in FeV2O4."

"Figure 1B shows the changes in lattice deformations across the successive phase transitions in FeV2O4(15–18). Hereafter, a, b, and c denote the lattice parameters in the standard crystallographic setting of the cubic phase, whereas a′, b′, and c′ are used as a unified notation, referenced to the cubic phase, for all temperature-dependent phases. As shown in Fig. 1B, upon cooling, the cubic lattice first transforms to tetragonal with I41/amd symmetry at Ts = 140 K, in which the lattice is contracted along the c′-axis. Below TN1 = 110 K, collinear ferrimagnetic (C-FM) order emerges, with Fe2+ and V3+ moments aligned antiparallel to each other [Fig. 1C], accompanied by a further lowering of lattice symmetry from tetragonal to orthorhombic, with space group Fddd. The ferrimagnetic moment appears along the longest a′-axis [Fig. 1B]."

"Below TN2 = 70 K, the structure reverts from orthorhombic to tetragonal with space group I41/amd, accompanied by a transition to a noncollinear ferrimagnetic (NC-FM) state(15–18). In contrast to the high-temperature tetragonal phase, the low-temperature tetragonal lattice is elongated along the fourfold axis, unlike in other spinel vanadates AV2O4 with orbitally inactive A ions(10–15). Note that Fig. 1B depicts the low-temperature tetragonal phase using a face-centered tetragonal cell with a′ as the unique axis. The Fe spins remain collinear, whereas the V spins are tilted by approximately 55 degrees Fig. 1C. In the two lower-temperature phases, the spin arrangements of Fe and V are strongly coupled to their orbital states(15, 16, 18)."

"In this study, we address this problem by applying external magnetic fields to a smaller crystal across the magnetic phase transitions. In the low-temperature phases, this approach strongly biases the domain populations, reducing the number of populated domains to two or fewer and thereby allowing us to obtain high-quality diffraction data suitable for precise VED analysis. This enables us to follow the temperature evolution of the VED across the successive phase transitions and to establish a comprehensive experimental picture of how the Fe and V orbital states are reconfigured through the successive phases. Comparison between the experimental VED and spin-polarized density-functional-theory (DFT) calculations further allows us to distinguish between nearly degenerate microscopic solutions for the ground-state orbital configuration and to clarify their relation to the spin structure [Figs. 1D and 1E]. More broadly, our results show that experimentally determined VED provides a practical route to identifying orbital states and resolving their coupling to spin and lattice degrees of freedom in strongly correlated electron systems."

"VED in the Cubic phase: Figure 2A shows the VED distribution at 160 K (Cubic phase), obtained from the CDFS analysis. As shown by yellow iso-density surfaces of 3.9e/Å3, the VED distribution appears isotropic around the O2- ions, which is consistent with the 2s22p6 electronic configuration. On the other hand, weak and pronounced anisotropies are observed around the Fe and V sites, respectively, as shown in the zoomed-in views [Figs. 2B and 2C]. Since the FeO4 forms a regular tetrahedron, the Fe 3d orbitals split into a lower-lying e doublet and a higher-lying t2 triplet. In the high-spin configuration for Fe2+ (3d6), the five majority-spin electrons fully occupy the 3d orbitals, while the remaining minority-spin electron resides in degenerate e orbitals (3z2–r2 and x2–y2) with equal probability [inset, dashed square in Fig. 2B]. Since x-ray diffraction probes orbitals as time- and space-averaged states, the observed VED anisotropy around the Fe site can be interpreted as a manifestation of this orbital disorder(35). A similar interpretation can be applied to the VED distribution around the V 3+ (3d2) ion. In a regular VO6 octahedron, the 3d orbitals split into a higher-lying eg doublet and a lower-lying t2g triplet. However, the VO6 octahedron is slightly elongated along the local three-fold axis, satisfying the.3̅ m symmetry at the V site. As a result, the t2g orbitals are further split into a1g and e'g orbitals. If considering only the crystal field from the adjacent six O2- ions forming the octahedron, the e’g orbitals would lie lower in energy than the a1g orbital. The influence of more distant ions using the Ewald method(39), however, reverses this hierarchy, stabilizing the a1g orbital relative to the e’g orbitals(35). The observed VED anisotropy around the V site confirms the latter situation, where one electron occupies the a1g orbital and the other is equally distributed to the e'g doublet [Fig. 2C]. Starting from this Cubic phase orbital state, we discuss how successive phase transitions lift the degeneracy of Fe and V orbitals and lead to the ground state."

"FeO4 deformation and VED around the Fe site associated with the phase transitions: Figure 3 shows the distortions of FeO4 and VO6 polyhedra associated with the phase transition, as determined by synchrotron x-ray diffraction. As shown in Fig. 3A, the FeO4 distortions contain two-dimensional normal modes of Q2,Fe and Q3,Fe, which cause splitting of the e orbitals [Fig. 2B]. First, we examine the successive change in FeO4 distortion. Figures 3B and 3C present the temperature dependence of Q2,Fe and Q3,Fe for the FeO4 tetrahedron. Orange iso-density surfaces in Fig. 3C show the minority-spin electron densities of Fe predicted from the local FeO4 distortions in each phase. In the Cubic phase, the Q2,Fe and Q3,Fe modes are absent. In the HT-tetra phase at 130 K, negative Q3,Fe coupled to the 3z2–r2 orbital emerges, followed by the emergence of Q2,Fe in the Ortho phase at 100 K. Here we consider the domain with Q2,Fe > 0. In the LT-tetra phase, the FeO4 tetrahedron is elongated along the x-axis and couples with the y2–z2 orbital. To directly elucidate Fe 3d orbital states, we examine the VED distribution. Figure 4A shows the temperature evolution of the VED around the Fe site. In the Cubic phase, high-density regions, shown by orange, are located on the x-, y-, and z-axes. In the HT-tetra phase, enhanced VED develops along the z-axis reflecting the splitting of the e orbitals and minority-spin occupation of 3z2–r2 orbital. Upon further cooling into the Ortho phase, the VED along the y-axis increases, described by a linear combination of the 3z2–r2 and y2–z2 orbitals. In the LT-tetra phase, high-density regions extending within the yz-plane are observed, corresponding to the stabilization of the y2–z2 orbital. Figure 4B presents the calculated VED distribution for the minority-spin of Fe, which reproduces the anisotropy observed in the high iso-density surfaces shown in orange [Fig. 4A]. These results are consistent with those expected from the distortion of the FeO4 tetrahedron [Fig. 3C]."

"The successive changes in Fe 3d orbital state can be understood in terms of the cooperative Jahn–Teller effect below TS and the spin–orbit coupling in the magnetic phases, which stabilize the 3z2–r2 and y2–z2 orbitals, respectively(15, 16, 18, 40). The effective Hamiltonians for the Jahn–Teller coupling and the spin-orbit coupling in the FeO4 tetrahedron are given by HJT=A(τxQ2,Fe−τzQ3,Fe), and HSOC = B/6[(3S2z − S2)τz −√3(S2x − S2y)τx] (41), respectively. Here, A (<0) is an electron-phonon coupling constant, τz and τx are pseudospin operators representing 3z2–r2 (τz = 1/2) and x2–y2 (τz = -1/2) orbitals. S, Sx, Sy, and Sz are spin operators, B (>0) represents a second-order perturbation of spin–orbit coupling. Firstly, we examine the HTtetra phase. Focusing on the anharmonic term of the lattice elastic energy, U = A3Q3cos(3θ), the compressed tetragonal distortion is stabilized. Here, A3 (>0) is the third-order elastic constant, Q = (Q2,Fe2 + Q3,Fe2)1/2, and θ is the polar angle of the distortion vector (Q3,Fe, Q2,Fe) in the Q2,Fe- Q3,Fe plane. This induces cooperative Jahn–Teller distortion, stabilizing the Fe 3z2-r2 orbitals. At lower temperatures where magnetic order emerges, the system undergoes a further structural transition. Under second-order spin-orbit coupling, the favored orbital depends on the spin operation: the y2−z2 orbital is stabilized when the spin moment is aligned along the x-axes. When Fe2+ magnetic order emerges at TN1, the orbital state is therefore modulated by spin–orbit coupling. As the spins align along the x-axis below TN1 in the Ortho phase, the y2−z2 orbital is stabilized by spin–orbit coupling, leading to a hybrid of 3z2–r2 and y2–z2 orbitals(15). Consequently, Q2,Fe becomes non-zero through Jahn–Teller coupling, resulting in an orthorhombic distortion. Here, Q2,Fe > 0, as we adopt the orthorhombicity of a > b > c. Indeed, in FeCr2O4, Cr3+ ions without orbital degrees of freedom occupy the octahedral sites, the ground state remains an orthorhombic phase with conical order of Fe and Cr spins(41, 42). This suggests that the cooperative Jahn–Teller effect and spin–orbit coupling at the Fe site account for the tetragonal-to-orthorhombic transition. However, the reentrant tetragonal symmetry in the ground state of FeV 2O4 cannot be explained within this Fe-only framework, indicating an essential role of the orbital degrees of freedom at the V sites."

"VO6 deformation and VED around the V site associated with the phase transitions: The symmetric distortion of a VO6 octahedron can be decomposed into Eg (Q2,V, Q3,V) and T2g (Q4,V, Q5,V, Q6,V) modes [Fig. 3D]. In the Cubic phase, Q4,V = Q5,V = Q6,V > 0 and their linear combination Qt1 v = (Q4 v + Q5 v + Q6 v)/√3 corresponds to elongation or compression of the VO6 octahedron along the [111] axis. This trigonal distortion leads to the splitting of the t2g into a1g and e’g orbitals. Figure 3E shows the temperature dependence of each Q mode of the VO6 octahedron. The Q2,V and Q3,V values, which directly couple to the splitting of the eg orbitals, remain nearly zero, whereas Q4 v, Q5 v, and Q6 v values exhibit pronounced temperature dependence, reflecting their coupling to the splitting of the t2g orbitals. Notably, the large Qt1 V value is nearly temperature-independent, indicating that one electron continues to occupy the a1g orbital down to the lowest temperatures. To further probe the splitting of the e’g orbitals, we examine Q(z)t2 V = (−Q4 v − Q5 v + 2Q6 v)/√6 and Q(z)t3 V = (−Q4 v + Q5 v)/√2 modes, which are orthogonal to Qt1 V [Fig. 3F]. In the Cubic phase at 160 K, both Q(z)t2 V and Q(z)t3 V values are zero, indicating that the e’g orbitals remain degenerate. In the HT-tetra phase at 130 K, a positive Q(z)t2 V mode emerges while Q(z)t3 V remains zero, corresponding to tetragonal distortion. Subsequently, in the Ortho phase at 100 K, a positive Q(z)t3 V mode appears. Upon entering the LT-tetra phase, Q(z)t3 V gradually increases. We emphasize that the value at 50 K lies along the Q(x)t2 V = (−Q5 v −Q6 v + 2Q4 v)/√6 axis when the Q4,V direction is taken as the principal axis, corresponding to Q(x)t3 V = (−Q5 v +Q6 v)/√2 = 0. The continuous evolution from positive Q(z)t2 V to negative Q(x)t2 V may reflect modifications of the 3d wave functions."

"Next, we turn to the VED distribution around the V site [Fig. 5A], which also continuously changes upon cooling. Here, we define t2g wavefunctions as Ψ1 = C1yz⟩+C2zx⟩+C3xy⟩, Ψ2 = D1yz⟩+D2zx⟩+D3xy>, and Ψ3 = E1yz⟩+E2zx⟩+E3xy>, with C1 2+C2 2+C3 2=1, D1 2+D2 2+D3 2=1, and E1 2+E2 2+E3 2=1, where the coefficients are chosen to satisfy the site-symmetry and mutual-orthogonality constraints (see SI Appendix section3). To elucidate the orbital states, the quantum parameters Ci, Di, and Ei (i = 1 to 3) are optimized by fitting the calculated VED distributions to reproduce the observed anisotropy. Assuming that the trigonal distortion is large enough for one electron to occupy Ψ1, the calculated 3d2 VED is expressed asrhocalc (r) = psi1 2 + ηpsi2 2 + (1 − η)psi3 2, where η represents the electron filling of Ψ2. The evaluation function s for the fitting is defined as s=∑rrhoobs (r) − κ 3 ρcalc (κr)/∑rrhoobs (r), where ρobs (r) is the observed VED around the V site, κ is a variable parameter representing contraction of the VED. The quantum parameters are optimized by minimizing the evaluation function in Eq. (2) within a radial range of 0.25 < r < 0.6 Å around the V site (see SI Appendix Fig. S6). The evaluation function is minimized at approximately η = 0.5 in all the phases."

"Spin-polarized DFT+U calculations in the LT-tetra phase: The 3d2 orbital states of V3+ in the low-temperature phases can be interpreted within two possible frameworks, as shown in Fig. 5D. One possible scenario is that the splitting of the e'g orbitals is sufficiently small that the system resides in a thermally accessible pseudo-degenerate state. Alternatively, the VED with η = 0.5 can also arise from wavefunctions stabilized in magnetically ordered states through relativistic spin–orbit interaction(27). The two scenarios are impossible to distinguish…

As a fastidious and diligent researcher, I have analyzed this paper, Real-space determination of orbital states driving successive phase transitions in FeV2O4. The core scientific contribution is establishing a direct link between experimentally determined valence electron density (VED) anisotropy and the microscopic orbital states that drive coupled structural and magnetic phase transitions in strongly correlated systems.

Here are the specific improvements to AI systems derived from this research, focusing on areas where current computational chemistry/materials science models struggle:

The improved AI system can achieve the following capabilities:

1. A more accurate and robust prediction of ground-state configurations for complex transition metal oxides (TMOs) by incorporating real-space experimental constraints.

2. Automated identification of hidden, nearly degenerate metastable states that are currently computationally inaccessible or ambiguous in purely first-principles calculations.

Here are the specific improvements to AI systems:

1. The system can be trained on the relationship between experimentally measured VED anisotropy maps (derived from synchrotron XRD/CDFS) and theoretical Density Functional Theory (DFT) solutions.

2. The system can perform real-space constraint matching on DFT output, effectively pruning thousands of theoretically possible metastable electronic configurations down to a few physically realizable ground states that match experimental orbital signatures.

3. The system can predict the precise sequence and nature (collinear vs. noncoplanar) of coupled spin/orbital/lattice transitions based on the evolution of the VED anisotropy across different temperature regimes, moving beyond simple energy minimization to capture cooperative effects like Kugel-Khomskii interactions and Jahn-Teller coupling.

The improved AI system can do the following specific tasks:

1. Predict the sequence of structural phase transitions (e.g., Cubic → HT-tetra → Ortho → LT-tetra) in FeV2O4 based on predicted changes in orbital occupation patterns derived from VED evolution.

2. Identify which specific orbital wavefunctions (e.g., stabilizing the 3z2−r2 vs y2−z2 orbitals) are responsible for driving a transition from one structural symmetry to another, linking it directly to spin-orbit coupling strengths and electron-phonon coupling constants derived from DFT+U+SOC calculations.

3. Distinguish between nearly degenerate theoretical solutions (e.g., collinear vs. noncoplanar spin arrangements in the LT-tetra phase) by matching their resulting VED distributions against experimental data, thereby resolving fundamental ambiguities in strong correlation physics that standard energy-based methods cannot resolve alone.

4. Generate a predictive framework for identifying hidden orbital orders in other complex Mott insulators by using experimentally determined VED as a decisive real-space filter to select the physically relevant ground state among competing theoretical solutions.