In this thesis we consider the hybrid quantum Monte Carlo method for simulations of the Hubbard and Su-Schrieffer-Heeger model. In the first instance, we discuss the hybrid quantum Monte Carlo method for the Hubbard model on a square lattice. We point out potential ergodicity issues and provide a way to circumvent them by a complexification of the method. Furthermore, we compare the efficiency of the hybrid quantum Monte Carlo method with a well established determinantal quantum Monte Carlo method for simulations of the half-filled Hubbard model on square lattices. One reason why the hybrid quantum Monte Carlo method loses the comparison is that we do not observe the desired sub-quadratic scaling of the numerical effort. Afterwards we present a formulation of the hybrid quantum Monte Carlo method for the Su-Schrieffer-Heeger model in two dimensions. Electron-phonon models like this are in general very hard to simulate using other Monte Carlo methods in more than one dimensions. It turns out that the hybrid quantum Monte Carlo method is much better suited for this model . We achieve favorable scaling properties and provide a proof of concept. Subsequently, we use the hybrid quantum Monte Carlo method to investigate the Su-Schrieffer-Heeger model in detail at half-filling in two dimensions. We present numerical data for staggered valence bond order at small phonon frequencies and an antiferromagnetic order at high frequencies. Due to an O(4) symmetry the antiferromagnetic order is connected to a superconducting charge density wave. Considering the Su-Schrieffer-Heeger model without tight-binding hopping reveals an additional unconstrained Z_2 gauge theory. In this case, we find indications for π-fluxes and a possible Z_2 Dirac deconfined phase as well as for a columnar valence bond ordered state at low phonon energies. In our investigations of the several phase transitions we discuss the different possibilities for the underlying mechanisms and reveal first insights into a rich phase diagram. / In der vorliegenden Arbeit betrachten wir die Hybrid-Quanten-Monte-Carlo-Methode für Simulationen des Hubbard- sowie des Su-Schrieffer-Heeger-Modells. Zunächst diskutieren wir die Hybrid-Quanten-Monte-Carlo-Methode am Beispiel des Hubbard-Modells auf dem Quadratgitter. Wir zeigen mögliche Ergodizitätsprobleme auf und präsentieren eine Möglichkeit, diese durch Verwendung komplexwertiger Hilfsfelder zu vermeiden. Für Simulationen des halbgefüllten Hubbard-Modells auf Quadratgittern vergleichen wir die Effizienz der Hybrid-Quanten-Monte-Carlo-Methode mit der einer weit verbreiteten und gebräuchlichen Determinanten-Quanten-Monte-Carlo-Methode. Ein Grund für die Niederlage der Hybrid-Quanten-Monte-Carlo-Methode in diesem Vergleich ist die Skalierung des benötigten Rechenaufwandes. Die erhoffte sub-quadratische Skalierung in Abhängigkeit von Systemgröße und inverser Temperatur wird nicht beobachtet. Anschließend präsentieren wir eine Formulierung der Hybrid-Quanten-Monte-Carlo-Methode zur Untersuchung des halbgefüllten Su-Schrieffer-Heeger-Modells in zwei Dimensionen. Elektron-Phonon-Modelle wie dieses sind in mehr als einer Dimension für gewöhnlich mit anderen Quanten-Monte-Carlo-Methoden nur schwer simulierbar. Es stellt sich heraus, dass sich die Hybrid-Quanten-Monte-Carlo-Methode deutlich besser zur Simulation dieses Modells eignet. Wir erreichen eine vorteilhafte Skalierung des Rechenaufwandes und präsentieren einen Machbarkeitsnachweis. Folglich verwenden wir die Hybrid-Quanten-Monte-Carlo-Methode für nähere Untersuchungen des Su-Schrieffer-Heeger-Modells. Wir zeigen numerische Resultate für eine gestaffelte Ordnung aus Valenzbindungen bei kleinen Phononfrequenzen und für eine antiferromagnetischen Ordnung bei hohen Frequenzen. Aufgrund einer O(4)-Symmetrie ist die antiferromagnetische Ordnung mit einer supraleitenden Ladungsdichtewelle verknüpft. Ohne Tight-Binding-Hüpfparameter offenbart das Su-Schrieffer-Heeger-Modell eine zusätzliche spezielle Z_2 -Eichsymmetrie, die nicht den Satz von Gauß erfüllt. In diesem Fall finden wir Hinweise für einen π-Flux-Zustand. Bei niedrigen Phononenergien gibt es außerdem Anzeichen für einen möglichen Z_2 Dirac deconfined Zustand sowie eine spaltenweise Ordnung von Valenzbindungen. Bei Untersuchungen der Phasenübergänge beleuchten wir die möglichen Mechanismen, die den Übergängen zugrunde liegen. Zum Abschluss diskutieren wir das vielfältige Phasendiagramm, in welches wir erste Einblicke ermöglichen.
| Identifer | oai:union.ndltd.org:uni-wuerzburg.de/oai:opus.bibliothek.uni-wuerzburg.de:19122 |
| Date | January 2020 |
| Creators | Beyl, Stefan |
| Source Sets | University of Würzburg |
| Language | English |
| Detected Language | German |
| Type | doctoralthesis, doc-type:doctoralThesis |
| Format | application/pdf |
| Rights | https://creativecommons.org/licenses/by-sa/4.0/deed.de, info:eu-repo/semantics/openAccess |
Page generated in 0.0074 seconds