Electron liquid state in the symmetric Anderson lattice

Using mean field approach, we provide analytical and numerical solution of the symmetric Anderson lattice for arbitrary dimension at half filling. The symmetric Anderson lattice is equivalent to the Kondo lattice, which makes it possible to study the behavior of an electron liquid in the Kondo lattice. We have shown that, due to hybridization (through an effective field due to localized electrons) of electrons with different spins and momenta \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{k} $$\end{document}k and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{k} +\overrightarrow{\pi }$$\end{document}k+π→, the gap in the electron spectrum opens at half filling. Such hybridization breaks the conservation of the total magnetic momentum of electrons, the spontaneous symmetry is broken. The state of electron liquid is characterized by a large Fermi surface. A gap in the spectrum is calculated depending on the magnitude of the on-site Coulomb repulsion and value of s–d hybridization for the chain, as well as for square and cubic lattices. Anomalous behavior of the heat capacity at low temperatures in the gapped state, which is realized in the symmetric Anderson lattice, was also found.

The processes of electron scattering with spin flip play a dominant role in the Kondo problem, they are determined by the exchange s-d interaction 1,2 . They should also be explicitly taken into account when considering the behavior of the Kondo lattice, since at half-filling the insulator state is also determined by this interaction. The one-particle modifications of the Kondo lattice (see for example 3 ), exact solvable many-particle models 4,5 , which do not taken into account the electron scattering with spin flip unable to describe the Kondo insulator state. Numerical calculations of the Anderson and Kondo lattices take into account clusters with a small number of particles, which significantly affects the adequacy of the results obtained [6][7][8] . This leads to the fact that the mechanism of the formation of the large Fermi surface and the gap in the electronic spectrum in the Kondo lattice at half filling is still unclear. In contrast to the single-impurity Anderson and Kondo models in their lattice version, the processes of electron scattering with spin flip cannot be calculated exactly.
In the weak s-d hybridization limit the Hamiltonian of the Anderson lattice reduces to that of the Kondo lattice. It should also be taken into account, that the local density of localized electrons must be equal to unity. In the weak s-d hybridization limit, the local density of d-electrons in the symmetric Anderson lattice is equal to unity. d-electrons determine the local spin-1 2 at the sites of the lattice, only in this case we can talk about the spin-1 2 Kondo lattice. Thus the symmetric Anderson model is similar to the Kondo lattice. In the strong coupling limit of the Hubbard model, when U ≫ t ( U and t are on-site repulsion and hopping integral), it is shown that the Hubbard model and the Kondo lattice model become identical 9 . It gives possibility to use the formalism, proposed for calculation of the Mott transition in the Hubbard model 10 , for solution of the Kondo lattice problem.
In the paper, we consider the solution of the symmetric Anderson model, using a mean field approach. An effective -field connects the states of d-electrons with the different spins and momenta k , k + − → π . Due to the s-d hybridization, the states of s-electrons with different spins and momenta k , k + − → π also hybridize, the gap opens in the electron spectrum at half filling. The value of the gap is determined by the magnitude of the -field (which in turn depends on the on-site repulsion) and the value of the s-d hybridization. We shall show that breaking spontaneous symmetry in the Anderson lattice Hamiltonian makes it possible to take into account the scattering of conduction electrons on d-electrons with a spin flip. They determine the ground state of symmetric Anderson and Kondo lattices, the value of the gap in the electron spectrum, cause formation of an insulator state at half-filling. The electron spectrum is symmetric about the zero energy, it corresponds to the symmetric Anderson model for a nontrivial solution of the -field. The gap in the electron spectrum has a many-particle nature.

