Open access peer-reviewed chapter

XAO Analysis - AO's and Their Populations in Crystal Fields

By Kiyoaki Tanaka and Yasuyuki Takenaka

Submitted: November 21st 2011Reviewed: June 7th 2012Published: September 19th 2012

DOI: 10.5772/50517

Downloaded: 1273

1. Introduction

There are two ways to the analysis of electron density distribution (EDD) measured by X-ray diffraction. Multipole refinement (Stewart, 1969, 1973; Hansen & Coppens, 1978) expresses EDD as a linear combination of spherical harmonics. Using the determined functions, EDD is analyzed with Bader topological analysis (Bader, 1994). The concept of critical points of Bader theory defines the characters of chemical bonds. However, it does not give atomic orbital (AO) and its electron population, except for orbitals in high symmetry crystal fields where AO's are defined by theory. The other is based on quantum-mechanical orbital-models and gives the physical quantities just mentioned. The objective of this chapter is to introduce one of the second methods, the X-ray atomic orbital analysis (XAO) (Tanaka et al., 2008; hereinafter referred to as I), and its application to metal and rare-earth complexes, where advantages of XAO are revealed quite clear.

XAO is closely related to the electron population analysis (Stewart, 1969; Coppens et al., 1971). However, it was abandoned since the orthonormal relationships between orbital functions caused severe parameter interactions in conventional least-squares refinements. Later, the method incorporating the ortho-normal relationships between AO or MO (molecular orbital) was formulated by employing Lagrange’s unknown multiplier method (Tanaka, 1988; hereinafter referred to as II).

The AO-based EDD refinement has been started from 3d-transition metal complexes with so high crystal symmetry that the AO's are known by the crystal field theory. Spin states of the metals in perovskites, KCoF3 and KMnF3 (Kijima et al., 1981, 1983) and KFeF3 (Miyata et al., 1983) were determined to be high spin. On the other hand, mixed orbitals of 3dx2y2and 3dz2in Jahn-Teller distorted KCuF3 were determined (Tanaka et al., 1979). The hybrid orbital of 3dz2and 4s of Cu+ which made O-Cu-O straight bond in CuAlO2 was also determined (Ishiguro et al., 1983). The investigations on KNiF3 and KCuF3 are, to our knowledge, the first ones that determined spin states and AO’s by X-ray diffraction, respectively. The Cu2+ ion in [Cu(daco)2](NO3)2 (daco: diazacyclooctane) is in the crystal field Ci. Five d-orbitals were determined with the least-squares method stated in II.

These investigations revealed that the anharmonic vibration (AHV) of atoms could not be ignored, since significant peaks still remained after removing the 3d-peaks on difference density maps. They were attributed to the AHV of the metal atoms (Tanaka & Marumo, 1982; Ishiguro et al. 1983). For these investigations, a method proposed by Dawson et al., (1967), which had been applied only to high symmetry crystals, was made applicable to atoms in a general crystal field (Tanaka & Marumo, 1983).

Since the ratio of the number of bonding electrons to those in the unit cell becomes smaller as the atomic number increases, very accurate structure factors are necessary for X-ray EDD investigations of rare-earth compounds. Thus, the EDD analysis based on chemical-bond theories had not been done when we started the study on CeB6. Actually, the ratio in CeB6 is 1/88, which demands us to measure structure factors with the accuracy less than 1 %. Rare-earth crystals are usually very hard and good resulting in enhanced extinction and multiple diffraction (MD). Therefore, MD was investigated using the method by Tanaka & Saito (1975) in which an effective way to detect MD and correct for it were proposed. It introduced the time-lag between the relevant reflections, which usually do not occur at exactly the same time, to the method by Moon and Shull (1964). The effect of MD was demonstrated through the study of PtP2 (Tanaka, et al., 1994) by measuring intensities avoiding MD and compared to those measured at the bisecting positions.

The frontier investigations aiming to measure 4f-EDD were done for CeB6 (Sato, 1985) at 100 and 298 K and for nonaaqualanthanoids (Ln:La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Yb and Lu) (Chatterjee at al. 1988). Significant peaks were found around these rare-earth atoms on the residual density maps. However, the aspherical features of 4f-EDD were not analyzed quantitatively by the X-ray scattering factors calculated with AO's in a crystal field. The first 4f-EDD analysis of CeB6 with the crystal field model but without spin-orbit interaction exhibited T1u state was occupied (Tanaka et al. 1997). The ground state of the 4f-state was found to be 4f(j=5/2)Γ8 by the inelastic neutron scattering (Zirngiebl et al., 1984). When the spin-orbit interaction is taken into account, the total sum J of the orbital and spin angular momentum L and S, J=L+S, becomes a good quantum number. In the following discussion J values of p, d, f orbitals are attached as a subscript on the right-hand side of each orbital. The experiment does not contradict to our results. It was, however, investigated again at 100, 165, 230 and 298 K (Tanaka & Ōnuki, 2002) by introducing the spin-orbit interaction. Our experiments ascertained that 4f electron occupied 4f5/2Γ8. orbitals. It revealed further that the electron population of them, n(4f5/2Γ8 ), decreased on lowering the temperature. 4f5/2Γ8 electrons are transferred to the B6 moiety below room temperature. In addition to it just significant amount of 4f5/2Γ7 electrons were found at 298 K, which was expected to be due to the thermal excitation. Thus, the electron population above room temperature became very interesting. Since the energy gap between the two 4f5/2 orbitals is about 500 K, the X-ray intensities were measured at 430 K expecting for 4f5/2Γ7 to be more populated (Makita et al., 2007). The results were surprising. The 4f5/2Γ7 is more populated than 4f5/2Γ8, and 5d5/2Γ8 is fully occupied at 430 K, which is expected to be at least a few ten thousands degrees higher than the 4f orbitals. The electron configurations of CeB6 at 340 and 535 K were further investigated to confirm the results (Makita et al., 2008). The EDD of SmB6 (4f5) was studied (Funahashi, et al., 2010) below room temperature. 4f5/2Γ8 and 4f5/2Γ7 orbitals are fully or partially occupied. It is noted that 5d5/2Γ8 occupation was also found. The physical properties such as electric resistivity are correlated to the electron configuration.

