Outline
Comptes Rendus

Account / Revue
Exchange and correlation effects in transition-metal oxides 3dn (n = 4, 5 and 6)
Comptes Rendus. Chimie, Volume 6 (2003) no. 1, pp. 135-145.

Abstracts

In order to clarify the competition of exchange and correlation energies vs crystal-field stabilisation energy, a simple approach is proposed. The exchange and correlation contributions are described on the base of Brandow and Kanamori’s U, U′ and JH parameters. The dependence of the crystal field effect on site distortion has simply been modelled using a linear interpolation between undistorted and fully distorted sites. This approach leads to establish phase diagrams for d4, d5 and d6 cations, allowing us to predict the spin-state stability range depending on exchange (JH), crystal field (Dq) and distortion (k) parameters. It can be used to complement the Extended Hückel Tight Binding calculations of the electronic structure of materials having a noticeable ionic character such as 3d transition-metal oxides, in order to interpret the electronic behaviour, and, particularly, discuss the insulating-vs-metallic character of these oxides.

Afin de mieux comprendre la compétition qui oppose les effets d’échange et de corrélation à ceux du champ cristallin, une approche simple est proposée. Les contributions d’échange et de corrélation sont décrites sur la base des paramètres U, U′ et JH de Brandow et Kanamori. L’évolution de la contribution du champ cristallin avec la distorsion du site a été simplement modélisée à l’aide d’une interpolation linéaire entre le site non distordu et complètement distordu. Cette approche conduit à l’établissement de diagrammes de phase pour les cations d4, d5 et d6, permettant de prévoir le domaine de stabilité d’un état de spin en fonction des paramètres d’échange (JH), de champ cristallin (Dq) et de distorsion (k). Elle peut être utilisée pour compléter les calculs de structure électronique de type liaison forte Hückel étendue dans le cas de matériaux ayant un caractère ionique notable, tels que les oxydes d’éléments de transition, dans le but d’interpréter le comportement électronique et, en particulier, de discuter le caractère isolant ou métallique de ces oxydes.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/S1631-0748(03)00014-6
Keywords: exchange, correlation, crystal field, transition metal oxides, electronic structure, metal–insulator transitions, spin-state phase diagrams
Keywords: échange, corrélation, champ cristallin, oxydes d’éléments de transition, structure électronique, transitions métal–isolant, diagrammes de phase d’états de spin

Michel Pouchard 1; Antoine Villesuzanne 1; Jean-Pierre Doumerc 1

1 ICMCB–CNRS, 87, avenue du Docteur-Albert-Schweitzer, 33608 Pessac, France
@article{CRCHIM_2003__6_1_135_0,
     author = {Michel Pouchard and Antoine Villesuzanne and Jean-Pierre Doumerc},
     title = {Exchange and correlation effects in transition-metal oxides 3d\protect\emph{\protect\textsuperscript{n}}~(\protect\emph{n} = 4, 5 and 6)},
     journal = {Comptes Rendus. Chimie},
     pages = {135--145},
     publisher = {Elsevier},
     volume = {6},
     number = {1},
     year = {2003},
     doi = {10.1016/S1631-0748(03)00014-6},
     language = {en},
}
TY  - JOUR
AU  - Michel Pouchard
AU  - Antoine Villesuzanne
AU  - Jean-Pierre Doumerc
TI  - Exchange and correlation effects in transition-metal oxides 3dn (n = 4, 5 and 6)
JO  - Comptes Rendus. Chimie
PY  - 2003
SP  - 135
EP  - 145
VL  - 6
IS  - 1
PB  - Elsevier
DO  - 10.1016/S1631-0748(03)00014-6
LA  - en
ID  - CRCHIM_2003__6_1_135_0
ER  - 
%0 Journal Article
%A Michel Pouchard
%A Antoine Villesuzanne
%A Jean-Pierre Doumerc
%T Exchange and correlation effects in transition-metal oxides 3dn (n = 4, 5 and 6)
%J Comptes Rendus. Chimie
%D 2003
%P 135-145
%V 6
%N 1
%I Elsevier
%R 10.1016/S1631-0748(03)00014-6
%G en
%F CRCHIM_2003__6_1_135_0
Michel Pouchard; Antoine Villesuzanne; Jean-Pierre Doumerc. Exchange and correlation effects in transition-metal oxides 3dn (n = 4, 5 and 6). Comptes Rendus. Chimie, Volume 6 (2003) no. 1, pp. 135-145. doi : 10.1016/S1631-0748(03)00014-6. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/S1631-0748(03)00014-6/

