U-perturbation treatment of heavy fermion systems

Heavy fermion systems are described by the periodic Anderson Model (PAM), i.e. a lattice of localized, highly correlatedf-electron states hybridized with the delocalized states of a conduction band. We treat the PAM within the second orderU perturbation theory around the non-magnetic Hartree-Fock solution (U on site Coulomb correlation between thef-electrons). This treatment has the advantage that Fermi liquid relations (Luttinger theorem) are automatically fulfilled. Thef-electron selfenergy and spectral function are calculated for different temperatures, and, for the symmetric PAM, we obtain single-particle peaks near toE _f andE _f +U and in addition many-particle (Kondo) resonance peaks near to the chemical potential (E _f baref-electron energy). The resonance peaks are strongly temperature dependent and vanish on a characteristic temperature scaleT _K. For the symmetric PAM and a constant on-site hybridization the Fermi energy falls into a hybridization gap. A second, smaller characteristic temperature scaleT _coh (coherence temperature), on which the hybridization gap vanishes, is observed within this approach. For the non-symmetric PAM (i.e.E _f andE _f +U not symmetric around the chemical potential) we obtain a similar behaviour, but the single-particle peaks are no longer at the correct positionsE _f andE _f +U. The proper behaviour for the symmetric PAM but less satisfactory behaviour for the non-symmetric PAM can be understood from the fact that only for the symmetric PAM the exactly solvable limit of a vanishing hybridization is reproduced within this approach.     

Coulomb repulsion versus Hubbard repulsion in a disordered chain

We study the difference between on site Hubbard and long range Coulomb repulsions for two interacting particles in a disordered chain. The system size L (in units of the lattice spacing) is of the order of the one particle localization length and the energies are taken near the band center. In the two cases, the limits of weak and strong interactions are characterized by uncorrelated energy levels and are separated by a crossover regime where the states are more extended and the spectra more rigid. U denoting the interaction strength and t the kinetic energy scale, the crossovers take place for interaction energy to kinetic energy ratios U/t and U/(2tL) of order one, for Hubbard and Coulomb repulsions respectively. While Hubbard repulsion can only yield weak critical chaos with intermediate spectral statistics, Coulomb repulsion can drive the two particle system to quantum chaos with Wigner-Dyson spectral statistics. The interaction matrix elements are studied to explain this difference. Received 21 March 2000 and Received in final form 5 February 2001     

CeO2的电子结构，光学和晶格动力学性质：第一性原理研究

基于考虑了Ce-4f电子间的库仑作用U和交换作用J的LDA+U方案,应用第一性原理计算系统研究了CeO₂的电子结构,光学和晶格动力学性质.电荷密度和电子局域函数的分布特征表明,CeO₂是属于共价键的绝缘体.介电常数、玻恩有效电荷张量和声子色散曲线的计算值和相应的实验结果符合得比较好. 关键词： 电子结构 光学性质 晶格动力学 第一性原理计算     

Density-matrix renormalization study of the Hubbard model [4]on a Bethe lattice

The half-filled Hubbard model on the Bethe lattice with coordination number z=3 is studied using the density-matrix renormalization group (DMRG) method. Ground-state properties such as the energy per site E, average local magnetization , its fluctuations and various spin correlation functions are determined as a function of the Coulomb interaction strength U/t. The local magnetic moments increase monotonically with increasing Coulomb repulsion U/t showing antiferromagnetic order between nearest neighbors []. At large U/t, is strongly reduced with respect to the saturation value 1/2 due to exchange fluctuations between nearest neighbors (NN) spins [ for ]. shows a maximum for U/t=2.4-2.9 that results from the interplay between the usual increase of with increasing U/t and the formation of important permanent moments at large U/t. While NN sites show antiferromagnetic spin correlations that increase with increasing Coulomb repulsion, the next NN sites are very weakly correlated over the whole range of U/t. The DMRG results are discussed and compared with tight-binding calculations for U=0, independent DMRG studies for the Heisenberg model and simple first-order perturbation estimates. Received 8 February 1999 and Received in final form 14 June 1999     

Magnetic impurities in gapless Fermi systems: perturbation theory