Model
The Hamiltonian of the Anderson lattice is the sum of two terms, the first of which is determined by energy of the bands of s-and d-electrons and hybridization between them, the second takes into account the on-site repulsion of d-electrons H = H 0 + H int where c † j,σ , c j,σ and d † j,σ , d j,σ ( σ =↑, ↓) are the fermion operators determined on a lattice site j, U is the value of the on-site Hubbard interaction determined by the density operator n j,σ = d † j,σ d j,σ , the band width of c-fermions is determined by the hopping integral equal to one, the energy of flat band of d-fermions equal to ε g , v defines the hybridizations of s-and d-electrons, N is the total number of atoms.
First of all, we will focus on the consideration of the model for single Anderson impurity, exact solution of which has been obtained by Wiegmann 1 . At ε g = − U 2 and U ≫ Ŵ ( Ŵ = v 2 ) trickly one electron is localized on the impurity, so we can talk about an impurity with spin-1 2 . This case corresponds to the symmetric Anderson model. According to 2 the behavior of impurity in the symmetric Anderson model is equivalent to the spin-1 2 Kondo impurity. In this context, in a non-magnetic (paramagnetic) state the behavior of an electron liquid in the symmetric Anderson lattice is similar to that in the spin-1 2 Kondo lattice. Should be notes also, when the hybridization is small the Hamiltonian (1) can be mapped into the Kondo lattice model with the exchange imtegral J ≃ − 2ŴU ε g (ε g +U) .
The ground-state -field breaks the spontaneous symmetry 11,12 , the magnetic moments of conduction and d-electrons are not conserved, the total number of electrons is conserved 10 . Only for q = − → π , the electron spectrum corresponds to the symmetric Anderson lattice. Thus, we are talking about the lattice with the doubled period, as in the case of the Peierls phase-transition. In this case the electron spectrum is symmetric with respect to zero and has the following form for the chain (see in Fig 1a) and square lattice (see in Fig 2a). Four branches of the spectrum ±E γ (k) ( γ = 1, 2 ) are determined by the following expression (see section Methods)) (1) For v = 0 , = 0 and an arbitrary dimension of the lattice the spectrum is gapped at half-filling. The gap is equal to zero at v = 0 or = 0 (see in Fig 3), v = = 0 corresponds to the atomic limit of the symmetric Anderson lattice with a bare on-site repulsion.
The dependence of the gap on v and has quite universal in nature, its value does not depend on the dimension of the model (see in Fig 3). As shown below, the value of is the solution of Eq (3), at given U its value is determined by v, depends on dimension of the lattice. In the chain the low energy spectrum is determined by two special points (see in Fig 1b). In the square lattice, they form a square in the k-plane (see in Fig 2b). The k-surface, formed by the points on the k-plane of the low-energy spectrum, is constructed as follows: each projection of the wave vector corresponds to two values of its other projections. The γ-branches of the spectrum are hybridized, the electron density is defined as γ =1,2 dk = n s + 1 , here n s is the density of s-electrons ( n s = 1 at half-filling). Thus, we can talk about a large Fermi surface in the Kondo lattice also, the behavior of which is similar to the symmetric Anderson lattice.
In the ground state the equation for follows from the form of the action (6) and the spectrum (2) where the value of the on-site repulsion corresponds to a condition U = −2ε g . The energy of the quasi-particle excitations depend on ε g (2), thus (3) is self-consistent equation. As a result, the state of the symmetric Anderson lattice can be realized for = 0.    Fig. 4, corresponds to stable insulator state at U = 2 , the gapped state at U = 1 , is unstable. Calculating the actions δS = S( ) β − S(0) β at T = 0K, which correspond to the states with the gap (for = 0 ) and gapless with a maximum magnetic moment (for = 0 ), we determine the stability of the insulator state at U = 2 . The numerical calculations are shown in Fig 5 for different dimensions of the lattice. The region of stability of the insulator state is highlighted in color, the solution for obtained for U = 2 (see Fig. 4a) is in the colored region. Note, that δS = 0 at = 0 and δS ∼ 2 at → ∞ , thus intermediate values of 0 < < 0.1 correspond to the stable insulator state (see in Fig 5).
Taking into account the spectrum of quasi particle excitations (2) in the action (6), we determine the low temperature behavior of the specific heat (its electronic part) for 3D system (cubic lattice). In the insulator state, the value of C exponentially depends on the value of the gap. Numerical calculations of the temperature dependence C T are shown in Fig. 6. The calculations are made for U = 2 and v = 0.43 , = 0.0075 , = 0.0023 (see Fig. 6a) and v = 0.5 , = 0.07 , = 0.0279 (see Fig.6b). Note that for these parameters, the gaps in the electron spectrum differ by an order of magnitude. The experimentally measured heat capacity C T SmB 6 14 shows the anomalous behavior in the low temperature range 2K < T < 10K , namely at T = 4K there is a minimum. According to numerical calculations (see insets in Fig. 6), a local minimum of C T is realized at low temperatures. The temperature range and the minimum value depend on the value of the gap.