Version originale du texte intégral

1 Introduction

Exchange and correlation phenomena describe the effects of electron–electron interactions beyond the one-electron mean-field approximation; they are necessary to build correct electronic configurations based on an orbital description for 3dn ions. A large part of the correlation energy, in the sense of electron–electron interactions, is included in mean-field approaches: for example, most one-electron schemes account for screening phenomena (i.e. Slater screening constant). In the Hartree–Fock framework, correlation energy is defined restrictively as the difference between the exact and the Hartree–Fock ground-state energy.

Exchange phenomena are more complex, including the symmetry of the wave functions and the Fermi hole, both giving rise to Hund’s rules [1]. In the case of partially occupied orbitals, for example 3d orbitals for transition metals, exchange will contribute to differentiate up- and down-spin populations.

In inorganic molecular chemistry, exchange and correlation are discussed on the basis of Racah A, B, C, or Slater parameters, which allow to quantify electron–electron interactions and result in the Tanabe–Sugano diagrams for octahedral complexes. What about solid-state chemistry?

Tight-binding methods as Extended Hückel (EHTB), which emphasise atomic orbital overlaps, do not include exchange and correlation effects. Moreover, in non self-consistent methods as EHTB, crystal or molecular orbital energies do not depend on the electron count.

Methods based on the Density Functional Theory (DFT) use approximate exchange-correlation functionals (so-called LDA, GGA, WDA...), based on the analytical expression for the exchange-correlation energy of an homogenous electron gas. Although closer to the exact energy than single-determinant Hartree–Fock, DFT remains at the mean-field stage. A step beyond correlation effects can be achieved through the LDA+U method, in which on-site electron–electron repulsion terms are added to the exchange-correlation functional [2].

Chemical bonding in transition metal oxides can be described on the basis of two types of atomic orbitals:

  • • (i) s and p orbitals, for which ζ Slater exponents take low values (ζ4s ∼ 1–1.5 [3]), corresponding to diffuse orbitals;
  • • d orbitals, more localised (ζ3d ∼ 2.5–4.5 [3]), describing strongly interacting electrons, for which exchange and correlation phenomena are often crucial.

As an example, a monoelectronic band description of the three following double oxides with rhombohedral perovskite-related structures would predict a metallic behaviour for LaFeO3 (d5) and LaNiO3 (d7), and an insulating one for LaCoO3 (d6) (Fig. 1). Whereas LaNiO3 indeed is a metal at any temperatures (within its stability range), LaFeO3 actually is always insulating and LaCoO3 undergoes an insulating to metal transition at high temperature. So a simple low-spin band description, neglecting exchange and correlation effects, does not correspond to experimental facts.

Fig. 1

Schematic monoelectronic band structure for a cubic perovskite LaMO3(M = d5, d6, and d7 cation) on the basis of EHTB method.

In LaFeO3, trivalent iron (3d5) is not in the low-spin t2g5eg0 configuration (favoured by the crystal field), but in the high-spin t2g3eg2 configuration (favoured by exchange and correlation). The competition between crystal field (which stabilises some energy levels) and exchange and correlation (which tend to keep electrons apart in space) is at the origin of spin equilibria and spin transition phenomena. At low temperature, the six 3d electrons of cobalt in LaCoO3 lie in the three low-lying t2g orbitals, giving an insulating state; at high temperature, exchange and correlation prevail and spread the six electrons in the five d orbitals. Then LaCoO3 becomes a metal, with the cobalt ion in an intermediate spin state.

