Angeregte Zustände
Temperature and bath size in exact diagonalization dynamical mean field theory
Dynamical mean field theory is one of the key methods to describe electronic properties of strongly correlated materials. If exact diagonalization is used as impurity solver, the size of the finite bath representing the infinite lattice is limited because of the rapid growth of the Hilbert space. In view of the increasing interest in the effect of multi-orbital and multi-site Coulomb correlations in transition metal oxides, high-Tc 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 bath size, (ii) 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 multi-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 to 10 meV two bath levels per correlated impurity orbital or site are usually adequate.
(A. Liebsch, H. Ishida) Topical Review: J. Phys. CM 24, 053201 (2012)
Correlated Dirac fermions on the honeycomb lattice studied within cluster dynamical mean field theory
The role of nonlocal Coulomb correlations in the honeycomb lattice is investigated within cluster dynamical mean field theory combined with finite-temperature exact diagonalization. The paramagnetic semimetal-to-insulator transition is found to be in excellent agreement with finite-size determinantal quantum Monte Carlo simulations and with cluster dynamical mean field calculations based on the continuous-time quantum Monte Carlo approach. As expected, the critical Coulomb energy is much lower than within a local or single-site formulation. Short-range correlations are shown to give rise to a pseudogap and concomitant non-Fermi-liquid behavior within a narrow range below the Mott transition.
(A. Liebsch) Phys. Rev. B 83, 035113 (2011)
High-energy pseudogap in iron pnictides induced via Hund coupling
The discovery of superconductivity in iron pnictides has stimulated intense discussions concerning the role of correlation effects, in particular, the importance of Hund exchange interactions. Optical data reveal a high-energy pseudogap not compatible with normal metal behavior. This pseudogap differs fundamentally from the low-energy gap in the antiferromagnetic spin-density wave phase. To analyze the effect of Hund coupling, the degenerate five-band Hubbard model is studied within dynamical mean field theory combined with exact diagonalization. A significant depletion of spectral weight is found above the Fermi level. It is shown that this pseudogap is associated with a collective mode in the self-energy caused by spin fluctuations. The pseudogap is remarkably stable over a wide range of Coulomb and exchange energies, but disappears for weak Hund coupling.
(A. Liebsch) Phys. Rev. B 84, 180505 (2011) (Rapid Communication)
Coulomb blockade and Kondo effect in the electronic structure of Hubbard molecules linked between metal leads
The electronic structure of small Hubbard molecules coupled between two non-interacting semi-infinite leads is studied in the low bias-voltage limit. To calculate the finite-temperature Green's function of the system, each semi-infinite lead is simulated by a small cluster, so that the problem is reduced to that of a finite-size system comprising the molecule and clusters on both sides. The Hamiltonian parameters of the cluster is chosen such that its embedding potential coincides with those of the semi-infinite leads on Matsubara points. Exact diagonalization based on the Arnoldi approach is used to evaluate the effect of Coulomb correlations on the electronic properties of the molecule at finite temperature. Depending on important Hamiltonian parameters, such as Coulomb repulsion and one-electron hopping within the molecule, and hybridization between molecule and leads, the self-energy of the molecule is shown to exhibit standard Fermi-liquid behavior or non-Fermi-liquid-like deviations giving rise to finite electronic scattering rates. The present method can also describe the formation of Kondo resonances inside correlation-induced pseudogaps, as far as the associated Kondo temperature is comparable with the temperature range studied by the exact diagonalization, whose lower boundary can be reduced by increasing the cluster size. These results demonstrate how the system can be tuned between the ballistic transport regime, the Coulomb blockade regime, and the Kondo regime.
(H. Ishida, A. Liebsch) Rhys. Rev. B 86, xxx (2012)
Spin excitations in nanostructures
Besides ground state properties, it is possible nowadays to grasp information on the magnetic excitations of nanostructures (adatoms, clusters and thin films) with state of the art inelastic scanning tunneling spectroscopy or electron energy loss spectroscopy. The theoretical investigation of these excitations hinges on the ability to access quantum mechanical processes behind the dynamics of a system. We use time dependent density functional theory combined with the Korringa-Kohn-Rostoker Green function method to devise a real-space scheme that enables the description of the main characteristics of magnetic excitations, i.e. their resonance frequency, their lifetime and their behavior upon application of external perturbations.
More details: Funsilab
Doppel- und Ladungstransferanregungen
Im Bestreben einer ab initio Beschreibung der elektronischen Struktur komplexer Systeme erlangen Doppel- und Ladungstransferanregungen wachsende Aufmerksamkeit aufgrund ihrer möglichen technologischen Relevanz. Erstere spielen in ultra-schnellen Prozessen, die gerade experimentell zugänglich geworden sind, eine Rolle, während letztere wahrscheinlich essentiell für die Beschreibung der komplexen Prozesse der Photosynthese sind. Die Herausforderungen Doppel- und Ladungstransferanregungen im Rahmen einer Dichtefunktionaltheorie zu beschreiben sind eng verknüpft, da beide ein Funktional voraussetzen, das sowohl im Raum als auch in der Zeit nicht-lokal ist. Insbesondere die zeitliche Nichtlokalität, d.h. eine Abhängigkeit von der Frequenz, fehlt in den momentan zur Verfügung stehenden Funktionalen.
Wir wollen ein frequenzabhängiges Dichtefunktional entwickeln, das uns die Beschreibung von Doppel- und Ladungstransferanregungen ermöglicht. Als alternative Methode werden wir reduzierte Dichtematrix-Funktionaltheorie einsetzen, die unter Beweis gestellt hat, viele hartnäckige Probleme der Dichtefunktionaltheorie lösen zu können. Wir werden zunächst die zeitabhängigen Version von reduzierter Dichtematrix-Funktionaltheorie formal korrekt herleiten und darauf aufbauend ein Funktional der reduzierten Dichtematrix ableiten, das für die Beschreibung von Doppel- und Ladunstransferanregungen geeignet ist. Die Eigenschaften aller Funktionale werden aus exakten Rechnungen für ein- bzw. zweidimensionale
Modellsysteme hergeleitet, für die eine Lösung der Schrödingergleichung ohne zusätzliche Näherungen für eine kleine Teilchenzahl möglich ist.
Vielteilchen-Störungstheorie und die GW-Näherung
Während die Dichtefunktionaltheorie eine effiziente Behandlung des elektronischen Grundzustands ermöglicht, lassen sich mit der Vielteilchen-Störungstheorie auch angeregte Zustände im Prinzip exakt beschreiben. Diese Methode ist daher besonders für die quantitative Berechnung von Anregungsenergien und spektroskopischen Größen geeignet. Sie basiert auf Greenschen Funktionen und der Entwicklung von Korrelationsfunktionen mit Hilfe von Feynman-Diagrammen. Die wichtigste praktische Realisierung ist die GW-Näherung für die elektronische Selbstenergie, die die Beschreibung von Quasiteilchen-Bandstrukturen wesentlich verbessert, aber auch Exziton- und Magnonresonanzen in der Dichte- und Spinsuszeptibilität lassen sich auf diese Weise untersuchen. Unsere Implementierung im Programm SPEX basiert ebenso wie FLEUR auf der FLAPW-Methode. Der Beitrag der Rumpfelektronen und relativistische Effekte werden genauso berücksichtigt wie die Spinpolarisierung von magnetischen Materialien. Es ist auch möglich, die Selbstenergie um die Elektron-Magnon-Streuung zu erweitern. Die kombinierte GWT-Selbstenergie kann zu starken Lebensdauereffekten sowie Band-Anomalien führen, wie in der obigen Abbildung am Beispiel von Eisen zu sehen.