Jahn-Teller distortion driven magnetic polarons in magnetite.

The first known magnetic mineral, magnetite, has unusual properties, which have fascinated mankind for centuries; it undergoes the Verwey transition around 120 K with an abrupt change in structure and electrical conductivity. The mechanism of the Verwey transition, however, remains contentious. Here we use resonant inelastic X-ray scattering over a wide temperature range across the Verwey transition to identify and separate out the magnetic excitations derived from nominal Fe2+ and Fe3+ states. Comparison of the experimental results with crystal-field multiplet calculations shows that the spin-orbital dd excitons of the Fe2+ sites arise from a tetragonal Jahn-Teller active polaronic distortion of the Fe2+O6 octahedra. These low-energy excitations, which get weakened for temperatures above 350 K but persist at least up to 550 K, are distinct from optical excitations and are best explained as magnetic polarons.

Since its first X-ray structural elucidation by Bragg1 a hundred years ago and the discovery of the Verwey transition23, magnetite (Fe3O4), has received much attention for decades. Even today, it attracts significant scientific and technological interest for its applications in ultrafast magnetic sensors4, palaeomagnetism5, nanomedicine carriers6, and so on. Fe3O4 becomes ferrimagnetic below Tc∼850 K, followed by an abrupt decrease in its electrical conductivity by two orders of magnitude as the temperature is cooled below TV. In this first known magnet to mankind, one-third of Fe sites, termed A-sites, are tetrahedrally (Td) coordinated with oxygens; the other two-thirds, termed B-sites, have octahedral (Oh) coordination. Verwey first suggested a Fe2+–Fe3+ charge-ordering occurring on the B-sites as the driving force of this transition. There are two major schools of interpretation: the first one interprets the Verwey transition as a transition driven by charge/orbital ordering7891011121314151617. The second one exploits the mechanism of a lattice distortion-driven electron–phonon coupling1819202122 enhanced by the on-site Coulomb interaction and thus opens a gap at the Fermi energy when the temperature is lowered below the Verwey transition temperature TV.

Although numerous investigations have been carried out to verify the charge localization on the B-sites, the charge-ordering pattern of magnetite is subtle and still elusive1920. While it is agreed that the charge disproportionation involves changes in the nominal Fe2+ and Fe3+ states associated with the B-sites, X-ray diffraction studies of the low-temperature phase of magnetite microcrystals1517 revealed that the t2g electrons of the B-sites are not fully localized in the form of Fe2+ states. Instead, the electrons are distributed over linear three-Fe-site units termed trimerons, which are coupled to the Td Jahn-Teller distortion of B-site Fe2+O6 octahedra, as illustrated in Fig. 1. To the first approximation, the B-site Fe3+O6 octahedra are Jahn-Teller-inactive. The tetragonal distortion of B-site Fe2+O6 octahedra removes the degeneracy of t2g orbitals, in going from Oh symmetry to D4h symmetry. In the absence of spin–orbit coupling, an effective energy separation between d and d/d is created if the four Fe–O bonds in the xy plane are elongated or contracted. The trimeron scenario then indicates that the Verwey transition is essentially due to an ordering of trimerons. Because previous results of optical conductivity23 and photoemission24252627 showed the pseudogap feature of magnetite above TV, and results of entropy analysis28, neutron/X-ray diffuse scattering29 and anomalous phonon broadening22 revealed the short-range order above TV, one important open question is whether trimeron correlations persist in the cubic phase at temperatures above TV. Combining these short-range correlations of polaronic characters with the spin degrees of freedom of t2g electrons, one can expect magnetic polarons in magnetite.

Figure 1

Trimeron scenario and t2 energy-level splitting.

Illustration of the orbital ordering of B-site Fe2+ in Fe3O4 and the corresponding t2g energy-level splitting for a Fe2+ ion in a negative crystal field. A trimeron is indicated with a dashed oval. The elongation of the four Fe–O bonds in the xy plane are indicated with arrows.