In the following, we will describe exchange and correlation effects for transition metal ions in various coordination polyhedra the symmetry of which strongly influences the energy splitting of the five 3d orbitals. We will extend the treatment of most-common octahedral (Oh) and tetrahedral (Td) sites to the square-based pyramidal (C4v), square planar (D4h), dumbbell (Dh), cubic (Oh) and antiprismatic (D4d) environments (Fig. 2).

Fig. 2

Various coordinations and symmetries of cationic sites in oxides.

We will restrict our discussion to d4, d5 and d6 ions, most interesting because they are to be found in three different spin states: S = 0, 1 and 2 for d4 and d6 ions; S = 1/2, 3/2 and 5/2 for the d5 ion. Our approach can then be applied to metal ions such as, for example, Fe4+ and Mn3+ (3d4), Fe3+ and Co4+ (3d5), Co3+ (3d6). Its extension to the 4d and 5d series should be handled cautiously because of the spin–orbit coupling.

The aim of this work is to guide inorganic chemists in the choice of appropriate crystal structures for stabilising a chosen electronic configuration for these 3d transition ions, and to help in the interpretation of their remarkable physical properties.

2 Expression of the exchange and correlation using U, U′ and JH parameters of Kanamori and Brandow

According to Kanamori [4] and Brandow [5], Coulomb repulsion energy between two electrons belonging to the same atomic orbital (U) or to two different orbitals (U′) is balanced by an exchange energy (JH) that stabilises parallel vs antiparallel spins. These parameters are expressed by the following integrals:

U=<μμVμμ>,U ' =<μμ ' Vμμ ' > and J H =<μμ ' Vμ ' μ>
where μ and μ' are different 3d orbitals from a single site and V is the Coulomb operator.

These energies are linked to Racah parameters through the following relationships [5]:

U=A+4B+3C,U ' =A-B+C and J H =5/2B+C

These relationships show that U, U′ and JH are not independent and immediately lead to UU′ = 2 JH, which will allow us to compare the energy of the various spin states using a single parameter, the Hund’s rule exchange energy, JH. For a 3d4 ion in an octahedral site and a HS configuration (t2g3eg1), the energy of a singly occupied (so) eg orbital can be written as:

E e g so =E core +3U ' -3J H +6 Dq
where Ecore stands for the interaction of the 3d electrons with the nucleus and the [Ar] core electrons, and Dq is the crystal field parameter. The three t2g electrons are lying at the energy:

E t 2g so =E core +3U ' -3J H -4 Dq

Summing the crystal field energy over the four electrons and taking half the sum of exchange and correlation energies (as each contribution is shared by two electrons), the ion overall energy (ET), including nucleus–electron and electron–electron interactions, becomes:

E T =4E core +3/23U ' -3J H +1/23U ' -3J H -3×4 Dq +6 Dq =4E core +6U ' -6J H -6 Dq

The so-called t2g4(S = 1) low spin (LS) state actually is an intermediate spin (IS) state, but in no case the S = 0 state can be the ground state in an octahedral site of Oh symmetry owing to the Hund’s rule. Depending on whether the site is singly (so) or doubly occupied (do) and whether the spin is a majority (denoted α) or minority (denoted β) one, the electron energies respectively are:

E t 2g α so =E core +3U ' -2J H -4 Dq

E t 2g α do =E core +U+2U ' -2J H -4 Dq

E t 2g ϐ do =E core +U+2U ' -4 Dq

In the same way as above, we obtain for the total energy:

E T =4E core +U+5U ' -3J H -16 Dq

It should be pointed out that in a Td tetrahedral site where the e orbitals are stabilised by –2.67 Dq and t2 orbitals are lifted by 1.78 Dq, the e4 LS (S = 0) configuration reversely becomes stable vs the S = 1 state. The energy of an electron belonging to one of the two doubly occupied degenerate e-orbitals (as there is the same number of spin up and spin down electrons) is:

