Research Article  Open Access
Electronic Structure of the Cubic Compounds ReGa_{3} (Re = Er, Tm, Yb, and Lu)
Abstract
The electronic structure of ErGa_{3} and its isostructural compounds with Tm, Yb, and Lu are investigated with a highly accurate band structure scheme in LDA and GGA and warped muffintin approximation. In contrast to other investigations, the 4f electrons of the constituent Re are also treated as part of the valence bands. The position of the corresponding 4f bands relative to the Fermi energy strongly depends on the nuclear charge of Re. In Lu, they lie almost by 0.5 Ryd below and are extremely narrow. In Er, both in LDA and GGA, the 4f bands are found to be very close to the Fermi level . Assuming most of the 4f electrons to be part of the core removes the disagreement almost completely but produces a Fermi surface with a topology markedly different from that proposed in previous investigations. The intersections of the Fermi surface with planes are strongly varying within the Brillouin zone, they do not well match with the sparse experimental results. Investigations using the scheme as well as investigations of the dielectric response function are sketched.
1. Introduction
Studies on the electronic structure of cubic ReGa_{3} (Re denotes the heavy earth metals Er, Tm, Yb, and Lu, resp.) compounds which crystallize in the AuCu_{3} structure (see Figure 1) are a challenge not only as it is supposed that their magnetic properties are controlled by peculiarities of their Fermi surfaces but also as the role of the localized electrons is not well understood.
In previous investigations of ErGa_{3} and of TmGa_{3} [1โ6], which have been done to get detailed information about the shape of the Fermi surface, the electrons are assumed to be part of the core contributing to the charge density without interacting with the other valence electrons. Using the linear muffintin orbital method (LMTO) in the atomic sphere approximation (ASA), Pluzhnikov et al. [2] and Petukhov et al. [3] found in the range of the high as well as the intermediate de Haasvan Alphen (dHvA) frequencies a quite good agreement between their ab initio calculations and the experimental data. In order to analyze measurements of the twodimensional angular correlation of the annihilation radiation spectra [5] which give more details about the FS topology more modern schemes have been used: the fullpotential LMTO and the fullpotential linear augmented plane waves (FLAPW) schemes have been applied [6]. Disappointingly, none of these codes was able to produce a satisfying description of some experimental features that were in satisfying agreement with the previous LMTOASA results [2, 3]. It was supposed [6] that the most probable reason for this difference lies in the choice of the linearization energy around which the expansion of the energydependent radial function is performed, in LMTO as well as in FLAPW.
It is a question of principle importance whether the possible breakdown of the LDA in the cases of the compounds ErGa_{3} and of TmGa_{3} is a failure of the LDA scheme or if it is caused by fact that in these compounds the bands fortuitously coincide with the Fermi level and thus yield a density of states unrealistically enhanced at . In the systematic analysis of the electronic structure of compounds consisting of rare earth metals neighboured in the periodic table done in present work, the second alternative turned out to be true.
The present investigations use the Modified Augmented Plane Wave (MAPW) method [7] scheme. It is a linearized version of Slater's APW method but differs from LAPW in that within the atomic spheres the radial part of the Bloch functions is a combination of radial functions , counting the energies chosen within a broad range, which distinctly differ from each other. Thus, the abovementioned problems of LAPW are avoided. In the context of the LDA DFT formalism, the MAPW scheme is almost free of any systematic restrictions and guarantees any desired accuracy by a suitable choice of the ansatz functions and of other intrinsic parameters.
This paper is organized as follows. Section 2 starts with a comparison of the essentials of the LAPW and the MAPW scheme. It is explained why the latter guarantees the high accuracy needed to properly describe the electronic structure of ErGa_{3}. The progress of selfconsistent calculations incorporating many valence electrons, especially in the case of narrow bands just below , sensitively depends on the choice of the starting potential. Therefore, in Section 3, our study of the series of RE compounds were started with the compound LuGa_{3} characterized by fully occupied bands. Thus, we could avoid that the iteration mechanism stops at an incidental minimum. By reducing the Re core charge in noninteger steps, we could show how the peak in the density of states curve approaches the Fermi level. This strategy actually makes selfconsistent calculations faster and more stable because we remain near an absolute minimum in the high dimensional potential landscape. Then, the band structure, the shape of the Fermi surfaces including extremal areas, and cyclotron masses are discussed. In Section 4, investigations of the dielectric response function are sketched.
2. Basic Theoretical Aspects
2.1. Relativistic Effects
The heavy Re rare earth atom makes a relativistic treatment necessary. In previous investigations dealing with this aspect, the core and the valence electrons have been treated by completely different schemes. The atomlike core states were obtained by use of the fourcomponent Dirac equation whereas the valence states were approximated by a twocomponent equation that has been derived, more or less, via successive applications of unitary, normconserving FoldyWouthuysen transformations. In the simplest case, the socalled scalar relativistic approach, a modified Pauli or radial Schrรถdinger equation containing the relativistic mass correction and the Darwin term, is solved within the spherical regions [2, 4]. More refined investigations [5, 6] make the expectation value of a Hamiltonian containing the spinorbit coupling stationary, with a Ritz ansatz obtained in the scalarrelativistic approximation. It is obvious that the orthogonality of the core states to the valence states is destroyed by the use of completely different Hamiltonians. Test calculations on Au [8] that go beyond the approximations just described (because they are based on nonrestricted twocomponent spinors) yielded the result that this improper treatment of the relativistic effects produces uncontrolled errors in the charge density and the total energy. Therefore, as long as a fully relativistic treatment of compounds similar to the treatment of solids containing only one atom in the cell [9] is outstanding, a nonrelativistic is regarded as more reliable.
2.2. LAPW versus MAPW: A Critical Comparison
It is common to both the LAPW and the MAPW [7] methods that the Bloch functions outside the APWspheres are approximated by superpositions of plane waves, and that within the APW spheres, the plane waves are suitably augmented, yielding a linear eigenvalue problem. This has been done in different ways. In LAPW, the wave function consists of a sum of products of a spherical harmonic and a linear combination of two radial functions. The first is a solution of the radial differential equation inside the atomic sphere with an averaged spherical potential and a linearization energy , around which the expansion of the energy dependent wave functions is performed. The second one is the energy derivative of the radial function. Up to the maximal value of the angular momentum , the linear combination at the surface of the APWsphere, joins continuously to the angular momentum expansion of the corresponding plane wave. As for all , no continuity exists, and the truncation error can be made small by choosing, large , for example, 10. As a consequence of this choice of the radial functions, the solutions of the eigenvalue problem are most reliable at energies around . The range of validity around can be quite small when the branches of the logarithmic derivatives follow in rapid succession as is the case in ErGa_{3}. Because the sets of eigenfunctions obtained with different values of are not strictly orthogonal to each other, a change of the relevant energy is not a good remedy.
In the MAPW scheme, a set of radial functions is generated in the spherical averaged potential, mostly by requiring that their logarithmic derivative is either or โ1. To avoid any truncation error, the augmentation of the plane waves is only performed for the leading angular numbers , whereas for , the spherical Bessel functions are kept. In contrast to all other versions of the APW scheme, the full wave function and its derivative are made exactly continuous on the surface of the APW spheres by the use of additional constraints. To account for the special role of the electrons, the basis of the radial functions within the Re spheres contains one radial function that is optimally localized by an appropriate choice of the radial energy, implying that the value of the radial wave function on the APW sphere is smaller by almost a factor of 100 than the others. Further radial functions are chosen according to the previously mentioned recipe, guaranteeing the continuity of the the corresponding partial wave.
2.3. Details of the MAPW Calculation
The following considerations are based on the assumption that the Re core has a Pdtype configuration consisting of 46 electrons. We have found that the states up to 4d are welllocalized; for example, their radial functions assume values less then on the surface of the Er sphere, and need not to be considered in the MAPW calculations. Nevertheless, they are numerically orthogonal to the wave functions of the valence states because they are not kept frozen during the selfconsistency cycles. In contrast, the best localized state of Er assumes the value 0.032413โa.u. on the Er sphere confirming our reservation against treating this state as part of the core.
The quantum number of the angular decomposition of the Bloch functions is restricted to be less than or equal to 2 within the Ga spheres and to be less than or equal to 3 within the Re spheres. Four radial functions are used for each of the . The proper choice of the radial functions needs special care as the corresponding logarithmic derivative is strongly varying in the energy region near the Fermi level.
The number of plane waves used to describe a Bloch state within the whole atomic polyhedron is restricted by the inequality where denotes the vector of a simple cubic reciprocal lattice, and a is the lattice constant. In contrast to the planewave cutoff mostly applied in LAPW it guarantees the full point symmetry of the Bloch energies at any point of the Brillouin zone (BZ). Depending on the specific value of the wave vector , this restriction yields a superposition consisting of between 160 and 181 plane waves. In total, each ansatz consists of more than 320 trial functions. It guarantees sufficient accuracy of the Fermi energy and of the total energy in the SCF cycles and leaves the topology of the band structure almost stationary, that is, no jumps are visible in the intersection of the Fermi surface with certain planes, as shown in Figures 7, 8, 9, 10, and 12.
The Brillouin zone integrations over all occupied states yielding the electron density , the Fermi energy , and the total energy are approximated by a sum over points of a simple cubic lattice of length within one octant of the BZ. Because the evaluation of requires those operations of the subgroup of the point group that leave the atoms of the basis invariant, a further restriction to the irreducible wedge of the fcc BZ is not possible, at least in the SCF cycles.
The characteristic groundstate properties are calculated by use of elaborate exchangecorrelation functionals in LDA [10] and GGA [11โ14]. Without any further truncation, exchange and correlation within the atomic spheres are obtained from the charge density along the special directions. Outside the spheres, the charge density is evaluated on a fine mesh of suitably chosen points [15] which allows to determine exchange and correlation in any accuracy wanted. Further details of the present calculation are described in previous publications [16, 17].
3. Results
3.1. Trends in the Electronic Structure of the Isostructural Compounds ErGa_{3}, TmGa_{3}, YbGa_{3}, and LuGa_{3} in LDA
The following investigations aim to provide an overall insight into the influence of the electrons. To avoid any influence of the lattice spacing, all investigations have been done with the lattice constant of ErGa_{3} at zero temperature, โร [4]. The narrow Re5s and 5p bands with energies around โ0.5764โRyd and 0.7636โRyd, respectively, in the case of ErGa_{3} are well separated from the complex of the other bands originating from the atomic Re 6s and 5dstates and are not to be considered.
In the course of these considerations, we have found that the electronic structure of the compound LuGa_{3} is standing out for its relative simplicity. As shown in Figure 2, the electrons occupy a band of extremely small width of 3โmRyd located by more than half of aโRydberg below the Fermi level. It is split according to the cubic symmetry and additionally hybridized with the single band starting at the lowest state. It is remarkable that it leaves the other bands almost unchanged with the favourable effect on the selfconsistency that a small number of cycles is sufficient. Thus, the final potential is best suited as a reliable starting point for the investigations of the other isostructural compounds.
Table 1 lists some characteristic quantities found by successive decrease of the nuclear charge of the Re constituent, , in noninteger steps, starting from the last line. It causes the energy of the best localized radial state, , continuously to rise whereas the Fermi energy slightly diminishes. Both energies get closer and closer but do not cross. This change is accompanied by a strong increase of the hybridization of the states with the other valence bands, although the wave function of the best localized radial state at the APW sphere only moderately increases from 0.0208โa.u. to 0.0323โa.u. in ErGa_{3}. Correspondingly, raises but does not tend to the unphysical values cited in [2]. There is no doubt that in the cases of TmGa_{3} and ErGa_{3} the corresponding values are still too large compared with experimental results. Figures 4, 5, and 6 illustrate how the Re nuclear charge influences the band structure along the lines of high symmetry in the Brillouin zone. In their overall look, all Figures are rather similar; however, significant differences are visible in the range of the bands and their location relative to the Fermi level.