Here we present measurements of resonant inelastic X-ray scattering (RIXS)3031 at the Fe L3-edge on magnetite to reveal the low-energy spin–orbital excitations of Fe2+ ions in both the monoclinic and cubic phases. To the best of our knowledge, the magnetic excitations derived from the local tetragonal distortion field of Fe2+ ions, that is, magnetic polarons, have not been reported to date. In combination with crystal-field multiplet calculations, we show the the existence of magnetic polarons in magnetite which is driven by the Jahn-Teller distortion.

Results

Fe L 3-edge RIXS

Figure 2a shows the Fe L-edge X-ray absorption spectrum of magnetite. By comparing with crystal-field multiplet calculations (see Supplementary Fig. 1), it is understood that the absorption-energy centroid of Fe2+ ions is lower than that of Fe3+ ions by ∼1–2 eV, consistent with earlier work323334. Accordingly, the features at X-ray energies of 706.0 and 707.5 eV originate from the absorption of octahedrally coordinated B-site Fe2+ states, while the maximum intensity feature at 708.8 eV is dominated by absorption from the the Fe3+ ions of both the B-site octahedral and A-site tetrahedral symmetries.

Figure 2

RIXS measurements of Fe3O4.

(a) Fe L-edge X-ray absorption spectrum (XAS) spectrum measured in the fluorescence yield mode through the summation of all inelastic X-ray intensities taken at room temperature T=300 K. The XAS is plotted with correction for self-absorption. The incident X-ray energy resolution was 0.5 eV. (b) Colour map of RIXS intensity after correction for self-absorption in the plane of incident photon energy versus energy loss recorded at T=80 K. (c–e) RIXS spectra plotted in terms of energy loss with a vertical offset for clarity. They were recorded by using π-polarized incident X-rays under the scattering geometry of the scattering angle ϕ=90° and the incident angle ϕ=20°. Panels c,d were measured at 80 K and e was measured at 550 K.

The colour map of RIXS intensity measured at 80 K in the plane of incident photon energy versus energy loss shown in Fig. 2b presents the evolution of the RIXS spectral profile associated with Fe2+ and Fe3+ ions as detailed in the following. When the incident X-ray energy was set to below 707.5 eV, we observed dd excitations of Fe2+ with energy losses at 2.8±0.05, 1.65±0.05 and 1.16±0.05 eV shown in Fig. 2c, and also a broad excitation centred at 200 meV shown in Fig. 2d. If the incident X-ray energy goes beyond 707.5 eV, the 1.16-eV dd excitation of Fe2+ begins to evolve into a fluorescence that has a constant X-ray emission energy independent of incident energy. With the incident X-ray energy set to 708.8 eV, RIXS excitations arise mostly from Fe3+ ions of octahedral or tetrahedral symmetry.

Figure 2d shows two RIXS features centred at 90 and 200 meV in a magnified plot of energy loss below 0.7 eV. Measurements carried out by varying the scattering angle suggested that these two low-energy excitations do not disperse in momentum space (see Supplementary Fig. 2). The 200-meV excitation has a full-width at half-maximum larger than the instrumental energy resolution. This broad RIXS feature resonates near the L3-edge of Fe2+ and almost disappears for incident energy above 708 eV, at which the other excitation centred at 90 meV emerges. The 90-meV excitation has a full-width at half-maximum nearly equal to the instrumental energy resolution and resonates at 708.4 eV. The distinct incident X-ray energies for these resonant excitations indicate that the 200- and 90-meV features arise from Fe2+ and Fe3+ states, respectively.

Many experimental353637 and theoretical3038394041 studies have shown that L-edge RIXS allows spin–flip processes that are not accessible with optical spectroscopy23. For example, if both the incident and scattered X-rays are π-polarized, the spin–flip excitation of d is allowed owing to the spin–orbit coupling in the 2p core state. In the present RIXS measurements with a 90°-scattering geometry (see Supplementary Fig. 2a), the intensity of elastic excitation with incident X-rays of π polarization is reduced in comparison with that of σ polarization, and spin–flip excitations are effectively revealed. In addition, the cross-section of L3-edge RIXS for a magnetic excitation is usually larger than that for a phonon excitation35. Because O K-edge RIXS-probe excitations derived from bimagnons42 and phonons43, we performed O K-edge RIXS measurements to probe the phonons of Fe3O4. Our data shown in Supplementary Fig. 4 reveal an excitation at 70 meV in the O K-edge RIXS, indicating that the observed 90-meV feature of the Fe L3-edge RIXS has a small contribution from phonon excitation4445.

