Mixtures of correlated bosons and fermions: Dynamical mean‐field theory for normal and condensed phases

We derive a dynamical mean‐field theory for mixtures of interacting bosons and fermions on a lattice (BF‐DMFT). The BF‐DMFT is a comprehensive, thermodynamically consistent framework for the theoretical investigation of Bose‐Fermi mixtures and is applicable for arbitrary values of the coupling parameters and temperatures. It becomes exact in the limit of high spatial dimensions d or coordination number Z of the lattice. In particular, the BF‐DMFT treats normal and condensed bosons on equal footing and thus includes the effects caused by their dynamic coupling. Using the BF‐DMFT we investigate two different interaction models of correlated lattice bosons and fermions, one where all particles are spinless (model I) and one where fermions carry a spin one‐half (model II). In model I the local, repulsive interaction between bosons and fermions can give rise to an attractive effective interaction between the bosons. In model II it can also lead to an attraction between the fermions.     

Electronic structure calculations using dynamical mean field theory

The calculation of the electronic properties of materials is an important task of solid-state theory, albeit particularly difficult if electronic correlations are strong, e.g., in transition metals, their oxides and in f-electron systems. The standard approach to material calculations, the density functional theory in its local density approximation (LDA), incorporates electronic correlations only very rudimentarily and fails if the correlations are strong. Encouraged by the success of dynamical mean field theory (DMFT) in dealing with strongly correlated model Hamiltonians, physicists from the bandstructure and the many-body communities have joined forces and developed a combined LDA + DMFT method recently. Depending on the strength of electronic correlations, this new approach yields a weakly correlated metal as in the LDA, a strongly correlated metal or a Mott insulator. This approach is widely regarded as a breakthrough for electronic structure calculations of strongly correlated materials. We review this LDA + DMFT method and also discuss alternative approaches to employ DMFT in electronic structure calculations, e.g., by replacing the LDA part with the so-called GW approximation. Different methods to solve the DMFT equations are introduced with a focus on those that are suitable for realistic calculations with many orbitals. An overview of the successful application of LDA + DMFT to a wide variety of materials, ranging from Pu and Ce, to Fe and Ni, to numerous transition metal oxides, is given.     

Consistent LDA’ + DMFT approach to the electronic structure of transition metal oxides: Charge transfer insulators and correlated metals

We discuss the recently proposed LDA’ + DMFT approach providing a consistent parameter-free treatment of the so-called double counting problem arising within the LDA + DMFT hybrid computational method for realistic strongly correlated materials. In this approach, the local exchange-correlation portion of the electron-electron interaction is excluded from self-consistent LDA calculations for strongly correlated electronic shells, e.g., d-states of transition metal compounds. Then, the corresponding double-counting term in the LDA’ + DMFT Hamiltonian is consistently set in the local Hartree (fully localized limit, FLL) form of the Hubbard model interaction term. We present the results of extensive LDA’ + DMFT calculations of densities of states, spectral densities, and optical conductivity for most typical representatives of two wide classes of strongly correlated systems in the paramagnetic phase: charge transfer insulators (MnO, CoO, and NiO) and strongly correlated metals (SrVO₃ and Sr₂RuO₄). It is shown that for NiO and CoO systems, the LDA’ + DMFT approach qualitatively improves the conventional LDA + DMFT results with the FLL type of double counting, where CoO and NiO were obtained to be metals. Our calculations also include transition-metal 4s-states located near the Fermi level, missed in previous LDA + DMFT studies of these monoxides. General agreement with optical and the X-ray experiments is obtained. For strongly correlated metals, the LDA’ + DMFT results agree well with the earlier LDA + DMFT calculations and existing experiments. However, in general, LDA’ + DMFT results give better quantitative agreement with experimental data for band gap sizes and oxygen-state positions compared to the conventional LDA + DMFT method.     

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.     

Computation of electronic structure and magnetic properties of strongly correlated materials with LDA+DMFT method

