Thema dieser Arbeit ist die Entwicklung und Kombination verschiedener numerischer Methoden, sowie deren Anwendung auf Probleme stark korrelierter Elektronensysteme. Solche Materialien zeigen viele interessante physikalische Eigenschaften, wie z.B. Supraleitung und magnetische Ordnung und spielen eine bedeutende Rolle in technischen Anwendungen. Es werden zwei verschiedene Modelle behandelt: das Hubbard-Modell und das Kondo-Gitter-Modell (KLM). In den letzten Jahrzehnten konnten bereits viele Erkenntnisse durch die numerische Lösung dieser Modelle gewonnen werden. Dennoch bleibt der physikalische Ursprung vieler Effekte verborgen. Grund dafür ist die Beschränkung aktueller Methoden auf bestimmte Parameterbereiche. Eine der stärksten Einschränkungen ist das Fehlen effizienter Algorithmen für tiefe Temperaturen.rnrnBasierend auf dem Blankenbecler-Scalapino-Sugar Quanten-Monte-Carlo (BSS-QMC) Algorithmus präsentieren wir eine numerisch exakte Methode, die das Hubbard-Modell und das KLM effizient bei sehr tiefen Temperaturen löst. Diese Methode wird auf den Mott-Übergang im zweidimensionalen Hubbard-Modell angewendet. Im Gegensatz zu früheren Studien können wir einen Mott-Übergang bei endlichen Temperaturen und endlichen Wechselwirkungen klar ausschließen.rnrnAuf der Basis dieses exakten BSS-QMC Algorithmus, haben wir einen Störstellenlöser für die dynamische Molekularfeld Theorie (DMFT) sowie ihre Cluster Erweiterungen (CDMFT) entwickelt. Die DMFT ist die vorherrschende Theorie stark korrelierter Systeme, bei denen übliche Bandstrukturrechnungen versagen. Eine Hauptlimitation ist dabei die Verfügbarkeit effizienter Störstellenlöser für das intrinsische Quantenproblem. Der in dieser Arbeit entwickelte Algorithmus hat das gleiche überlegene Skalierungsverhalten mit der inversen Temperatur wie BSS-QMC. Wir untersuchen den Mott-Übergang im Rahmen der DMFT und analysieren den Einfluss von systematischen Fehlern auf diesen Übergang.rnrnEin weiteres prominentes Thema ist die Vernachlässigung von nicht-lokalen Wechselwirkungen in der DMFT. Hierzu kombinieren wir direkte BSS-QMC Gitterrechnungen mit CDMFT für das halb gefüllte zweidimensionale anisotrope Hubbard Modell, das dotierte Hubbard Modell und das KLM. Die Ergebnisse für die verschiedenen Modelle unterscheiden sich stark: während nicht-lokale Korrelationen eine wichtige Rolle im zweidimensionalen (anisotropen) Modell spielen, ist in der paramagnetischen Phase die Impulsabhängigkeit der Selbstenergie für stark dotierte Systeme und für das KLM deutlich schwächer. Eine bemerkenswerte Erkenntnis ist, dass die Selbstenergie sich durch die nicht-wechselwirkende Dispersion parametrisieren lässt. Die spezielle Struktur der Selbstenergie im Impulsraum kann sehr nützlich für die Klassifizierung von elektronischen Korrelationseffekten sein und öffnet den Weg für die Entwicklung neuer Schemata über die Grenzen der DMFT hinaus.
Thema dieser Arbeit ist die Entwicklung und Kombination verschiedener numerischer Methoden, sowie deren Anwendung auf Probleme stark korrelierter Elektronensysteme. Solche Materialien zeigen viele interessante physikalische Eigenschaften, wie z.B. Supraleitung und magnetische Ordnung und spielen eine bedeutende Rolle in technischen Anwendungen. Es werden zwei verschiedene Modelle behandelt: das Hubbard-Modell und das Kondo-Gitter-Modell (KLM). In den letzten Jahrzehnten konnten bereits viele Erkenntnisse durch die numerische Lösung dieser Modelle gewonnen werden. Dennoch bleibt der physikalische Ursprung vieler Effekte verborgen. Grund dafür ist die Beschränkung aktueller Methoden auf bestimmte Parameterbereiche. Eine der stärksten Einschränkungen ist das Fehlen effizienter Algorithmen für tiefe Temperaturen.rnrnBasierend auf dem Blankenbecler-Scalapino-Sugar Quanten-Monte-Carlo (BSS-QMC) Algorithmus präsentieren wir eine numerisch exakte Methode, die das Hubbard-Modell und das KLM effizient bei sehr tiefen Temperaturen löst. Diese Methode wird auf den Mott-Übergang im zweidimensionalen Hubbard-Modell angewendet. Im Gegensatz zu früheren Studien können wir einen Mott-Übergang bei endlichen Temperaturen und endlichen Wechselwirkungen klar ausschließen.rnrnAuf der Basis dieses exakten BSS-QMC Algorithmus, haben wir einen Störstellenlöser für die dynamische Molekularfeld Theorie (DMFT) sowie ihre Cluster Erweiterungen (CDMFT) entwickelt. Die DMFT ist die vorherrschende Theorie stark korrelierter Systeme, bei denen übliche Bandstrukturrechnungen versagen. Eine Hauptlimitation ist dabei die Verfügbarkeit effizienter Störstellenlöser für das intrinsische Quantenproblem. Der in dieser Arbeit entwickelte Algorithmus hat das gleiche überlegene Skalierungsverhalten mit der inversen Temperatur wie BSS-QMC. Wir untersuchen den Mott-Übergang im Rahmen der DMFT und analysieren den Einfluss von systematischen Fehlern auf diesen Übergang.rnrnEin weiteres prominentes Thema ist die Vernachlässigung von nicht-lokalen Wechselwirkungen in der DMFT. Hierzu kombinieren wir direkte BSS-QMC Gitterrechnungen mit CDMFT für das halb gefüllte zweidimensionale anisotrope Hubbard Modell, das dotierte Hubbard Modell und das KLM. Die Ergebnisse für die verschiedenen Modelle unterscheiden sich stark: während nicht-lokale Korrelationen eine wichtige Rolle im zweidimensionalen (anisotropen) Modell spielen, ist in der paramagnetischen Phase die Impulsabhängigkeit der Selbstenergie für stark dotierte Systeme und für das KLM deutlich schwächer. Eine bemerkenswerte Erkenntnis ist, dass die Selbstenergie sich durch die nicht-wechselwirkende Dispersion parametrisieren lässt. Die spezielle Struktur der Selbstenergie im Impulsraum kann sehr nützlich für die Klassifizierung von elektronischen Korrelationseffekten sein und öffnet den Weg für die Entwicklung neuer Schemata über die Grenzen der DMFT hinaus.
Jüngste experimentelle Fortschritte erlauben die Beobachtung von elektronischen Relaxationsprozessen in Festkörpern in Echtzeit. Nach der optischen Anregung eines Systems relaxiert dieses in einen Gleichgewichtszustand. Die Dynamik des Systems während des Relaxationsprozesses ist von den vorherrschenden mikroskopischen Wechselwirkungen zwischen den Konstituenten des Systems abhängig. Die Wechselwirkung zwischen Elektronen und Gitterschwingungen - den Phononen - ist allgegenwärtig in Festkörpern, weshalb sie einer der wichtigsten Relaxationskanäle ist. Eine analytische Beschreibung der Relaxationsdynamik ist schwierig, und nur wenige exakte Lösungen existieren, selbst für die Beschreibung von Gleichgewichtsproblemen. Prinzipiell sind numerische Methoden in der Lage, eine Lösung sowohl für Gleichgewichts- als auch für Nichtgleichgewichtssituationen zu finden. Methoden wie die exakte Diagonalisierung oder die Dichtematrix-Renormierungsgruppenmethode, die mit Wellenfunktionen arbeiten, skalieren jedoch ungünstig mit der Dimension des lokalen Hilbertraumes. Dies ist besonders problematisch bei der Studie von Elektron-Phonon gekoppelten Systemen, bei der im allgemeinen besonders große lokale Basen benötigt werden. Für Grundzustandsprobleme gibt es zwei unabhängige Vorgehensweisen, um dem entgegenzuwirken: die Ein-Gitterplatz Dichtematrix Renormierungsgruppenmethode, die linear in der lokalen Dimension skaliert, sowie die Implementierung eines Schemas, das die Dimension der lokalen Basis verringert, indem es im Eigenraum der lokalen reduzierten Dichtematrix trunkiert was auf sogenannte optimale Moden führt. In dieser Dissertation kombinieren wir diese Strategien zu einer Methode, die linear in der Dimension der optimierten lokalen Basis anstelle der vollen lokalen Dimension skaliert. Wir demonstrieren das verbesserte Skalierungsverhalten der Methode anhand des Holstein-Polaron und des halb gefüllten Hubbard-Holstein Modells. Zudem beschreiben wir die Implementierung einer optimalen Phononbasis in den Time-Evolving Block Decimation Algorithmus, um das Skalierungsverhalten dieser Methode mit der lokalen Dimension zu verbessern. Für die Simulation des Polaron-Problems auf einer unendlich langen Kette hat sich die Krylovraum Zeitentwicklung in einem limitierten funktionalen Raum, als effizient herausgestellt. Wir adaptieren diese Methode für periodische Randbedingungen und endliche Systemgröße. Außerdem betrachten wir die Eigenschaften der lokalen reduzierten Dichtematrix als Funktion der Modellparameter und in Nichtgleichgewichtssituationen in drei Modellen: dem Bose-Bose Resonanz-Modell, dem Holstein Modell und dem Hubbard-Holstein Modell. Für fermionische und Spin-Modelle weiß man, dass die lokale von Neumann Entropie ein Indikator für Phasenübergänge ist. Für das Bose-Bose Resonanzmodell finden wir, dass die beiden Größen, lokale von Neumann Entropie und Eigenzustände der lokalen Dichtematrix, besondere Eigenschaften an Phasenübergängen zeigen. Wir finden außerdem, dass die Struktur der optimalen Moden sich als Funktion der Zeit in der Nichtgleichgewichtsdynamik ändert. Desweiteren beschäftigen wir uns mit der Relaxationsdynamik eines einzigen Elektrons gekoppelt an Holstein-Phononen. Im adiabatischen Fall finden wir einen Netto-Energiefluss von dem Elektron zu den Phononen und wir leiten einen analytischen Ausdruck für die Relaxationszeit des Elektrons im schwach wechselwirkenden, adiabatischen Fall her. Ein weiteres Hauptthema dieser Dissertation ist die Thermalisierung in geschlossenen Quanten- Vielteilchensystemen. Unser erstes Beispiel ist der Zerfall der Néel Ordnung unter Zeitentwicklung im eindimensionalen Fermi-Hubbard Hamiltonian. Wir finden, dass die Relaxation spin-verwandter Observablen von Spin-Anregungen kontrolliert werden und dass die Doppelbesetzung nicht im Sinne der Eigenzustands-Thermalisierungs-Hypothese thermalisiert, da das System nicht ergodisch ist. Als zweites Beispiel behandeln wir Vielteilchenlokalisierung in einem eindimensionalen System von attraktiv wechselwirkenden Fermionen, das im Grundzustandsphasendiagramm auch mit Unordnung eine delokalisierte Phase zeigt. Wir untersuchen diesen Phasenübergang mit modernen Hilfsmitteln, die auf den Verschränkungseigenschaften und der Existenz von Quasiteilchen in lokalisierten Phasen beruhen.
Die theoretische und experimentelle Untersuchung stark korrelierter Quantensysteme im und ausserhalb des Gleichgewichts ist in den letzten Jahren zu einem zentralen Forschungsschwerpunkt in der Physik der kondensierten Materie geworden. In dieser Habilitationsschrift werden anhand eigener Beiträge aktuelle Entwicklungen diskutiert und ein Ausblick auf zukünftige Entwicklungen des Feldes gegeben. Es wird dabei ein direkter Bezug von theoretischen Zugängen zu aktuellen Experimenten mit ultrakalten Gasen und mit stark korrelierten Materialien genommen. Insbesondere werden mikroskopische Modelle diskutiert, die zur Modellierung stark korrelierter Materialien (wie z.B. Quantenmagnete) herangezogen werden. Das zweite Thema der Schrift behandelt die Realisierung und Untersuchung von Quantenvielteilcheneffekten in aktuellen Experimenten mit ultrakalten Gasen durch sogenannte Quantensimulatoren. Ziel ist es, aufgrund theoretischer Überlegungen Vorhersagen zur Emulation stark korrelierter Systeme zu treffen und, basierend auf numerisch exakten Simulationen, Vorhersagen für die Experimente zu treffen. Da solche Experimente insbesondere auch ausserhalb des Gleichgewichts durchgeführt werden können, knüpft dies direkt an das dritte Thema der Schrift an, in dem die Physik von Quantenvielteilchensystemen ausserhalb des Gleichgewichts untersucht wird. Dazu gehören die Behandlung des Relaxationsverhaltens nach sogenannten Quantum Quenches und Fragestellungen zum Transport, die hier anhand von Situationen auf optischen Gittern diskutiert werden. Die drei Themenbereiche werden überspannt von der Fragestellung, wie man neuartige Quantenzustände realisieren kann, und wie man diese charakterisieren kann. Ausserhalb des Gleichgewichts erwächst die Fragestellung, ob durch das Relaxationsverhalten ungewöhnliche Zustände erzeugt werden können. Wie in dieser Arbeit diskutiert, ergeben sich in der Tat unerwartete Ergebnisse, wie die Bildung negativer Temperaturen, oder das Anwachsen der Teilchendichte verursacht durch repulsive Wechselwirkungen. Im Gleichgewicht kann das Vorhandensein starker Magnetfelder zu interessantem Verhalten führen, wie anhand verschiedener Systeme diskutiert wird. Hervorzuheben sind hier die Bildung unkonventioneller Mott-Isolatoren, von Phasen mit 'spin-nematischer' Ordnung, und unkonventionelle Quantenphasenübergänge. Eine aktuelle Entwicklung der aktuellen Forschung behandelt die Bildung sogenannter topologischer Phasen, die nicht durch einen lokalen Ordnungsparameter beschrieben werden und sich somit der Charakterisierung im Rahmen der Landau-Theorie von Phasen und Phasenübergängen entziehen. Aktuelle Fragestellungen sind die Realisierung solcher Phasen im Vorhandensein von Wechselwirkungen, auch langreichweitiger Natur, wie sie in Quantensimulatoren mit polaren Molekülen in Erscheinung treten. In dieser Schrift wird diskutiert, wie diese Fragestellungen durch einen numerischen Zugang - die Dichtematrixrenormierungsgruppe - behandelt werden können und wie diese sehr flexible Methode zu interessanten Einsichten in die facettenreiche Physik stark korrelierter Quantensysteme führt.
Numerische Simulationen sind ein unverzichtbares Werkzeug der modernen Physik geworden, insbesondere im Kontext stark korrelierter Quantensysteme, die häufig nicht analytisch behandelt werden können. Viele Physiker betrachten die numerische Physik mittlerweile als eigene Fachrichtung neben der experimentellen und der theoretischen Physik. In den Veröffentlichungen, aus denen diese kumulative Dissertation besteht, werden stark korrelierte Quantensysteme studiert. Dafür werden eine Vielzahl moderner numerischer Methoden verwendet, unter anderem exakte Diagonalisierung, Matrix-Produktzustände, Quanten-Monte-Carlo und neuronale Quantenzustände. In dieser Dissertation werden drei Klassen von stark korrelierten Quantensystemen betrachtet: 1) Vom Hubbard-Modell abgeleitete frustrierte Quantenmagnete, 2) Toric Codes und 3) Z2 Gittereichtheorien. Seit der Entdeckung der Hochtemperatursupraleitung in Cupraten in den 1980er-Jahren sind die Niedertemperaturphasen von Cupraten von zentralem Interesse für experimentelle und theoretische Studien. Modelle wie das Hubbard-Modell wurden eingeführt, um die effektive zweidimensionale Physik in CuO2-Schichten zu modellieren. Das Verschwinden der antiferromagnetischen Ordnung als Konsequenz von Störstellen und deren Zusammenhang mit supraleitenden Phasen ist bedeutsam, um den Bindungsmechanismus von Elektronen in supraleitenden Phasen zu entschlüsseln. Dieser Bindungsmechanismus ist bis heute Gegenstand aktueller Forschung. Mikroskopische Studien des t-J-Modells, eines Grenzfalls des Hubbard-Modells bei starken Interaktionen, legen Szenarien nahe, in denen Elektronen in zwei Partons fraktionalisieren: Ein Spinon (trägt den Spin) und ein Holon (trägt die Ladung). Mit Hilfe von Matrix-Produktzuständen wird in dieser Dissertation für einen paradigmatischen dotierten, frustrierten Quantenmagneten - das Majumdar-Ghosh-Modell - gezeigt, dass die Existenz eines gebundenen Parton-Zustands eng mit einer Feshbach-Resonanz verknüpft ist. Frustrierte Quantenmagnete werden zudem als Plattform vorgeschlagen, um Systeme weniger Teilchen zu untersuchen. Der Toric Code ist ein Grundpfeiler der modernen Physik der kondensierten Materie und ein grundlegendes Modell für die Felder Topologie und Quantum Error Correction. Der Grundzustand des Toric Codes ist eine Z2 Quantenspinflüssigkeit und ist eines der einfachsten Beispiele für eine topologische Phase. Es wurde ein C++-Softwarepaket (ParaToric) entwickelt, in dem ein Continuous-Time-Quanten-Monte-Carlo-Verfahrem implementiert wird, um den Toric Code in einem parallel angelegten Feld zu simulieren. Der Algorithmus von Wu, Deng und Prokof’ev wurde um weitere Zustandsaktualisierungen erweitert, die Simulationen in bislang unzulänglichen Parameterbereichen ermöglichen. ParaToric erlaubt bislang unerreichte Systemgrößen und bildet die Grundlage für Publikationen mit den bisher größten Systemen auf dreieckigen, quadratischen, hexagonalen und kubischen Gittern. Es werden Schnittstellen für C, C++ und Python bereitgestellt, die eine nahtlose Integration von ParaToric in Projekte Dritter ermöglichen. In dieser Dissertation werden neue geometrische Ordnungsparameter für Z2-Gittereichtheorien auf Basis von Perkolationstheorie eingeführt und getestet. Diese Ordnungsparameter können direkt aus experimentellen Messungen von Quantensimulatoren berechnet werden und sind konkurrierenden Ordnungsparametern in vielen Bereichen überlegen. Die Anwendung auf den Toric Code offenbart, dass selbst beim absoluten Nullpunkt Deconfined Phasen existieren, die nicht topologisch sind. Des Weiteren wird eine neue Netzwerkarchitektur für neuronale Quantenzustände eingeführt, die ungefähre Symmetrien in das Netzwerk kodiert. Diese Netzwerkarchitektur wird verwendet, um die Phasendiagramme nicht-stoquastischer Varianten des Toric Codes zu berechnen. Zu guter Letzt werden symmetriegeschützte topologische Phasen vorgeschlagen, die auch bei einer endlichen Temperatur existieren.
In stark korrelierte Elektronensystemen lassen sich eindrucksvolle Phänomene der Physik der kondensierten Materie beobachten: Hier treten neben dem berühmten Mott-Hubbard Übergang von einem Metall zu einem Isolator auch klassische Phasenübergänge wie Magnetismus, Ladungsordnung und Supraleitung auf. Darüber hinaus zeigen diese Systeme Quantenphasenübergänge, die Gegenstand intensiver modernster Forschung sind. Was diese Übergänge so bemerkenswert macht, ist die Tatsache, dass sie im Gegensatz zu klassischen Übergängen am absoluten Temperaturnullpunkt stattfinden. Dies macht Quantenphasenübergänge zum Einen sehr faszinierend, zum Anderen theoretisch sehr schwer zu beschreiben: Bis dato existiert keine umfassende Theorie der Quantenkritikalität. Das Ziel dieser Dissertation ist die theoretische Beschreibung von klassischen und Quantenphasenübergängen durch die Anwendung von hochaktuellen quantenfeldtheoretischen Methoden auf das fundamentalste Modell für elektronische Korrelationen, das Hubbard-Modell. Die hierfür verwendeten Methoden sind die dynamische Molekularfeldtheorie (DMFT), welche in der Lage ist, lokale Korrelationen zu beschreiben, aber räumliche Korrelationen vernachlässigt, und die dynamische Vertexapproximation (D-A), welche zusätzlich räumliche Fluktuationen auf beliebigen Längenskalen berücksichtigt. Zuerst werden Vorboten von Phasenübergängen näher untersucht. Beim Mott-Hubbard Übergang werden diese durch divergierende irreduzible Vertexfunktionen und bei klassischen Phasenübergängen durch starke Ladungs-, Spin- oder Paarfluktuationen repräsentiert. Danach wird der Einfluss von Dimensionalität und Nähe zu Phasenübergängen auf das Spektrum des Systems untersucht. Spezifische Relationen für die Selbstenergie in verschiedenen Dimensionen werden herausgearbeitet. Im nächsten Schritt wird das Schicksal des Mott-Hubbard Metall-Isolator-Übergangs in zwei Dimensionen bestimmt, wo die DMFT zu dessen Beschreibung nicht ausreicht, weil hier räumliche Korrelationen sehr stark werden. Schließlich wird das magnetische Phasendiagramm für das lochdotierte Hubbard-Modell in drei Dimensionen berechnet und analysiert. Die D-A ist hierfür prädestiniert, weil sie die simultane Behandlung von zeitlichen und räumlichen Korrelationen erlaubt, was insbesondere wichtig für die Beschreibung eines quantenkritischen Punktes bei starker Wechselwirkung ist. Die (quanten)kritischen Exponenten der magnetischen Suszeptibilität und Korrelationslänge werden mittels der D-A bestimmt. Diese stehen im Widerspruch zur konventionellen Hertz-Millis-Moriya Theorie.
Schwerpunkt der vorliegenden Arbeit ist die Bestimmung und Analyse von Einund Zwei-Teilchen-Funktionen stark korrelierter Elektronensysteme. Gerade die Analyse der Zwei-Teilchen-Funktionen legt die Einfuhrung neuer bosonischer Quasiteilchen, den -tonen, nahe. Dies geschieht durch Anwendung von Modellsystemen, konkret, dem Hubbard, erweiterten Hubbard und Pariser-Parr-Pople Modell, sowie auf Basis unterschiedlicher Methoden. Mit der Determinanten-Quanten-Monte-Carlo-Simulation wird die Ein-Teilchen Vertexfunktion, die Selbstenergie , bestimmt. Dabei wird eine Darstellung ausgearbeitet, welche den zweidimensionalen Impulsvektor k auf die eindimensionale Dispersionsrelation k abbildet, d.h. k k . Dies erm öglicht eine vollständige, intuitive Darstellung und somit Rückführung der zugrunde liegenden physikalischen Prozesse. Die Analyse von Zwei-Teilchen-Funktionen erfolgt durch den Parquetformalismus in der Parquetnäherung. Ferner, wird die Wechselwirkung zwischen den Elektronen erweitert, sodass z.B. das eindimensionale Benzolmolekül in der Pariser-Parr-Pople N äherung oder das zweidimensionale erweiterte Hubbard Modell simuliert werden kann. Bei immer stärker werdenden Ladungsdichtefluktuationen wird hierbei ein verstärkter Einfluss der Zwei-Teilchen-Vertexfunktion beobachtet,während die Ein-Teilchen-Vertexfunktion weniger stark korreliert. Die Zwei-Teilchen-Vertexfunktion wird mittels der optischen Leitfähigkeit weiter analysiert. Diese Vertexkorrekturen sind dabei von transversalen Teilchen-Loch Anregungen verursacht. Da diese mit einem Impulsvektor q = (, ..) verbunden sind, koppelt das einund ausgehende Lichtquantum an je zwei Elektron-Loch-Paare. Für stark korrelierte Elektronensysteme scheint dies von allgemeiner Gültigkeit zu sein und motiviert deshalb die Einführung der entsprechenden Polaritonen als -tonen.
In den letzten Jahren haben Quantenpunkte viel Aufmerksamkeit erhalten, da ihre magnetischen Eigenschaften für Spintronik Geräte benutzt oder ein in einem Quantenpunkt eingesperrter Spin als Qubit für einen Quantencomputer verwendet werden könnte. In dieser Arbeit untersuchen wir die Gleichgewichts- und Nichtgleichgewichtseigenschaften von verschiedenen Quantenstörstellensystemen (QIS). Als erstes wird das Kondo-Modell betrachtet, in dem das lokale magnetische Moment einer Störstelle über eine Heisenberg-Wechselwirkung an ein fermionisches Bad gekoppelt ist. Wir untersuchen, wie die Kondokorrelationen aufgebaut werden und wie sie sich durch das Leitungsband ausbreiten. Die Nichtgleichgewichtsergebnisse zeigen, dass die meisten Korrelationen sich innerhalb eines Lichtkegels ausbreiten, der bestimmt ist durch die Fermi-Geschwindigkeit. Interessanterweise sehen wir auch Korrelationen außerhalb des Lichtkegels. Mittels einer Störungstheorie in zweiter Ordnung in der Kopplung, können wir zeigen, dass diese Korrelationen mit den intrinsischen Korrelationen des Fermi-Sees zusammenhängen. Weiterhin stellen wir ein Experiment vor, in dem Metall-Molekül-Komplexe durch Reaktion von einem PTCDA Molekül, das auf einer Au(111) Oberfläche adsorbiert, mit einem einzelnen Au Atom gebildet werden. Dies führt zur Entstehung von Radikalen. Der beobachtete Kondoeffekt in den Rastertunnelmikroskopie (STM) Spektren ist ein eindeutiger Beweis für die Bildung eines lokalen Momentes. Die physikalisch korrekte Beschreibung der Bildung des lokalen Momentes und des Kondoeffekt ist jedoch mit der Dichtefunktionaltheorie (DFT) nicht möglich. Daher verwenden wir einen kombinierten DFT + Numerische Renormierungsgruppen (NRG) Ansatz, in dem wir die Ergebnisse der DFT als “ab initio” Eingangs-Parameter für unsere NRG verwenden. Wir zeigen, dass die Ergebnisse des Experiments und der NRG ohne Anpassungen der Modellparameter sehr gut übereinstimmen. Wenn zwei Au Atome auf benachbarten PTCDA Molekülen adsorbieren, weist das Spektrum in manchen Fällen eine Lücke um die Fermi-Energie, anstatt einer Kondoresonanz auf. Indem wir das System auf ein Zwei-Störstellen-Anderson-Modell (TIAM) abbilden, können wir zeigen, dass diese Lücke durch die starke nichtmagnetische chemische Wechselwirkung zwischen den beiden Komplexen entsteht. Mittels unserer NRG Ergebnisse liefern wir eine ausführliche Analyse des auftretenden Quantenphasenübergangs (QPT). Am Schluss betrachten wir die Gleichgewichts- und Nichtgleichgewichtseigenschaften des Zwei-Störstellen-Kondo-Modells (TIKM), welches einen QPT aufweist, wenn man die Energieabhängigkeit der Kopplungen an das Leitungsband vernachlässigt. Wir zeigen, dass es möglich ist, diesen QPT in dem voll energieabhängigen Modell mittels passenden Potentialstreutermen wiederherzustellen. Danach wird die Zeitdynamik nach einem Quench in der direkten Spin-Spin-Wechselwirkung zwischen den beiden Störstellenspins untersucht. Es wird gezeigt, dass eine ausreichend große antiferromagnetische Wechselwirkung die Thermalisierung der Spin-Spin-Korrelationsfunktion verhindert, da in diesem Fall die Störstellen vom Leitungsband entkoppeln.