Mott-hubbard metal-insulator transition in paramagnetic V2O3: an LDA+DMFT(QMC) study

The electronic properties of paramagnetic V2O3 are investigated by the computational scheme LDA+DMFT(QMC). This approach merges the local density approximation (LDA) with dynamical mean-field theory (DMFT) and uses quantum Monte Carlo simulations (QMC) to solve the effective Anderson impurity model of DMFT. Starting with the crystal structure of metallic V2O3 and insulating (V0.962Cr0.038)2O3 we find a Mott-Hubbard transition at a Coulomb interaction U approximately 5 eV. The calculated spectrum is in very good agreement with experiment. Furthermore, the orbital occupation and the spin state S = 1 determined by us agree with recent polarization dependent x-ray-absorption experiments.     

Nonlocal Coulomb interactions and metal-insulator transition in Ti2O3: a cluster LDA+DMFT approach

We present an ab initio quantum theory of the metal-insulator transition in Ti2O3. The recently developed cluster LDA+DMFT scheme is applied to describe the many-body features of this compound. The conventional single site DMFT cannot reproduce a low temperature insulating phase for any reasonable values of the Coulomb interaction. We show that the nonlocal Coulomb interactions and the strong chemical bonding within the Ti-Ti pair is the origin of the small gap insulating ground state of Ti2O3.     

Calculation of photoemission spectra of the doped Mott insulator using LDA+DMFT(QMC)

The spectral properties of La_1–xSr_xTiO₃, a doped Mott insulator with strong Coulomb correlations, are calculated with the ab initio computational scheme LDA+DMFT(QMC). It starts from the non-interacting electronic band structure as calculated by the local density approximation (LDA), and introduces the missing correlations by the dynamical mean-field theory (DMFT), using numerically exact quantum Monte-Carlo (QMC) techniques to solve the resulting self-consistent multi-band single-impurity problem. The results of the LDA+DMFT(QMC) approach for the photoemission spectra of La_1–xSr_xTiO₃ are in good agreement with experiment and represent a considerable qualitative and quantitative improvement on standard LDA calculations. Received 20 May 2000 and Received in final form 27 July 2000     

Consistent LDA’ + DMFT—an unambiguous way to avoid double counting problem: NiO test

We present a consistent way of treating a double counting problem unavoidably arising within the LDA + DMFT combined approach to realistic calculations of electronic structure of strongly correlated systems. The main obstacle here is the absence of systematic (e.g., diagrammatic) way to express LDA (local density approximation) contribution to exchange correlation energy appearing in the density functional theory. It is not clear then, which part of interaction entering DMFT (dynamical mean-field theory) is already taken into account through LDA calculations. Because of that, up to now there is no accepted unique expression for the double counting correction in LDA + DMFT. To avoid this problem we propose here the consistent LDA’ + DMFT approach, where LDA exchange correlation contribution is explicitly excluded for correlated states (bands) during self-consistent band structure calculations. What is left out of Coulomb interaction for those strongly correlated states (bands) is its non-local part, which is not included in DMFT, and the local Hartreelike contribution. Then the double counting correction is uniquely reduced to the local Hartree contribution. Correlations for strongly correlated states are then directly accounted for via the standard DMFT. We further test the consistent LDA’ + DMFT scheme and compare it with conventional LDA + DMFT calculating the electronic structure of NiO. Opposite to the conventional LDA + DMFT our consistent LDA’ + DMFT approach unambiguously produces the insulating band structure in agreement with experiments.     

Phase diagram of electron-hole systems: Interplay between exciton Mott transition and quantum pair condensation

