Antiferromagnetic Ising models on frustrated lattices can realize classical spin liquids, with highly degenerate ground states and, possibly, fractionalized excitations and emergent gauge fields. Motivated by the recent interest in many-body systems in negatively curved space, we study hyperbolic frustrated Ising models. Specifically, we consider nearest-neighbor Ising models on tesselations with odd-length loops in two-dimensional hyperbolic space. For finite systems with open boundaries we determine the ground-state degeneracy exactly, and we perform extensive finite-temperature Monte Carlo simulations to obtain thermodynamic data as well as correlation functions. We show that the shape of the boundary, constituting an extensive part of the system, can be used to control low-energy states: Depending on the boundary, we find ordered or disordered ground states. Our results demonstrate how geometric frustration acts in curved space to produce classical spin liquids.
Antiferromagnetic Ising models on frustrated lattices can realize classical spin liquids, with highly degenerate ground states and, possibly, fractionalized excitations and emergent gauge fields. Motivated by the recent interest in many-body systems in negatively curved space, we study hyperbolic frustrated Ising models. Specifically, we consider nearest-neighbor Ising models on tesselations with odd-length loops in two-dimensional hyperbolic space. For finite systems with open boundaries we determine the ground-state degeneracy exactly, and we perform extensive finite-temperature Monte Carlo simulations to obtain thermodynamic data as well as correlation functions. We show that the shape of the boundary, constituting an extensive part of the system, can be used to control low-energy states: Depending on the boundary, we find ordered or disordered ground states. Our results demonstrate how geometric frustration acts in curved space to produce classical spin liquids.
In the field of frustrated magnetism, Kitaev models provide prototypical examples of exactly solvable quantum spin liquids, in which the spin degrees of freedom fractionalize into Majorana fermions coupled to an emergent Z2 gauge field. While the ground states of these models are already well understood as exhibiting different kinds of Majorana (semi)metals, it has still been an open task to achieve a detailed understanding of the thermodynamics, and, namely, the low-temperature ordering of the gauge field into different Z2 flux configurations. In this thesis, we investigate the thermodynamics of Kitaev systems in two and three spatial dimensions with sign-problem-free quantum Monte Carlo simulations. In a first study of elementary 3D Kitaev models, we verify that the ground state Z2 flux sectors of these systems are entirely determined by their elementary plaquette length – a result which shows the validity of Lieb’s theorem for lattice geometries which lack the geometric requirements for its proof. We closely investigate the low-temperature phase transition associated with gauge-ordering, which is a particular realization of an inverted Ising phase transition, as it occurs in general lattice Z2 gauge theories. Our results corroborate the understanding of this transition as separating different regimes in terms of vison-loop excitations, by showing a clear correlation between the critical temperature and the vison gap. We also introduce the concept of “gauge frustration”, which, for a particular 3D Kitaev model, leads to the suppression of the phase transition and a more complex physical behavior at low temperatures. A second quantum Monte Carlo study focuses on a generalized 2D Kitaev model on a five-coordinated lattice system. Here, the gauge ordering is accompanied by the spontaneous breaking of time-reversal symmetry, a scenario which not only allows for a variety of topological ground states, but also for the occurrence of a phase transition in two spatial dimensions. We show that this model exhibits such a transition at a particularly high temperature scale. Moreover, the system possesses a number of phases where the Z2 fluxes are only partially ordered.
The Lieb-Schultz-Mattis (LSM) theorem and its higher-dimensional extensions forbid the existence of a unique, symmetric, and gapped ground state at fractional fillings in quantum many-body systems with a conserved particle number (or spin angular momentum) and the conventional translation symmetry of Euclidean lattices. In this work, we propose a generalization of the LSM theorem to quantum many-body systems on hyperbolic lattices, i.e., regular tessellations of two-dimensional negatively curved space. By leveraging concepts from hyperbolic band theory in a many-body setting, we adapt Oshikawa's flux-threading argument to periodic hyperbolic lattices with a non-Euclidean (Fuchsian) translation symmetry and compute a lower-bound to the ground-state degeneracy as a function of filling and lattice geometry. We explore the consequences of LSM constraints for gapped phases of hyperbolic quantum matter and suggest frustrated spin models on hyperbolic analogs of the square and triangular lattices as promising platforms for realizing symmetric spin liquids in hyperbolic space.
We introduce a geometrically frustrated classical Ising model, dubbed the “spin-vorticity model,” whose ground-state manifold is a classical spin liquid: a 2-form Coulomb phase. We study the thermodynamics of this model both analytically and numerically, exposing the presence of algebraically decaying correlations and an extensive ground-state entropy, and give a comprehensive account of its ground-state properties and excitations. Each classical ground state may be decomposed into collections of closed two-dimensional membranes, supporting fractionalized string excitations attached to the edges of open membranes. At finite temperature, the model can then be described as a gas of closed strings in a background of fluctuating membranes. The emergent gauge structure of this spin liquid is naturally placed in the language of 2-form electrodynamics, which describes one-dimensional charged strings coupled to a rank-2 antisymmetric gauge field. After establishing the classical spin-vorticity model, we consider perturbing it with quantum exchange interactions, deriving an effective “membrane exchange” model of the quantum dynamics, analogous to ring exchange in quantum spin ice. We demonstrate the existence of a Rokhsar-Kivelson point where the quantum ground state is an equal-weight superposition of all classical ground-state configurations, i.e., a quantum spin liquid. The quantum aspects of this spin liquid are exposed by mapping the membrane exchange model to a strongly coupled frustrated 2-form U(1) lattice gauge theory. We further demonstrate how to quantize the string excitations by coupling a 1-form string field to the 2-form U(1) gauge field, thus mapping a quantum spin model to a 2-form U(1) gauge-Higgs model. We discuss the stability of the gapless deconfined phase of this gauge theory and the possibility of realizing a class of quantum phases of matter: 2-form U(1) quantum spin liquids.
In spin systems, geometrical frustration describes the impossibility of minimizing simultaneously all the interactions in a Hamiltonian, often giving rise to macroscopic ground-state degeneracies and emergent low-temperature physics. In this thesis, combining tensor network (TN) methods to Monte Carlo (MC) methods and ground-state energy lower bound approaches, we study two-dimensional frustrated classical Ising models. In particular, we focus on the determination of the residual entropy in the presence of farther-neighbor interactions in kagome lattice Ising antiferromagnets (KIAFM). In general, using MC to determine the residual entropy is a significant challenge requiring ad-hoc updates, a precise evaluation of the energy at all temperatures to allow for thermodynamic integration, and a good control of the finite-size scaling behavior. As an alternative, we turn to TNs; however, we argue that, in the presence of frustration and macroscopic ground-state degeneracy, standard algorithms fail to converge at low temperatures on the usual TN formulation of partition functions. Inspired by methods for constructing ground-state energy lower bounds, we propose a systematic way to find the ground-state local rule using linear programming. Characterizing the rules as tiles that can be tessellated to form ground states of the model gives rise to a natural contractible TN formulation of the partition function. This method provides a direct access to the ground-state properties of frustrated models and, in particular, allows an extremely precise determination of their residual entropy. We then study two models inspired by artificial spin systems on the kagome lattice with out-of-plane (OOP) anisotropy. The first model is motivated by experiments on an array of chirally coupled nanomagnets. We argue that the farther-neighbor to nearest-neighbor couplings ratios in this system are much smaller than in the dipolar case, J2/J1 being of the order of 2%. A comparison of the experimental correlations with the results of extensive TN and MC simulations shows that (1) the experimental second- and third-neighbor correlations are inverted as compared to those of a pure nearest-neighbor model at equilibrium (even with a magnetic field), and (2) second-neighbor couplings as small as 1% of the nearest-neighbor couplings will affect the spin-spin correlations even at fairly high temperatures. Motivated by dipolar coupled artificial spin systems, we turn to the progressive lifting of the ground-state degeneracy of the KIAFM. We provide a detailed study of the ground-state phases of this model with up to third neighbor interactions, for arbitrary J2, J3 such that J1 >> J2, J3, obtaining exact results for the ground-state energies. When all couplings are antiferromagnetic, we exhibit three macroscopically degenerate ground-state phases and establish their residual entropy using our TN approach. Furthermore, in the phase corresponding to the dipolar KIAFM truncated to third neighbors, we use the ground-state tiles to establish the existence of a mapping to the ground-state manifold of the triangular Ising antiferromagnet.
In this work the ordering processes in classical frustrated Ising models and the ground-state behavior of frustrated quantum spin models are investigated. In the first part the phase transition from a high-temperature paramagnetic phase to a antiferromagnetic ground state is analyzed for the square-lattice Ising model with competing antiferromagnetic interactions between nearest- and next-nearest-neighbor spins. In the case of spatially isotropic interactions the phase transition to a collinear state is of special interest. Using Monte-Carlo simulations on finite lattices as well as field-theoretical methods for the continuous model it is shown that depending on the frustration three different sectors can be distinguished that show (i) Ising-universal behavior (ii) Ashkin-Teller-like non-universal behavior (varying critical exponents) and (iii) non-critical first-order behavior. For spatially anisotropic interactions the emergence of an intermediate incommensurate phase is verified by means of Monte-Carlo simulations which has been predicted by earlier mean-field calculations. The nature of this phase is studied analyzing spin-spin correlation functions and its ordering wave-vector reveals a strong temperature and size dependence. In the second part the influence of quantum fluctuations on classical antiferromagnetic ground states are investigated in anisotropic quantum spin-1/2 Heisenberg models with competing interactions. Using quantum Monte-Carlo simulations, exact diagonalization and perturbation theory the ground-state phase diagrams are calculated for a square lattice with nearest- and next-nearest-neighbor interactions and a honeycomb lattice with up to third-nearest-neighbor interactions whereby the fluctuations are chosen to be non-frustrating. For both cases the stability regions of classical magnetic states are identified and an intermediate phase is detected that shows no sign of long-range magnetic order. Further calculations also rule out valence-bond-crystal or topological order for this state, thus a quantum disordered ground state is stabilized in these models.
In the field of frustrated magnetism, Kitaev models provide a unique framework to study the phenomena of spin fractionalization and emergent lattice gauge theories in two and three spatial dimensions. Their ground states are quantum spin liquids, which can typically be described in terms of a Majorana band structure and an ordering of the underlying Z2 gauge structure. Here we provide a comprehensive classification of the ``gauge physics'' of a family of elementary three-dimensional Kitaev models, discussing how their thermodynamics and ground state order depends on the underlying lattice geometry. Using large-scale, sign-free quantum Monte Carlo simulations we show that the ground-state gauge order can generally be understood in terms of the length of elementary plaquettes---a result which extends the applicability of Lieb's theorem to lattice geometries beyond its original scope. At finite temperatures, the proliferation of (gapped) vison excitations destroys the gauge order at a critical temperature scale, which we show to correlate with the size of vison gap for the family of three-dimensional Kitaev models. We also discuss two notable exceptions where the lattice structure gives rise to ``gauge frustration'' or intertwines the gauge ordering with time-reversal symmetry breaking. In a more general context, the thermodynamic gauge transitions in such 3D Kitaev models are one of the most natural settings for phase transitions beyond the standard Landau-Ginzburg-Wilson paradigm.