Multiplet RIXS calculations

In order to characterize the origin of the observed excitations, we undertook crystal-field multiplet calculations for the B-site Fe3+ and Fe2+ ionic configurations. See the Methods section, Supplementary Note 2 and Supplementary Figs 5 and 6 for calculation details.

Multiplet calculations carried out for the B-site Fe3+ ions under an exchange molecular field of 90 meV, as shown in Supplementary Fig. 5d, explain the observed excitation energy of 90 meV well. This is consistent with the 100-meV Zeeman splitting induced by the molecular field deduced from the Curie temperature of magnetite and the exchange coupling constants46. This spin–flip energy also agrees with the energy of the nearly dispersionless mode at 80–85 meV observed in inelastic neutron scattering4647. For an individual Fe3+ site, the intensity of the 90-meV RIXS feature can change under a spin reorientation as magnetic RIXS is sensitive to the spin direction with respect to the incident polarization. However, RIXS measurements of B-sites Fe3+ reflect an average of eight non-equivalent Fe3+. Each of these non-equivalent Fe3+ B-sites can make different contributions to the intensity of the 90-meV feature, and our measured RIXS spectra suggest that these changes are beyond our experimental sensitivity. Although we cannot totally rule out phonon contributions, the 90-meV L3-edge RIXS excitation is best explained as a result of spin–flip excitations of Fe3+ ions, like magnetic excitations observed in the RIXS of Fe pnictide superconductors36, cuprates303740 and nickelates3841.

The ground state of the octahedral Fe2+ ion is a high-spin 5T2g state with S=2. According to Hund’s rule, out of the six 3d electrons of the Fe2+ ion, five 3d electrons occupy spin-up states ; the remaining one electron occupies one of the three spin-down orbitals . When the spin–orbit effect of 3d electrons couple a pseudo-orbital angular momentum to S=2, the 5T2g state splits into three manifolds of pseudo-angular momenta and 3. That is, there are effectively 15 separate states from Fe2+, as the 5T2g ground state is split by the combination of these interactions. For the broad 200-meV RIXS feature associated with the octahedral Fe2+ states, the excitation energy is too large to be explained in terms of spin–flip excitations only. We carefully checked the effect of the local Jahn-Teller distortion to explain the energy of the observed excitations to understand the nature of this feature. From an extensive set of RIXS calculations of Fe2+ with varied tetragonal distortions as shown in Supplementary Fig. 6, we found that the average RIXS spectrum calculated using Hex=90 meV and =−22, −26 and −30 meV explains the measured spectrum most satisfactorily, as demonstrated in Supplementary Fig. 5f. Figure 3a shows the calculated low-energy RIXS excitations of Fe2+ in the form of incident photon energy versus energy-loss maps. The calculated RIXS obtained by including the tetragonal distortion, exchange interaction and 3d spin–orbit coupling matches fairly well with the experimental data.

Figure 3

Calculated RIXS of Fe2+ in comparison with measurements.

(a) Calculated RIXS intensity map of B-site Fe2+ by using Hex=90 meV,  meV and the 3d spin–orbit coupling  meV. The resonance photon energy of Fe2+ is set to the experimental resonance energy 707.5 eV. The core–hole lifetime width is set to 200 meV, and the final-state lifetime width is set to 10 meV. This calculated intensity map presents the average RIXS intensity for the magnetic easy axis along the [100], [010] and [001] directions, and are plotted after Gaussian broadening of width 500 and 80 meV for the incident photon energy and the energy loss, respectively. (b) A magnified intensity map of RIXS measurements of single-crystal Fe3O4 extracted from Fig. 2b. (c) Comparison of measured (expt) and calculated (calc) RIXS spectra. Open circles are measurements with incident X-rays of 707 eV at 80 K; the solid line presents the calculated RIXS spectra of incident X-ray energy 707.5 eV.