2. Basic formalism of XAO

2.1. AO's in crystal fields

The i-th AO,ψα,i(r)(i=1,2,...,M) of the α-th atom are expressed with N basisfunctions, k(r)

ψα,i(r)==kNcikψk(r), i=1,2,...,M,E1

where cik is a constant to be determined in the XAO analysis with the least-squares method incorporating orthonormal relationship between AO's (Tanaka, 1988). Matrices and vectors are written in upper-case letters and bold lower-case letters, and superscript ‘ means a transposed matrix or a row vector, respectively. The basis functions, ψk(r)are listed in Table 1 of I. They are expressed using polar coordinate (r, θ, ϕ) as a linear combination of a product of a radial function Rnl(r) and a spherical harmonic Ylml(θ,ϕ)as,

ψk(r)=ml=lldk,mlRnl(r)Ylml(θ,ϕ),E2

where n, l and ml are principal, azimuthal and magnetic quantum numbers. The non-relativistic radial functions, Rnl(r) (Mann, 1968) and relativistic functions calculated by HEX (Liberman et al., 1971) were mainly used for the XAO analysis. dk,ml's are listed in Table 1 of I. Since approximate cik's were necessary at the start of the non-linear least-squares calculation, they were calculated by taking the crystal field, which was calculated by placing a proper point charge on each atom, as a perturbation to the system (Kamimura, et al., 1969). For details, see eqs. (52)-(56) of I. The adjustable variables cik's in the least-squares calculation are listed for all point group symmetries in Tables 4 and 5 in I.

2.2. Least-squares method incorporating ortho-normal condition

In a non-linear least-squares refinement, mathematical and physical relationships between parameters should be taken into account to avoid parameter interaction. The conventional one used in X-ray crystallography was improved to obtain AO/MO by taking into account the orthonormal relationship between wave functions (Tanaka, 1988),

mNnNcim*cjnsmn=δijE3

where smn=ψmψndrand N is the number of basis functions. δij is Kronecker’s δ. Assuming cim's are real, it is rewritten with (M,N) matrix C={cim} and (N,N) matrix S={smn},

CSC'=I,E4

where I is the unit matrix of order M. Using Lagrange’s unknown multiplier λij,, the value Q to be minimized in the least-squares method becomes,

Q=v'Mf1viMjMλij(mNnNcimcjnsmnδij)E5

For definition of v see eqs. (10a) to (10d) of II. Putting cim as the sum of the true value cim0and its small shift Δcim, Q is expressed by ignoring the terms higher than second order of Δcim,

Q=v'Mf1viMjMλij(mNnN(Δcimcjn0+cim0cjn)smnδij),E6

