We investigate the quantum dynamics of monopolelike excitations in quantum square ice, as captured by the strongly anisotropic spin-1/2 XXZ model on the checkerboard lattice. We obtain exact results for excitation dynamics in both analytically solvable effective models and a fully interacting model of quantum square ice on finite clusters. We find that the dispersive lower bound of the dynamic response of freely propagating spinons is recovered in the dynamic structure factor of the interacting system, yielding a marked fingerprint of coherent spinon dispersion. Our results provide unbiased evidence for the formation of coherent quasiparticles propagating freely in the correlated ``vacuum'' of quantum square ice.
We study the quantum spin dynamics of a frustrated XXZ model on a pyrochlore lattice by using large-scale quantum Monte Carlo simulation and stochastic analytic continuation. In the low-temperature quantum spin ice regime, we observe signatures of coherent photon and spinon excitations in the dynamic spin structure factor. As the temperature rises to the classical spin ice regime, the photon disappears from the dynamic spin structure factor, whereas the dynamics of the spinon remain coherent in a broad temperature window. Our results provide experimentally relevant, quantitative information for the ongoing pursuit of quantum spin ice materials.
The antiferromagnetic Ising model on a checkerboard lattice has an ice-like ground state manifold with extensive degeneracy. and, to leading order in
Jxy, deconfined spinon excitations. We explore the role of cyclic exchange arising at order
Jxy2∕Jz on the ice states and their associated spinon excitations. By mapping the original problem onto an equivalent quantum six-vertex model, we identify three different phases as a function of the chemical potential for flippable plaquettes---a phase with long range N\'eel order and confined spinon excitations, a nonmagnetic state of resonating square plaquettes, and a quasicollinear phase with gapped but deconfined spinon excitations. The relevance of the results to the square-lattice quantum dimer model is also discussed.
Frustrated spin-ice systems support emergent gauge fields and fractionalized quasiparticles that act as magnetic monopoles. Although artificial platforms have enabled their direct visualization, access to their quantum-coherent dynamics has remained limited. Here we realize a programmable dipolar square spin-ice model using a superconducting-qubit quantum annealer, providing access to a previously unexplored quantum-coherent regime of artificial spin ice. By implementing a direct one-to-one mapping between lattice spins and physical qubits, together with engineered extended couplings, we realize effective dipolar interactions on frustrated lattices comprising more than 400 vertices. Tuning transverse-field fluctuations enables us to probe the real-time dynamics of Dirac-string defects and interacting monopole plasmas. We observe super-diffusive monopole transport, with scaling exponents intermediate between classical diffusion and ballistic motion, indicating dynamics beyond classical stochastic relaxation and consistent with coherent propagation within an emergent gauge manifold. These results establish programmable quantum spin ice as a scalable platform for investigating fractionalized excitations and emergent gauge dynamics in engineered quantum matter.
Three dimensional frustrated magnets realized in the rare-earth pyrochlores R_2B_2O_7, R=rare-earth ion, B transition metal ion) have attracted much interest since the discovery of spin ice in Ho_2Ti_2O_7 by Harris et al. Since then, a number of interesting phenomena in pyrochlore magnets have been observed and the interest in these systems, from an experimental and theoretical point of view, has nothing but increased. Perhaps the most appealing and unique possibility is that these systems may realize a quantum disordered phase; a so-called quantum spin ice. Theoretically, it has been predicted by Hermele et al. that a U(1) quantum spin liquid is realized in the XXZ quantum spin ice (with global U(1)xZ_2 symmetry) in the perturbative regime and more recently Savary and Balents have introduced a quantum rotor model with bosonic spinons coupled to
U(1) gauge fields to describe non-peturbatively the deconfined phase predicted. They studied an anisotropic quantum spin ice (QSI) model which corresponds to the XXZ QSI with an additional anisotropic coupling. They solved the quantum rotor model using a mean-field approximation combined with a large-N approximation and predicted the possibility of a U(1) liquid in addition to an exotic ferromagnetic phase. Further work by Lee et al. has addressed the more complicated fully frustrated XXZ quantum spin ice regime, where all components of the model are frustrated, by using some insights provided by the previous perturbative results. \n \nIn this thesis, we will address the non-perturbative regime of the fully frustrated XXZ QSI beyond the simplest perturbative treatment. We will begin by mapping the quantum spin ice model to a compact U(1) gauge theory with deconfined electric and magnetic monopoles. Neglecting the interaction between the gapless photon excitation and the electric charges, we will study the matter sector with only electric charges, the so-called spinons. Building up on insights provided by high order degenerate perturbation theory calculations, we will construct some symmetric spin liquid ansatz for the static part of the gauge fields. The two novel spin liquids that we construct break time reversal symmetry and thus have a non-trivial chirality. In order to solve the matter sector of the theory and determine the stability of the deconfined phases, we will introduce a novel representation of a XY quantum rotor in terms of ``exclusive" bosons, a representation which is dynamically implemented in the dilute limit. We will benchmark this mapping in different ways by comparing our results with Quantum Monte Carlo simulations and exact diagonalization in some limiting cases. Using this mapping we will compute the region of stability of the U(1) liquid phase of XXZ QSI. Within this mapping we will reproduce qualitatively the results of previous work by finding only the so-called zero and \\pi-flux U(1) liquids in the vicinity of the classical spin ice point, while improving greatly on the estimation of the exchange coupling at which the bosonic spinons condense to form a superfluid phase. In the fully frustrated XXZ QSI, we will provide evidence, using exact diagonalization of the quantum rotor model, that one of the chiral spin liquid may be realized in the strongly correlated and frustrated QSI regime, where all previous approximations break down.
We study the excitation spectra of two-dimensional spin-
1/2 Heisenberg model with checkerboard structures using stochastic analytic continuation of imaginary-time correlation function obtained from quantum Monte Carlo simulation. The checkerboard models have two different antiferromagnetic nearest-neighbor interactions
J1 and
J2, and the tuning parameter
g is defined as
J2/J1. The dynamic structure factors are systematically calculated in all phases of the models as well as at the critical points. To give a full understanding of the dynamic spectra, wave theory is employed to explain some features of numerical results, especially for the low-energy part. When
g is close to
1, the features of the excitation spectra of each checkerboard model are roughly the same as those of the original square lattice antiferromagnetic Heisenberg model, and the high-energy continuum among them are discussed. In contrast to the other checkerboard structures investigated in this paper, the
3×3 checkerboard model has distinctive excitation features, like a gap between a low-energy gapless branch and a gapped high-energy part exists when
g is small. The gapless branch in this case can be regarded as wave in N
eˊel order formed by block spin in each
3×3 plaquette with effective exchange interaction originated from renormalization. One finding is that the continuum also appears in this low-energy branch, which may originated from nearly deconfined process in block spin.