E e α do =E e ϐ do =E core +U+2U ' -J H -2.67 Dq
and the total energy is:

E T =4E core +2U+4U ' -2J H -10.68 Dq

Making equal the ET energies of the HS (S = 2) and IS (S = 1) states of a d4 ion in Oh symmetry leads to JH/Dq = 2. For high values of JH, the HS state is stable and for high values of Dq the IS state on the contrary becomes the ground state.

In Td symmetry, the condition for passing from S = 2 to S = 0 becomes JH/Dq = 1.1125.

For the d5and d6 cations, the IS state (S = 3/2 or S = 1, respectively) is unstable in Oh symmetry from the following conditions: in the d5 case, the LS→IS transition requires JH/Dq > 5/2, but the IS→HS transition requires only JH/Dq > 5/3. Therefore, one predicts a direct LS→HS transition which is actually expected for JH/Dq > 2. In the d6 case, the LS→IS transition requires JH/Dq > 10/3, the IS→HS transition only JH/Dq > 2. The direct LS→HS transition occurs at JH/Dq > 2.5.

3 Influence of a structural distortion on spin transitions of a d4 ion

Let us first consider the case in which one of the six oxygen atoms moves away from the central ion, the coordination progressively changing from octahedral of Oh symmetry to a (5 + 1) coordination of C4v symmetry until the limiting case of the square based pyramid is reached. In this type of coordination, a simple crystal field approach gives the following energy values in Dq units for the five d orbitals [6]:

E z 2 C 4v =+0.86,E x 2 -y 2 C 4v =+9.14,E xy C 4v =-0.86,E xz C 4v =E xz C 4v =-4.57

In octahedral (Oh) coordination:

E z 2 O h =E x 2 -y 2 O h =+6,E xy O h =E xz O h =E yz O h =-4

We assume that the magnitude of the distortion can be measured by a coefficient k (0 ≤ k ≤ 1) and that the energy of the orbital i is a function of k:

E i k=1-kE i 0+kE i
where Ei(0) is the energy of the orbital i in the undistorted site and Ei(1) the energy in the fully distorted site.

For instance, in the present case, with E i 0=E i O h and E i 1=E i C 4v , we obtain:

E z 2 =6-5.14k Dq ,E x 2 -y 2 =6+3.14k DqE xy =-4+3.14k Dq ,E xz =E yz =-4-0.57k Dq

Using these energy values Ei(k,Dq) to take into account the crystal field effect in the computation of the atomic orbital energies including the exchange correlation effects and then in the computation of the overall energy of the ion for each configuration, the following relations are obtained for the three possible spin equilibria:

S=2S=1J H / Dq =2-0.914kS=1S=0J H / Dq =1.237kS=2S=0J H / Dq =1.25-0.1075k

It is now possible to delimit in a JH/Dq-vs-k diagram the existence domains of each electronic configuration, S = 2, S = 1 and S = 0 (Fig. 3).

Fig. 3

Spin-state phase diagram for a d4 ion in C4v symmetry depending on exchange (JH), crystal field (Dq) and distortion (k) parameters.