We consider a symmetric Anderson impurity model with a soft-gap hybridization vanishing at the Fermi level, with r>0. Three facets of the problem are examined. First the non-interacting limit, which despite its simplicity contains much physics relevant to the U>0case: it exhibits both strong coupling (SC) states (for r<1) and local moment states (for r>1), with characteristic signatures in both spectral properties and thermodynamic functions. Second, we establish general conditions upon the interaction self-energy for the occurence of a SC state for U>0. This leads to a pinning theorem, whereby the modified spectral function is pinned at the Fermi level for any U where a SC state obtains; it generalizes to arbitrary r the pinning condition upon familiar in the normal r=0 Anderson model. Finally, we consider explicitly spectral functions at the simplest level: second order perturbation theory in U, which we conclude is applicable for and r>1 but not for . Characteristic spectral features observed in numerical renormalization group calculations are thereby recovered, for both SC and LM phases; and for the SC state the modified spectral functions are found to contain a generalized Abrikosov-Suhl resonance exhibiting a characteristic low-energy Kondo scale with increasing interaction strength. Received 26 August 1999     

Phase transitions in the generalized spin-one-half Falicov–Kimball model in two dimensions

The extrapolation of small-cluster exact-diagonalisation calculations and the Monte-Carlo method is used to study the spin-one-half Falicov–Kimball model extended by the spin-dependent Coulomb interaction (J) between the localized f and itinerant d electrons as well as the on-site Coulomb interaction (U _ff ) between the localized f electrons. It is shown that in the symmetric case the ground-state phase diagram of the model has an extremely simple structure that consists of only two phases, and namely, the charge-density-wave (CDW) phase and the spin-density-wave (SDW) phase. The nonzero temperature studies showed that these phases persist also at finite temperatures. The critical temperature T _c for a transition from the low-temperature ordered phases to the high-temperature disordered phase is calculated numerically for various values of J and U _ff .     

Luttinger liquids with boundaries: Power-laws and energy scales

We present a study of the one-particle spectral properties for a variety of models of Luttinger liquids with open boundaries. We first consider the Tomonaga-Luttinger model using bosonization. For weak interactions the boundary exponent of the power-law suppression of the spectral weight close to the chemical potential is dominated by a term linear in the interaction. This motivates us to study the spectral properties also within the Hartree-Fock approximation. It already gives power-law behavior and qualitative agreement with the exact spectral function. For the lattice model of spinless fermions and the Hubbard model we present numerically exact results obtained using the density-matrix renormalization-group algorithm. We show that many aspects of the behavior of the spectral function close to the boundary can again be understood within the Hartree-Fock approximation. For the repulsive Hubbard model with interaction U the spectral weight is enhanced in a large energy range around the chemical potential. At smaller energies a power-law suppression, as predicted by bosonization, sets in. We present an analytical discussion of the crossover and show that for small U it occurs at energies exponentially (in -1/U) close to the chemical potential, i.e. that bosonization only holds on exponentially small energy scales. We show that such a crossover can also be found in other models. Received 8 February 2000 and Received in final form 25 April 2000     

Triplet <Emphasis Type="Italic">p</Emphasis>-wave superconductivity in the low-density extended hubbard model with Coulomb repulsion

We analyze superconducting instabilities in 3D and 2D extended Hubbard model with Coulomb repulsion between electrons on neighboring sites in the limit of low electron density (n _el → 0) on simple cubic (square) lattice. We show that in a realistic strong-coupling case U ≫ V ≫ W (U and V are the onsite and the intersite Coulomb repulsions, respectively, and W the bandwidth) the main SC instability corresponds to the p-wave pairing and in the leading order is correctly described by the equations obtained earlier in the absence of the intersite Coulomb interaction V = 0.     

Electronic structures and mechanical properties of uranium monocarbide from first-principles and calculations

The electronic structure, elastic constants, Poisson's ratio, and phonon dispersion curves of UC have been systematically investigated from the first-principles calculations by the projector-augmented-wave (PAW) method. In order to describe precisely the strong on-site Coulomb repulsion among the localized U 5f electrons, we adopt the local density approximation (LDA)+U and generalized gradient approximation (GGA)+U formalisms for the exchange correlation term. We systematically study how the electronic properties and elastic constants of UC are affected by the different choice of U as well as the exchange-correlation potential. We show that by choosing an appropriate Hubbard U parameter within the GGA+U approach, most of our calculated results are in good agreement with the experimental data. Therefore, the results obtained by the GGA+U with effective Hubbard parameter U chosen around 3 eV for UC are considered to be reasonable.     

