Abstract In recent years, an increasing number of systems displaying exotic quantum states like unconventional superconductivity, quantum spin‐liquids, or topological states were experimentally found. Here we summarize findings in quantum Monte Carlo simulations of correlated electrons on a honeycomb lattice, the structure of graphene, that revealed an unexpected spin‐liquid emerging between a state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. Moreover, we found that this quantum‐disordered state is a resonating valence‐bond (RVB) liquid, akin to the one proposed for high temperature superconductors. This was the first unbiased determination of a RVB‐liquid in an electronic system.
Abstract In recent years, an increasing number of systems displaying exotic quantum states like unconventional superconductivity, quantum spin‐liquids, or topological states were experimentally found. Here we summarize findings in quantum Monte Carlo simulations of correlated electrons on a honeycomb lattice, the structure of graphene, that revealed an unexpected spin‐liquid emerging between a state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. Moreover, we found that this quantum‐disordered state is a resonating valence‐bond (RVB) liquid, akin to the one proposed for high temperature superconductors. This was the first unbiased determination of a RVB‐liquid in an electronic system.
The Hubbard model encapsulates the physics of strongly correlated quantum systems in its most basic form. It has been studied intensively in the context of the high-temperature superconductivity. A number of novel phases were recently proposed for Hubbard-like models on the honeycomb lattice, the structure of graphene. We analyzed the Hubbard model of spin- 21 fermions on the honeycomb lattice at half-filling using large-scale quantum Monte Carlo simulations. We find that the weak coupling semimetal and the antiferromagnetic Mott insulator at strong interaction are separated by an extended gapped phase in an intermediate coupling regime. Exploring excitation gaps, various correlation functions as well as probing for flux quantization, we conclude that a quantum spin liquid, lacking any conventional order, emerges with local charge and spin correlations, best described by a resonating valence bonds state.
Resonating valence bond (RVB) states are fundamental for understanding quantum spin liquids in two-dimensional (2D) systems. The RVB state is a collective phenomenon in which spins are uncoupled. 2D lattices such as triangular, honeycomb, and dice lattices were investigated using the Hubbard model and exact diagonalization method. We analyzed the total spin, spin-spin correlation functions, local magnetic moments, and spin and charge gaps as a function of on-site Coulomb repulsion, electron concentration, and electronic hopping parameters. Phase diagrams showed that RVB states can live in half-filled and hole-doped anisotropic triangular lattices. We found two types of RVB states: one in the honeycomb sublattice and the other in the centered hexagons in the triangular lattices. Owing to the novel discovery of exotic magnetic ordering in triangular moiré patterns in twisted bilayer graphene and transition metal dichalcogenide systems, our results provide physical insights into the onset of magnetism and possible spin liquid states in these layered materials.
Three models of strongly correlated electron systems have been studied in this thesis: the Shastry-Sutherland quantum antiferromagnet, the quantum spin liquid emerging out of correlated Dirac fermions on the honeycomb lattice, and the edge-state magnetism in the graphene nanoribbons. The methods applied are quantum Monte Carlo simulations, namely, stochastic series expansion quantum Monte Carlo for bosonic (spin) systems and projector auxiliary field quantum Monte Carlo for fermionic systems. The stochastic series expansion quantum Monte Carlo method is applied to study the ground state properties as well as the magnetization process of the spin-1/2 easy-axis Heisenberg model on the Shastry-Sutherland lattice. The quantum phase transitions between a N\'eel state, a dimer triplet state and a ferromagnetic XY phase have been characterized by finite size scaling analysis. The crystalline structures of the triplon magnetic excitations in the magnetic plateaus have been identified. The projector auxiliary field quantum Monte Carlo is applied to explore the ground state phase diagram of the SU(2) Hubbard model on the honeycomb lattice. The interplay between the relativistic Dirac dispersion, the vanishing density of states at half-filling and the strong quantum fluctuations, drives the system from a semimetal at weak coupling to an antiferromagnetic Mott insulator in the strong coupling regime. In between these phases there emerges an unexpected quantum spin liquid with no symmetry breaking in any channel. After careful finite size scaling analysis of various correlation functions, the quantum spin liquid phase is characterized as a short-range resonating valence bond state with a finite spin gap to all magnetic excitations. The study of edge-state magnetism in graphene nanoribbons is more experimentally motivated. The projector auxiliary field quantum Monte Carlo method is used to measure the static spin-spin correlation function, the single-particle spectral function, as well as the dynamic spin structure factor along the zigzag edge of graphene nanoribbons. A dominant low energy peak generated by the Coulomb interaction in the single-particle spectral function is identified as the dynamical signature of edge-state magnetism. The magnetic excitations along the zigzag edge show a linear dispersion at low energies due to the effective antiferromagnetic coupling of the edges across the ribbon. The magnetic excitations are damped by scattering processes into the particle-hole continuum at high energies.
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.
In the past few decades, tremendous efforts have been made toward understanding the exotic physics emerging from competition between various ordering tendencies in strongly correlated systems. Employing state-of-the-art quantum Monte Carlo simulation, we investigate an interacting SU(N) fermionic model with varying interaction strength and value of N, and we unveil the ground-state phase diagram of the model exhibiting a plethora of exotic phases. For small values of N-namely, N=2, 3-the ground state is an antiferromagnetic (AFM) phase, whereas in the large-N limit, a staggered valence bond solid (VBS) order is dominant. For intermediate values of N such as N=4, 5, remarkably, our study reveals that distinct VBS orders appear in the weak and strong coupling regimes. More fantastically, the competition between staggered and columnar VBS ordering tendencies gives rise to a Mott insulating phase without spontaneous symmetry breaking (SSB), existing in a large interacting parameter regime, which is consistent with a gapped quantum spin liquid. Our study not only provides a platform to investigate the fundamental physics of quantum many-body systems-it also offers a novel route toward searching for exotic states of matter such as quantum spin liquid in realistic quantum materials.
At sufficiently low temperatures, interacting electron systems tend to develop orders. Exceptions are quantum critical point (QCP) and quantum spin liquid (QSL), where fluctuations prevent the highly entangled quantum matter to an ordered state down to the lowest temperature. While the ramification of these states may have appeared in high-temperature superconductors, ultracold atoms, frustrated magnets, and quantum moir\'e materials, their unbiased presence remains elusive in microscopic two-dimensional lattice models. Here, we show, by means of large-scale quantum Monte Carlo simulations of correlated electrons on the $\ensuremath{\pi}$-flux square lattice subjected to plaquette Hubbard interaction, that a Gross-Neveu QCP separating massless Dirac fermions and a columnar valence bond solid at finite interaction and a possible Dirac QSL at the infinite yet tractable interaction limit emerge in a coherent sequence. These unexpected quantum states reside in this simple-looking model, unifying ingredients including emergent symmetry, deconfined fractionalization, and the dynamic coupling between emergent matter and gauge fields and will have profound implications both in quantum many-body theory and understanding of the aforementioned experimental systems.