To emphasize the role of the Re bands as possible sinks of electrons, we have calculated the f orbital projected density. The number of felectrons per rareearth atom, , listed in Table 1 in the last but one column, increases up to almost in the same measure as the virtual core charge and reversely the Fermi energy decreases. Above this value, all like states are occupied and the additional electrons occupy the spd like states which are less and less hybridized with the states. The Fermi energy starts to rise and approaches moderate values. This change occurs at a noninteger value of : the compound YbGa_{3} stands on one side and the compound LuGa_{3} on the other, and consequently their crystal potentials are significantly different. According to Table 1, the projected density substantially contributes to the large value of especially in the case of ErGa_{3} and TmGa_{3}.
These considerations are further illustrated by the plots of the density of states (DOS) displayed in the right panels, respectively. The sharp peak caused by the band approaches the Fermi level with decreasing values of . In ErGa_{3}, the energy distance amounts up to 4โmRyd. The non part of the DOS is rather similar in the four compounds and is quite smooth. From these results, we learn that the large value of , especially in the case of ErGa_{3}, in contradiction to experimental facts is not a systematic deficit of LDA, but it is almost accidental as it locates the position of the bands near to the Fermi level. Only a small shift of the position of the bands by about 5โmRyd will improve the agreement with the experimental data considerably. On the other hand, these results do not support at all the concept of assuming 11 or 12 electrons as part of the core [2, 4].
3.2. Fermi Surface in Strict LDA and GGA
It is evident that the position of the bands relative to the Fermi level has significant influence on the shape of the Fermi surface. To tie on to previous investigations [2, 4โ6], we first consider intersections of the Fermi surface with high symmetry planes, the panels ฮXMX, RMXM, and ฮXRM.
In the case of LuGa_{3}, displayed in Figure 7, we find a certain resemblance to the FLAPW results derived for ErGa_{3} [6], especially in the lower panel showing the intersection with the (011) plane. The intersection with the plane displayed in the upper panel is a bounded contour of almost elliptic shape whereas in FLAPW the central contour merges with the contours encircling the occupied states around the Xpoints. The LMTOASA results of ErGa_{3} discloses a minor similarity as the contours of the occupied states around the M points extend far in the direction of the ฮ point. In YbGa_{3}, the intersections look quite similar.
The contour lines of ErGa_{3} shown in Figure 8 disclose a completely different shape of the Fermi surface. This already follows from the band structure displayed in Figure 6 which shows that twofold degenerated sheets of the Fermi surface cross the ฮ and the Tlines producing cusps on the respective boundary lines. Contours around the M and Rpoints encircle occupied states whereas the surrounding of the ฮ and Xpoints are not occupied. The lines encircling the points M, R, and X are the tracks of electron sheets, whereas the line around ฮ is the track of a hole sheet. It is remarkable that the intersections of the Fermi surface for TmGa_{3}, displayed in Figure 9, are quite similar to that of ErGa_{3} although their densities of states differ by more than a factor 3.
An extensive analysis in 3 dimensions which is straightforward but quite arduous had the result that the Fermi surface of ErGa_{3} consists of four sheets: three, more or less, spherically shaped surfaces denoted by b_{1}, b_{2}, and b_{3} all having the full cubic point symmetry . The hole FS b_{1} is centered at the point whereas the electron surfaces b_{2} and b_{3} are located at the point. Finally, the FS b_{4} is multiplyconnected and extends over the whole BZ and is best described by its intersections with planes . displayed in Figure 10 in steps of . For completeness, the crosssections with other surfaces are also shown. From this Figure, we learn that the FSs b_{1}, b_{2}, and b_{3} are monotonically shrinking with increasing value of , indicating that the crosssections are extremal at the surface of the first BZ. According to Figure 10, the FS b_{1} is strongly indented in the (111) direction. The FS b_{4} continuously changes its shape. Up to it is outside of the cylinder around , then rapid changes of its shape occur and beyond , it is again outside of the cylinder around . It is evident that the intersections with planes oriented in (110) and (111) with the FS b_{4} also reflect this great variety of orbits. The investigations in Section 3.3 will have similar results. Intersections of the type b_{1}, b_{2}, or b_{3} which are easy to comprehend, and those of type b_{4} which show a great diversity of intersections. These findings distinctly demonstrate that the customary plots [2, 4โ6] of the panels ฮXMX, RMXM and ฮXRM do not capture the great variety of orbits which change rapidly mostly inside the Brillouin zone. Our results also cast doubts on the methods usually applied to reconstruct the Fermi surfaces of such compounds like ErGa_{3} directly from experiments (de Haasvan Alphen, cyclotron resonance, electronpositron annihilation measurements) which are mostly based on the assumption that the FS has some welldefined crosssections [19].
The GGA using the functionals [11โ14] almost rigidly shifts the band structure of ErGa to lower energies by 50.0 mRyd but has only minor influence on the shape of the Fermi surface. The energy of the bestlocalized radial states is even closer to the Fermi energy than in LDA, that is, 1.8 mRyd, with the fatal consequence that the density of states at the Fermi level assumes the unphysical value of 195.2 electrons cell^{1}โRyd^{1}. Consequently, sections of the Fermi surface with the high symmetry planes of the Brillouin zone look quite similar to those displayed in Figure 8. Hence, GGA does not improve agreement with the experimental data at all.
3.3. Exclusion of the Hybridization of the Valence Bands by the Electrons: Following Previous Tracks
Our investigations have the result that, based on experimental evidence [4, 6], the threedimensional mapping of the Fermi surface of ErGa_{3} could not be explained by highaccurate LDAMAPW calculations if the electrons are considered to be constituents of the valence bands. The strong influence of the bands on the topology of the Fermi surface of ErGa_{3} and TmGa_{3} is avoided by treating them as part of the core as has been done in previous investigations [2, 4โ6]. Then, in the MAPW scheme, the quantum number of the angular decomposition of the Bloch functions may also be restricted to be less than or equal 2 within the Er or Tm spheres. Selfconsistent calculations yielded the result that the energy of the orbitals is below the Fermi level by half aโRydberg in the case of the Er core whereas it is distinctly above the Fermi level when an additional electron is assumed as belonging to the core. The band structure of the Er configuration, displayed in Figure 11, as well as that of the Tm configuration shows a great similarity with that of the LuGa_{3}, but they differ markedly from the band structure found when the electrons are assumed to be part of the valence electrons (see Figure 6). Consequently, the Fermi topology is completely changed. This is illustrated in Figure 12, showing the contour lines of the Fermi surface. Compared with FLAPW results obtained with the same concept [6], the sections with the (110) plane displayed in the XMRM face are quite similar, but some sections with the (100) plane are completely different. A coarse resemblance to the LMTOASA [4] result also exists. In contrast, the contour plots of the Er , not shown here, disclose a completely different Fermi surface topology for which we have no experimental evidence. The plot of the density of states displayed in the right panel of Figure 11 shows some vanHove singularities with modest heights.
The values of are in a reasonable range, for example, 10.34 and 9.60 electrons cell^{1}โRyd^{1} in the Er and the Tm configurations, respectively, and are even smaller by 30% if the number of the electrons in the core is increased by one. From these results, we learn that within a heuristic approach, we could achieve closer agreement with the low temperature specific heat data by choosing the nonintegral parameter of the Er core configuration to be near the value 11.0.
The Fermi surface in the Er core configuration consists of four sheets: B_{1} is almost spherical in shape around the point , slightly dented in (011), with a mean radius . B_{2} is a prolate spheroid around the point with the longest axis in the (001) direction. The lengths of its axis are approximately and . B_{3} is almost ellipsoidally shaped around the point (0.2730,0.2730,0) with the longest axis oriented in (110) directions. The lengths of its axes are approximately , , and . Finally, B_{4} is multiplyconnected and again is best described by its intersections with planes , as displayed in Figure 13 together with those of the other Fermi surfaces. Similar to the FS b_{4} shown in Figure 10, its shape exhibits rapid variations within the BZ. According to Figure 12, besides the contours of B_{2} and B_{3}, two contours of B_{4} exist near the face ฮXMX. One is almost circular the centre being over the point M whereas the other has fourfold symmetry around the point . With increasing value of the component , the contours of B_{2} and B_{3} shrink whereas those of B_{4} grow, merge for , and then split into two other parts with contours centred above X points. Being well separated in the range , both curves merge again at and form considerably larger contours centred below the R point for . For still larger values of , the size of the contours, again centred around the line gradually shrinks. The intersections of the sheet with planes normal to the (110) and (111) directions even show a still greater diversity. Besides closed contours, now in certain ranges of the projection of on the normal direction, contours exist that extend over the whole reciprocal lattice space, yielding the socalled open orbits. In general, the shape of the contours shows rapid changes, especially near .
It is surprising that even in the present case, considering the electrons as part of the core possible crosssections of FS can rapidly change inside the BZ, and they are not at all described by the usual plot of the contour line only showing the faces , , and . We suspect that a thorough analysis using the LAPW scheme will produce a similar behaviour provided the plane waves satisfy a criterion analogous to that quoted in Section 2.3.
3.4. Extremal Areas and Cyclotron Masses of ErGa_{3}
The extremal crosssectional areas and the corresponding cyclotron masses for the closed Fermi sheets are listed in Table 2 (the areas of the crosssections are related to the magnetic field B, in gauss, by the relation 10^{โ9} B in the case of the lattice constant cited by Pluzhnikov et al. [4]). Because of their spherical shapes, the areas of the Fermi sheet b_{1} are slightly orientationally dependent, within a margin of 10%. Both FS b_{2} and b_{3} have, as a consequence of their dimension, such small areas that they need no further consideration. As a consequence of the strong change of the shape of the FS within the BZ, the areas and the cyclotron masses of b_{4} listed in Table 3 cover a wide spectrum of values strongly depending on the vector pointing to the centre of the corresponding contour. It is questionable whether the usual approach which simply connects the de Haasvan Alphen frequencies with extremal areas of the FS, is applicable [19] when the crosssections show the strong variations found in the previous sections.