Quasi-thermal-equilibrium states of electron-hole (e-h) systems in photoexcited insulators are studied from a theoretical viewpoint, stressing the exciton Bose-Einstein condensation (BEC), the e-h BCS-type pair-condensed state, and the exciton Mott transition between an insulating exciton/biexciton gas phase and a metallic e-h plasma phase. We determine the quasi-equilibrium phase diagram of the e-h system at zero and finite temperatures with applying the dynamical mean-field theory (DMFT) to the e-h Hubbard model with both repulsive and attractive on-site interactions. Effects of inter-site interactions on the exciton Mott transition are also clarified with applying the extended DMFT to the extended e-h Hubbard model.     

Quasiparticle properties of strongly correlated electron systems with itinerant metamagnetic behavior

A brief account of the zero temperature magnetic response of a system of strongly correlated electrons in strong magnetic field is given in terms of its quasiparticle properties. The scenario is based on the paramagnetic phase of the half-filled Hubbard model, and the calculations are carried out with the dynamical mean field theory (DMFT) together with the numerical renormalization group (NRG). As well known, in a certain parameter regime one finds a magnetic susceptibility which increases with the field strength. Here, we analyze this metamagnetic response based on Fermi liquid parameters, which can be calculated within the DMFT-NRG procedure. The results indicate that the metamagnetic response can be driven by field-induced effective mass enhancement. However, also the contribution due to quasiparticle interactions can play a significant role. We put our results in context with experimental studies of itinerant metamagnetic materials.     

Normal phase and superconducting instability in the attractive Hubbard model: a DMFT(NRG) study

We study the normal (nonsuperconducting) phase of the attractive Hubbard model within the dynamical mean field theory (DMFT) using the numerical renormalization group (NRG) as an impurity solver. A wide range of attractive potentials U is considered, from the weak-coupling limit, where superconducting instability is well described by the BCS approximation, to the strong-coupling region, where the superconducting transition is described by Bose condensation of compact Cooper pairs, which are formed at temperatures much exceeding the superconducting transition temperature. We calculate the density of states, the spectral density, and the optical conductivity in the normal phase for this wide range of U, including the disorder effects. We also present the results on superconducting instability of the normal state dependence on the attraction strength U and the degree of disorder. The disorder influence on the critical temperature T _c is rather weak, suggesting in fact the validity of Anderson’s theorem, with the account of the general widening of the conduction band due to disorder.     

Mott-Hubbard transition and Anderson localization: A generalized dynamical mean-field theory approach

The DOS, the dynamic (optical) conductivity, and the phase diagram of a strongly correlated and strongly disordered paramagnetic Anderson-Hubbard model are analyzed within the generalized dynamical mean field theory (DMFT + Σ approximation). Strong correlations are taken into account by the DMFT, and disorder is taken into account via an appropriate generalization of the self-consistent theory of localization. The DMFT effective single-impurity problem is solved by a numerical renormalization group (NRG); we consider the three-dimensional system with a semielliptic DOS. The correlated metal, Mott insulator, and correlated Anderson insulator phases are identified via the evolution of the DOS and dynamic conductivity, demonstrating both the Mott-Hubbard and Anderson metal-insulator transition and allowing the construction of the complete zero-temperature phase diagram of the Anderson-Hubbard model. Rather unusual is the possibility of a disorder-induced Mott insulator-to-metal transition. The text was submitted by the authors in English.     

Momentum Space Quantum Monte Carlo on Twisted Bilayer Graphene

We report an implementation of the momentum space quantum Monte Carlo(QMC) method on the interaction model for the twisted bilayer graphene(TBG). The long-range Coulomb repulsion is treated exactly with the flat bands, spin and valley degrees of freedom of electrons taking into account. We prove the absence of the minus sign problem for QMC simulation when either the two valleys or the two spin degrees of freedom are considered.By taking the realistic parameters of the twist angle and interlayer tunnelings into the simulation, we benchmark the QMC data with the exact band gap obtained at the chiral limit, to reveal the insulating ground states at the charge neutrality point(CNP). Then, with the exact Green's functions from QMC, we perform stochastic analytic continuation to obtain the first set of single-particle spectral function for the TBG model at CNP. Our momentum space QMC scheme therefore offers the controlled computation pathway for systematic investigation of the electronic states in realistic TBG model at various electron fillings.     