The three straight lines intercept at a triple point t of abscissa kt where the three spin states are simultaneously stable (kt = 0.923 and JH/Dq = 1.142). As usual, if we look for the more stable state, only three half-line starting from t correspond to an actual equilibrium. As already pointed out, S = 0 cannot be stable in Oh symmetry. On the other hand, S = 1 cannot exist beyond the k value of 0.923. We can extend the calculations of the crystal field contribution to the five 3d atomic orbitals to the seven types of distortions illustrated in Fig. 2 , which encompasses most of the crystal sites usually observed in transition metal oxides. In cases (4) and (7), k represents a rotation from 0 to 60° and from 0 to 45°, respectively. Case (5) corresponds to a hypothetic case progressively passing from an octahedron (CN = 6) to a tetrahedron (CN = 4). The orbital energies are given in Table 1. The JH/Dq-vs-k diagrams are given for each spin state in Figs. 4–6 , for 3d4, 3d5, and 3d6 ions, respectively. These diagrams show that, whereas in Oh symmetry no value of JH/Dq allows to stabilise the three possible spin states of the d4, d5 and d6 ions – at least for narrow d levels –, most of the distortions (C4v, D4h, D h, D3 h, D4d) allow it, with one triple point t or even two triple points t1 and t2. For a d4 ion in D4h symmetry (i.e. an octahedron distorting into a square) kt1 = 0.465 and kt2 = 0.735, and for a d5 ion in D3 h symmetry (i.e. an octahedron distorting into a triangle based prism), kt1 = 0.175 and kt2 = 0.728. Some diagrams exhibit particular points of abscissa kc generally corresponding to the energy crossover of some atomic orbitals as k varies. Let us finally notice that predicted spin transitions are occurring for JH/Dq values close to 2 (±0.5). As according to literature, depending on the covalence of the bond, JH values are ranging from 0.6 to 0.9, such spin transitions are expected for Dq values ranging from 1.2 to 1.8 eV (i.e. from 2400 to 3600 cm–1 ± 25%).

Table 1

Energies (in Dq units) of the five d orbitals for the site distortions illustrated in Fig. 2.

dz2dx2y2dxydxz,yz
C4v6 – 5.14 k6 + 3.14 k–4 + 3.14 k–4 – 0.57 k
D4h6 – 10.28 k6 + 6.28 k–4 + 6.28 k–4 – 1.14 k
Dh6 + 4.28 k6 – 12.28 k–4 – 2.28 k–4 + 5.14 k
D3h6 – 5.04 k6 – 11.84 k–4 – 1.84 k–4 + 9.36 k
Td (hyp.)6 – 8.67 k6 – 8.67 k–4 + 5.78 k–4 + 5.78 k
Oh (CN8)–2.67 – 2.67 k–2.67 – 2.67 k1.78 + 1.78 k1.78 + 1.78 k
D4d–5.34–5.34 + 4.45 k3.56 – 4.45 k3.56
Fig. 4

Spin-state phase diagram for a d4 ion in the various symmetry sites shown in Fig. 2.

Fig. 5

Spin-state phase diagram for a d5 ion in the various symmetry sites shown in Fig. 2.

Fig. 6

Spin-state phase diagram for a d6ion in the various symmetry sites shown in Fig. 2.

4 Energies of atomic orbitals for dn (n = 4, 5, 6) ions at the triple point taking into account exchange and correlation

In Fig. 7 we give the orbital splitting for dn (n = 4, 5, 6) ions at some triple points (t) of the above diagrams (D4d: d4; C4v: d5 and d6) for each of the three spin states.

Fig. 7

Atomic orbital energies at the triple point determined in Figs. 4–6 for d4 in D4d, d5 and d6 in C4v symmetry. The origin for the energy scale is E0 = 3 U’ – 3 JH, E 0 = 4 U’ – 4 JH, E0 = 5 U’ – 4 JH for d4, d5 and d6, respectively. Dotted lines and dashed lines indicate JH and 3 JH energy differences, respectively. In some cases, (′): singly occupied orbitals, (′′): doubly occupied orbitals.

In each case, we can point out that various orbitals with majority α spin and minority β spin are degenerate. For instance, for a d6 ion in C4v symmetry, the dαdoxy or d 𝑥𝑦 ϐ do orbital from S = 0 is degenerate with dαz2 from S = 1. In the same way, d 𝑥𝑧,𝑦𝑧 ϐ do from S = 1 is degenerate, with dαx2-y2 from S = 2.

For a band diagram, in the vicinity of the triple point, bands with different symmetry and spin directions will thus overlap.