The Fermi surface and the cyclotron masses of ErGa_{3} were investigated by Pluzhnikov et al. [4] by de Haasvan Alphen measurements in sufficiently strong magnetic fields destroying the antiferromagnetic phase. They found that the measured dHvA frequencies arranged in four branches that are orientationally independent to within a margin of at most 2%. The highest branch denoted by can be uniquely attributed to the FS b_{1}; the measured values 0.4211, 0.4100, and 0,4210, in units of in (001), (110), and (111), respectively, agree satisfactorily with the corresponding values in Table 2. The following two branches, denoted by , , and , are due to the FS b_{4}. However, the values listed in Table 3 cover a far broader range than the experimental values, all in the same units, 0.1763, 0.1298, 0.0519 in (001), 0.0678, 0.0649, 0.0513 in (110), and 0.1522 in the (111)direction, and are far more orientationally independent. Finally, the FS b_{2}, b_{3}, and some of the contour lines of b_{4} produce the lowfrequency branch, . The good agreement between the experimental data with the ab initio calculations (carried out using the linear muffintin orbital method in the atomic sphere approximation with the Er electrons as part of the core) performed by Pluzhnikov et al. [4] seems to be rather fortuitous in the light of our more accurate investigations, according to which only the orbit around R is expected to be orientationally independent and is characterized by a high dHvA frequency. Finally, the large discrepancies between the measured and theoretical values of the cyclotron masses found in both investigations can hardly be explained by manybody enhancement, as has been done by the previous authors [4], but raises the suspicion that the high magnetic fields applied to destroy the antiferromagnetic phase of ErGa_{3} have a nonnegligible influence on the electronic structure.
3.5. LDA + U Investigations [20]
A simple way to describe the correlation between the localized Er electrons beyond the usual LDA is provided by the socalled โLDA + Uโ totalenergy functional [21โ23]. We have found that the MAPW scheme is especially suited for the implementation of this functional. Due to the lack of detailed information about the screened interaction, we assumed that the effective onsite Coulomb interaction has the form of a screened Coulomb potential with the screening length q as a free parameter. As desired, the reasonable value of โBohr shifts the narrow bands by โ 0.3โRy below the Fermi level and leaves the other bands almost unchanged with the consequence that is found to be in quite good agreement with the experimental data. Details concerning the band structure are described in [20] giving evidence that the band structure is quite similar to those in LDA displayed in Figure 6 apart form the line MX. But it shows no similarity at all with the frozen core results described in Section 3.3.
4. A Sketch of the Investigations of the Dielectric Response
The response to a scalar electric field is fully described in the framework of the timedependent densityfunctional theory by the the KohnSham response function [24โ27]. In the case of a slowly varying perturbation, it uniquely splits into an intraband and interband contribution. The first can be expressed by a sum over the Fermi surface, Einstein summation convention over and . This relation is physically interesting: the mean value of the inner product of the velocities of the Bloch electrons at the Fermi surface essentially determine the dc behaviour. A highly simplified treatment of the dcconductivity based on the assumption of a relaxation time describing the scattering of the electrons by lattice defects or phonons has the result volume of the elementary cell. This expression is based on more general assumptions than those leading to the wellknown Drude formula: (i) the solution of the Boltzmann equation by a dependent relaxation time and (ii) Kohler's variational principle [28, 29] using the trial function grad .
Provided the relaxation can be assumed to be constant the cofactor of in (2) gives an estimate of the static conductivity: it roughly depends on the product of the inverse mass and the density of electrons.
In Table 4, a diagonal element of the inverse mass averaged over the occupied states, the mean value of grad over the Fermi surface and the plasma frequencies and are listed. A crude estimate of the dc conductivity of these compounds has been obtained by using a relaxation time which appropriately explains the room temperature conductivity of polycrystalline Cu, say โ(ฮฉcm)^{(1)}. The last column lists the ratio of these conductivities. From these results, we learn that asymptotic behaviour of is dominated by the interband contribution as the intraband contribution is smaller up to one order. The effective inverse mass and the inner product of grad show a stronger variation; in the compounds with Lu and Yb, the corresponding values are comparable with those of Cu, but in the compounds with Tm and Er the inner product of grad is up to three order of magnitudes smaller with the consequence that in this context ErGa_{3} turns out to be a semimetal. The close similarity of the intersections with the Fermi surface displayed in Figures 7 and 11 arouses the suspicion that the model considering the electrons as part of the core will produce a dcconductivity comparable with Cu.