Discussion

In comparison with the magnified intensity map of RIXS measurements shown in Fig. 3b, calculations using a molecular field Hex=90 meV and =−26±4 meV reproduce the energy-loss features suitably. For the 200-meV excitation, the experimental resonance starts at an energy lower than that of the resonant quasi-elastic scattering of Fe2+ and its energy range is broad, while the calculated resonance starts at a higher energy with a narrow range. This discrepancy is attributed to differences in the dynamics of 3d orbitals due to core–hole effects in intermediate states. These effects do not affect the energy loss of excitation spectra because the core holes are filled in the RIXS final state. As is typical of RIXS calculations reported in the literature, our calculations do not include such core–hole effects, and hence do not reproduce the incident energy dependence perfectly, but our calculations correctly reproduce the energy-loss features.

Figure 3c presents calculated RIXS spectrum in comparison with measurements of the incident X-ray energy set to 707 eV, at which the 200-meV RIXS feature is most pronounced. The negative value of signifies that the energy of d is lower than that of d/d, that is, tetragonally distorted Fe2+O6 octahedra with elongated Fe–O bonds in the xy plane. This shows that the tetragonal distortion is directly related to a polaronic distortion of the Fe2+O6 octahedra, which in turn couple to the neighbouring Fe3+O6 octahedra constituting the trimerons, although, as mentioned earlier, they are Jahn-Teller-inactive in the first approximation. Our results are consistent with the locally distorted structure of the FeO6 octahedra and the short-range order above TV observed by X-ray absorption48 and diffuse scattering29. These short-range correlations are polaronic in nature2326272949. Because a local molecular field and a tetragonal Jahn-Teller polaronic distortion are required to correctly simulate the excitation energy, the observed spin–orbital excitations are, indeed, magnetic polarons.

The magnitude of obtained is comparable with the 3d spin–orbit coupling strength, and thus confirms the observation of the unquenched orbital moment50, which is known from work on Fe2+ impurities in MgO thin films51. These results are also consistent with conclusions of band-structure calculations using the monoclinic P2/c crystal structure of magnetite9 (see Supplementary Note 3 and Supplementary Fig. 7), which give an energy splitting ∼50 meV between minority-spin d and d/d bands at the Γ point, conforming to the deduced .

We also performed RIXS measurements above the Verwey transition and found that the spin–orbital excitations driven by polaronic distortion do exist in the cubic phase of Fe3O4 at high temperatures as shown in Fig. 2e and Supplementary Fig. 8. Figure 4 plots the temperature-dependent RIXS spectra with the incident X-ray energy set to the pre-edge absorption at 706 eV, an incident X-ray energy at which the elastic component is weak and the RIXS arises predominantly from octahedral Fe2+ ions with a negligible contribution from Fe3+. The temperature-dependent results show that, when the temperature is varied across TV, the spin–orbital excitation of 200 meV does not abruptly change its intensity and persists at least up to 550 K, albeit with a gradual decrease above 350 K. We interpret this as a gradual weakening of the polarons. RIXS results shown here serve as a fast probe to snapshot the dynamic lattice–spin–orbital excitations of Fe3O4. These temperature-dependent RIXS results indicate that the FeO6 octahedra are already locally distorted in the cubic phase of magnetite, in good agreement with the existence of the short-range correlations in the lattice structure above TV. The temperature dependence of these distortions follows that of the magnetization of magnetite, suggesting short-range ordering of the Jahn-Teller distortion, which gets weakened as the temperature approaches the Curie temperature TC, and providing further evidence for magnetic polarons. These observations suggest that the local distortion in the cubic phase could be attributed to the precursor of the monoclinic phase across the Verwey transition.