In this review, we describe general ideas of the LDA+DMFT method which merges dynamical mean-field theory (DMFT) and density functional theory (in particular the local density approximation (LDA)). Nowadays, the LDA+DMFT computational scheme is the most powerful numerical tool for studying physical properties of real materials and chemical compounds. It incorporates the advantage of DMFT to treat the full range of local dynamical Coulomb correlations and the ability of band methods to describe material-specific band dispersion caused by the lattice periodicity. We briefly discuss underlying physical ideas of LDA+DMFT and its mathematical implementation. Then different algorithms applied to solution of the DMFT impurity problem are briefly described. We then give examples of successful applications of the LDA+DMFT method to study spectral and magnetic properties of recently synthesized compounds like pnictide superconductors as well as classic charge-transfer systems NiO and MnO.     

Dynamical mean-field approach to materials with strong electronic correlations

We review recent results on the properties of materials with correlated electrons obtained within the LDA+DMFT approach, a combination of a conventional band structure approach based on the local density approximation (LDA) and the dynamical mean-field theory (DMFT). The application to four outstanding problems in this field is discussed: (i) we compute the full valence band structure of the charge-transfer insulator NiO by explicitly including the p-d hybridization, (ii) we explain the origin for the simultaneously occuring metal-insulator transition and collapse of the magnetic moment in MnO and Fe₂O₃, (iii) we describe a novel GGA+DMFT scheme in terms of plane-wave pseudopotentials which allows us to compute the orbital order and cooperative Jahn-Teller distortion in KCuF₃ and LaMnO₃, and (iv) we provide a general explanation for the appearance of kinks in the effective dispersion of correlated electrons in systems with a pronounced three-peak spectral function without having to resort to the coupling of electrons to bosonic excitations. These results provide a considerable progress in the fully microscopic investigations of correlated electron materials.     

Temperature and bath size in exact diagonalization dynamical mean field theory

Dynamical mean field theory (DMFT), combined with finite-temperature exact diagonalization, is one of the methods used to describe electronic properties of strongly correlated materials. Because of the rapid growth of the Hilbert space, the size of the finite bath used to represent the infinite lattice is severely limited. In view of the increasing interest in the effect of multi-orbital and multi-site Coulomb correlations in transition metal oxides, high-T(c) cuprates, iron-based pnictides, organic crystals, etc, it is appropriate to explore the range of temperatures and bath sizes in which exact diagonalization provides accurate results for various system properties. On the one hand, the bath must be large enough to achieve a sufficiently dense level spacing, so that useful spectral information can be derived, especially close to the Fermi level. On the other hand, for an adequate projection of the lattice Green's function onto a finite bath, the choice of the temperature is crucial. The role of these two key ingredients in exact diagonalization DMFT is discussed for a wide variety of systems in order to establish the domain of applicability of this approach. Three criteria are used to illustrate the accuracy of the results: (i) the convergence of the self-energy with the bath size, (ii) the quality of the discretization of the bath Green's function, and (iii) comparisons with complementary results obtained via continuous-time quantum Monte Carlo DMFT. The materials comprise a variety of three-orbital and five-orbital systems, as well as single-band Hubbard models for two-dimensional triangular, square and honeycomb lattices, where non-local Coulomb correlations are important. The main conclusion from these examples is that a larger number of correlated orbitals or sites requires a smaller number of bath levels. Down to temperatures of 5-10?meV (for typical bandwidths W?≈?2?eV) two bath levels per correlated impurity orbital or site are usually adequate.     

Implementation of LDA+DMFT with the pseudo-potential-plane-wave method