In the longwave limit, the interband contribution of the response function reduces to the sum are the matrix elements of the KohnSham Bloch functions and , the corresponding eigenvalues and occupation numbers. The factor 2 in front accounts for the spin degeneracy. The sum over the wave vector runs over the BZ whereas over n and extends over all Bloch states obtained in the KS scheme for a fixed value of , occupied or nonoccupied, in present case up to 240 states with energies up to 250โRyd above the Fermi level. This requirement has often been overseen in previous investigations. As usual this expression may be decomposed in a real and in an imaginary part. By using the identity which corresponds to the ThomasReicheKuhn sum rule in atom physics [30] at large values of the real part of the interband contributions reduces to in the longwave limit. The abbreviation is the plasma frequency of a free electron gas, which only depends on the density of the valence electrons . Thanks to the completeness of the Bloch functions (see (5)), all specific information for example, the band structure or the matrix elements of the momentum operator has no influence on the asymptotic behaviour of the interband contributions at long wavelength.
Partial integration of (2) has the result that has the same asymptotic behaviour with plasma frequencies analogously defined as in (7) with the sole difference that is substituted by the mean value of the tensor of the inverse mass. Both objects are only abbreviations which have nothing to do with collective excitations of the valence electrons. As consequence of the causality the infinite integral over the imaginary part at long wave length again gives the square of the free plasma frequency, It can be considered as a sum rule of a suitably defined conductivity [31].
The numerical work is largely analogous to previous investigations [32, 33]. The eigenvalues , their gradients, and the matrixelements of the momentum operator were evaluated using the MAPW results without any further approximation. As a check, the sum rule (5) was used. Figure 14 shows the frequency dependence of the conductivity in the longwave limit in the region โRydberg, on the right side the real part and on the left side the real part of the reciprocal conductivity, respectively, formally defined by . This expression is sensitive to collective excitations of the valence electrons similar to the energy loss function in the conventional theory.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
At low values of , the series of compounds shows a margin difference: ErGa_{3} and TmGa_{3} have a structure with two wellseparated peaks whereas the other two compounds have a singular peak only. Above 0.25โRyd, the real parts of the conductivity of all four compounds look quite similar. Up to 1โRyd, the reciprocal conductivity is without any structures. The peak structures near 1.6โRyd, especially pronounced in LuGa_{3}, indicated the possibility of collective excitations of the valence electrons of plasmon type.
5. Summary
The electronic structure of the compounds ReGa_{3} is found to be very sensitive to the constituent Re due to the location of the bands relative to the Fermi level. In the heaviest compound, Lu, the bands are quite small and far from the Fermi level . The remaining valence bands have the characteristic features of an spd complex and the sections of the Fermi surface with high symmetry planes of the BZ are quite similar to the results obtained for ErGa_{3} [4, 5] and TmGa_{3} [2] when the bands are assumed to be part of the cores. With decreasing core charge these bands approach and strongly interact with the other valence bands. Thus the topological structure of the Fermi surface is drastically changed and has no similarity at all with the momentum density obtained by deconvoluting the measurements of the angular correlation of the electronpositron annihilation radiation. In ErGa_{3} and TmGa_{3}, the closeness of the bands to the Fermi level produces values of which are not realistic. GGA investigations enforce the disagreement with the experimental results.
Analogous to previous investigations, the unrealistic high value of is avoided by considering 11 electron as part of the core. A detailed analysis covering the whole Brillouin zone shows that the Fermi surfaces of the compounds ErGa_{3} and TmGa_{3} consist of different sheets allowing a great variety of orbits in magnetic fields. One of them is multiconnected, and the usual plots showing intersections with high symmetry planes do not all describe its complex shape. Therefore, it is questionable whether the Fermi surface may be solely reconstructed from measurements. It is suspected that highprecision investigations of other compounds with a nonsimple basis will yield similar results.
With regard to the ground state properties and the response to an external scalar potential, the compounds of this series behave quite similarly up to one distinct exception: in ErGa_{3} and TmGa_{3}, the flat bands are found so close to the Fermi level that they strongly influence the shape of the Fermi surface and the inverse band mass. In DFT, both compounds are semimetals. This obvious defect is cured by use of the LDA + U scheme which shifts the bands slightly below the Fermi level [20].
Acknowledgments
Dr. R. Bader, LRZ Mรผnchen, and Dr. H. Stรถhr have critically read the paper. The generous hospitality of Professor J. van Delft and Professor U. Schollwรถck is gratefully acknowledged.
References
 V. B. Pluzhnikov, A. Czopnik, and I. V. Svechkarev, โde Haasvan Alphen effect in ScGa_{3}, LuGa_{3} and YIn_{3},โ Physica B, vol. 212, no. 4, pp. 375โ378, 1995. View at: Google Scholar
 V. B. Pluzhnikov, A. Czopnik, G. E. Grechnev, N. V. Savchenko, and W. Suski, โBand structure and Fermi surface of TmGa_{3},โ Physical Review B, vol. 59, no. 12, pp. 7893โ7900, 1999. View at: Google Scholar
 A. G. Petukhov, W. R. L. Lambrecht, and B. Segall, โElectronic structure of rareearth pnictides,โ Physical Review B, vol. 53, no. 8, pp. 4324โ4329, 1996. View at: Google Scholar
 V. B. Pluzhnikov, A. Czopnik, and G. E. Grechnev, โThe Fermi surface of ErGa_{3},โ Journal of Physics Condensed Matter, vol. 11, no. 23, pp. 4507โ4516, 1999. View at: Google Scholar
 M. Biasini, G. Ferro, G. KontrymSznajd, and A. Czopnik, โFermi surface nesting and magnetic structure of ErGa_{3},โ Physical Review B, vol. 66, no. 7, Article ID 075126, pp. 751261โ751268, 2002. View at: Google Scholar
 G. KontrymSznajd, M. SamselCzekała, G. E. Grechnev, and H. Sormann, โFermi surface of ErGa_{3},โ Physica Status Solidi C, vol. 4, no. 10, pp. 3879โ3882, 2007. View at: Publisher Site  Google Scholar
 H. Bross, G. Bohn, G. Meister, W. Schubö, and H. Stöhr, โNew version of the modified augmentedplanewave method,โ Physical Review B, vol. 2, no. 8, pp. 3098โ3103, 1970. View at: Publisher Site  Google Scholar
 H. Reinisch and H. Bross, โRelativistic density functional calculation of the total energy and Fermi surface of gold,โ Zeitschrift für Physik B, vol. 95, no. 2, pp. 145โ150, 1994. View at: Publisher Site  Google Scholar
 V. Theileis and H. Bross, โRelativistic modified augmented plane wave method and its application to the electronic structure of gold and platinum,โ Physical Review B, vol. 62, no. 20, pp. 13338โ13346, 2000. View at: Publisher Site  Google Scholar
 S. H. Vosko, L. Wilk, and M. Nusair, โAccurate spindependent electron liquid correlation energies for local spindensity calculations—a critical analysis,โ Canadian Journal of Physics, vol. 58, no. 8, pp. 1200โ1211, 1980. View at: Google Scholar
 J. P. Perdew, J. A. Chevary, S. H. Vosko et al., โAtoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation,โ Physical Review B, vol. 46, no. 11, pp. 6671โ6687, 1992. View at: Publisher Site  Google Scholar
 J. P. Perdew, K. Burke, and M. Ernzerhof, โGeneralized gradient approximation made simple,โ Physical Review Letters, vol. 77, no. 18, pp. 3865โ3868, 1996. View at: Google Scholar
 J. P. Perdew, โDensityfunctional approximation for the correlation energy of the inhomogeneous electron gas,โ Physical Review B, vol. 33, no. 12, pp. 8822โ8824, 1986. View at: Publisher Site  Google Scholar
 J. P. Perdew, โErratum: densityfunctional approximation for the correlation energy of the inhomogeneous electron gas,โ Physical Review B, vol. 34, no. 10, p. 7406, 1986. View at: Publisher Site  Google Scholar
 H. Bross, โThe local density approximation limit of the momentum density and the Compton profiles of Al,โ Journal of Physics Condensed Matter, vol. 16, no. 41, pp. 7363โ7378, 2004. View at: Publisher Site  Google Scholar
 H. Bross and R. Eder, โSelfconsistent MAPW calculation with a warped muffintin potential. I. The electronic structure of Al and its pressure dependence,โ Physica Status Solidi B, vol. 144, no. 1, pp. 175โ193, 1987. View at: Google Scholar
 H. Bross and R. Stryczek, โSelfconsistent MAPW calculation with a warped muffintin potential. II. The electronic structure of Li and its pressure dependence,โ Physica Status Solidi B, vol. 144, no. 2, pp. 675โ692, 1987. View at: Google Scholar
 H. Jones, The Theory of Brillouin Zones and Electronic States in Crystals, NorthHolland, Amsterdam, The Netherlands, 1962.
 I. M. Lifshitz and A. M. Kosevich, โQuantum Theory of the Electrical Conductivity of Metals in a Magnetic Field,โ Journal of Experimental and Theoretical Physics, vol. 2, p. 636, 1956. View at: Google Scholar
 to be published.
 V. I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein, โFirstprinciples calculations of the electronic structure and spectra of strongly correlated systems: the LDA + U method,โ Journal of Physics Condensed Matter, vol. 9, no. 4, pp. 767โ808, 1997. View at: Google Scholar
 A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, โDensityfunctional theory and strong interactions: orbital ordering in MottHubbard insulators,โ Physical Review B, vol. 52, no. 8, pp. R5467โR5470, 1995. View at: Publisher Site  Google Scholar
 A. B. Shick, A. I. Liechtenstein, and W. E. Pickett, โImplementation of the LDA+U method using the fullpotential linearized augmented planewave basis,โ Physical Review B, vol. 60, no. 15, pp. 10763โ10769, 1999. View at: Google Scholar
 E. K. U. Gross and W. Kohn, โLocal densityfunctional theory of frequencydependent linear response,โ Physical Review Letters, vol. 55, no. 26, pp. 2850โ2852, 1985. View at: Publisher Site  Google Scholar
 E. K. U. Gross, โErratum: local densityfunctional theory of frequencydependent linear response,โ Physical Review Letters, vol. 57, no. 7, p. 923, 1986. View at: Publisher Site  Google Scholar
 E. K. U. Gross and W. Kohn, โDensity Functional Theory of ManyFermionSystems,โ in Advances in Quantum Chemistry, S. B. Trickey, Ed., vol. 21, p. 255, Academic Press, San Diego, Calif, USA, 1990. View at: Google Scholar
 E. K. U. Gross, J. F. Dobson, and M. Petersilka, โDensityfunctional theory of timedependent phenomena,โ in Topics in Current Chemistry, R. F. Nalewajski, Ed., vol. 181, pp. 81โ172, Springer, Berlin, Germany, month 1996. View at: Google Scholar
 M. Kohler, โBehandlung von Nichtgleichgewichtsvorgängen mit Hilfe eines Extremalprinzips,โ Zeitschrift für Physik, vol. 124, no. 712, pp. 772โ789, 1948. View at: Publisher Site  Google Scholar
 J. M. Ziman, โElectrons,โ in The Physics of Metals, J. M. Ziman, Ed., chapter 5, Cambridge University Press, Cambridge, UK, 1969. View at: Google Scholar
 W. Thomas, โÜber die Zahl der Dispersionselektronen, die einem stationären Zustande zugeordnet sind. (Vorläufige Mitteilung),โ Die Naturwissenschaften, vol. 13, no. 28, p. 627, 1925. View at: Publisher Site  Google Scholar
 D. Pines and P. Nozieres, โNormal Fermi Liquids,โ The Theory of Quantum Liquids, vol. 1, W. A. Benjamin, Menlo Park, Calif, USA, 1966. View at: Google Scholar
 H. Bross, O. Belhachemi, B. Mekki, and A. Seoud, โLocal field corrections and the influence of exchange and correlation on the microscopic longitudinal dielectric function of copper,โ Journal of Physics, vol. 2, no. 17, pp. 3919โ3939, 1990. View at: Publisher Site  Google Scholar
 H. Bross and M. Ehrensperger, โCalculation of the Dielectric Function of Li Metals,โ Zeitschrift für Physik, vol. 97, p. 17, 1995. View at: Google Scholar
 P. T. Landsberg, โSimple Perturbation Methods in Band Theory,โ in Solid State Theory, P. T. Landsberg, Ed., chapter B17, WileyInterscience, London, UK, 1969. View at: Google Scholar
Copyright
Copyright © 2011 Helmut Bross. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.