Renormalized quasiparticles in antiferromagnetic states of the Hubbard model

We analyze the properties of the quasiparticle excitations of metallic antiferromagnetic states in a strongly correlated electron system. The study is based on dynamical mean field theory (DMFT) for the infinite dimensional Hubbard model with antiferromagnetic symmetry breaking. Self-consistent solutions of the DMFT equations are calculated using the numerical renormalization group (NRG). The low energy behavior in these results is then analyzed in terms of renormalized quasiparticles. The parameters for these quasiparticles are calculated directly from the NRG derived self-energy, and also from the low energy fixed point of the effective impurity model. From these the quasiparticle weight and the effective mass are deduced. We show that the main low energy features of the k-resolved spectral density can be understood in terms of the quasiparticle picture. We also find that Luttinger's theorem is satisfied for the total electron number in the doped antiferromagnetic state.     

Gap-anisotropic model for the narrow-Gap kondo insulators

A theory is presented which accounts for the dynamical generation of a hybridization gap with nodes in the Kondo insulating materials CeNiSn and CeRhSb. We show that Hund's interactions acting on virtual 4f(2) configurations of the cerium ion can act to dynamically select the shape of the cerium ion by generating a Weiss field which couples to the shape of the ion. In low symmetry crystals where the external crystal fields are negligible, this process selects a nodal Kondo semimetal state as the lowest energy configuration.     

Janus-faced influence of Hund's rule coupling in strongly correlated materials

We show that in multiband metals the correlations are strongly affected by Hund's rule coupling, which depending on the filling promotes metallic, insulating or bad-metallic behavior. The quasiparticle coherence and the proximity to a Mott insulator are influenced distinctly and, away from single- and half-filling, in opposite ways. A strongly correlated bad metal far from a Mott phase is found there. We propose a concise classification of 3d and 4d transition-metal oxides within which the ubiquitous occurrence of strong correlations in Ru- and Cr-based oxides, as well as the recently measured high Néel temperatures in Tc-based perovskites are naturally explained.     

Spin order due to orbital fluctuations: cubic vanadates

We investigate the highly frustrated spin and orbital superexchange interactions in cubic vanadates. The fluctuations of t(2g) orbitals trigger a novel mechanism of ferromagnetic interactions between spins S = 1 of V3+ ions along one of the cubic directions which operates already in the absence of Hund's rule exchange J(H), and leads to the C-type antiferromagnetic phase in LaVO3. The Jahn-Teller effect can stabilize the orbital ordering and the G-type antiferromagnetic phase at low temperatures, but large entropy due to orbital fluctuations favors again the C phase at higher temperatures, as observed in YVO (3).     

Construction and solution of a Wannier-functions based Hamiltonian in the pseudopotential plane-wave framework for strongly correlated materials

Ab initio determination of model Hamiltonian parameters for strongly correlated materials is a key issue in applying many-particle theoretical tools to real narrow-band materials. We propose a selfcontained calculation scheme to construct, with an ab initio approach, and solve such a Hamiltonian. The scheme uses a Wannier-function-basis set, with the Coulomb interaction parameter U obtained specifically for theseWannier functions via constrained Density functional theory (DFT) calculations. The Hamiltonian is solved by Dynamical Mean-Field Theory (DMFT) with the effective impurity problem treated by the Quantum Monte Carlo (QMC) method. Our scheme is based on the pseudopotential plane-wave method, which makes it suitable for developments addressing the challenging problem of crystal structural relaxations and transformations due to correlation effects. We have applied our scheme to the “charge transfer insulator” material nickel oxide and demonstrate a good agreement with the experimental photoemission spectra.     

Quantum phase transition in the two-band hubbard model