Figure 4

Temperature-dependent spin–orbital excitations of Fe3O4.

(a,b) RIXS spectra after the subtraction of the elastic component at selected temperatures. The spectra were recorded with the incident X-ray energy set to 706 eV. (c) Plot of the integrated 200-meV RIXS intensity versus temperature. The data were analysed by subtracting the elastic component from the measured spectrum and normalizing to the intensity of the dd excitation feature at 2.8 eV (see Supplementary Fig. 9). The dashed line indicates the Verwey transition temperature TV. The RIXS data comprise an average of four runs of experimental results. The error bars are deduced from the average value of variations in the spectral intensity of four sets of measurements.

To summarize, our results demonstrate the usefulness of RIXS to unravel the local electronic structure of a mixed-valence compound by selecting the energy and polarization of incident X-rays. We revealed dd excitons in magnetite that have an energy centroid 200 meV and arise from polaronic distortion-driven spin–orbital excitations, which are best explained as magnetic polarons. We also applied crystal-field multiplet calculations to obtain the t2 crystal field =−26±4 meV induced by the tetragonal Jahn-Teller distortion. These results are consistent with the mechanism of ordering trimerons for the Verwey transition. It would be interesting to carry out RIXS experiments with an improved energy resolution to study the change of spin–orbital excitations across the Verwey transition.

Methods

RIXS measurements

Using the AGM–AGS spectrometer at beamline 05A1 of the Taiwan Light Source31, we measured RIXS on a single-crystal Fe3O4(001) at incident photon energies set to specific energies about the L3(2p3/2→3d) absorption edge of Fe. See Supplementary Fig. 2a for the scattering geometry. Both the scattering angle ϕ defined as the angle between the incident and the scattered X-rays, and the incident angle θ from the crystal ab plane, were variable. The polarization of the incident X-ray was switchable between π and σ polarizations, that is, the polarization within and perpendicular to the scattering plane, respectively, and the polarization of scattered X-rays was not analysed. The energy bandwidth of the incident X-rays was 500 meV, and the total RIXS energy resolution was ∼80 meV because the energy compensation method was used to ensure a high-resolution measurement in the energy-loss scheme31. The beam diameter of incident X-ray at the sample is ∼0.5 mm.

Sample preparation

Single-crystal growth of magnetite was carried out in an infrared image furnace in high-purity argon gas (99.999% purity) atmosphere. Measurements of the temperature-dependent specific heat and resistivity of the synthesized magnetite crystal showed that it exhibits a sharp first-order Verwey transition at TV=122 K. The synthesized single crystal has a chemical composition of Fe3(1−O4 with |δ|≤0.00018, indicative of a nearly ideal chemical stoichiometry. See Supplementary Figs 10 and 11 for the sample characterization.

Multiplet calculations

We undertook crystal-field multiplet RIXS calculations of B-site Fe2+ and Fe3+ using CTM4RIXS52 and MISSING (Dallera and Gusmeroli http://www.esrf.eu/computing/scientic/MISSING/) with the scattering angle 90° and the magnetization axis perpendicular to the scattering plane or in the scattering plane with angles 20° or 70° to the incident beam. The polarization of incident X-rays was selected to be π-polarized. Since the polarization of scattered X-rays was not analysed in the measurements, we summed calculated RIXS spectra of scattered X-rays with σ and π polarizations. We used a Lorentzian broadening 0.2 eV for the lifetime width of the intermediate states. The calculated spectra are obtained as an average of the spectra calculated for magnetic domains with the easy axis along the [100], [010] and [001] directions. The crystal-field parameter 10Dq was set to 1.13 eV, and the Slater integrals were reduced to 79% of their atomic values for accurately reproducing the dd excitation energies.

Data availability

The data that support the findings of this study are available from the corresponding authors on request.

Additional information

How to cite this article: Huang, H. Y. et al. Jahn-Teller distortion driven magnetic polarons in magnetite. Nat. Commun. 8, 15929 doi: 10.1038/ncomms15929 (2017).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.