Hence, starting from S = 1, a spin equilibrium will be described as:

  • • (i) the excitation of an electron from d 𝑥𝑧,𝑦𝑧 ϐ do into (empty) dαx2-y2 with a spin flip; all the α spin states will be stabilised by –JH and the energy of the remaining d 𝑥𝑧,𝑦𝑧 ϐ electron will increase by + JH; in this way we have got the S = 2 (HS) state of Fig. 7;
  • • (ii) the excitation of an electron from dαz2 into dαsoxy with a spin flip; it will increase the stability of the β spin states (by –JH); we shall then obtain the S = 0 (LS) state ofFig. 7.

In the case of the d4 ion, we even observe a degeneracy of the HOMO – including both spin directions – of all the three (S = 0, 1 and 2) spin states simultaneously.

For S = 1, the HOMO itself contain both majority d 𝑥𝑦 α and minority d 𝑧 2 ϐ spin states.

5 Discussion

5.1 Structural distortions and stabilisation of unusual electronic configurations

5.1.1 d4 ion

Fig. 4 shows that any distortion of the coordination polyhedron with respect to the Oh (CN = 6) symmetry – with the exception of the prismatic D3h distortion – favours the HS configuration (S = 2). This is the well-known Jahn–Teller effect. On the other hand, a very unusual LS (S = 0) e4 configuration is favoured in cubic (CN = 8) coordination as well as in triangular prismatic (CN = 6) coordination as a result of the strong stabilisation of the dx2-y2 and dxy orbitals.

5.1.2 d5 ion

The usual configurations are HS (S = 5/2) and LS (S = 1/2). An IS (S = 3/2) configuration is possible, especially in C4 v and D4 h symmetry, and for realistic (close to 2) values of JH/Dq. This IS configuration could also be associated with a metallic state.

5.1.3 d6 ion

This is the most favourable number of d electrons for stabilising a LS (S = 0) configuration in an oxygen surrounding as the transition is occurring at JH/Dq = 2.5, corresponding to the lowest Dq value.

Obviously, the square planar D4h symmetry is the most appropriate for the IS (S = 1) state (see, for instance, [7] ) as well as the much more unusual Dh (CN = 2) one.

5.2 Exchange, correlation and covalence

U, U′ and JH first are atomic parameters measuring the exchange and correlation effects. As for Racah parameters and particularly B parameter, their values decrease with the hybridisation of metal 3d orbitals with oxygen 2p orbitals, i.e. the covalence of the M–O bond. As we have pointed out earlier [7], the energies of the electrons with majority α spin state are lowered by a value that can reach several times (up to 6 times in some cases) that of JH, which therefore favours the covalence of the M–O bond, whereas the contribution of electrons with minority β spins is more ionic – the covalent stabilisation is proportional to L ¯O - or even L ¯ ¯O 0 where βMO is the resonance integral of the bond and Hii (X) the Coulomb integral of the atom X [8]. This model provides a simple explanation for the high rate of metal-to-ligand hole transfer generally denoted L ¯(O) or even t 2g α3 e g α2 and observed by XAS for ions such as Fe4+, Co4+, Ni3+ or Cu3+.

5.3 Exchange, correlation and insulating vs metallic character

The metallic-vs-insulating character of a transition-metal oxide depends on the relative magnitude of the transfer energy βMO or of the bandwidth W and of the so-called Coulomb intra-atomic repulsion UHubbard, corresponding, in an ionic model, to the charge transfer Mn+ + Mn+ → M(n–1)+ + M(n+1)+. An upper limit of U Hubbard is then given by the difference between the (n + 1)th and nth ionisation energies. Actually UHubbard and W are not independent and UHubbard decreases as W increases. The Mott–Hubbard metal–insulator transition is thus governed by the magnitude of the U Hubbard/W ratio.

