Spin-state order/disorder and metal–insulator transition in GdBaCo2O5.5: experimental determination of the underlying electronic structure

We have investigated the electronic structure of GdBaCo2O5.5 across the metal–insulator transition (MIT) using soft x-ray absorption and photoelectron spectroscopy. For the low-temperature insulating phase, we find that half of the Co3+ ions at the octahedral sites are in the low spin (LS) and the other half in the high spin (HS) state, while the Co3+ ions at the pyramidal sites are in the HS configuration. Upon increasing the temperature across the MIT, part of the LS octahedral Co3+ undergoes a spin-state transition into the HS configuration. We infer that this destroys the spin-state ordering and thus explains the decrease in resistivity. We observed that the band gap is reduced but not closed in the high-temperature phase.

The layered cobaltates RBaCo 2 O 5.5 (R = rare earth) have attracted considerable attention of the scientific community in recent years because of their peculiar magnetic and transport properties, showing giant magneto-resistance as well as metal-insulator (MIT) and antiferro-ferro-paramagnetic transition phenomena . Moreover, these compounds possess the so-called spinstate degree of freedom, i.e. the possibility of the Co 3+ ions to be in a low spin (LS, S = 0), high spin (HS, S = 2) and even intermediate spin (IS, S = 1) state [28,29]. This spin-state degree of freedom is generally considered to play a key role for the unconventional transport properties of the cobaltates, and a spin-blockade mechanism [13,30] has recently been proposed [13,16] as the driving force for the MIT in RBaCo 2 O 5. 5 .
The crystal lattice of RBaCo 2 O 5.5 is composed of an equal number of CoO 6 octahedra and CoO 5 pyramids and the valence of all the Co ions is 3+. These materials exhibit a series of magnetic orders, namely, a ferromagnetic-like phase below T C ∼ 290 K, an antiferromagnetic phase below T N 1 ∼ 260 K which is accompanied by the appearance of strong magneto-resistive effects and a second antiferromagnetic phase below T N 2 ∼ 200 K. Despite a large number of neutron diffraction studies, the nature of these magnetic orderings is still unresolved [7,[18][19][20]22]. One main reason for the disagreement is the confusion about the spin state of the Co 3+ ions in both octahedral (Co 3+ oct ) and pyramidal (Co 3+ pyr ) coordination. This is particularly evident in the discussion about the sharp drop in the resistivity, the so-called MIT, occurring in the paramagnetic phase at T MI ∼ 360 K. This MIT is commonly attributed to a sudden spin-state switch of Co 3+ oct , but its modality is controversial and intensely debated. Contradictory scenarios have been proposed, including full or partial LS→IS or LS→HS state transitions for the Co 3+ oct and IS or HS configurations for the Co 3+ pyr [1-6-26]. It is noteworthy that at present there is, in fact, no direct evidence for any of the claimed spin-state transitions. There is therefore a pressing need to address the issue of the spin state of the Co ions and its evolution as a function of temperature. We have performed experiments utilizing soft x-ray absorption spectroscopy (XAS) at the Co-L 2,3 and O-K edges, a method proven to be extremely sensitive to the spin state of transition metal ions in general [31,32] and 3d 6 systems in particular [30,[33][34][35]. We have also carried out valence band photoelectron spectroscopy  (PES) measurements in order to determine the changes of the band width and band gap in relationship to the temperature evolution of the spin state. Twinned single crystals of GdBaCo 2 O 5.5 were grown from a high-temperature flux melt. The purity and cation composition of the crystals were checked by x-ray diffraction and fluorescence analysis, respectively. The oxygen content was determined by iodometric titration to be 5.48 (2). Co-L 2,3 spectra in the total-electron-yield mode and O-K XAS spectra in the totalfluorescence-yield mode were measured at the Dragon beamline of the NSRRC in Taiwan with a photon-energy resolution of 0.45 and 0.25 eV, respectively. XAS spectra of CoO were measured simultaneously as energy calibration for the Co-L 2,3 and O-K edges, respectively. The O-K and Co-L 2,3 XAS spectra of two reference compounds LaCoO 3 and Sr 2 CoO 3 Cl were collected under the same conditions. The PES spectra were recorded at the ID08 beamline of the ESRF in Grenoble, using a photon energy of 700 eV and a Scienta SES200 electron analyser. The overall experimental resolution was 110 meV and the Fermi level was calibrated using polycrystalline silver. Figure 1(a) displays the isotropic O-K XAS measured on GdBaCo 2 O 5.5 at different temperatures. The spectral structures from 528 to 532 eV are due to transitions from the O 1s core level to the O 2p orbitals that are mixed with the unoccupied Co 3d t 2g and e g states. In the XAS of GdBaCo 2 O 5.5 we observe a very modest spectral weight transfer from the broad feature at 530 eV to the dominant peak at 528.4 eV. These spectra are the sum of two contributions related to the Co 3+ pyr and Co 3+ oct sites. To explain their line shapes we will make use of Sr 2 CoO 3 Cl as a reference compound containing Co 3+ pyr ions [34,37], and of LaCoO 3 as a reference material having Co 3+ oct ions which undergo a spin-state transition as a function of temperature [35,38]. The O-K XAS of Sr 2 CoO 3 Cl is reported in figures 1(b) and (c) as the curve of blue dots. The Co 3+ pyr ions here have the HS configuration [34,37]. The O-K spectra of LaCoO 3 are displayed in figure 1(d) for T = 20 K (black), 300 K (blue) and 650 K (red). Figure 1(d) shows a significant spectral-weight transfer from the higher energy structure at around 529.4 eV to the lower energy one at around 528.4 eV when going from 20 to 300 K and further to 650 K. The LaCoO 3 data demonstrate that O-K XAS is also very sensitive to the spin-state transition of Co 3+ ions, which can be explained as follows. At low temperatures, the LS Co 3+ oct has its t 2g shell completely occupied [35] and only transitions to the higher lying empty e g states are possible. At higher temperatures, with part of the Co ions in the HS state [35], the t 2g states become partially unoccupied; then transitions to lower lying t 2g states are also allowed, leading to the feature at around 528.4 eV. As the temperature increases this lower energy feature gains gradually spectral weight at the expense of the higher energy feature. From the Co-L 2,3 XAS, it was estimated that at 650 K, 50% of the Co ions are in the HS state and 50% in the LS [35].
In determining the spin-state configurations for GdBaCo 2 O 5.5 at 78 K, we assume that its Co 3+ pyr ions are in the HS state like in Sr 2 CoO 3 Cl and we start by investigating the LS scenario for the Co 3+ oct ions as proposed in several studies [12,13,16]. To this end we sum the Sr 2 CoO 3 Cl spectrum (blue) and the LaCoO 3 at 20 K (black). The result is displayed in figure 1(c) (magenta). One can clearly observe that this is very different from the GdBaCo 2 O 5.5 spectrum in figure 1(a). The LS Co 3+ oct scenario can thus be safely ruled out. Next, we take the LaCoO 3 at 650 K (red) as ansatz for the Co 3+ oct part. The sum is shown in figures 1(b) and (a) (green). Surprisingly, this sum reproduces the GdBaCo 2 O 5.5 spectrum at 78 K very well. This strongly suggests that GdBaCo 2 O 5.5 in the low-temperature phase has 50% of its Co 3+ oct ions in the LS state and 50% in the HS. All of its Co 3+ pyr ions are HS. Remarkable is also the very modest temperature dependence in the GdBaCo 2 O 5.5 spectra, i.e. very much unlike the LaCoO 3 case, suggesting that only a small part of the Co oct ions participate in a spin-state transition across the MIT.
More support for the above results can be found from Co-L 2,3 spectroscopy. Figure 2(a) depicts the isotropic Co-L 2,3 XAS spectrum of GdBaCo 2 O 5.5 at 78 K. The Co 2p core-hole spin-orbit coupling splits the spectrum roughly into two parts, namely the L 3 (at 780.5 eV) and L 2 (at 795 eV) white line regions. The strong absorption at 5 eV above the Co-L 2,3 white lines is due to the Ba-M 4,5 signal (shaded area) as shown by the BaFeO 3 spectrum in figure 2(b). We now use this BaFeO 3 spectrum to remove the Ba-M 4,5 signal from the GdBaCo 2 O 5.5 Co-L 2,3 spectra. The result is given in figure 3(a), which displays the net Co-L 2,3 spectra of GdBaCo 2 O 5.5 taken at 78, 285 and 400 K. The set is at first sight similar to those reported earlier [23,24], but it is in fact essentially different in detail. Firstly, our spectra do not show a shoulder at 3 eV below the main peak at the Co-L 3 edge, which would represent a significant amount of Co 2+ in the sample. This indicates that our spectra are free from Co 2+ impurity signal [30]. Secondly, our set includes the very important low-temperature (78 K) spectrum. To analyse the Co-L 2,3 spectra, we have performed theoretical simulations based on the multiplet cluster method [34,35] using the XTLS 8.3 code [36], taking into account the real local environment of the Co 3+ ions in GdBaCo 2 O 5.5 , which we assume to be the same as in TbBaCo 2 O 5.5 [18]. The simulations for the octahedral HS (orange), octahedral LS (black) and pyramidal HS (light blue) Co 3+ clusters are shown in figure 3(e). We first make a sum of these contributions in a 1:1:2 ratio and we also add a standard edge jump (dashed curve) to take into account transitions to continuum states 12 . Figure 3(b) clearly shows that this sum agrees excellently with the experimental 78 K spectrum. This confirms our O-K XAS finding, namely that the low-temperature phase of GdBaCo 2 O 5.5 consists of 50% LS-50% HS Co oct . 13 To double check, we also make a sum in a 0:2:2 ratio to mimic the commonly accepted scenario for the low-temperature phase in which all Co oct are believed to be in the LS state. Figure 3(c) reveals unambiguously significant discrepancies from the experiment as marked by ellipses. We can therefore safely reject this scenario. Figure 3(a) shows also that there is only a very small temperature dependence in the Co-L 2,3 spectra of GdBaCo 2 O 5.5 , in strong contrast to, for example, LaCoO 3 [35]. This observation indicates again a very modest and only partial spin-state transition of the Co oct ions in GdBaCo 2 O 5.5 across the MIT, fully in agreement with the findings of the O-K XAS study above. Our data therefore also rule out the commonly accepted scenario of a sharp and massive LS-HS transition across T MI . For completeness, we have investigated what the spectral lineshape of an all HS Co oct high-temperature phase could look like, and we can see from figure 3(d) that the 2:0:2 sum strongly disagrees with the experimental 400 K spectrum.
Recent diffraction data showed that in going from high to low temperature, the MIT is accompanied by a structural transition towards a distorted structure with two inequivalent Co 3+ oct sites, one with the average Co-O bond length of 1.928 Å, called Co 3+ oct1 , and the other with 12 See the appendix for details of the cluster calculations. 13 Earlier work proposed the 100% LS Co 3+ oct scenario, claiming that this is necessary to reproduce the weak feature C (see figure 2) [23]. We find that feature C stems from the 3 P 1 final state in the Ba-M 4,5 spectrum and not from the Co L 2,3 signal, see the BaFeO 3 spectrum in figure 2.  1.978 Å, called Co 3+ oct2 [18,22]. This transition was interpreted as the effect of an orbital [17,23] or a spin-state order transition [22]. Our results showing a 50:50 mixture of LS/HS Co 3+ oct ions provide direct evidence that these low temperature distortion and superstructure are caused by LS/HS spin-state ordering of the Co 3+ ions at the octahedral sublattice. Also a direct link can then be made between the differences in Co-O bond lengths with differences in the spin state [35]. In LaCoO 3 , for example, the bond length is 1.925, 1.949 and 1.962 Å at 5, 650 and 1000 K, respectively [39]. The Co 3+ oct1 ion in GdBaCo 2 O 5.5 has a bond length close to LaCoO 3 at 5 K, i.e. LS. The Co 3+ oct2 ion has a bond length even larger than LaCoO 3 at 1000 K, i.e. HS, again fully consistent with the above finding.  We now investigate the effect of the spin-state configurations on the band width and band gap of GdBaCo 2 O 5.5 . Figure 4 shows the valence band PES at temperatures below and above T MI . Surprisingly, and in contradiction with a previous report 14 , the valence band PES of GdBaCo 2 O 5.5 reveals that the spectral weight at the Fermi level is negligible at 78 K (black) and 285 K (blue), and still very weak at 400 K (red), so that we have to conclude that actually the material is an insulator or semiconductor even above T MI . What we found is that the band gap is indeed reduced when the temperature is increased: the top of the valence band (figure 4, left panel) moves up by about 60 meV in going from 285 to 400 K, while the bottom of the conduction band (figure 4, right panel) shifts downwards by 70 meV from 285 to 400 K. So, in total, the band gap is reduced by about 130 meV across T MI . Using ρ ∝ e E g /k B T MI with k B T MI ∼ 30 meV, this gap narrowing accounts well for the observed reduction of resistivity by two orders of magnitude at T MI [5,9,11,14,16]. From figure 4 one can see that there is no pronounced LS peak at 0.8 eV as found in the spectrum (magenta) of LaCoO 3 at 65 K, where the Co 3+ ions are basically in the LS state [35]. However, we see that the weak feature at 0.6 eV in GdBaCo 2 O 5.5 corresponding to 50% LS Co 3+ oct ions loses spectral weight with increasing temperature.
The following picture now emerges from our electronic structure measurements concerning the MIT in this system. In the low temperature phase, the superstructure associated with the 50% LS-50% HS spin-state ordering of Co 3+ oct ions stabilizes the insulating state of the GdBaCo 2 O 5.5 system. On increasing the temperature across T MI , part of the LS Co 3+ oct ions undergo a transition to the HS state, thereby breaking down the superstructure as is shown by x-ray diffraction [17]. This loss of order causes a reduction of the band gap. In converting an HS-LS neighbouring pair into an antiferromagnetic HS-HS pair by a rise in temperature, one can transfer an e g electron both ways instead of only from HS to LS, thereby increasing the amount of charge fluctuations. However, the band gap does not vanish since the hopping of an extra hole or electron is energetically quite difficult in such an antiferromagnetic situation. Reversely, the observation that the changes in widths of the valence (PES) and conduction (O-K XAS) bands are very modest reaffirms that the IS state scenario does not become active, since otherwise the presence of IS Co 3+ oct ions would allow for free propagation of an extra hole or electron in a ferromagnetically aligned e g band.
To summarize, we have utilized O-K and Co-L 2,3 XAS spectroscopy to reveal that the lowtemperature insulating phase of GdBaCo 2 O 5.5 is characterized by a LS-HS spin-state ordering in the Co 3+ oct sublattice. We found a very modest spin-state transition in the Co 3+ oct with increasing temperature across the insulator-metal transition. Our photoemission experiments show that the band gap is reduced but does not close, in quantitative agreement with the change in resistivity. The MIT should be modelled in terms of a spin-state order-disorder transition with a modest gain in band width.

A.2. Cobalt L-edge data analysis and model calculation (figures 2 and 3)
The total electron yield was divided by I 0 (the intensity of the incoming beam). Then a (constant) pre-edge background was subtracted. As normalization we divided each background-subtracted spectrum by the average of its spectral intensity in the energy range 807-815 eV.
The energy of each GdBaCo 2 O 5.5 spectrum was calibrated using a (simultaneously recorded) reference CoO spectrum identifying the first peak of CoO to be at 778.1 eV [40]. For each temperature we averaged the E ab and E c spectra with a weighting ratio 2:1 in order to obtain the isotropic spectrum. In order to subtract the Ba M-edge signal we measured a polycrystalline BaFeO 3 spectrum. This BaFeO 3 spectrum has a background with a slope; the pre-edge background was fitted using a linear function (for the fit we considered the first 40 points) and we subtracted this linear background from the spectrum. Then the BaFeO 3 spectrum was subtracted from the GdBaCo 2 O 5.5 spectra-the scaling factor and energy position of the BaFeO 3 spectrum were adjusted such that the irregularities in the difference spectrum were minimized.
A theoretical calculation of the x-ray absorption spectra was performed with the program XTLS 8.3 [42]. The following parameters have been used (HF refers to Hartree-Fock calculations for a single Co ion): • Slater integrals were reduced to 80% of the HF values 15 .
• For the 3d spin-orbit coupling HF values were used.
• Hybridization parameters of Co d-electrons and oxygen ligand p-electrons were estimated according to the distance to the power of 7/2 law by Harrison [44], starting from V pdσ = −1.7 eV for a Co-O distance of 1.925 Å [43]. • The ionic crystal field was taken into account by a purely octahedral splitting parameter 10Dq.
• For the final state we used the same values of the hybridization parameters and reduced the ionic crystal field parameter 10Dq by 30% with respect to the initial state value.
• Ionic CF parameter of the initial state: 10Dq = 0.5 eV.
• Ligand positions [45] and the corresponding hybridization parameters V pdσ , V pdπ : The results of the model calculation are compared with the experimental data in figure 3 of the main text. The sum of oct-HS, oct-LS and pyr-HS contributions with the 1:1:2 spectral weight agrees best with the experimental data and supports our conclusion about the presence of a LS-HS spin-state ordering in the Co 3+ oct sublattice of GdBaCo 2 O 5.5 .

A.3. PES data analysis (figure 4)
The valence band PES spectra of GdBaCo 2 O 5.5 were calibrated in energy by fitting a Gaussian broadened Fermi function to Ag reference spectra taken under the same conditions immediately after and before. The Fermi energy is shown by a vertical line in figure 4. A constant background has been subtracted from each spectrum. For normalization, we divided each backgroundsubtracted spectrum by the average of its spectral intensity in the energy range 2-3 eV.