@phdthesis{Alboteanu2007, author = {Alboteanu, Ana Maria}, title = {The Noncommutative Standard Model : Construction Beyond Leading Order in Theta and Collider Phenomenology}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-24334}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2007}, abstract = {Trotz seiner pr{\"a}zisen {\"U}bereinstimmung mit dem Experiment ist die G{\"u}ltigkeit des Standardmodells (SM) der Elementarteilchenphysik bislang nur bis zu einer Energieskala von einigen hundert GeV gesichert. Abgesehen davon erweist sich schon das Einbinden der Gravitation in einer einheitlichen Beschreibung aller fundamentalen Wechselwirkungen als ein durch gew{\"o}hnliche Quantenfeldtheorie nicht zu l{\"o}sendes Problem. Das Interesse an Quantenfeldtheorien auf einer nichtkommutativen Raumzeit wurde durch deren Vorhersage als niederenergetischer Limes von Stringtheorien erweckt. Unabh{\"a}ngig davon, kann die Nichtlokalit{\"a}t einer solchen Theorie den Rahmen zur Einbeziehung der Gravitation in eine vereinheitlichende Theorie liefern. Die Hoffnung besteht, dass die Energieskala Lambda_NC, ab der solche Effekte sichtbar werden k{\"o}nnen und f{\"u}r die es einerlei theoretischen Vorhersagen gibt, schon bei der n{\"a}chsten Generation von Beschleunigern erreicht wird. Auf dieser Annahme beruht auch die vorliegende Arbeit, im Rahmen deren eine m{\"o}gliche Realisierung von Quantenfeldtheorien auf nichtkommutativer Raumzeit auf ihre ph{\"a}nomenologischen Konsequenzen hin untersucht wurde. Diese Arbeit ist durch fehlende LHC (Large Hadron Collider) Studien f{\"u}r nichkommutative Quantenfeldtheorien motiviert. Im ersten Teil des Vorhabens wurde der hadronische Prozess pp-> Z gamma -> l+l- gamma am LHC sowie die Elektron-Positron Paarvernichtung in ein Z-Boson und ein Photon am ILC (International Linear Collider) auf nichtkommutative Signale hin untersucht. Die ph{\"a}nomenlogischen Untersuchungen wurden im Rahmen dieses Modells in erster Ordnung des nichtkommutativen Parameters Theta durchgef{\"u}hrt. Eine nichtkommutative Raumzeit f{\"u}hrt zur Brechung der Rotationsinvarianz bez{\"u}glich der Strahlrichtung der einlaufenden Teilchen. Im differentiellen Wirkungsquerschnitt f{\"u}r Streuprozesse {\"a}ussert sich dieses als eine azimuthale Abh{\"a}ngigkeit, die weder im SM noch in anderen Modellen jenseits des SM auftritt. Diese klare, f\"ur nichtkommutative Theorien typische Signatur kann benutzt werden, um nichtkommutative Modelle von anderen Modellen, die neue Physik beschreiben, zu unterscheiden. Auch hat es sich erwiesen, dass die azimuthale Abh{\"a}ngigkeit des Wirkungsquerschnittes am besten daf\"ur geeignet ist, um die Sensitivit{\"a}t des LHC und des ILC auf der nichtkommutativen Skala \$\Lnc\$ zu bestimmen. Im ph{\"a}nomenologischen Teil der Arbeit wurde herausgefunden, dass Messungen am LHC f{\"u}r den Prozess pp-> Z gamma-> l+l- gamma nur in bestimmten F{\"a}llen auf nichtkommutative Effekte sensitiv sind. F{\"u}r diese F{\"a}lle wurde f{\"u}r die nichtkommutative Energieskala Lambda_NC eine Grenze von Lambda_NC > 1.2 TeV bestimmt. Diese ist um eine Gr{\"o}ßenordnung h{\"o}her als die Grenzen, die von bisherigen Beschleunigerexperimenten hergeleitet wurden. Bei einem zuk{\"u}nftigen Linearbeschleuniger, dem ILC, wird die Grenze auf Lambda_NC im Prozess e^+e^- -> Z gamma -> l^+ l^- gamma wesentlich erh{\"o}ht (bis zu 6 TeV). Abgesehen davon ist dem ILC gerade der f{\"u}r den LHC kaum zug{\"a}ngliche Parameterbereich der nichtkommutativen Theorie erschlossen, was die Komplementarit{\"a}t der beiden Beschleunigerexperimente hinsichtlich der nichtkommutativen Parameter zeigt. Der zweite Teil der Arbeit entwickelte sich aus der Notwendigkeit heraus, den G{\"u}ltigkeitsbereich der Theorie zu h{\"o}heren Energien hin zu erweitern. Daf{\"u}r haben wir den neutralen Sektor des nichtkommutativen SM um die n{\"a}chste Ordnung in Theta erg{\"a}nzt. Es stellte sich wider Erwarten heraus, dass die Theorie dabei um einige freie Parameter erweitert werden muss. Die zus{\"a}tzlichen Parameter sind durch die homogenen L{\"o}sungen der Eich{\"a}quivalenzbedingungen gegeben, welche Ambiguit\"aten der Seiberg-Witten Abbildungen darstellen. Die allgemeine Erwartung war, dass die Ambiguit{\"a}ten Feldredefinitionen entsprechen und daher in den Streumatrixelementen verschwinden m\"ussen. In dieser Arbeit wurde jedoch gezeigt, dass dies ab der zweiten Ordnung in Theta nicht der Fall ist und dass die Nichteindeutigkeit der Seiberg-Witten Abbildungen sich durchaus in Observablen niederschl{\"a}gt. Die Vermutung besteht, dass jede neue Ordnung in Theta neue Parameter in die Theorie einf{\"u}hrt. Wie weit und in welche Richtung die Theorie auf nichtkommutativer Raumzeit entwickelt werden muss, kann jedoch nur das Experiment entscheiden.}, subject = {Feldtheorie}, language = {en} } @phdthesis{Parragh2013, author = {Parragh, Nicolaus}, title = {Strongly Correlated Multi-Orbital Systems : A Continuous-Time Quantum Monte Carlo Analysis}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-85253}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2013}, abstract = {In this thesis I present results concerning realistic calculations of correlated fermionic many-body systems. One of the main objectives of this work was the implementation of a hybridization expansion continuous-time quantum Monte Carlo (CT-HYB) algorithm and of a flexible self-consistency loop based on the dynamical mean-field theory (DMFT). DMFT enables us to treat strongly correlated electron systems numerically. After the implementation and extensive testing of the program we investigated different problems to answer open questions concerning correlated systems and their numerical treatment.}, subject = {Monte-Carlo-Simulation}, language = {en} } @phdthesis{Walther2008, author = {Walther, Markus}, title = {Simulation of strain-induced and defect-controlled self-organization of nanostructures}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-27931}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2008}, abstract = {In this PhD thesis, the effect of strain on heteroepitaxial growth is investigated by means of Kinetic Monte Carlo simulations. In this context the lattice misfit, arising from the different lattice constants of the adsorbate and the substrate material, is of particular interest. As a consequence, this lattice misfit leads to long-range elastic strain effects having strong influence on the entire growing crystal and its resulting surface morphology. The main focus of this work is the investigation of different strain relaxation mechanisms and their controlling parameters, revealing interesting consequences on the subsequent growth. Since epitaxial growth is carried out under conditions far away from thermodynamic equilibrium, it is strongly determined by surface kinetics. At this point the relevant kinetic microscopic processes are described, followed by theoretical considerations of heteroepitaxial growth disclosing an overview over several independent methodological streams, used to model epitaxy in different time and length scales, as well as the characterization of misfit dislocations and the classification of epitaxial growth modes based on thermodynamic considerations. The epitaxial growth is performed by means of Kinetic Monte Carlo simulations which allows for the consideration of long range effects in systems with lateral extension of few hundred atoms. By using an off-lattice simulation model the particles are able to leave their predefined lattice sites, which is an indispensable condition for simulating strain relaxation mechanisms. The main idea of our used model is calculating the activation energy of all relevant thermally activated processes by using simple pair potentials and then realizing the dynamics by performing each event according to its probability by means of a rejection-free algorithm method. In addition, the crystal relaxation procedure, the grid-based particle access method, which accelerates the simulation enormously, and the efficient implementation of the algorithm are discussed. To study the influence of long range elastic strain effects, the main part of this work was realized on the two dimensional triangular lattice, which can be treated as a cross section of the real three dimensional case. Chapter 4 deals with the formation of misfit dislocations as a strain relaxation mechanism and the resulting consequences on the subsequent heteroepitaxial growth. We can distinguish between two principally different dislocation formation mechanisms, depending strongly on the sign as well as on the magnitude of the misfit, but also the surface kinetics need to be taken into account. Additionally, the dislocations affect the lattice spacings of the crystal whose observed progression is in qualitative good agreement with experimental results. Furthermore, the dislocations influence the subsequent growth of the adsorbate film, since the potential energy of an adatom is modulated by buried dislocations. A clear correlation between the lateral positions of buried dislocations and the positions of mounds grown on the surface can be observed. In chapter 5, an alternative strain relaxation mechanism is studied: the formation of three dimensional islands enables the particles to approach their preferred lattice spacing. We demonstrate that it is possible to adjust within our simulation model each of the three epitaxial growth modes: Volmer-Weber, Frank-van der Merve or layer-by-layer, and Stranski-Krastanov growth mode. Moreover, we can show that the emerging growth mode depends in principle on two parameters: on the one hand the interaction strength of adsorbate particles with each other, compared to the interaction of adsorbate with substrate particles, and on the other hand the lattice misfit between adsorbate and substrate particles. A sensible choice of these two parameters allows the realization of each growth mode within the simulations. In conclusion, the formation of nanostructures controlled by an underlying dislocation network can be applied in the concept of self-organized pattern formation as well as by the tendency to form ordered arrays of strain-induced three dimensional grown islands. In chapter 6, we extend our model to three dimensions and investigate the effect of strain on growth on bcc(100) surfaces. We introduce an anisotropic potential yielding a stable bcc lattice structure within the off-lattice representation. We can show that the strain built up in submonolayer islands is mainly released at the island edges and the lattice misfit has strong influence on the diffusion process on the plane surface as well as on the situation at island edges with eminent consequences on the appearance of submonolayer islands.}, subject = {Kristallwachstum}, language = {en} } @phdthesis{Weber2008, author = {Weber, Sebastian}, title = {Simulation of self-assembled nanopatterns in binary alloys on the fcc(111) surface}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-27914}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2008}, abstract = {In this PhD thesis, we study the heteroepitaxial crystal growth by means of Monte Carlo simulations. Of particular interest in this work is the influence of the lattice mismatch of the adsorbates relative to the substrate on surface structures. In the framework of an off-lattice model, we consider one monolayer of adsorbate and investigate the emerging nanopatterns in equilibrium and their formation during growth. In chapter 1, a brief introduction is given, which describes the role of computer simulations in the field of the physics of condensed matter. Chapter 2 is devoted to some technical basics of experimental methods of molecular beam epitaxy and the theoretical description. Before a model for the simulation can be designed, it is necessary to make some considerations of the single processes which occur during epitaxial growth. For that purpose we look at an experimental setup and extract the main microscopic processes. Afterwards a brief overview of different theoretical concepts describing that physical procedures is given. In chapter 3, the model used in the simulations is presented. The aim is to investigate the growth of an fcc crystal in the [111] direction. In order to keep the simulation times within a feasible limit a simple pair potential, the Lennard-Jones potential, with continuous particle positions is used, which are necessary to describe effects resulting from the atomic mismatch in the crystal. Furthermore the detailed algorithm is introduced which is based on the idea to calculate the barrier of each diffusion event and to use the barriers in a rejection-free method. Chapter 4 is attended to the simulation of equilibrium. The influence of different parameters on the emerging structures in the first monolayer upon the surface, which is completely covered with two adsorbate materials, is studied. Especially the competition between binding energy and strain leads to very interesting pattern formations like islands or stripes. In chapter 5 the results of growth simulations are presented. At first, we introduce a model in order to realize off-lattice Kinetic Monte Carlo simulations. Since the costs in simulation time are enormous, some simplifications in the calculation of diffusion barriers are necessary and therefore the previous model is supplemented with some elements from the so-called ball and spring model. The next point is devoted to the calculation of energy barriers followed by the presentation of the growth simulations. Binary systems with only one sort of adsorbate are investigated as well as ternary systems with two different adsorbates. Finally, a comparison to the equilibrium simulations is drawn. Chapter 6 contains some concluding remarks and gives an outlook to possible further investigations.}, subject = {Kristallwachstum}, language = {en} } @phdthesis{Lang2010, author = {Lang, Thomas C.}, title = {Quantum Monte Carlo methods and strongly correlated electrons on honeycomb structures}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-53506}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2010}, abstract = {In this thesis we apply recently developed, as well as sophisticated quantum Monte Carlo methods to numerically investigate models of strongly correlated electron systems on honeycomb structures. The latter are of particular interest owing to their unique properties when simulating electrons on them, like the relativistic dispersion, strong quantum fluctuations and their resistance against instabilities. This work covers several projects including the advancement of the weak-coupling continuous time quantum Monte Carlo and its application to zero temperature and phonons, quantum phase transitions of valence bond solids in spin-1/2 Heisenberg systems using projector quantum Monte Carlo in the valence bond basis, and the magnetic field induced transition to a canted antiferromagnet of the Hubbard model on the honeycomb lattice. The emphasis lies on two projects investigating the phase diagram of the SU(2) and the SU(N)-symmetric Hubbard model on the hexagonal lattice. At sufficiently low temperatures, condensed-matter systems tend to develop order. An exception are quantum spin-liquids, where fluctuations prevent a transition to an ordered state down to the lowest temperatures. Previously elusive in experimentally relevant microscopic two-dimensional models, we show by means of large-scale quantum Monte Carlo simulations of the SU(2) Hubbard model on the honeycomb lattice, that a quantum spin-liquid emerges between the state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. This unexpected quantum-disordered state is found to be a short-range resonating valence bond liquid, akin to the one proposed for high temperature superconductors. Inspired by the rich phase diagrams of SU(N) models we study the SU(N)-symmetric Hubbard Heisenberg quantum antiferromagnet on the honeycomb lattice to investigate the reliability of 1/N corrections to large-N results by means of numerically exact QMC simulations. We study the melting of phases as correlations increase with decreasing N and determine whether the quantum spin liquid found in the SU(2) Hubbard model at intermediate coupling is a specific feature, or also exists in the unconstrained t-J model and higher symmetries.}, subject = {Monte-Carlo-Simulation}, language = {en} } @phdthesis{Hofmann2020, author = {Hofmann, Johannes Stephan}, title = {On the interplay of topology and interaction: A quantum Monte Carlo study}, doi = {10.25972/OPUS-20507}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-205071}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2020}, abstract = {Adding interactions to topological (non-)trivial free fermion systems can in general have four different effects: (i) In symmetry protected topological band insulators, the correlations may lead to the spontaneous breaking of some protecting symmetries by long-range order that gaps the topological boundary modes. (ii) In free fermion (semi-)metal, the interaction could vice versa also generate long-range order that in turn induces a topological mass term and thus generates non-trivial phases dynamically. (iii) Correlation might reduce the topological classification of free fermion systems by allowing adiabatic deformations between states of formerly distinct phases. (iv) Interaction can generate long-range entangled topological order in states such as quantum spin liquids or fractional quantum Hall states that cannot be represented by non-interacting systems. During the course of this thesis, we use numerically exact quantum Monte Carlo algorithms to study various model systems that (potentially) represent one of the four scenarios, respectively. First, we investigate a two-dimensional \$d_{xy}\$-wave, spin-singlet superconductor, which is relevant for high-\$T_c\$ materials such as the cuprates. This model represents nodal topological superconductors and exhibits chiral flat-band edge states that are protected by time-reversal and translational invariance. We introduce the conventional Hubbard interaction along the edge in order to study their stability with respect to correlations and find ferromagnetic order in case of repulsive interaction as well as charge-density-wave order and/or additional \$i\$s-wave pairing for attractive couplings. A mean-field analysis that, for the first time, is formulated in terms of the Majorana edge modes suggests that any order has normal and superconducting contributions. For example, the ferromagnetic order appears in linear superposition with triplet pairing. This finding is well confirmed by the numerically exact quantum Monte Carlo investigation. Second, we consider spinless electrons on a two-dimensional Lieb lattice that are subject to nearest-neighbor Coulomb repulsion. The low energy modes of the free fermion part constitute a spin-\$1\$ Dirac cone that might be gapped by several mass terms. One option breaks time-reversal symmetry and generates a topological Chern insulator, which mainly motivated this study. We employ two flavors of quantum Monte Carlo methods and find instead the formation of charge-density-wave order that breaks particle-hole symmetry. Additionally, due to sublattices of unequal size in Lieb lattices, this induces a finite chemical potential that drives the system away from half-filling. We argue that this mechanism potentially extends the range of solvable models with finite doping by coupling the Lieb lattice to the target system of interest. Third, we construct a system with four layers of a topological insulators and interlayer correlation that respects one independent time-reversal and a unitary \$\mathbb{Z}_2\$ symmetry. Previous studies claim a reduced topological classification from \$\mathbb{Z}\$ to \$\mathbb{Z}_4\$, for example by gapping out degenerate zero modes in topological defects once the correlation term is designed properly. Our interaction is chosen according to this analysis such that there should exist an adiabatic deformation between states whose topological invariant differs by \$\Delta w=\pm4\$ in the free fermion classification. We use a projective quantum Monte Carlo algorithm to determine the ground-state phase diagram and find a symmetry breaking regime, in addition to the non-interacting semi-metal, that separates the free fermion insulators. Frustration reduces the size of the long-range ordered region until it is replaced by a first order phase transition. Within the investigated range of parameters, there is no adiabatic path deforming the formerly distinct free fermion states into each other. We conclude that the prescribed reduction rules, which often use the bulk-boundary correspondence, are necessary but not sufficient and require a more careful investigation. Fourth, we study conduction electron on a honeycomb lattice that form a Dirac semi-metal Kondo coupled to spin-1/2 degrees of freedom on a Kagome lattice. The local moments are described by a variant of the Balents-Fisher-Girvin model that has been shown to host a ferromagnetic phase and a \$\mathbb{Z}_2\$ spin liquid at strong frustration. Here, we report the first numerical exact quantum Monte Carlo simulation of the Kondo-coupled system that does not exhibit the negative-sign problem. When the local moments form a ferromagnet, the Kondo coupling induces an anti-ferromagnetic mass term in the conduction-electron system. At large frustration, the Dirac cone remains massless and the spin system forms a \$\mathbb{Z}_2\$ spin liquid. Owing to the odd number of spins per unit cell, this constitutes a non-Fermi liquid that violates Luttinger's theorem which relates the Fermi volume to the particle density in a Fermi liquid. This phase is a specific realization of the so called 'fractional Fermi liquid` as it has been first introduced in the context of heavy fermion models.}, subject = {Monte-Carlo-Simulation}, language = {en} } @phdthesis{Bercx2014, author = {Bercx, Martin Helmut}, title = {Numerical studies of heavy-fermion systems and correlated topological insulators}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-116138}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2014}, abstract = {In this thesis, we investigate aspects of the physics of heavy-fermion systems and correlated topological insulators. We numerically solve the interacting Hamiltonians that model the physical systems using quantum Monte Carlo algorithms to access both ground-state and finite-temperature observables. Initially, we focus on the metamagnetic transition in the Kondo lattice model for heavy fermions. On the basis of the dynamical mean-field theory and the dynamical cluster approximation, our calculations point towards a continuous transition, where the signatures of metamagnetism are linked to a Lifshitz transition of heavy-fermion bands. In the second part of the thesis, we study various aspects of magnetic pi fluxes in the Kane-Mele-Hubbard model of a correlated topological insulator. We describe a numerical measurement of the topological index, based on the localized mid-gap states that are provided by pi flux insertions. Furthermore, we take advantage of the intrinsic spin degree of freedom of a pi flux to devise instances of interacting quantum spin systems. In the third part of the thesis, we introduce and characterize the Kane-Mele-Hubbard model on the pi flux honeycomb lattice. We place particular emphasis on the correlations effects along the one-dimensional boundary of the lattice and compare results from a bosonization study with finite-size quantum Monte Carlo simulations.}, subject = {Schwere-Fermionen-System}, language = {en} } @phdthesis{Volkmann2004, author = {Volkmann, Thorsten}, title = {Lattice gas models and simulations of epitaxial growth}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-13812}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2004}, abstract = {In this PhD thesis, we develop models for the numerical simulation of epitaxial crystal growth, as realized, e.g., in molecular beam epitaxy (MBE). The basic idea is to use a discrete lattice gas representation of the crystal structure, and to apply kinetic Monte Carlo (KMC) simulations for the description of the growth dynamics. The main advantage of the KMC approach is the possibility to account for atomistic details and at the same time cover MBE relevant time scales in the simulation. In chapter 1, we describe the principles of MBE, pointing out relevant physical processes and the influence of experimental control parameters. We discuss various methods used in the theoretical description of epitaxial growth. Subsequently, the underlying concepts of the KMC method and the lattice gas approach are presented. Important aspects concerning the design of a lattice gas model are considered, e.g. the solid-on-solid approximation or the choice of an appropriate lattice topology. A key element of any KMC simulation is the selection of allowed events and the evaluation of Arrhenius rates for thermally activated processes. We discuss simplifying schemes that are used to approximate the corresponding energy barriers if detailed knowledge about the barriers is not available. Finally, the efficient implementation of the MC kinetics using a rejection-free algorithm is described. In chapter 2, we present a solid-on-solid lattice gas model which aims at the description of II-VI(001) semiconductor surfaces like CdTe(001). The model accounts for the zincblende structure and the relevant surface reconstructions of Cd- and Te-terminated surfaces. Particles at the surface interact via anisotropic nearest and next nearest neighbor interactions, whereas interactions in the bulk are isotropic. The anisotropic surface interactions reflect known properties of CdTe(001) like the small energy difference between the c(2x2) and (2x1) vacancy structures of Cd-terminated surfaces. A key element of the model is the presence of additional Te atoms in a weakly bound Te* state, which is motivated by experimental observations of Te coverages exceeding one monolayer at low temperatures and high Te fluxes. The true mechanism of binding excess Te to the surface is still unclear. Here, we use a mean-field approach assuming a Te* reservoir with limited occupation. In chapter 3, we perform KMC simulations of atomic layer epitaxy (ALE) of CdTe(001). We study the self-regulation of the ALE growth rate and demonstrate how the interplay of the Te* reservoir occupation with the surface kinetics results in two different regimes: at high temperatures the growth rate is limited to one half layer of CdTe per ALE cycle, whereas at low enough temperatures each cycle adds a complete layer. The temperature where the transition between the two regimes occurs depends mainly on the particle fluxes. The temperature dependence of the growth rate and the flux dependence of the transition temperature are in good qualitative agreement with experimental results. Comparing the macroscopic activation energy for Te* desorption in our model with experimental values we find semiquantitative agreement. In chapter 4, we study the formation of nanostructures with alternating stripes during submonolayer heteroepitaxy of two different adsorbate species on a given substrate. We evaluate the influence of two mechanisms: kinetic segregation due to chemically induced diffusion barriers, and strain relaxation by alternating arrangement of the adsorbate species. KMC simulations of a simple cubic lattice gas with weak inter-species binding energy show that kinetic effects are sufficient to account for stripe formation during growth. The dependence of the stripe width on control parameters is investigated. We find an Arrhenius temperature dependence, in agreement with experimental investigations of phase separation in binary or ternary material systems. Canonical MC simulations show that the observed stripes are not stable under equilibrium conditions: the adsorbate species separate into very large domains. Off-lattice simulations which account for the lattice misfit of the involved particle species show that, under equilibrium conditions, the competition between binding and strain energy results in regular stripe patterns with a well-defined width depending on both misfit and binding energies. In KMC simulations, the stripe-formation and the experimentally reported ramification of adsorbate islands are reproduced. To clarify the origin of the island ramification, we investigate an enhanced lattice gas model whose parameters are fitted to match characteristic off-lattice diffusion barriers. The simulation results show that a satisfactory explanation of experimental observations within the lattice gas framework requires a detailed incorporation of long-range elastic interactions. In the appendix we discuss supplementary topics related to the lattice gas simulations in chapter 4.}, subject = {Kristallwachstum}, language = {en} } @phdthesis{Beyl2020, author = {Beyl, Stefan}, title = {Hybrid Quantum Monte Carlo for Condensed Matter Models}, doi = {10.25972/OPUS-19122}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-191225}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2020}, abstract = {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.}, subject = {Monte-Carlo-Simulation}, language = {en} } @phdthesis{Hausoel2022, author = {Hausoel, Andreas}, title = {Electronic magnetism in correlated systems: from quantum materials down to Earth's core}, doi = {10.25972/OPUS-25444}, url = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-254444}, school = {Universit{\"a}t W{\"u}rzburg}, year = {2022}, abstract = {In the last decade continuous-time quantum Monte Carlo in the hybridization expansion (CTHYB) was one of the most successful Monte Carlo techniques to describe correlated quantum phenomena in conjunction with dynamical mean field theory (DMFT). The first part of the thesis consists of algorithmical developments regarding CTHYB and DMFT. I provide a complete derivation and an extensive discussion of the expansion formula. We generalized it to treat spin-orbit coupling, and invented the superstate sampling algorithm to make it efficient enough for describing systems with general interactions, crystal fields and spin-orbit coupling at low temperatures. But CTHYB is known to fail in the standard implementation for equal-time correlators, certain higher-order Green's functions and the atomic limit; we discovered that its estimator for the Greens function is also inconsistent for Anderson impurities with finite, discrete baths. I focus then on further improvements of CTHYB that we have conceived and worked on, in particular for f-orbitals and for taking physical symmetries into account in the calculation of the Monte Carlo observables. The second part of the thesis presents selected physical applications of these methods. I show DMFT calculations of highest accuracy for elemental iron and nickel and discover a new mechanism of magnetic ordering in nickel: the ordering of band structure-induced local moments. Then we analyze the stability of this phenomenon under pressure and temperatures, that characterize in the Earth's core. We find, that the mechanism survives these conditions and may give a significant contribution to the generation of the Earth's magnetic field. The next topic is the stability of double Dirac fermions against electronic correlations. We find, that the Coulomb interaction in the corresponding material Bi2 CuO4 are strong enough to destroy the double Dirac cone, and substantial uniform pressure is necessary to restore them. In the last chapter I derive the properties of Higgs and Goldstone bosons from Ginzburg-Landau theory, and identify these excitations in a model of an excitonic magnet.}, subject = {Monte-Carlo-Simulation}, language = {en} }