For a given dn configuration, it is then important to evaluate the orbital energy of the d n +1 configuration for the various possible spin states as a function of U or U′, and the type of orbital overlap that controls the value of W and thus the stabilisation due to the transfer. For instance, one sees that for a HS (S = 5/2) d5 ion in C4v symmetry only the HS d6 configuration involving the π-type dxz,yz orbitals is concerned (Fig. 7), whereas a IS (S = 3/2) ion can lead to two possible states S = 2 or S = 1 of the d6 configuration. Thus, the insulating-vs-metallic character of a dn transition metal oxide results from the magnitude of the energy gap separating the dn from the dn+1 configurations.

Confining us to the Oh symmetry and to the above d5/d6 case, the following electronic transitions can be considered, in which the type of bonding involved in the charge transfer is given between parentheses:

  • • (i) d5HS+ d5HS → d4HS + d6HS (σ, π)
  • • (ii) d5IS+ d5IS → d4IS + d6HS (σ)
  • • (iii) d5IS+ d5IS → d4HS + d6IS (π)
  • • (iv) d5LS+ d5LS → d4IS + d6LS (π)

for each of them we can calculate the energy gaps using the total energies including exchange, correlation and crystal field contributions for all the valence electrons (Table 2).

Table 2

Exchange and correlation energies for valence electrons of various d4, d5 and d6 configurations.

HSISLS
d4U' – 6 JH – 6 DqU + 5 U' – 3 JH – 16 Dq
d510 U' – 10 JHU + 9 U' – 6 JH – 10 DqU + 8 U' – 4 JH – 20 Dq
d6U + 14 U' – 10 JH – 4 DqU + 13 U' – 7 JH – 14 DqU + 12 U' – 6 JH – 24 Dq

With the exception of (i), all the transfers correspond to a minimal energy of U′ – JH. For the first case, the gap (U′ + JH, assuming JH/Dq = 2) is less favourable to such a transfer.

However, only (ii) involves σ-type bands that are broader than π-type ones. The ratio (U′ – JH)/Wσ is then the most favourable for the formation of a metallic state from an IS state.

In the case of a lower symmetry favouring the HS↔IS spin equilibrium (Fig. 5), it should be noticed that case (i) can be split into two steps:

  • • (i′) d5HS + d5HS → d5IS + d5HS
  • • (i′′) d5IS + d5HS → d4HS + d6HS (σ)

the gap corresponding to (i′′) is equal to U′.

This HS/IS configuration mixing could also lead to a metallic state, as it implies a σ-type transfer.

This simple ionic-type approach of exchange and correlation explains the insulating behaviour of LaFeO3and of LaCoO3 at low temperature as well as the transition to a metallic behaviour at high temperature for the latter.

Indeed, for Fe3+ in LaFeO3, the JH/Dq ratio, much larger than 2.5, favours the exchange energy gain with respect to the crystal field and leads to a HS t 2g α3 e g α2 t 2g ϐ1 state. The upper Hubbard band corresponding to the HS t 2g 6 e g 0 d6 configuration lies at a much higher energy than the lower Hubbard band and is well separated from it, which accounts for the insulating character.

On the other hand, in LaCoO3, Co3+ is in a LS

M 5+ + Fe 3+ M 6+ + Fe 2+
state at low temperature and thus insulating, because the crystal field energy overcomes the exchange contribution. As the temperature increases, the Co–O distances become longer, leading to a decrease of the covalence and of Dq, and the exchange and correlation parameters slightly increase, which leads to a change in the electronic configuration. LaCoO3 was first described as undergoing a LS(S = 0)↔HS(S = 2) transition [9]. Recent works now describe the transition as LS(S = 0)↔IS(S = 1), which better agrees with our description of the insulator to metal transition given above [10, 11]. Double perovskites synthesised several decades ago such as Sr2FeMoO6 (M = Mo, W, Re...) and characterised by very different physical properties (either ferromagnetic metal or antiferromagnetic insulator) depending on the nature of M, are now revisited, since giant magnetoresistance has been observed for these materials [12]. If we consider the possibility of a charge transfer between M5+ and Fe3+ ions such as:

