Strongly correlated systems on geometrically frustrated lattices can stabilize a large number of interesting phases that includes a wide array of novel Mott insulators in both bosonic and electronic systems. Charge fluctuations in a Mott insulator are suppressed due to strong mutual interaction among the particles. The presence of frustration is of particular importance as the physics it offers is often rich, unexpectedly complicated, and continues to raise many open questions. The thesis elucidates some of these issues on a kagome lattice where strong interactions among the particles in the Mott phase impose non-trivial local constraints depending on the filling fraction on the lattice. These Mott insulators, in addition to featuring unusual magnetic and/or charge ordering, can also harbor topologically ordered states of quantum matter, e.g., resonating valence bond liquids realized in certain quantum dimer models on non-bipartite lattices. The dimer models can be regarded as low-energy effective theories for different types of bosonic models in the strong-coupling limit. Exploring this connection is a central theme of this thesis with the aim of realizing novel strongly correlated ground states. Past studies of these models have revealed the existence of various ordered and disordered phases with distinct signatures. Among these low-energy phases, the presence of a stable topological liquid at a particular point, known as Rokhsar-Kivelson point, in the phase diagram is notable. The classical versions of the dimer model are also known to have garnered a vast interest in various fields ranging from problems of pure mathematical origin to ones in physical chemistry as well as statistical physics. Pioneered by Kasteleyn, several analytical works came forward to exactly calculate the partition function of the problem from which other physical observables can be derived. Classical numerical methods are extensively applied to these models to verify the analytical predictions. We introduce a new classical algorithm here to compute the correlation functions of a classical dimer model on a square (bipartite) and a triangular (non-bipartite) lattice based on a tensor network construction. The method, called tensor network renormalization group, turns out to be a powerful tool for simulating short-ranged gapped systems as inferred from our results benchmarked against the classical Monte-Carlo technique and compared with past analytical studies. One should note that the quantum dimer model at the Rokhsar-Kivelson point can also be described as an infinite temperature canonical ensemble of classical dimers because of the particular structure of the ground state which is an equal weight superposition in the configuration manifold. The geometry of the lattice plays a pivotal role in deciding the nature of the phases that arise in the dimer models. Many physical properties of the dimer liquid phase can be extracted in the simple classical setting which certainly allows for a deep understanding of the classical models to be developed. The liquid phase is gapped on non-bipartite lattices and gapless on bipartite lattices, which is reflected in the decay of correlation functions with spatial distances. In general on non-bipartite lattices, the topological nature of the dimer liquid is characterized by a Z2 topological order which survives even when the model is perturbed away from the Rokhsar-Kivelson point. Stability of this liquid phase not only depends on the lattice geometries but notably on dimer concentrations also. In this context, we focus on a particular variant of the dimer model on a triangular lattice which is known as the quantum fully packed loop model. The model is composed of nonintersecting closed loops made of dimers and governed by the same Hamiltonian as the quantum dimer model. The loop model provides an effective low-energy description of a strongly correlated bosonic system at 1/3 filling on the kagome lattice. The corresponding Bose-Hubbard Hamiltonian consists of nearest-neighbor hopping and all possible repulsive interactions within a hexagonal plaquette. Conspicuous features of the zero-temperature phase diagram for this model include (i) presence of a stable Z2 liquid even without any Rokhsar-Kivelson potential term (in distinction to the standard quantum dimer model), and (ii) an unconventional phase transition from the liquid phase to a novel crystalline phase that has nematic order (dubbed lattice nematic). For a deeper understanding of the physics, a mapping to an Ising gauge theory is presented. The gauge theoretic description provides a useful way to predict the nature of the quantum phase transition to lie in the O(3) universality class. Finally a fermionic model at the same 1/3 filling is considered in which the ground state exhibits a number of exotic local orderings resulting from the spin-charge interplay of electrons. The Hamiltonian comprises nearest-neighbor hopping, strong on-site Coulomb interaction, and repulsive interaction terms only between nearest-neighbors. In the strong correlation limit, this fermionic problem maps to a two-color fully packed loop model – a model in which the loop segments carry an additional quantum number as color on a honeycomb lattice. The effective theory is governed by coherent three-particle ring exchanges and nearest-neighbor antiferromagnetic spin exchanges. The competition between these two leads to a phase diagram composed of a novel plaquette ordered state (known as the plaquette phase) that undergoes phase transition to a new kind of charge ordered state which we call a short loop phase. From our numerical analysis, we conclude that the plaquette phase features an unusual antiferromagnetic order with gapless spin excitations while the charge-ordered state is subjugated by spin fluctuations of localized electrons arranged in small hexagonal loops on the kagome lattice.
Strongly correlated systems on geometrically frustrated lattices can stabilize a large number of interesting phases that includes a wide array of novel Mott insulators in both bosonic and electronic systems. Charge fluctuations in a Mott insulator are suppressed due to strong mutual interaction among the particles. The presence of frustration is of particular importance as the physics it offers is often rich, unexpectedly complicated, and continues to raise many open questions. The thesis elucidates some of these issues on a kagome lattice where strong interactions among the particles in the Mott phase impose non-trivial local constraints depending on the filling fraction on the lattice. These Mott insulators, in addition to featuring unusual magnetic and/or charge ordering, can also harbor topologically ordered states of quantum matter, e.g., resonating valence bond liquids realized in certain quantum dimer models on non-bipartite lattices. The dimer models can be regarded as low-energy effective theories for different types of bosonic models in the strong-coupling limit. Exploring this connection is a central theme of this thesis with the aim of realizing novel strongly correlated ground states. Past studies of these models have revealed the existence of various ordered and disordered phases with distinct signatures. Among these low-energy phases, the presence of a stable topological liquid at a particular point, known as Rokhsar-Kivelson point, in the phase diagram is notable. The classical versions of the dimer model are also known to have garnered a vast interest in various fields ranging from problems of pure mathematical origin to ones in physical chemistry as well as statistical physics. Pioneered by Kasteleyn, several analytical works came forward to exactly calculate the partition function of the problem from which other physical observables can be derived. Classical numerical methods are extensively applied to these models to verify the analytical predictions. We introduce a new classical algorithm here to compute the correlation functions of a classical dimer model on a square (bipartite) and a triangular (non-bipartite) lattice based on a tensor network construction. The method, called tensor network renormalization group, turns out to be a powerful tool for simulating short-ranged gapped systems as inferred from our results benchmarked against the classical Monte-Carlo technique and compared with past analytical studies. One should note that the quantum dimer model at the Rokhsar-Kivelson point can also be described as an infinite temperature canonical ensemble of classical dimers because of the particular structure of the ground state which is an equal weight superposition in the configuration manifold. The geometry of the lattice plays a pivotal role in deciding the nature of the phases that arise in the dimer models. Many physical properties of the dimer liquid phase can be extracted in the simple classical setting which certainly allows for a deep understanding of the classical models to be developed. The liquid phase is gapped on non-bipartite lattices and gapless on bipartite lattices, which is reflected in the decay of correlation functions with spatial distances. In general on non-bipartite lattices, the topological nature of the dimer liquid is characterized by a Z2 topological order which survives even when the model is perturbed away from the Rokhsar-Kivelson point. Stability of this liquid phase not only depends on the lattice geometries but notably on dimer concentrations also. In this context, we focus on a particular variant of the dimer model on a triangular lattice which is known as the quantum fully packed loop model. The model is composed of nonintersecting closed loops made of dimers and governed by the same Hamiltonian as the quantum dimer model. The loop model provides an effective low-energy description of a strongly correlated bosonic system at 1/3 filling on the kagome lattice. The corresponding Bose-Hubbard Hamiltonian consists of nearest-neighbor hopping and all possible repulsive interactions within a hexagonal plaquette. Conspicuous features of the zero-temperature phase diagram for this model include (i) presence of a stable Z2 liquid even without any Rokhsar-Kivelson potential term (in distinction to the standard quantum dimer model), and (ii) an unconventional phase transition from the liquid phase to a novel crystalline phase that has nematic order (dubbed lattice nematic). For a deeper understanding of the physics, a mapping to an Ising gauge theory is presented. The gauge theoretic description provides a useful way to predict the nature of the quantum phase transition to lie in the O(3) universality class. Finally a fermionic model at the same 1/3 filling is considered in which the ground state exhibits a number of exotic local orderings resulting from the spin-charge interplay of electrons. The Hamiltonian comprises nearest-neighbor hopping, strong on-site Coulomb interaction, and repulsive interaction terms only between nearest-neighbors. In the strong correlation limit, this fermionic problem maps to a two-color fully packed loop model – a model in which the loop segments carry an additional quantum number as color on a honeycomb lattice. The effective theory is governed by coherent three-particle ring exchanges and nearest-neighbor antiferromagnetic spin exchanges. The competition between these two leads to a phase diagram composed of a novel plaquette ordered state (known as the plaquette phase) that undergoes phase transition to a new kind of charge ordered state which we call a short loop phase. From our numerical analysis, we conclude that the plaquette phase features an unusual antiferromagnetic order with gapless spin excitations while the charge-ordered state is subjugated by spin fluctuations of localized electrons arranged in small hexagonal loops on the kagome lattice.
The quantum nature of electrons emerges in various exotic ways when they are put together with quantum many-body interactions. In Mott insulating systems with strong electron correlation, frustration and strong quantum fluctuations promote localized electron spins to form nontrivial quantum paramagnetic phases such as quantum spin liquid and valence bond solid preventing the development of conventional magnetic orders. In this thesis, we study such exotic quantum phases in various different frustrated quantum magnets. First, encouraged by the spin liquid phenomenology observed in Ba3XSb2O9 (X=Cu, Ni), we investigate Z2 spin liquids on a three dimensional triangular lattice by employing the Schwinger boson mean-field theory with projective symmetry group analysis. It is found that a nematic spin liquid is energetically favorable on this lattice structure rather than symmetric spin liquids. Next, we study various valence bond solid phases that can appear in kagome lattice antiferromagnets by generalizing the recently proposed phenomenological Z2 gauge theory for the topological Z2 spin liquid state found in the nearest-neighbor kagome lattice Heisenberg model. This work provides a bridge between the topological spin liquid and valence bond solid phases previously proposed for kagome lattice antiferromagnets via vison condensation in the gauge theory. Then, we consider optical conductivity in Mott insulators as an experimental probe to detect quantum paramagnetic phases like valence bond solid. Applying this idea to deformed kagome lattice antiferromagnet Rb2Cu3SnF12 with a 12-site valence bond solid state, we estimate the size of optical conductivity in this quantum magnet. After this, another quantum paramagnet BiCu2PO6 is carefully investigated to clarify its puzzling ground state nature and quasiparticle excitations. Utilizing recent high resolution neutron scattering data, we construct a microscopic spin model and also an effective theory for the low energy spin excitations in the compound. We find significant influences of magnetic anisotropies on the ground state and quasiparticle dynamics. Lastly, we investigate emergent quantum phases in the thin film geometries of the pyrochlore iridates, where a number of exotic quantum ground states are proposed to occur in bulk materials as a result of the interplay between electron correlation and strong spin-orbit coupling. The fate of these bulk phases, as well as the novel quantum states that may arise only in a thin film setting, are studied via a theoretical model that allows layer-dependent magnetic structures.
Recent work [Kundu, Damle, Phys. Rev. X 15, 011018 (2025)] has demonstrated that a class of S=1 kagome lattice antiferromagnets with competing anisotropies can exhibit unusual low-temperature behavior that is best described in terms of the statistical mechanics of an interesting interpolation between the fully-packed dimer model and the fully-packed loop model on the bipartite honeycomb lattice. On this and other bipartite lattices in two and three dimension, such fully-packed dimer-loop models generically have two distinct power-law correlated Coulomb phases, both of which admit a low-energy description in terms of a fluctuating divergence-free polarization vector field [Phys. Rev. X 15, 011018 (2025); arXiv:2510.11607]. The transition between these phases corresponds to a confinement transition for the Z2 part of the corresponding U(1) flux [Phys. Rev. X 15, 011018 (2025); arXiv:2510.11607]. Here we demonstrate that a variant of such a topological confinement transition can also exist in dimer-loop models on non-bipartite lattices such as the kagome and triangular lattices although the low-energy description does not feature any conserved U(1) flux. In these non-bipartite dimer-loop models, both phases are short-range correlated liquids, and the topological distinction between them is in terms of a Z2-valued winding number instead of the flux of a divergence-free polarization field. In spite of this key difference, this transition appears to be in the same Ising universality class as the analogous transition of the bipartite dimer-loop models. These kagome and triangular lattice dimer-loop models correspond respectively to the low-energy limit of a class of S=1 antiferromagnets with competing anisotropies on the ruby lattice comprising a corner-sharing arrangement of fully-frustrated squares and a lattice of corner-sharing fully-frustrated hexagons.
We derive exact results for close-packed dimers on the triangular kagome lattice (TKL), formed by inserting triangles into the triangles of the kagome lattice. Because the TKL is a nonbipartite lattice, dimer-dimer correlations are short ranged so that the ground state at the Rokhsar-Kivelson (RK) point of the corresponding quantum dimer model on the same lattice is a short-ranged spin liquid. Using the Pfaffian method, we derive an exact form for the free energy, and we find that the entropy is 31ln2 per site, regardless of the weights of the bonds. The occupation probability of every bond is 41 in the case of equal weights on every bond. Similar to the case of lattices formed by corner-sharing triangles (such as the kagome and squagome lattices), we find that the dimer-dimer correlation function is identically zero beyond a certain (short) distance. We find in addition that monomers are deconfined on the TKL, indicating that there is a short-ranged spin liquid phase at the RK point. We also find exact results for the ground-state energy of the classical Heisenberg model. The ground state can be ferromagnetic, ferrimagnetic, locally coplanar, or locally canted, depending on the couplings. From the dimer model and the classical spin model, we derive upper bounds on the ground-state energy of the quantum Heisenberg model on the TKL.
We analyze the antiferromagnetic spin-1/2 XXZ model on the kagome lattice at finite external magnetic field with the help of a non-perturbative zero-temperature renormalization group (RG) technique. The exact nature of the ground and excited state properties (e.g. gapped or gapless spectrum etc) of this system are still debated. Approximate methods have typically been adopted towards understanding the low-energy spectrum. Following the work of Kumar et al (2014 Phys. Rev. B 90 174409), we use a Jordan–Wigner transformation to map the spin problem into one of spinless fermions (spinons) in the presence of a statistical gauge field, and with nearest-neighbor interactions. While the work of Kumar et al was confined mostly to the plateau at 1/3-filling (magnetization per site) in the XY regime, we analyze the role of inter-spinon interactions in shaping the phases around this plateau in the entire XXZ model. The RG phase diagram obtained contains three spin liquid phases whose position is determined as a function of the exchange anisotropy and the energy scale for fluctuations arising from spinon scattering. Two of these spins liquids are topologically ordered states of matter with gapped, degenerate states on the torus. The gap for one of these phases corresponds to the one-spinon band gap of the Azbel–Hofstadter spectrum for the XY part of the Hamiltonian, while the other arises from two-spinon interactions. The Heisenberg point of this problem is found to lie within the interaction gapped spin liquid phase, in broad agreement with a recent experimental finding. The third phase is an algebraic spin liquid with a gapless Dirac spectrum for spinon excitations, and possess properties that show departures from the Fermi liquid paradigm. The three phase boundaries correspond to critical theories, and meet at a SU (2)-symmetric multicritical point. This special critical point agrees well with the gap-closing transition point predicted by Kumar et al . We discuss the relevance of our findings to various recent experiments, as well as results obtained from other theoretical analyses.
The spin-21 Heisenberg model on the kagome lattice, which is closely realized in layered Mott insulators such as ZnCu3(OH)6Cl2, is one of the oldest and most enigmatic spin-21 lattice models. While the numerical evidence has accumulated in favor of a quantum spin liquid, the debate is still open as to whether it is a Z2 spin liquid with very short-range correlations (some kind of resonating valence bond spin liquid), or an algebraic spin liquid with power-law correlations. To address this issue, we have pushed the program started by Rokhsar and Kivelson in their derivation of the effective quantum dimer model description of Heisenberg models to unprecedented accuracy for the spin-21 kagome, by including all the most important virtual singlet contributions on top of the orthogonalization of the nearest-neighbor valence bond singlet basis. Quite remarkably, the resulting picture is a competition between a Z2 spin liquid and a diamond valence bond crystal with a 12-site unit cell, as in the density-matrix renormalization group simulations of Yan et al. Furthermore, we found that, on cylinders of finite diameter d, there is a transition between the Z2 spin liquid at small d and the diamond valence bond crystal at large d, the prediction of the present microscopic description for the two-dimensional lattice. These results show that, if the ground state of the spin-21 kagome antiferromagnet can be described by nearest-neighbor singlet dimers, it is a diamond valence bond crystal, and, a contrario, that, if the system is a quantum spin liquid, it has to involve long-range singlets, consistent with the algebraic spin liquid scenario.
Ever since the introduction of resonating valence bond states by P.\\ W.\\ Anderson in 1970s, the search for exotic many--body quantum behavior in peculiar states of matter known as quantum spin liquids has been a central challenge in the field of condensed matter physics. An important searching ground for such \\emph{quantum--ness} is the low energy physics of frustrated magnets. In a quantum system, frustration---the inability of a physical system to achieve a global state that consistently minimizes its energy locally---enables the unavoidable quantum fluctuations to hybridize an extensive number of classical states into highly entangled quantum states at low temperatures. It can be argued that this is the most practical approach to achieve genuine many--body quantum behavior at large scales. The fundamental source of frustration in many physical systems, is the existence of local non--commuting terms. In the case of antiferromagnet spin theories the simplest examples are triangles of nearest--neighboring sites. Consequently, lattices with triangles as building blocks have long been studied in hopes of discovering spin liquids physics. In particular the kagome lattice, which is composed of corner--sharing triangles, has been the center of at least four decades of spin liquid research. This is due to the underconstrained nature of system made out of corner--sharing triangles. After the first three chapters, introduction in chapter~\\ref{chap:intro}, preliminary physics in chapter~\\ref{chap:prelim}, and a brief introduction to spin--liquids in chapter~\\ref{chap:QSLs}, in chapter~\\ref{chap:kagome} I tackle a simpler version of the long-standing problem of spin--1/2 Heisenberg antiferromagnets by considering a quasi--1D lattice consisting entirely of corner--sharing triangles, \\emph{kagome strip}. I illustrate that the standard Heisenberg antiferromagnet Hamiltonian over kagome strip is an extended gapless quantum phase, that is well characterized by two fermionic/bosonic gapless modes and power--law decaying spin and bond--energy correlations. I also demonstrate that the correlation functions oscillate at tunably incommensurate wave vectors. It turns out that this phase can be identified by a particular marginal instability of a two-band spinon Fermi surface coupled to an emergent U(1) gauge field. This interpretation is supported by analytic Abelian bosonization and with extensive numerical large--scale density matrix renormalization group study as well as variational Monte Carlo calculations on Gutzwiller ans\\"atz wave functions. This intriguing result is the first numerical demonstration of emergent fermionic spinons in a simple SU(2) invariant nearest-neighbor Heisenberg model beyond the strictly 1D (Bethe chain) limit. The unexpected success of the fermionic spinons in describing the low energy physics of kagome antiferromagnet, reintroduces the questions about the validity of spinon physics as an effective theory of quantum spin liquids. However, numerical methods are still far behind the theoretical advances. And still even today, an accurate description of the quantum spin liquid states using tensor network methods is notoriously challenging. It is known that for large quasi-1D systems, the density matrix renormalization group and related methods usually require significant computational resources and sometimes fail to converge to a satisfactory state. On the other hand, variational wavefunctions acquired from the Gutzwiller projection of gaussian fermionic theories have long served as both a theoretical starting point for the construction of such spin liquid states and as an inspiration for numerical variational Monte Carlo (VMC) to calculate observables of interest. Noting this observation I examine a different method by exploring the possibility of constructing a matrix product state (MPS) representation for a Gutzwiller--projected state from two given MPS representations of gaussian fermionic theories in the \\ref{chap:gutz} chapter. I investigate the complexity of different approaches to achieve Gutzwiller projection for MPSs and introduce the novel algorithm which we call the Guzwiller zipper method. The performance of the algorithm is tested against two copies of a single half-filled band of spin-1/2 fermionic spinons. In a successful attempt to describe the nature of spin liquid states on quasi-1D strips of triangular and kagome-like lattices, I apply this method to two MPS of multi-band fermionic spinon theories and compare with the complexity of the traditional VMC approach. In particular, we methodically disprove the conjecture of a spinon fermi surface spin liquid for the triangular lattice. Finally, I conclude this thesis by laying out a bird's eye view of the current and future of spin liquid physics in kagome and triangular lattices and pointing out numerous possible applications of the novel Gutzwiller zipper method. In addition, I would also discuss its possible extensions to more complicated tensor networks as well as versions for higher spatial dimensions.