Spectral densities of response functions for the O(3) symmetric Anderson and two channel Kondo models

The O(3) symmetric Anderson model is an example of a system which has a stable low energy marginal Fermi liquid fixed point for a certain choice of parameters. It is also exactly equivalent, in the large U limit, to a localized model which describes the spin degrees of freedom of the linear dispersion two channel Kondo model. We first use an argument based on conformal field theory to establish this precise equivalence with the two channel model. We then use the numerical renormalization group (NRG) approach to calculate both one-electron and two-electron response functions for a range of values of the interaction strength U. We compare the behaviours about the marginal Fermi liquid and Fermi liquid fixed points and interpret the results in terms of a renormalized Majorana fermion picture of the elementary excitations. In the marginal Fermi liquid case the spectral densities of all the Majorana fermion modes display a dependence on the lowest energy scale, and in addition the zero Majorana mode has a delta function contribution. The weight of this delta function is studied as a function of the interaction U and is found to decrease exponentially with U for large U. Using the equivalence with the two channel Kondo model in the large U limit, we deduce the dynamical spin susceptibility of the two channel Kondo model over the full frequency range. We use renormalized perturbation theory to interpret the results and to calculate the coefficient of the ln divergence found in the low frequency behaviour of the T=0 dynamic susceptibility. Received 29 January 1999     

Nonlocal Coulomb interaction in the two-dimensional spin-<Emphasis Type="Bold">1/2</Emphasis> Falicov–Kimball model

The two-dimensional (2D) extended Falicov–Kimball model has been studied to observe the role of nonlocal Coulomb interaction (U _nc) using an exact diagonalization technique. The f-state occupation (n ^f ), the f–d intersite correlation function (c _fd), the specific heat (C), entropy (S) and the specific heat coefficient (γ) have been examined. Nonlocal Coulomb interaction-induced discontinuous insulator-to-metal transition occurs at a critical f-level energy. More ordered state is obtained with the increase of U _nc. In the specific heat curves, two-peak structure as well as a single-peak structure appears. At low-temperature region, a sharp rise in the specific heat coefficient γ is observed. The peak value of γ shifts to the higher temperature region with U _nc.     

Destruction of the Fermi surface due to pseudogap fluctuations in strongly correlated systems

We generalize the dynamical-mean field theory (DMFT) by including into the DMFT equations dependence on the correlation length of the pseudogap fluctuations via the additional (momentum dependent) self-energy Σ_k. This self-energy describes nonlocal dynamical correlations induced by short-ranged collective SDW-like antiferromagnetic spin (or CDW-like charge) fluctuations. At high enough temperatures, these fluctuations can be viewed as a quenched Gaussian random field with finite correlation length. This generalized DMFT + Σ_k approach is used for the numerical solution of the weakly doped one-band Hubbard model with repulsive Coulomb interaction on a square lattice with nearest and next nearest neighbor hopping. The effective single impurity problem is solved by using a numerical renormalization group (NRG). Both types of strongly correlated metals, namely, (i) doped Mott insulator and (ii) the case of the bandwidth W ? U (U-value of local Coulomb interaction) are considered. By calculating profiles of the spectral densities for different parameters of the model, we demonstrate the qualitative picture of Fermi surface destruction and formation of Fermi arcs due to pseudogap fluctuations in qualitative agreement with the ARPES experiments. Blurring of the Fermi surface is enhanced with the growth of the Coulomb interaction.     

Superconductivity in weakly correlated electron systems

We study the influence of the short-ranged Hubbard correlation U between the conduction electrons on the Cooper pair formation in normal (s-wave) superconductors. The Coulomb correlation is considered within the standard second order perturbation theory, which becomes exact in the weak coupling limit but goes beyond the simple Hartree-Fock treatment by yielding a finite lifetime of the quasiparticles at finite temperature. An attractive pairing interaction V, which may be mediated by the standard electron-phonon mechanism, is considered between nearest neighbor sites. A critical value for the attractive interaction is required to obtain a superconducting state. For finite temperature a gapless superconductivity is obtained due to the finite lifetime of the quasiparticles, i.e. the Coulomb correlation has a pair-breaking influence. The energy gap and depend very sensitively on U, V and band filling n and develop a maximum away from half filling as function of n. The ratio varies with n, being higher than the BCS value near half filling and reaching the BCS value for lower n. Received 17 February 1999     

