This thesis presents results covering several topics in correlated many fermion systems. A Monte Carlo technique (CT-INT) that has been implemented, used and extended by the author is discussed in great detail in chapter 3. The following chapter discusses how CT-INT can be used to calculate the two particle Green’s function and explains how exact frequency summations can be obtained. A benchmark against exact diagonalization is presented. The link to the dynamical cluster approximation is made in the end of chapter 4, where these techniques are of immense importance. In chapter 5 an extensive CT-INT study of a strongly correlated Josephson junction is shown. In particular, the signature of the first order quantum phase transition between a Kondo and a local moment regime in the Josephson current is discussed. The connection to an experimental system is made with great care by developing a parameter extraction strategy. As a final result, we show that it is possible to reproduce experimental data from a numerically exact CT-INT model-calculation. The last topic is a study of graphene edge magnetism. We introduce a general effective model for the edge states, incorporating a complicated interaction Hamiltonian and perform an exact diagonalization study for different parameter regimes. This yields a strong argument for the importance of forbidden umklapp processes and of the strongly momentum dependent interaction vertex for the formation of edge magnetism. Additional fragments concerning the use of a Legendre polynomial basis for the representation of the two particle Green’s function, the analytic continuation of the self energy for the Anderson Kane Mele Model, as well as the generation of test data with a given covariance matrix are documented in the appendix. A final appendix provides some very important matrix identities that are used for the discussion of technical details of CT-INT.
This thesis presents results covering several topics in correlated many fermion systems. A Monte Carlo technique (CT-INT) that has been implemented, used and extended by the author is discussed in great detail in chapter 3. The following chapter discusses how CT-INT can be used to calculate the two particle Green’s function and explains how exact frequency summations can be obtained. A benchmark against exact diagonalization is presented. The link to the dynamical cluster approximation is made in the end of chapter 4, where these techniques are of immense importance. In chapter 5 an extensive CT-INT study of a strongly correlated Josephson junction is shown. In particular, the signature of the first order quantum phase transition between a Kondo and a local moment regime in the Josephson current is discussed. The connection to an experimental system is made with great care by developing a parameter extraction strategy. As a final result, we show that it is possible to reproduce experimental data from a numerically exact CT-INT model-calculation. The last topic is a study of graphene edge magnetism. We introduce a general effective model for the edge states, incorporating a complicated interaction Hamiltonian and perform an exact diagonalization study for different parameter regimes. This yields a strong argument for the importance of forbidden umklapp processes and of the strongly momentum dependent interaction vertex for the formation of edge magnetism. Additional fragments concerning the use of a Legendre polynomial basis for the representation of the two particle Green’s function, the analytic continuation of the self energy for the Anderson Kane Mele Model, as well as the generation of test data with a given covariance matrix are documented in the appendix. A final appendix provides some very important matrix identities that are used for the discussion of technical details of CT-INT.
One of the most intriguing areas of physics is the study of strongly correlated many-body systems. In these kinds of systems, a large number of particles interacts via interactions that are strong enough to play a major role in determining the properties of the system, which is particularly interesting when multiple interaction terms compete and cannot be satisfied simultaneously. Examples of phases they realize range from simple magnetism to puzzling phenomena such as superconductivity or topological order. Of special interest in condensed matter are fermionic models which have an extra layer of complexity added by the possibly intricate sign structure of the wave function. \nHowever, exactly solvable models that exhibit the sought-after phenomena are scarce and Hamiltonians that are supposed to model experiments often elude analytical solutions precisely due to the strong interactions. Numerical methods are thus essential to gain insight into these systems and provide an important link between theoretical models and real world experiments. There are many different methods that in principle allow solving a problem numerically exact, most prominently exact diagonalization, tensor network based methods such as DMRG and (quantum) Monte Carlo. Exact diagonalization is certainly the most straightforward technique which allows measuring any observable of interest but its application to interacting models is severely limited by the exponential growth of the Hilbert space with the system size which only allows studying rather small systems. Tensor network methods excel in the study of gapped, one-dimensional systems and while they have shown some impressive results for two-dimensional systems, they need a lot of fine tuning and careful analysis to be considered reliable in this very important domain. Quantum Monte Carlo, on the other hand, is in principle neither limited by the size of the Hilbert space nor by the strength of the interactions and has thus been a tremendously important tool for studying condensed matter systems which is why it is the focus of this thesis. \nIt is necessary to continuously advance existing and develop new techniques to keep up with theoretical as well as with experimental progress. In this thesis, new methods to identify and characterize conventional and novel phases of matter within quantum Monte Carlo are developed and tested on archetypical models of condensed matter. \nIn the first part of this thesis, a method to study entanglement properties of strongly interacting fermionic models in quantum Monte Carlo is presented. Measurements of this kind are needed because some characteristics of a system remain hidden from conventional approaches based on the calculation of correlation functions. Prime examples are systems that realize so-called topological order entailing long-range entanglement that can be positively diagnosed using entanglement techniques. Realizing such measurements in a concrete algorithm is not straightforward for fermions. A large section of this part is thus devoted to developing a solution to the problem of implementing entanglement measurements and benchmarking the results against known data. \nThe second part is concerned with a truly novel and recent approach to the many-body problem, namely the symbiosis of machine learning and quantum Monte Carlo techniques. Out of the many different possibilities to combine the two, this part is concerne with exploring ways in which machine learning can help to effortlessly explore phase diagrams of hitherto unknown Hamiltonians. \nDespite the power and versatility of quantum Monte Carlo techniques, it does suffer from one significant weakness called the fermion sign problem. At its core related to the peculiar exchange statistics of fermions, it is in principle merely a technical problem of the Monte Carlo procedure but one that actually prohibits the use of Monte Carlo methods for a large number of interesting problems. The sign problem has been known for a long time but remains unsolved to this day. Thus, it is time to approach the problem from a different perspective with the hope of learning something new that can either help to better understand the problem itself or the properties of the models affected by it. Both of the two major research areas, entanglement measures and machine learning, are capable of providing such new perspectives and the sign problem will be revisited in more detail in the context of each of them.
This thesis divides very naturally into three chapters, reflecting the three separate areas I have worked in throughout my PhD studies. During my PhD, I published two papers, one relating to the work in the first chapter, and one to that of the second. At the time of submission, the project which the third chapter relates to was still ongoing, with plans for a future publication. The first chapter discusses work done on the real space renormalisation group (RG) for Ising and Potts models in 2 dimensions. Taking inspiration from Hasenbusch's work, a new framework for carrying out the RG is developed, its computational implementation discussed in some detail, as well as the results for a variety of systems and the implications of these results. The numerical scaling of the procedure and the consequences of this for future work are also covered. The results documented in chapter two are purely theoretical and presented in closed form. The research conducted is to do with properties of free interfaces. The second chapter is by and large critical of previous assumptions which have been made about the so-called wave vector dependent, in order to attempt to experimentally measure it and use these measurements to make further predictions. Using some simple toy models, many of these assumptions are shown to be false. In a sense, the goal of the research presented in the chapter two, is not to motivate further research, but to dissuade research in a direction we consider to be misguided, due to the faulty assumptions it is based on. The third chapter covers a small subset of a project concerned with studying correlations within the so-called Abelian Manna Model. The majority of the project involves (computational) Monte Carlo simulations of the dynamics of such systems, but as these results are not ready to present at the time of writing, the chapter is mainly concerned with some analytical results which were derived in order to validate our models for small systems, explain certain quirky phenomena arising from our simulations, and help quantify errors. Finally, there is an appendix which expands upon various topics from the first two chapters.
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.
Strong electronic correlations are the root of many fascinating phenomena of condensed matter physics. However, they also make calculations notoriously difficult. This thesis is concerned with the development and improvement of computational methods for strongly correlated electron systems. After a brief recapitulation of the Green’s function formalism,which is the theoretical basis of most such methods, we discuss the analytical continuation from imaginary to real frequencies. Especially the maximum entropy method is explained in detail and an implementation there of is presented.The remainder of this thesis deals with methods to compute Green’s functions on thei maginary frequency axis. First, the Anderson impurity model (AIM) is considered. It can be solved numerically – a task for which several methods exist. We focus on continuous-time quantum Monte Carlo with worm sampling, through which various kinds of correlation functions can be calculated. On this basis we derive and implement symmetric improved estimators that allow for the computation of the self-energy and vertex functions where the level of stochastic noise is largely reduced with respect to conventional estimators.As a next step the dynamical mean-field theory (DMFT) is introduced as a method to capture local correlations in the Hubbard model (HM) by self-consistently mapping it onto the AIM. We study the convergence behavior of this iterative procedure and the improvement of convergence due to mixing. Then, we compare electron scattering rates from DMFT in the weak-coupling regime of the square-lattice HM to results from Boltzmann scattering theory.While local correlations account for many interesting phenomena, such as the Mott-Hubbard metal-insulator transition, they are not sufficient to correctly describe, e.g., the physics of phase transitions in two- or three-dimensional systems. It is thus the final step in this thesis to also include nonlocal correlations by means of the dynamical vertex approximation (DΓA). For this we develop a new self-consistency scheme, which is then applied to several systems: We study the square-lattice HM, where we compare to benchmark methods in the weak-coupling regime, but also present results for intermediate and strong coupling. Furthermore, a two-band model with doping and a three-band Wannier model for strontium vanadate are studied. As an application of our new self-consistent DΓA, we study the effects of geometrical frustration in the kagome lattice, where we also compare our results to determinant quantum Monte Carlo.
An essential role for the description of correlated quantum many-particle systems is played by impurity problems. These consist of a small number of localized orbitals, where the electrons are subject to a Coulomb-interaction, coupled to noninteracting bath degrees of freedom. Impurity models are used for the description of magnetic impurities in metals, of quantum points in nanostructures, and in the context of the dynamical mean field theory (DMFT). For their solution a large number of numerical methods exists. We mention here quantum Monte Carlo methods, the numerical renormalization group, or the exact diagonalization method. In this thesis we introduce a new approach to impurity problems, which is based on the functional renormalization group (fRG). In contrast to the methods mentioned above this approach is approximative, but it can be applied with comparatively lower numerical effort. The aim of this thesis is to examine the advantages and limits of this method. An important role in this connection is played by the single impurity Anderson model (SIAM), which is already understood very well with other methods and is seen as a prototype for models with strong local correlations. The spin fluctuations of this model are governed by the so-called Kondo energy scale and a central question of this thesis is whether this scale can be reproduced with our approach. We begin with a short introduction to the Green's function formalism for correlated quantum many particle systems. Then the fRG flow equations for the one-particle irreducible vertex functions are deduced, and several approximation schemes are discussed. In chapter 4 we introduce the fRG scheme used in this thesis. At first the SIAM is mapped to a semi-infinite chain, in which the interacting orbital is given by the first site of this chain. The system is then subdivided into two parts. The first part, which is called 'core', consists of the interacting orbital and the first L bath sites. The remaining bath sites form the second part. At the beginning both parts are decoupled from each other such that the core can be solved exactly. Starting with this exact solution the coupling between the core and the remaining bath sites is switched on slowly in a renormalization group flow. This way it is possible to calculate the one-particle and two-particle correlation functions on the correlated site. We call this flow scheme 'hybridization flow'. The flow equations are formulated in an effective theory on the first bath site outside the core, which turns out to be advantageous compared to other implementations of the hybridization flow. In chapter 5 this hybridization flow method is applied to the SIAM with semi-elliptic bath density of states. We discuss the differences which arise for different core-sizes (L = 0,1,2,3) and for different truncations of the fRG flow equations. In the cases L=0 and L=2 the local Fermi liquid properties of the SIAM and the associated Kondo scale can not be reproduced in both approximation schemes. For the core sizes L=1 and L=3 the ground state of the core is already a spin singlet state in the beginning of the flow and the Fermi liquid properties are reproduced. The Kondo energy scale, which determines the width of the central quasi-particle peak and the size of the local spin-susceptibility, is however not accurately resolved. This is partly due to the finite temperature used in our implementation such that for larger interaction strength the temperature is above the Kondo energy scale. The application of the hybridization flow method to quantum impurity problems in the context of the DMFT is discussed in chapter 6. In this approach lattice models with a local interaction, like the Hubbard model, are mapped to an effective impurity problem for a single lattice site embedded in a dynamical mean field representing the influence of the other electrons. This scheme is exact in infinitely many dimensions and a nonperturbative approximation method for finite-dimensional systems. In particular the Mott metal-insulator transition in the Hubbard model is successfully described by this method. One class of quantities, which is easily accessible from the hybridization flow scheme, is given by the local two-particle correlation functions. These play an essential role in non-local extensions of the DMFT. Therefore we mainly focus on these quantities. At first we apply our flow scheme to the Hubbard model on the Bethe lattice in infinite dimensions and calculate the local one-particle irreducible vertex functions in the insulating and the metallic phase in good agreement with previous calculations, which use the exact diagonalization method as impurity solver. After this we extend our flow scheme to a cluster DMFT method, which includes short-range antiferromagnetic correlations and calculate again the local and non-local vertex functions for the two-dimensional Hubbard model. We compare these to the results of the single-site DMFT.
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.