We propose an efficient implementation of combining dynamical mean field theory(DMFT) with electronic structural calculation based on the local density approximation(LDA).The pseudo-potential-plane-wave method is used in the LDA part,which enables it to be applied to large systems.The full loop self consistency of the charge density has been reached in our implementation,which allows us to compute the total energy related properties.The procedure of LDA+DMFT is introduced in detail with a complete flow chart.We have also applied our code to study the electronic structure of several typical strong correlated materials,including cerium,americium and NiO.Our results fit quite well with both the experimental data and previous studies.     

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.     

Quantum many‐body phenomena in coupled cavity arrays

The increasing level of experimental control over atomic and optical systems gained in recent years has paved the way for the exploration of new physical regimes in quantum optics and atomic physics, characterised by the appearance of quantum many‐body phenomena, originally encountered only in condensed‐matter physics, and the possibility of experimentally accessing them in a more controlled manner. In this review article we survey recent theoretical studies concerning the use of cavity quantum electrodynamics to create quantum many‐body systems. Based on recent experimental progress in the fabrication of arrays of interacting micro‐cavities and on their coupling to atomic‐like structures in several different physical architectures, we review proposals on the realisation of paradigmatic many‐body models in such systems, such as the Bose‐Hubbard and the anisotropic Heisenberg models. Such arrays of coupled cavities offer interesting properties as simulators of quantum many‐body physics, including the full addressability of individual sites and the accessibility of inhomogeneous models.     

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.     

Fast multi-orbital equation of motion impurity solver for dynamical mean field theory

We propose a fast multi-orbital impurity solver for dynamical mean field theory (DMFT). Our DMFT solver is based on the equations of motion (EOMs) for local Green's functions and is constructed by generalizing from the single-orbital case to the multi-orbital case with the inclusion of the inter-orbital hybridizations and applying a mean field approximation to the inter-orbital Coulomb interactions. The two-orbital Hubbard model is studied using this impurity solver within a large range of parameters. The Mott metal-insulator transition and the quasiparticle peak are well described. A comparison of the EOM method with the quantum Monte Carlo method is made for the two-orbital Hubbard model and good agreement is obtained. The developed method hence holds promise as a fast DMFT impurity solver in studies of strongly correlated systems.     

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.     

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.     

Identifying Particle Correlations in Quantum Hall Regime

We introduce a method that allows the disclosure of correlations between particle positions in an arbitrary many‐body system. The method is based on a well‐known simulated annealing algorithm and the proposed artificial distribution technique. Additionally, we investigate correlations in quantum Hall liquids (we consider many‐body wave functions that have been recently determined via the cyclotron subgroup model) and present three‐dimensional plots of configuration probability distributions that have been established from numerical simulations. We demonstrate that the preferred simultaneous positions of particles (configurations of positions, which correspond to large values of a system's probability distribution, ) tend to form complicated geometric structures, which are equivalent to classical Wigner crystals only for Laughlin states. Furthermore, we claim that quantum Hall liquids attributed to non‐Laughlin fillings are correlated on subdomains rather than on a whole particle domain (due to a quantizing magnetic field, which modifies the topology of a system's dynamics). Finally, we characterize Hall‐like internal orders in terms of statistical correlations (one‐dimensional unitary representations of cyclotron subgroups). Our conclusions concerning the stability of many‐body states agree with transport measurements and various numerical studies.     

Pseudogap behavior in Bi<Subscript>2</Subscript>Ca<Subscript>2</Subscript>SrCu<Subscript>2</Subscript>O<Subscript>8</Subscript>: Results of the generalized dynamical mean-field approach