Effective electron Hamiltonian for a quasi one-dimensional system. The (TMTSF)2X and the (TMTTF)2X systems

In this paper we have introduced a variational approach to investigate the ground state of a model which includes both the Holstein electron-phonon interaction and the extended Hubbard electron-electron interaction. We have considered a variational state for the phonon subsystem which generalizes the previous used forms. This state allows to take into account the possibility of extended phonon mediated correlations. The effective electron Hamiltonian, which we have obtained, includes first and second neighbor electron-electron interaction terms. We have treated exactly, through a Lanczos method, this effective model in the one-dimensional case. We have applied our method to two Bechgaard salts and in these cases we have estimated the correlation parameters. We have shown that the introduction of electron-phonon interaction allows an estimate of the on site U and nearest-neighbor V Coulomb repulsion, which are in agreement with the experimental optical spectra of the above mentioned two compounds. Received: 30 October 1997 / Revised: 28 January 1998 / Accepted: 10 April 1998     

Thermoelectric power of high-T c superconductors using extended Hubbard model

The thermoelectric power (TEP) for a one dimensional lattice has been studied by using the extended Hubbard model in the limitU≠∞, whereU is the on-site Coulomb interaction. A new expression for TEP, derived in this study, has been found to successfully reproduce the experimental results of the insulating samples of different hole- and electron-doped high temperature superconducting systems.     

Stability of Ground States of 2d Strongly Asymmetric Correlated-Hopping Hubbard Model

The asymmetric correlated-hopping Hubbard model is analysed perturbatively for large values of the Coulomb interaction U. An effective Hamiltonian is obtained up to terms of the order U ^–3. For d=2 and in the limit of the strong asymmetry, the orderings of the ground states are found (confirming earlier nonrigorous results). Their thermal and quantum stability is proved. These results have been obtained by an application of the quantum Pirogov–Sinai theory in the variant developed by Datta, Fernandez, Fröhlich, and Rey-Bellet.     

Conductance and shot noise of inelastic tunneling through a quantum dot in the Kondo regime

We investigate the inelastic transport properties of a quantum dot connected to two leads, based on the combination of a recently developed nonperturbative technique and slave-boson methods involving the approximate mapping of the many-body electron–phonon coupling problem onto a multichannel scattering problem in the Kondo regime. The nonequilibrium Green's function method is adopted in calculations for the inelastic transport processes of electrons in the limit of large Coulomb interaction U→∞$U\to \infty$ under nonequilibrium conditions. The electron–phonon interactions, which are the main source of the inelasticity, are taken into account. For a single quantum dot, we find that the differential conductance and the shot noise exhibit new structures of peaks and dips which are absent in the case without electron–phonon interactions.     

欠掺杂高温超导体中的涡旋电荷结构相变

利用平均场t-t′-U-V-V_c模型,通过自洽求解Bogoliubov-de Gennes方程,研究了高温超导体中涡旋结构的相变.发现增大原位排斥势U,自旋密度波、电荷密度波以及d波序参量由棋盘结构转变为条纹结构.模型哈密顿量中引入合适强度的长程库仑势后,欠掺杂高温超导体样品中也可以出现二维或者棋盘结构,结果与文献报道的扫描隧道显微镜实验结果一致. 关键词： 高温超导 涡旋结构 长程库仑势     

Role of rotational symmetry in magnetism of multiband models

The role of the complete Coulomb interaction matrix in the determination of the magnetic susceptibility of the two-band Hubbard model has been investigated in the dynamical mean-field theory. The off-diagonal component of the U-matrix has been found to reduce the Néel temperature, which agrees with the available experimental data. The magnitude of the effect upon a change in the dimension of the U-matrix is discussed. The results have been compared with the thermodynamic characteristics of an isolated Hubbard atom.     

Kondo lattice state within the slave boson approach: spin-split masses and effective interaction among heavy quasiparticles

We provide a brief overview of our Hamiltonian approach to the the Kondo-lattice state at both finite Coulomb interaction U and nonzero applied magnetic field. In the mean-field approximation we introduce the spindependent effective masses, as well as calculate the effective Kondo (hybridization) temperature as a function of U. The effective interaction coming from the slave-boson exchange is shown to contain both the local (intraband) pairing and the Kondo (interband) interaction.