Putting λij’=(1+δijij/2, Q is differentiated to give,

δQ=2δx'(A'Mf1AxA'Mf1f)2iMjMλij'mNnNδ(Δcimcjn0smn)=0,E7

where xrepresents least-squares estimates of parameters. The first term is obtained following the usual procedure. For x, A and Mf, see Hamilton (1964). After alligning Δcim's linearly and adding at the end of x for simplicity, λij is expressed in terms of cim by using the orthonormal relationship in eq. (3). Then the final secular equation is obtained,

(IR)A'Mf1Ax=(IR)A'Mf1f,E8

where I is unit matrix of order P, total number of variables. The (P,P) matrix R is expressed in terms of cik’s. It makes the inverse matrix of (IR)A'Mf1Ain (8) calculable. The explicit form of R={rij} under the limiting conditions (10) for i, j, k and l is,

rij={SC'C}kl,E9
PMN+(m1)N+1i,jPMN+mN, k=i{PMN+(m1)N}, l=j{PMN+(m1)N},E10

where m runs from 1 to M. The other elements rij not defined here are zero.

2.3. Electron density and structure factors

The electron density ρα(r) of the αth atom centered at rαatomis divided into that of the core and valence orbitals,ρa,valence(r). It is further expressed from eq. (1),

ρα,valence(r)=inα,iρα,i(κi(rrαatom))=inα,iΨα,i(κi(rrαatom))Ψα,i(κi(rrαatom)),E11

where nα,iis the number of electrons of Ψα,iexpresses the expansion (κi>1)or contraction (κi>1)of i-th orbital (Coppens et al., 1979). Atom α in the asymmetric unit is translated to rασatomby the σ-th crystal symmetry operation. Then the structure factor is,

F(k)=αωασfασ(k)exp[ikrασatom]Tασ(k),E12

where fασ(k)is an atomic scattering factor and k is a scattering vector. The αth atom at rασatomis displaced by atomic vibration torασatom+u. It perturbs the periodicity of the crystal and reduces diffracted intensity. Then the atomic displacement parameter (ADP),

Tασatom(k)=exp(iku),E13

is calculated as an average in space and time of exp(ik∙u) by assuming each atom vibrates independently from each other. fασis the sum of the scattering factors of core and valence electrons, fασcore(k)andfασvalence(k). fασvalence(k)is expressed as,

fασvalence(k)=inifα,σ,ivalence(k),E14

where ni is the number of electrons of the ith valence orbitals. fα,σ,ivalence(k)is a Fourier transform ofρα,σ,ivalence(r). Final expression of fα,σ,ivalencebecomes,

fα,σ,ivalence=kk'cikcikK=|ll'|l+l'<jK>M=KKiK2(2K+1)ΘKM(βσ)exp(iMγσ)×ml=llml'=l'l'dk,mldk',ml'cK(lml,l'ml'),E15

where K+l+l' is even. M=mlml'and|ll'|Kl+l'.jKis expressed as,

jK=Rnl(r)Rn'l'(r)jK(kr)r2drE16

jK is Bessel function of order K. For σσ), see Appendix B of I.

2.4. Anharmonic vibration (AHV)

Gaussian probability density functions (p.d.f.) of atoms are Fourier transformed to give harmonic ADP. Gram-Charier (G-C) formalism is widely used to introduce AHV. However, it expresses aspherical EDDs both from orbital functions and AHV, as shown clearly by Mallinson et al. (1988). Thus, it may not be adequate for accurate EDD researches. The method based on Boltzmann statistics (Dawson et al. 1967, Willis, 1969) was employed,

Tασatom(k)=exp(V(u)/kBT)exp(iku)du/exp(V(u)/kBT)du,E17
(17)

where kB is Boltzmann constant, and T is absolute temperature. Assuming each atom vibrates independently and defining u(u1,u2,u3) on the Cartesian coordinates parallel to the principal axes of the thermal ellipsoid of harmonic ADP, potential V(u) is expressed as,

V(u)=V0+12i=13biui2+i,j,k'cijkuiujuk+i,j,k,l'qijkluiujukul+.....E18

where i, j, k, l run from 1 to 3 and ‘ means they are permutable. There are 10 cubic and 15 quartic terms. Potential expansion terminates at the quadratic terms in the harmonic ADP.

By assuming the sum of the third and the fourth terms in (18) is much smaller than kBT, bi, cijk and qijkl are introduced into structure factor formalism using (12), (17) and (18).

3. Experimental

3.1. Multiple diffraction (MD)

MD occurs when two or more reflections are on the Ewald sphere. When the incident beam with wave length λenters parallel to the unit vector j, and a primary reflection, the intensity of which is measured, is at the reflecting position, the conditions for a secondary reflection, h1h2h3, to occur become as follows (Tanaka & Saito, 1975),

|e3(j/λ+h3b3|1/λE19
|e2(j/λ+h2b2+h3b3|[1/λ2{e3(j/λ+h3b3}2]1/2r3E20
e1(j/λ+h1b1+h2b2+h3b3)=±[r32{e2(j/λ+h2b2+h3b3)}2]E21

b1, b2, b3 are reciprocal-lattice vectors when the primary reflection exactly fulfills the reflecting condition and e1, e2, e3 are the unit vectors defined as,

e3=b1×b2|b1×b2|,e2=e3×b1|e3×b1|,e1=b1|b1|.E22

Secondary reflections only around the surface of the Ewald sphere are searched using eqs. (19) to (21). Integers h3 are selected from (19) and integers h2 are calculated from (20) for each h3. Then the value of h1 is evaluated from eq. (21) for each (h2, h3). The point (h1,h2, h3) is judged to be a reciprocal lattice point if h1 is an integer. The Ewald sphere has width due to the divergence of the incident beam, finite size and mosaicity of the crystal, and wavelength spread. It permits some margin from an integer to h1.

The change of intensity by MD was first formulated by Moon & Shull (1964). Primary and secondary beams usually come on the Ewald sphere with an angle-lag Δε. Using Δε The formalism was modified for a spherical crystal with a radius r (Tanaka & Saito, 1975),

ΔI1I1=πr28Q01i[g01;0i(Q0iQ01)g01;1i(Q1iQ01)+g01;i1(Q0iQ01)(Qi1Q01)].E23

Subscripts ‘0’, ‘1’ and ‘i’ stand for the incident, primary and secondary beams. Qmn is the integrated reflectivity per unit volume of a crystal for the m-th incident beam and n-th diffracted beam (m-n process). It is proportional to the squares of the structure factor. In i-1 process, for example, a secondary beam i acts as an incident beam and is diffracted along the direction of the primary beam 1, which enhances the intensity of the primary beam. gij;mn is a function of Δε, Lorentz and polarization factors of the simultaneous i-j and m-n processes. The polarization factor for a general multi-diffracted process was recently formulated (Tanaka et al., 2010). For gij;mn see Tanaka & Saito (1975). Since ΔI1/I1 depends on the ratio of Qmn's and Δε of sharp primary and secondary reflections, Δ I1 fluctuates sharply. Eq. (23) was applied to diformohydrazide (Tanaka, 1978) and two peaks on the N-N bond in the difference density map became one peak at the middle of the bond.

The EDD analyses of compounds with metals and rare-earth elements need accurate measurements of structure factors, however MD affects the measured intensity seriously since Qmn is large. However correction for it does not seem fruitful, since each primary reflection has many secondary reflections and MD depends sharply on crystal orientation and half-width of each Bragg peak. Thus it is better to avoid MD by rotating a crystal around the reciprocal vector of the primary reflection (ψrotation) so that secondary reflections move apart from the Ewald sphere and measure at where ΔI1/I1 is minimum. For actinoid compounds, MD cannot be avoided and the correction for it becomes necessary.

4. XAO analysis for crystals in the Oh crystal field

4.1. 3d, 4f and 5d orbitals in the Oh crystal field

In the present chapter EDD investigations on the first-transition metal complexes and rare earth compounds in Oh crystal field are mainly stated. Thus it is worth while to represent the energy level splitting of d- and f-states in Figs.1(a) and (b) (Funahashi, 2010). The spin-orbit splitting is neglected in a strong field model. For (3d)n systems strong field model is often employed. In (4f)n and (5d)n systems the two splittings are of the same order of magnitude. Weak field models are employed for CeB6 and SmB6. Their quantization axes point to the face-centre of the cubic unit cell, where no B atom exists. Note that the order of the levels in Fig. 1(a) is inverted, that is, t2g orbitals lies lower than eg, when ligands are on the quantization axes like KCoF3 (section 4). In weak field 5d- and 4f-states split into j=3/2 and j=5/2 states, and into j=5/2 and j=7/2 states by the spin-orbit interaction. The Oh crystal field splits them further. Lobes of 4f5/2Γ8 extend along <100> directions while those of 4f5/2Γ7 orbitals along <110> and <111> as illustrated in Fig. 7 of Makita et al. (2007).

Figure 1.

Energy-level splitting of (a) d- and (b) f-states in the Oh crystal field with the degeneracy of each state in parentheses. Each orbital of a strong field can have two electrons.

4.2. 3d-EDDand spin states in KCoF3, KMnF3 and KFeF3

The peak due to 3d electrons was first reported for [Co(NH3)6][Co(CN)6] (Iwata & Saito, 1975) around the Co atom in 3¯crystal field. Later electron distribution at 80 K was reported (Iwata, 1977). 3d-EDD in γ−Co2SiO4 (Marumo et al., 1977) and CoAl2O4 (Toriumi, et al., 1978) were also observed. The electron density around Co in KCoF3 was first analysed quantitatively based on the 3d-orbitals in Oh crystal field (Kijima et al., 1981). A crystal was shaped into a sphere of diameter 0.13 mm, and intensities were measured up to 2θ=130°with a four-circle diffractometer using MoKα radiation without avoiding MD. cik's in (1) are fixed constants because of the high symmetry as listed in Tables 1(b) and 5(b) of I. κi in eq. (11) was not introduced to the investigations in sections 4 5.1 and 5.2. Scattering factors of dx2y2,dz2,dyz,dzxanddxyorbitals are calculated with eq. (15). When electrons occupy the five d-orbitals equally, that is (t2g)4.2(eg)2.8 for Co2+, atoms have spherical EDD. Hereinafter the refinement with the spherical scattering factors and without assuming AHV is called as spherical-atom refinement, in which type I or II anisotropic extinction (Becker & Coppens, 1974a, b, 1975) was assumed. The difference density map assuming the spherical-atom, high and low spins, (t2g)5(eg)2 and (t2g)6(eg)1, are shown in Figs. 2(a) to (c). Peaks remained after spherical-atom refinement were reduced and enhanced after high- and low-spin refinements, respectively. Large and almost no peaks in Fig. 2(c) and Fig. 2(b) demonstrate that Co2+ is in the high-spin state. Spin-state was first determined by X-ray diffraction in this study. In the similar way, the spin-states of 3d-transition metals in KMnF3 (Kijima et al., 1983) and KFeF3 (Miyata et al., 1983) were determined to be high-spin.

Figure 2.

Difference density after refinements of (a) spherical-atom, (b) high-spin, (c) low- spin models. Contours at 0.2 eÅ-3. Negative and zero in broken and dashed-dotted lines.

5. Experimental AO determination

When symmetries of crystal fields are low, cik's in (1) become unknown variables, and it is necessary to determine them by the least-squares method keeping the orthonormal relationship between AO's. Independent variables, cik's, and the relation between them are listed for 32 point group symmetries in Tables 4 and 5 of I.

5.1. Jahn-Teller distortion in KCuF3 and mixed dx2y2and dz2orbital

In KCuF3 each F- ion between Cu2+ ions (3d9) shifts from the centre by Jahn-Teller effect. It makes short Cu-Fs, medium Cu-Fm, parallel to c-axis, and long Cu-Fl bonds, resulting in mmm point group symmetry for Cu2+. Difference density after spherical-atom refinement exhibits non-equivalent four holes in Fig. 3 (b) (Tanaka et al. 1979). Table 5(b) of I indicates a mixing of 3dz2and3dx2y2in the crystal field mmm,

ψG=cos(φ/2)dz2+sin(φ/2)dx2y2,ψE=sin(φ/2)dz2cos(φ/2)dx2y2.E24

The peaks on Cu-Fl in Fig. 3(a) correspond to the lone pair ofψG. Putting two and one electrons to ψGandψE, cos(/2) became 0.964(18). 3d-peaks in Fig. 3 (a) and (b) were deleted in (c) and (d). It confirms orbitals in eq. (24) are quite reasonable. The orbital functions except the phase factor were determined for the first time from X-ray diffraction experiment. However, significant peaks still remained. The positive and negative peaks along <100> and <110> in Fig.3(c) indicate preferential and inhibitive vibration of Cu2+, respectively. Therefore, AHV of Cu2+ was analysed (Tanaka & Marumo, 1982). The obtained potential in (18) corresponds to the peaks and removes them as seen in the residual density in Fig. 4. Our study on AHV started from this investigation. Since the AHV peaks appear in Fig. 3(c) and (d) after the d-electron peaks were removed, the aspherical peaks due to electron configuration and vibration of Cu2+ seem to be separated well.

Figure 3.

Difference density after spherical-atom refinement on (a) (001) and (b) (220) planes passing Cu2+ at (1/2,0,0). That after d-orbital analysis on (c) (001) and (d) (220). Contours are the same as those in Fig. 2.

Figure 4.

Difference density after AHV analysis. Contours are the same as those in Fig.2.

5.2. 3d orbitals in 1¯crystal field of Cu(daco)2(NO3)2

The Cu2+ ion (3d9) of Cu(daco)2(NO3)2(daco: 1,5-diazacyclooctane) is on the centre of symmetry and forms a coordination plane with the four N atoms in daco (Hoshino, et al., 1989). Difference densities on and perpendicular to the Cu-N4 plane after the spherical-atom refinement are shown in Fig. 5(a) and (b). Negative four peaks in Fig. 5 (a) near Cu2+ correspond to a 3d-hole. However, the peaks do not locate on the Cu-N bonds. Since the crystal field of the Cu atom is Ci, d-orbitals of Cu2+ are expressed as a linear combination of dx2y2,dz2,dyz,dzxanddxyorbitals and all the 25 coefficients cik were refined with the method in II. Since approximate cik's are necessary for the non-linear least-squares refinement, they are calculated putting a point charge on each atom (Kamimura et al., 1969). jK3dof Cu2+ were taken from International Tables for X-ray crystallography (1974, vol. IV). After each cycle of refinement, new set of coefficients were orthonormalized by the Löwdin’s method (Löwdin, 1950). Refinement was converged and the coefficients cik's, population ni and κi are listed in Table 1. The number of significant cik’s is seven among 25 coefficients. Orbital 1 in Table 1 is the hole orbital. It is not a pure dx2y2but a mixed orbital of dx2y2and dzxorbitals, c14=0.57(59) is not significant. The hole and large negative and positive peaks in Fig. 5 were removed in Fig. 6. Further AHV refinement removed the peaks in Fig.6, but that prior to the d-orbital analysis did not improve Fig. 5.

Figure 5.

Difference density after spherical-atom refinement on (a) Cu-N4 plane and (b) the plane perpendicular to it. Contours are the same as those in Fig. 2.

Figure 6.

Difference density after 3d-orbital analysis. Contours are the same as in Fig. 2

When G-C formalism was applied to an iron complex, it removed the peaks equally well as the multipole refinement did (Mallinson, et al., 1988). Our AHV analysis is based on the classical thermodynamics and does not have such a problem, though the p.d.f. function exp[-V(u)/kBT] is not integrable (Scheringer, 1977) at a place far from the nucleus. However, we can apply it safely as long as it is applied within the region where the sum of the third and fourth terms in (18) is much smaller than kBT. The condition is usually fulfilled in the investigations of EDD.

iniκidx2y2dz2dyzdzxdxy11.00.97(11)0.75(34)0.21(59)0.23(216)0.57(59)0.10(79)22.01.37(7)0.08(28)0.12(17)0.19(19)0.10(26)0.97(6)32.00.81(5)0.63(22)0.15(52)0.44(44)0.61(34)0.06(22)42.00.71(4)0.04(32)0.93(13)0.10(29)0.34(41)0.06(10)52.00.90(5)0.15(76)0.22(45)0.84(43)0.42(85)0.22(15)

Table 1.

3d-orbital parameters cik’s of Cu2+. Significant ones in thick lines

6. 4f-EDD in rare-earth hexa-borides

4f-EDD analysis has become more and more important since many rare-earth compounds with interesting physical properties, such as high-temperature super conductors, have been found. The EDD of rare-earth hexa-borides were investigated since they are famous for their properties related to the Kondo effect.

6.1. XAO analysis of CeB6

Ce and B atoms are at the body-centre and on the edges of the cubic unit cell with the B6 regular octahedra at the corners as shown in Fig. 7. Ce is in an Oh crystal field and B atoms are on 4-fold axes. CeB6 is a typical dense Kondo material with Kondo-temperature Tk=2 K. The magnetic moment decreases with that of temperature and vanishes below TK. However, Kondo effect starts above room temperature (Ōnuki et al. 1984). Thus, change of EDD below room temperature is interesting although our diffractometer did not permit to measure below 2K. Intensities were measured at 100, 165, 230 and 298 K aiming to separate the two kinds of asphericity of EDD stated in 5.1 and 5.2 (Tanaka & Ōnuki, 2002). MD was avoided by ψ−scan with the program IUANGL (Tanaka et al., 1994). Crystal was shaped into a sphere with a radius 36 μm to reduce absorption, extinction and MD.

Figure 7.

Crystal structure of CeB6

Spin-orbit interaction was introduced into the XAO analysis. Four-fold degenerate 4f5/2Γ8 orbitals were taken from Table 4 (b) of I as,

Ψ14f=56ψ1+16ψ5,Ψ24f=16ψ2+56ψ6,Ψ34f=ψ3,Ψ44f=ψ4,E25

where ψi(i=1 to 6) are basis functions in Table 1(c) of I. In the similar way 5p1/2, 5p3/2 and 4f5/2Γ7 and 4f7/2Γ678 orbitals were introduced into the XAO analysis.

Figure 8.

Difference density on (100) plane around Ce at (a) 100 K, (b) 165 K, (c) 230 K, (d) 298 K. Contours are the same as those in Fig. 2.

Temp.Cen(Γ8)*κ(Γ8)κ(5p3/2)q1111q1122Bn(2px)n(2pz)κ(2px)κ(2pz)CeB100K13(1)179(8)95(2)47(19)131(48)27(6)103(8)93(31)92(12)3.0416(3)165K15(2)158(12)98(2)14(6)39(15)52(4)52(8)85(16)108(11)3.0440(4)230K19(2)157(11)100(2)0(4)1(11)69(4)16(8)51(19)130(18)3.0459(7)298K21(2)160(10)100(2)6(2)13(5)67(4)19(8)54(17)130(13)3.0439(4)

Table 2.

Orbital and AHV parameters. n(X), and κ(X) are a number of electrons and κ of the orbital X. n(x10-2), κ (x10-2), qiijj (x10-19 JÅ-4). 5p3/2,1/2 and 2s are fully occupied. n(2px)=n(2py) and κ(2px)=κ(2py). Unit of Ce-B is Å. *:n(Γ7) is 4.2(3.6)x10-2 only at 298 K.

6.1.1. XAO analysis of CeB6 below room temperature

After spherical-atom refinement for ions, Ce3+ and B0.5-, the populations of the orbitals were refined keeping the sum of them equal to that of nuclear charges. When populations of them exceeded one/two or became negative, they were fixed to one/two or zero. In our program QNTAO (Tanaka, 2000) each sub-shell is treated as an pseudo-atom sharing atomic parameters with the other sub-shells. It enables us to analyze the EDD of non-stoichiometric atoms, since valence electrons can be treated independently from the core electrons. Basis functions in Table 1 of I are automatically assigned to each AO labelled with 2l+1 and the number of basis functions. Since high symmetry fixes cik's in (1), the other orbital parameters and atomic parameters, including AHV were refined. Electron density around Ce on (001) after spherical-atom refinement is shown in Fig. 8. The peak heights around Ce increase from 0.6, 1.2 to 2.0 eÅ-3on lowering the temperature to 165 K and reduce to 0.6 eÅ-3 at 100 K. After XAO analysis they were removed, exhibiting they were due to the 4f-electrons. Parameters of AO's and AHV are listed in Table 2. No electrons were found in 4f7/2 orbitals but B-2s and Ce-5p are always fully occupied. n(4f5/2Γ7) has value 0.042(36) only at 298 K but 4f5/2Γ8 is occupied at the four temperatures. It agrees with the previous experiments (Zirngiebl et al., 1984;Sato et al., 1984). Shorter Ce-B at 298 K than that at 230 K seems to be correlated to 4f5/2Γ7. It extends along 111or to the centre of B6 octahedron making them extend to reduce electrostatic repulsion. Expanded 2px still has electrons. Total number of 4f-elctrons mainly composed of 4f5/2Γ8 are 0.98(11), 0.77(8), 0.61(7) and 0.52(6) at 298, 230, 165 and 100 K. They decrease with temperature. Table 2 tells 2px,y and 4f electrons are transferred to 2pz, main contributor to the shortest covalent B-B bonds between B6 octahedra (called as B-Bout).

Electron accumulation at B-Bout enhances q1111 and q1122 of Ce at 165 K. They change signs and become more larger at 100 K. Negative and positive q1111 at 165 K and 100K indicates that vibration of Ce at body-centre to the face-centre is favorable at 165 K, and becomes unfavorable at 100 K. This is the reason for the enhanced peak in Fig. 8(b) and reduced one in (a). Negative q1122 at 100 K indicates the attractive force to the edge centre increased by the accumulated 2pz electrons. Since sum of the mean radii of B-2p and Ce-5p in free space, 2.2048+1.7947=3.9995 Å (Mann, 1968) is longer than Ce-B distance in Table 2, a slight expansion/contraction of them affects the potential seriously. Since the 5p orbitals are fully occupied, EDD's of them are spherical. B-2px is closer than 2pz to the spherically distributed Ce-5p3/2 electrons. κ(5p3/2) exhibits slight expansion on lowering the temperature. Expansion and contraction of 5p3/2and 4f5/2Γ7 in Table 2 reduces the potential of 4f orbitals lying closer to the nucleus. However, it contradicts to the decrease of 4f population. It indicates that the system resists against losing 4f electrons from Ce by making the potential of Ce more stable. It reminds us of the Le Chatelier’s law: the system resists against the change. Then why are 4f electrons allowed to flow out of Ce in spite of the reduced potential energy of the Γ8 state? The enhanced AHV at 165 and 100 K produced new ways of vibration along 100and 110directions increasing the ways to attain the energy of the system. Therefore, the enhanced AHV at lower temperatures means an increase in entropy. Since the electron transfer itself from Ce to (B-B)out increases entropy, it cannot be stopped.

The XAO analysis of CeB6 was applied to the weak-field model. The crystal structure, EDD, electron populations, expansion/contraction parameters and AHV parameters were quite consistent with each other at different temperatures. Therefore, the 4f-EDD is concluded to be measured and analyzed successfully.6.1.2. 4f population inversion and fully occupied 5d states at 430 K

The electron population at 298 K exhibited slightly occupied 4f5/2Γ7. Since the energy gap between 4f5/2Γ8 and Γ7 was reported to be 530-560 K (Loewenhaupt et al., 1985; Zirngiebl et al., 1984), XAO analysis at 430 K was performed to observe more electrons in Γ7 (Makita et al., 2007). For details of the high-temperature diffraction equipment, see PhD thesis of Makita (Makita, 2008). Scattering factors of Ce were evaluated from relativistic AO's calculated by the program GRASP (Dyall et al., 1989). Difference density around Ce on (001) after the spherical-atom refinement is shown in Fig. 9(a). The 4f-peaks closest to Ce elongate along <110> in contrast to those below room temperature in Fig. 8. They are surrounded by the four peaks of 0.65 eÅ-3, which extend along <100>, and disappeared after the analysis of the 4f peaks. However, the population of 4f5/2Γ8 and Γ7 were 0.06(3) and 0.36(11). Higher temperature reverses the populations of 4f5/2Γ8 and Γ7. The peaks outside the 4f-peaks, called 5d peaks, remained almost unchanged. Accordingly 5d5/2Γ7 and Γ8 orbitals were further refined. The population of 5d5/2Γ8, n(5d5/2Γ8), exceeded 1.0 and was fixed. n(4f5/2Γ8, Γ7) and n(5d5/2Γ8) were 0.06(2), 0.37(1) and 1.0, while 5d5/2Γ7 was vacant. R factor was reduced from 1.25 to 1.16 %. The peaks outside the 4f peaks reduced slightly from 0.65 to 0.59 eÅ-3 in Fig. 10(b). Since 5d orbitals extend in a large area, in contrast to 4f orbitals, slight change of EDD is significant. Thus electron densities around 5d–peaks in Figs. 9 (a) and (b) were numerically integrated to give 3.41 and 1.57 electrons. XAO analysis explained 54 % of the 5d-peaks. Therefore d5/2Γ8 is concluded to be occupied, though peaks still remained. It may be ascribed to the inaccurate 5d-AO's (Claiser et al., 2004).

Figure 9.

Differnece density around Ce at 430 K on (100) after (a) spherical-atom refinement and (b) final residual density. Contours at 0.2 eÅ-3. Negative and zero contours in dotted and thick lines.

Why are 5d orbitals occupied? The energy level of the Ce-5d and B-2p calculated by Liberman et al., (1971) and Mann (1967) are -0.63 and -0.62 a.u. They are very close compared to the Ce-4f level at -0.75 a.u. Thus B-2p electrons are transferred to 5d levels predominantly according to the first-order perturbation theory. The radial distribution functions of the relevant AO's are illustrated in Fig. 10. Putting the origin of B at r=3.045 Å, 2p was drawn to the reverse side. 5d and 2p overlap well, and it is expected that 2p electrons which distribute all over the crystal through the network of B-B covalent bonds can be transferred most probably to 5d orbitals. Electrons traveling in the crystal are expected to stay at 5d5/2Γ orbitals when they come close to Ce.

Figure 10.

Radial distribution functions of Ca and B.

Figure 11.

red) and 5d5/2Γ8 obritals.

Figure 12.

and 4f-energy levels expected from the electron populations.

Why does 4f5/2Γ87 inversion occur? 5d5/2Γ8 orbitals are located outside of 4f5/2Γ having exactly the same symmetry (Fig. 11). Since 5d5/2Γ8 is fully occupied, the potential of the 4f5/2Γ8 orbital is enhanced. It is the reason for the population inversion of 4f5/2Γ8 and Γ7. The energy level diagram expected from the electron populations obtained by the XAO analysis is shown in Fig.12. Since the quantization axes are defined parallel to <100>, T2g states locate higher than Eg. 5d3/2Γ8 orbital closely related to T2g is located highest among the 5d orbitals as shown in Fig. 12. From the populations obtained, 2p electrons seem to be transferred to 5d5/2Γ8 orbitals, directly or first to 5d3/2Γ8 and then to 5d5/2Γ8.

6.1.3. Electron configuration at 340 and 535 K

In order to confirm the 5d-occupation in CeB6 at 430 K, XAO analyses were done at 340 and 535 K. The difference densities, as well as electron configurations, are different from each other as shown in Fig. 13. Since the electrons are continually transferred between 2p and 5d orbitals, the crystal field seems to change with temperature resulting in different directions of the 4f and 5d peaks in Fig. 13. However, it is confirmed that 5d orbitals in CeB6 are occupied above room temperature (Makita et al., 2008).

Figure 13.

Bird-eye view of difference densities around Ce at (a) 340, (b) 430 and (c) 535 K. Contours as in Fig. 9 except the interval of 0.1 eÅ-3. Frames are parallel to <100>axes.

6.2. XAO analysis of SmB6

SmB6 formally has five 4f electrons. In order to extend XAO analysis to multi 4f-electron system, EDD of SmB6 was measured at 100, 165, 230 and 298 K (Funahashi, et al., 2010). It is interesting to see how physical properties as a Kondo insulator are explained by the XAO analysis. The shadow in the difference density at 230 K in Fig. 14 specifies roughly the area of 5d peaks. 5p3/2, 4f5/2Γ87 and AHV were refined reducing the 4f peaks in Fig. 14(a). However, since 4f-peaks still remained along <100>, the populations of 4f7/2Γ678 were refined but only Γ6 orbital which extends along <100> more sharply than 4f5/2Γ8 survived. Since peaks (peak A) remained in the 5d-area, 5d5/2Γ87 and 5d3/2Γ8 were added in the refinement and only 5d5/2Γ8 orbital, which stems from eg orbital of the storing field model had population. Fig. 14(b) exhibits that peak A is reduced from 0.43 to 0.17 eÅ-3 and 4f-peaks are almost deleted. Final parameters are listed in Table 3.

The electron configuration in Table 3 is correlated to the physical properties of SmB6 as follows: (a) SmB6 is a Kondo insulator. Its electric resistivity increases gradually like semiconductors below room temperature and begins to increase steeply below 30 K with a decrease in temperature. It also begins to increase like metals above room temperature (Ueda & Onuki, 1998). 2p and 5d3/2Γ8 orbitals consist of the conduction band (Kimura et al., 1990). Since B-2pz extends along the edge, it does not overlap with 5d3/2Γ8 effectively and does not seem to contribute to the band. The population of 2px in Table 3 as well as the resistivity steadily increases on lowering temperature. n(5d5/2Γ8) also increases from 230 K. The increase of populations indicates that of localized electrons. It may be correlated to the increase of the resistivity. (b) 4f7/2Γ6 are vacant only at 100 K, although they are occupied at the other three temperatures. It may be correlated to the band gap between the 4f states, which is reported to start developing between 150 and 100 K (Souma, et al., 2002). (c) Among the 5d orbitals only the 5d5/2Γ8 are occupied. Since 5d5/2Γ8 orbitals correspond to eg in the strong field model as illustrated in Fig. 1, it agrees with the band calculation of LaB6 by Harima (1988), reporting that 5d-eg and 2p of B consist of the conduction band.

Figure 14.

Difference density at 230 K after (a)spherical-atom refinement and (b) XAO analysis. Contours as in Fig. 13.

n(4f5/2Γ8)n(4f5/2Γ7)n(4f7/2Γ6)n(5d5/2Γ8)n(2px)n(2pz)100K0.96(9)0.48(12)0.01.00.21(8)0..45(17)165K1.00.67(13)0.44(12)0.75(24)0.16(13)0.55(39)230K1.00.61(13)0.59(13)0.56(16)0.13(7)0.66(28)298K0.77(3)0.59(7)0.47(6)0.78(6)0.07(3)0.81(12)

Table 3.

Temperature dependence of electron populations of Sm 4f, 5d and B 2p orbitals.

7. Bright future for X-ray crystallography

EDD investigation was limited up to 3d-transition-metal complexes. However, XAO analysis made EDD investigations of rare-earth compounds as well as non-stoichiometric ones possible. Its application to organic compounds can be attained when it is developed to X-ray molecular orbital analysis (XMO). Since the least-squares method stated in 2.2 was formulated for MO models, XMO analysis will be accomplished in a near future. When actinoid compounds become our targets, however, MD-effect is expected to be so much that the avoidance of it by ψrotation is impossible and the correction for it will be inevitable.

CeB6 is a possible quantum-material to emit UV light when electrons in 5d5/2Γ8 could be transferred to 4f5/2Γ8, as the investigation at 430 K revealed. The d-f transition is a permitted one by quantum mechanics. Since the 5d-occupation is found in the ground state of CeB6, some external force is necessary to make the transition occur. A electron populations in CeB6 and SmB6 found by the XAO analysis demostrate the importance of the EDD analysis based on quantum-mechanical orbitals.

As discussed in 5.3, the aspherical properties of EDD and AHV are separated better by the method based on classical Boltzmann statistics than the G-C method. However, recent development of neutron diffraction will make it possible to get intrinsic ADP's and use them as known parameters in X-ray EDD analyses. It will improve XAO analysis of rare-earth complexes and makes the XMO analysis surer and easier.

The accuracy of X-ray structure-factor measurements has been improved so much that every crystallographer will investigate EDD easily as a part of their X-ray structure analysis. The top-up operation with constant incident beam intensity at SR facilities has improved the accuracy of the structure factor measurements from 1 % to 0.1 %. Future of X-ray diffraction is bright.

© 2012 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Kiyoaki Tanaka and Yasuyuki Takenaka (September 19th 2012). XAO Analysis - AO's and Their Populations in Crystal Fields, Recent Advances in Crystallography, Jason B. Benedict, IntechOpen, DOI: 10.5772/50517. Available from:

chapter statistics

1273total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Computational Potential Energy Minimization Studies on the Prion AGAAAAGA Amyloid Fibril Molecular Structures

By Jiapu Zhang

Related Book

First chapter

Introductory Chapter: An Outline of Carbon Dioxide Chemistry, Uses and Technology

By Janah Shaya, Hassan Srour and Iyad Karamé

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us