Pseudogap phenomena are observed for the normal underdoped phase of different high-T _c cuprates. Among others, the Bi₂Sr₂CaCu₂O_{8 − δ} (Bi2212) compound is one of the most studied experimentally. To describe the pseudogap regime in Bi2212, we use a novel generalized ab initio LDA + DMFT + Σ_k hybrid scheme. This scheme is based on the strategy of one of the most powerful computational tools for real correlated materials: the local density approximation (LDA) + dynamical mean-field theory (DMFT). Conventional LDA + DMFT equations are here supplied with an additional (momentum-dependent) self-energy Σ_k in the spirit of our recently proposed DMFT + Σ_k approach taking into account pseudogap fluctuations. In the present model, Σ_k describes nonlocal correlations induced by short-range collective Heisenberg-like antiferromagnetic spin fluctuations. The effective single-impurity problem of the DMFT is solved by the numerical renormalization group (NRG) method. Material-specific model parameters for the effective x ² − y ² orbital of Cu-3d shell of the Bi2212 compound, e.g., the values of intra-and interlayer hopping integrals between different Cu sites, the local Coulomb interaction U, and the pseudogap potential Δ were obtained within the LDA and LDA + DMFT schemes. Here, we report on the theoretical LDA + DMFT + Σ_k quasiparticle band dispersion and damping, Fermi surface renormalization, momentum anisotropy of (quasi)static scattering, densities of states, spectral densities, and angular-resolved photoemission (ARPES) spectra, taking into account pseudogap and bilayer splitting effects for normal (slightly) underdoped Bi2212 (δ = 0.15). We show that LDA + DMFT + Σ_k successfully describes strong (pseudogap) scattering close to Brillouin zone boundaries. Our calculated LDA + DMFT + Σ_k Fermi surfaces and ARPES spectra in the presence of pseudogap fluctuations are almost insensitive to the bilayer splitting strength. However, our LDA-calculated value of bilayer splitting is rather small to describe the experimentally observed peak-dip-hump structure. The results obtained are in good semiquantitative agreement with various recent ARPES experiments. The article was submitted by the authors in English.     

Electron correlations and many-body techniques in magnetism

Recent progress in the theory of magnetism and electron correlations is reviewed to clarify the theories developed in the last decade and their mutual relations. A historical development of the theory of magnetism is outlined, and the dynamical coherent potential approximation (CPA) which completely takes account of the dynamical spin and charge fluctuations within the single-site approximation is introduced. Both the dynamical effects on various magnetic properties and the many-body band structure are shown to be explained on the same footing. It is shown that the dynamical CPA is equivalent to the other single-site theories of strongly correlated electrons: the many-body CPA, the dynamical mean-field theory (DMFT), and the projection operator method CPA (PM-CPA). These theories are elucidated with use of a common concept of effective medium or coherent potential. The effects of orbital degeneracy and the realistic calculation scheme are discussed with an emphasis on Hund’s rule coupling. Non-local theories of magnetism and electron correlations which go beyond the single-site approximation are presented. They include the molecular dynamics approach to the magnetic short range order, the dynamical cluster methods as a direct extension of the DMFT, and the self-consistent projection operator approach as an extension of the PM-CPA with use of the incremental cluster expansion. The current problems of their approaches and their future perspective are discussed.     

Semi‐Classical Model of Strongly Correlated Coulomb Systems in Weak Magnetic Field

The integer and fractional quantum Hall effects are two remarkable macroscopic quantum phenomena occurring in two‐dimensional strongly correlated electronic systems at high magnetic fields and low temperatures. Quantization of Hall resistivity in the very high magnetic field regime at partial filling of the lowest Landau level indicates the stabilization of an electronic liquid quantum Hall phase of matter. Other interesting phases that differ from the quantum Hall phases take prominence in weaker magnetic fields when many more Landau levels are filled. These states manifest anisotropic magneto‐transport properties and, under certain conditions, appear to mimic charge density waves and/or liquid crystalline phases. One way to understand such a behavior has been in terms of effective interaction potentials confined to the highest Landau level partially filled with electrons. In this work we show that, for weak magnetic fields, such a quantum treatment of these strongly correlated Coulomb systems resembles a semi‐classical model of rotating electrons in which the time‐averaged interaction potential can be expressed solely in terms of guiding center coordinates. We discuss how the features of this semi‐classical effective potential may affect the stability of various strongly correlated electronic phases in the weak magnetic field regime (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)