Methods
Canonical functional of the symmetric Anderson model. Let us introduce the operator χ † j = d † j,↑ d j,↓ and redefine the term H int (1) is the following form H int = −U j χ † j χ j 10 . The Hubbard-Stratonovich transformation converts the interacting problem into a non-interacting one in a stochastic external field (hereinafter we will define it as the -field). We define the interaction term, taking into account the action S 0 The canonical functional is determined as is the wave function, k is the wave vector of an electron. We expect that j does not depend on τ , since the translation invariance is conserved in an electron liquid state.
At the on-site hybridization between d-states of electrons with different spins and due to translation invariance of the model Hamiltonian, only the phases of j depends on j, a namely j = exp(iqj) 10 , where q is an unknown wave vector. The phase can also fluctuate at the lattice site; subsequent averaging over local fluctuations restores translational invariance 13 . The task is reduced to moving fermions in a static inhomogeneous -field that sets the form of H eff (k) The spectrum of non-interacting s-electrons is given by ε(k) = −2 D i=1 cos k i , here D is dimension of the model. We can integrate out fermions to obtain the following action S per an atom where ω n = T(2n + 1)π are the Matsubara frequencies, k , q are the momenta of electrons, four quasiparticle excitations E γ (k, q) ( γ = 1, ..., 4 ) determine the electron states in the -field. In the saddle point approximation the canonical functional Z will be dominated by the minimal action S (6), that satisfies the following equation ∂S/∂ = 0 . In the Kondo problem, the processes of electron scattering from spin flip dominate, in our case, the -field connects the states of s-electrons with opposite spins and different momenta k and k + q . Due to the on-site repulsion between electrons in the Hubbard model 10 , an effective field connects states of electrons with different spins and different momenta k , k + − → π . This state is stable at half-filling, a gap opens in the electron spectrum 10 . In this case, the same mechanism of coupling between s-electrons is realized indirectly through an intermediate subsystem of d-electrons. The gap in the electron spectrum will naturally be smaller.

Conclusion
Using mean field approximation we have considered the solution of the symmetric Anderson lattice at half-filling for different dimensions of the lattice. It is shown that an effective field, which binds the states of s-electrons with different spins and momenta, leads to the gap in the electron spectrum. The electron spectrum is mirror symmetric (with respect to zero energy), has the type of the Majorana spectrum. The states of s-and d-electrons are hybridized, therefore the Fermi surface is determined by the total density of electrons. In the symmetric Anderson lattice model the gaped state of electron liquid corresponds to breaking spontaneous symmetry. This made it possible to take into account the processes of spin flip scattering of conduction electrons on localized electrons. This is important because it is these processes that lead to the appearance of the Abrikosov -Suhl resonance in the Kondo problem and an insulator state in the Kondo lattice. The gapped state is formed at finite values of the s-d hybridization and on-site repulsion. Numerical calculations of the low-temperature behavior of the heat capacity are carried out in the gapped state. The results explained the presence of a local minimum of the heat capacity C T at low temperatures observed in SMB 6 . In the symmetric Anderson lattice, a local spin-1 2 is realized at lattice site, therefore its behavior is similar to the spin-1 2 Kondo lattice. The proposed approach allows us to describe the symmetric Anderson and Kondo lattices in one formalism. One can speak of a large Fermi surface and the gapped state of an electron liquid in the Kondo lattice at half filling.