The interaction between itinerant and Mott localized electronic states in strongly correlated materials is studied within dynamical mean field theory in combination with the numerical renormalization group method. A novel nonmagnetic zero temperature quantum phase transition is found in the bad-metallic orbital-selective Mott phase of the two-band Hubbard model, for values of the Hund's exchange which are relevant to typical transition metal oxides.     

Dimerization versus orbital-moment ordering in a Mott insulator YVO3

We use exact diagonalization combined with mean-field theory to investigate the phase diagram of the spin-orbital model for cubic vanadates. The spin-orbit coupling competes with Hund's exchange and triggers a novel phase, with the ordering of t(2g) orbital magnetic moments stabilized by the tilting of VO6 octahedra. It explains qualitatively spin canting and reduction of magnetization observed in YVO3. At finite temperature, an orbital instability in the C-type antiferromagnetic phase induces modulation of magnetic exchange constants even in the absence of lattice distortions. The calculated spin structure factor shows a magnon splitting at q-->=(0,0,pi / 2) due to the orbital dimerization.     

Double counting in LDA + DMFT—The example of NiO

An intrinsic issue of the LDA + DMFT approach is the so called double counting of interaction terms. How to choose the double-counting potential in a manner that is both physically sound and consistent is unknown. We have conducted an extensive study of the charge-transfer system NiO in the LDA + DMFT framework using quantum Monte Carlo and exact diagonalization as impurity solvers. By explicitly treating the double-counting correction as an adjustable parameter we systematically investigated the effects of different choices for the double counting on the spectral function. Different methods for fixing the double counting can drive the result from Mott insulating to almost metallic. We propose a reasonable scheme for the determination of double-counting corrections for insulating systems.     

Two-dimensional Anderson-Hubbard model in the DMFT + Σ approximation

The density of states, the dynamic (optical) conductivity, and the phase diagram of the paramagnetic two-dimensional Anderson-Hubbard model with strong correlations and disorder are analyzed within the generalized dynamical mean field theory (DMFT + Σ approximation). Strong correlations are accounted by the DMFT, while disorder is taken into account via the appropriate generalization of the self-consistent theory of localization. We consider the two-dimensional system with the rectangular “bare” density of states (DOS). The DMFT effective single-impurity problem is solved by numerical renormalization group (NRG). The “correlated metal,” Mott insulator, and correlated Anderson insulator phases are identified from the evolution of the density of states, optical conductivity, and localization length, demonstrating both Mott-Hubbard and Anderson metal-insulator transitions in two-dimensional systems of finite size, allowing us to construct the complete zero-temperature phase diagram of the paramagnetic Anderson-Hubbard model. The localization length in our approximation is practically independent of the strength of Hubbard correlations. But the divergence of the localization length in a finite-size two-dimensional system at small disorder signifies the existence of an effective Anderson transition.     

Pseudogap behavior in the emery model for the electron-doped superconductor Nd2 − x Ce x CuO4: Multiband LDA + DMFT + Σk approach

We propose a generalization of the LDA + DMFT + Σ _k approach to the multiband case, in which correlated and uncorrelated states are present in the model simultaneously. Using the multiband version of the LDA + DMFT + Σ _k approach, we calculate the density of states and spectral functions for the Emery model in a wide energy interval around the Fermi level. We also obtain the Fermi surfaces for the electron-doped high-temperature superconductor Nd_{2 ? x} Ce _x CuO₄ in the pseudogap phase. The self-energy part Σ _k introduced additionally to take into account pseudogap fluctuations describes the nonlocal interaction of correlated electrons with collective Heisenberg short-range spin fluctuations. To solve the effective impurity model, the numerical renorm-group (NRG) method is used for the DMFT equations. Good qualitative agreement of the Fermi surfaces calculated using the LDA + DMFT + Σ _k approach and experimental angle-resolved photoemission spectroscopic data is attained. The stability of the dielectric solution with charge transfer in the Emery model with correction for double counting is analyzed in the Appendix.