t 2g ϐ do we have to compare the energy of iron d orbitals (EFe:3d = –16.5 eV) modified by crystal field, exchange and correlation effects (tαso2g of Fe2+ is shifted by 6 JH above π*t 2g ϐ ) with that of M5+ ions orbitals (Mo5+: 4d1, W5+: 5d1, Re5+: 5d2), the energies of which depend on the main quantum number n and atomic number Z (EMo:4d = –11.5 eV, EW:5d = –11.0 eV, ERe:5d = –12.3 eV) (Fig. 8).

Fig. 8

M5+/Fe3+charge transfer in Sr2FeMO6 (M = Mo, W, Re) perovskites.

We now understand why a partial transfer is possible in the molybdenum and rhenium cases: a ϐ MO 2 /H ii M-H ii O, band of the Fe–Mo(Re)–O type involving mixed valences FeIII/FeII and MoV/MoVI gives rise to a metallic character and strong ferromagnetic interactions between localised iron moments via the itinerant π*β electrons accounting for a high Tc value.

The highest energy of tungsten 5d orbitals leads to a complete charge transfer, giving single valence Fe2+(3d6) and W6+ (5d0) ions, resulting in an insulating character and weak antiferromagnetic interactions (TN ≈ 20 K) between Fe2+ ions.

6 Conclusions

In order to clarify the competition of exchange and correlation energies vs crystal field stabilisation energy, a simple approach is proposed. The exchange and correlation contributions have been described on the base of Brandow and Kanamori’s U, U′ and JH parameters. The dependence of the crystal field on site distortion has simply been modelled using a linear interpolation between undistorted and fully distorted site. This approach led to establish phase diagrams for d4, d5 and d6 cations, allowing us to predict the spin-state stability range depending on exchange (JH), crystal field (Dq) and distortion (k) parameters. It can be used to complement the Extended Hückel Tight Binding calculations of the electronic structure of materials with a noticeable ionic character such as 3d transition-metal oxides, in order to interpret their electronic behaviour and, particularly, discuss their insulating-vs-metallic character, as we recently did [13] .


References

[1] C. Cohen-Tannoudji; B. Diu; F. Laloë Mécanique quantique, tome II, Hermann, Paris, 1998, p. 1410; J.D. Morgan III, W. Kutzelnigg, J. Phys. Chem., Volume 97 (1993), p. 2425

[2] V.I. Anisimov; J. Zaanen; O.K. Andersen Phys. Rev. B, 44 (1991), p. 943

[3] E. Clementi; D.L. Raimondi J. Chem. Phys., 38 (1963), p. 2686

[4] J. Kanamori J. Phys. Chem. Solids, 10 (1963) no. 87

[5] J. Brandow Adv. Phys., 26 (1977), p. 651

[6] R. Krishnamurthy; W.B. Shaap J. Chem. Educ., 46 (1969), p. 799

[7] M. Pouchard; A. Villesuzanne; J.P. Doumerc J. Solid State Chem., 162 (2001), p. 282

[8] W.E. Pickett; D.J. Singh Phys. Rev. B, 53 (1996), p. 1146

[9] J.B. Goodenough Progr. Solid State Chem., 5 (1971), p. 145

[10] M.A. Señaris-Rodriguez; J.B. Goodenough J. Solid State Chem., 116 (1995), p. 224

[11] M.A. Korotin; S.Y. Ezhov; I.V. Solovyev; V.I. Anisimov; D.I. Khomskii; G.A. Sawatzky Phys. Rev. B, 54 (1996), p. 5309

[12] S. Ray; A. Kumar; S. Majumdar; E.V. Sampathkumaran; D.D. Sarma J. Phys.: Condens. Mat., 13 (2001), p. 607

[13] M.H. Whangbo; H.J. Koo; A. Villesuzanne; M. Pouchard Inorg. Chem., 41 (2002), p. 1920


Comments - Policy