Spin-liquid Model Research Articles

Dynamics of K2Ni2(SO4)3 governed by proximity to a 3D spin liquid model

Quantum spin liquids (QSLs) have become a key area of research in magnetism due to their remarkable properties, such as long-range entanglement, fractional excitations, and topologically protected phenomena. Recently, the search for QSLs has expanded into the three-dimensional world, despite the suppression of quantum fluctuations due to high dimensionality. A new candidate material, K2Ni2(SO4)3, belongs to the langbeinite family and consists of two interconnected trillium lattices. Although magnetically ordered, it exhibits a highly dynamical and correlated state. In this work, we combine inelastic neutron scattering measurements with density functional theory (DFT), pseudo-fermion functional renormalization group (PFFRG), and classical Monte Carlo (cMC) calculations to study the magnetic properties of K2Ni2(SO4)3, revealing a high level of agreement between experiment and theory. We further reveal the origin of the dynamical state in K2Ni2(SO4)3 to be centred around a magnetic network composed of tetrahedra on a trillium lattice.

Generalized Kitaev spin liquid model and emergent twist defect

The Kitaev spin liquid model on a honeycomb lattice offers an intriguing feature that encapsulates both Abelian and non-Abelian phases (Alexei Kitaev, 2006). Recent studies suggest that the comprehensive phase diagram of the possible generalized Kitaev model largely depends on the specific details of the discrete lattice, which somewhat deviates from the traditional understanding of “topological” phases. In this paper, we propose an adapted version of the Kitaev spin liquid model on arbitrary planar lattices. Our revised model recovers the surface code model under certain parameter selections within the Hamiltonian terms. Changes in parameters can initiate the emergence of holes, domain walls, or twist defects. Notably, the twist defect, which presents as a lattice dislocation defect, exhibits non-Abelian braiding statistics upon tuning the coefficients of the Hamiltonian. Additionally, we illustrate that the creation, movement, and fusion of these defects can be accomplished through natural time evolution by linearly interpolating the static Hamiltonian. These defects demonstrate the Ising anyon fusion rule as anticipated. Our findings hint at possible implementation in actual physical materials owing to a more realistically achievable two-body interaction.

Amorphous and Polycrystalline Routes toward a Chiral Spin Liquid.

We show that a chiral spin liquid spontaneously emerges in partially amorphous, polycrystalline, or ion-irradiated Kitaev materials. In these systems, time-reversal symmetry is broken spontaneously due to a nonzero density of plaquettes with an odd number of edges n_{odd}. This mechanism opens a sizable gap, at small n_{odd} compatible with that of typical amorphous materials and polycrystals, and which can alternatively be induced by ion irradiation. We find that the gap is proportional to n_{odd}, saturating at n_{odd}∼40%. Using exact diagonalization, we find that the chiral spin liquid is approximately as stable to Heisenberg interactions as Kitaev's honeycomb spin-liquid model. Our results open up a significant number of noncrystalline systems where chiral spin liquids can emerge without external magnetic fields.

Single-shot quantum measurements sketch quantum many-body states

Quantum measurements allow us to observe quantum many-body systems consisting of a multitude of microscopic degrees of freedom. However, the quantum uncertainty and the exponentially large Hilbert space pose natural barriers to simple interpretations of the quantum measurement outcomes. We propose a nonlinear ``measurement energy'' based upon the measurement outcomes and a general approach akin to quantum machine learning to extract the most probable states (maximum likelihood estimates), naturally reconciling noncommuting observables and getting more out of the quantum measurements. Compatible with established quantum many-body ansatzes and efficient optimization, our strategy offers state-of-art capacity with control and full information. We showcase the versatility and accuracy of our strategy on random long-range fermion and Kitaev quantum spin-liquid models, where smoking-gun signatures were lacking.

Topological Defect Network Representations of Fracton Stabilizer Codes

A topological defect network (TDN) is formed by a network of topological defects embedded within a topological quantum field theory (TQFT). TDNs were introduced recently for the purpose of describing fracton topological phases of matter using the framework of defect TQFT. Their effectiveness has been demonstrated through numerous examples, yet a systematic construction was lacking. Here we solve this problem by formulating a method to construct TDNs for a wide range of lattice Hamiltonians. Our method takes a lattice Hamiltonian as input, applies an ungauging procedure, then creates a refined lattice within each unit cell, followed by regauging the system to produce a TDN as output. For topological Calderbank-Shor-Steane (CSS) Pauli stabilizer models, this procedure is guaranteed to produce a phase equivalent TDN. This provides TDN representations of canonical fracton models for which no such construction was previously known, including Haah's cubic code and Yoshida's infinite family of fractal spin liquid models. We demonstrate the applicability of our method beyond CSS stabilizer models by constructing TDNs for non-CSS models including Chamon's model and the semionic X-cube model.

Higher-dimensional Jordan-Wigner transformation and auxiliary Majorana fermions

We discuss a scheme for performing Jordan-Wigner transformation for various lattice fermion systems in two and three dimensions which keeps internal and spatial symmetries manifest. The correspondence between fermionic and bosonic operators is established with the help of auxiliary Majorana fermions. The current construction is applicable to general lattices with even coordination numbers and an arbitrary number of fermion flavors. The approach is demonstrated on the single-orbital square, triangular, and cubic lattices for spin-$1/2$ fermions. We also discuss the relation to some quantum spin liquid models.

Spin dynamics, entanglement, and the nature of the spin liquid state in YbZnGaO4

Electron spin dynamics was studied down to 80 mK in the triangular-lattice quantum spin-liquid candidate ${\mathrm{YbZnGaO}}_{4}$ using muon spin relaxation, finding no evidence for freezing or ordering of the Yb spins. The muon spin relaxation rate can be represented by the sum of two contributions, one dependent on longitudinal magnetic field and the other independent of field. The field-dependent term follows the form expected for two-dimensional diffusion of mobile spin excitations. The spin-diffusion rate obtained for these excitations in the high temperature paramagnetic regime is comparable with the exchange coupling frequency $J/h$, reducing significantly in the low temperature quantum regime. This slowdown is assigned to the effect of quantum entanglement. The exchange coupling $J$ is estimated to be 2.0(2) K from the crossover between the two regimes. The field-independent term is only weakly dependent on temperature, and at 15 K its absolute value is consistent with dipolar coupling of the muon to the three Yb moments closest to the muon site, where the spin dynamics of these moments is determined by exchange fluctuations. The temperature-dependent properties in the quantum regime are compared against the three possible U(1) spin-liquid models that have been obtained for the strongly spin-orbit coupled triangular lattice by Y.-D. Li, Y.-M. Lu, and G. Chen [Phys. Rev. B 96, 054445 (2017)]. The comparison with theory takes published specific heat and thermal conductivity data into account, along with the spin-diffusion rate obtained from the muons. It is found that the nodal spin-liquid model U1A11 containing both linear and quadratic nodes provides better agreement with experiment than either the U1A00 spinon Fermi surface (FS) model or the U1A01 model that contains only linear nodes.

Measuring Rényi entanglement entropy with high efficiency and precision in quantum Monte Carlo simulations

We develop a nonequilibrium increment method in quantum Monte Carlo simulations to obtain the Rényi entanglement entropy of various quantum many-body systems with high efficiency and precision. To demonstrate its power, we show the results on a few important yet difficult (2 + 1)d quantum lattice models, ranging from the Heisenberg quantum antiferromagnet with spontaneous symmetry breaking, the quantum critical point with O(3) conformal field theory (CFT) to the toric code {{mathbb{Z}}}_{2} topological ordered state and the Kagome {{mathbb{Z}}}_{2} quantum spin liquid model with frustration and multi-spin interactions. In all these cases, our method either reveals the precise CFT data from the logarithmic correction or extracts the quantum dimension in topological order, from the dominant area law in finite-size scaling, with very large system sizes, controlled errorbars, and minimal computational costs. Our method, therefore, establishes a controlled and practical computation paradigm to obtain the difficult yet important universal properties in highly entangled quantum matter.

Floquet engineering of Kitaev quantum magnets

In recent years, there has been an intense search for materials realizing the Kitaev quantum spin liquid model. A number of edge-shared compounds with strong spin-orbit coupling, such as RuCl3 and iridates, have been proposed to realize this model. Nevertheless, an effective spin Hamiltonian derived from the microscopic model relevant to these compounds generally contains terms that are antagonistic toward the quantum spin liquid. This is consistent with the fact that the zero magnetic field ground state of these materials is generally magnetically ordered. It is a pressing issue to identify protocols to drive the system to the limit of the Kitaev quantum spin model. In this work, we propose Floquet engineering of these Kitaev quantum magnets by coupling materials to a circularly polarized laser. We demonstrate that all the magnetic interactions can be tuned in situ by the amplitude and frequency of the laser, hence providing a route to stabilize the Kitaev quantum spin liquid phase.

Low-energy structure of spiral spin liquids

In this paper we identify a previously unexplored type of topological defect in spiral spin liquids---the momentum vortex---and reveal its dominant role in shaping the low-energy physics of such systems. Spiral spin liquids are a class of classical spin liquids featuring subextensively degenerate ground states. They are distinct from spin liquids on geometrically frustrated lattices, in which the ground-state degeneracy is extensive and connected by local spin flips. Despite a handful of experimental realizations and many theoretical studies, a concrete physical picture of their spin liquidity has not been established so far. In this paper, we study a 2D spiral spin liquid model to answer this question. We find that the local momentum vector field can carry topological defects in the form of vortices, which, however, have very different properties from the commonly known spin vortices. The fluctuations of such vortices lead the system into a liquid phase at intermediate temperatures. Furthermore, the effective low-energy theory of such vortices indicates their equivalence to quadrupoles of fractons in a rank-2 U(1) gauge theory or, alternatively, to quadrupoles of disclinations in elasticity theory. At very low temperatures, the system freezes into a glassy state in which these vortices form a rigid network with straight-line domain walls. Our paper sheds light on the nature of spiral spin liquids and also paves the way toward understanding their quantum limit.

Bifurcating subsystem symmetric entanglement renormalization in two dimensions

We introduce the subsystem symmetry-preserving real-space entanglement renormalization group and apply it to study bifurcating flows generated by linear and fractal subsystem symmetry-protected topological phases in two spatial dimensions. We classify all bifurcating fixed points that are given by subsystem symmetric cluster states with two qubits per unit cell. In particular, we find that the square lattice cluster state is a quotient-bifurcating fixed point, while the cluster states derived from Yoshida's first order fractal spin liquid models are self-bifurcating fixed points. We discuss the relevance of bifurcating subsystem symmetry-preserving renormalization group fixed points for the classification and equivalence of subsystem symmetry-protected topological phases.

Vestigial anyon condensation in kagome quantum spin liquids

We construct a lattice model of topological order (kagome quantum spin liquids) and solve it with unbiased quantum Monte Carlo simulations. A three-stage anyon condensation with two transitions from a $\mathbb Z_2\boxtimes\mathbb Z_2$ topological order to a $\mathbb Z_2$ topological order and eventually to a trivial symmetric phase is revealed. These results provide concrete examples of phase transitions between topological orders in quantum magnets. The designed quantum spin liquid model and its numerical solution offer a playground for further investigations on vestigial anyon condensation.

Heisenberg-Kitaev model in a magnetic field: 1/S expansion

The exact solution of Kitaev's spin-$1/2$ honeycomb spin-liquid model has sparked an intense search for Mott insulators hosting bond-dependent Kitaev interactions, of which $\mathrm{Na}_{2}\mathrm{IrO}_{3}$ and $\alpha-\mathrm{RuCl}_{3}$ are prime examples. Subsequently, it has been proposed that also spin-$1$ and spin-$3/2$ analogs of Kitaev interactions may occur in materials with strong spin-orbit coupling. As a minimal model to describe these Kitaev materials, we study the Heisenberg-Kitaev Hamiltonian in a consistent $1/S$ expansion, with $S$ being the spin size. We present a comprehensive study of this model in the presence of an external magnetic field applied along two different directions, [001] and [111], for which an intricate classical phase diagram has been reported. In both settings, we employ spin-wave theory in a number of ordered phases to compute phase boundaries at the next-to-leading order in $1/S$ and show that quantum corrections substantially modify the classical phase diagram. More broadly, our work presents a consistent route to investigate the leading quantum corrections in spin models that break spin-rotational symmetry.

Probing quantum spin liquids in equilibrium using the inverse spin Hall effect

We propose an experimental method utilizing a strongly spin-orbit coupled metal to quantum magnet bilayer that will probe quantum magnets lacking long range magnetic order, e.g., quantum spin liquids, via examination of the voltage noise spectrum in the metal layer. The bilayer is held in thermal and chemical equilibrium, and spin fluctuations arising across the single interface are converted into voltage fluctuations in the metal as a result of the inverse spin Hall effect. We elucidate the theoretical workings of the proposed bilayer system, and provide precise predictions for the frequency characteristics of the enhancement to the ac electrical resistance measured in the metal layer for three candidate quantum spin liquid models. Application to the Heisenberg spin-$1/2$ kagom{\'e} lattice model should allow for the extraction of any spinon gap present. A quantum spin liquid consisting of fermionic spinons coupled to a $U(1)$ gauge field should cause subdominant $\W^{4/3}$ scaling of the resistance of the coupled metal. Finally, if the magnet is well-captured by the Kitaev model in the gapless spin liquid phase, then the proposed bilayer can extract the two-flux gap which arises in spite of the gapless spectrum of the fermions. We therefore show that spectral analysis of the ac resistance in the metal in a single interface, equilibrium bilayer can test the relevance of a quantum spin liquid model to a given candidate material.

Vanishing Wilson ratio as the hallmark of quantum spin-liquid models

We present numerical results for finite-temperature $T>0$ thermodynamic quantities, entropy $s(T)$, uniform susceptibility $\chi_0(T)$ and the Wilson ratio $R(T)$, for several isotropic $S=1/2$ extended Heisenberg models which are prototype models for planar quantum spin liquids. We consider in this context the frustrated $J_1$-$J_2$ model on kagome, triangular, and square lattice, as well as the Heisenberg model on triangular lattice with the ring exchange. Our analysis reveals that typically in the spin-liquid parameter regimes the low-temperature $s(T)$ remains considerable, while $\chi_0(T)$ is reduced consistent mostly with a triplet gap. This leads to vanishing $R(T \to 0)$, being the indication of macroscopic number of singlets lying below triplet excitations. This is in contrast to $J_1$-$J_2$ Heisenberg chain, where $R(T \to 0)$ either remains finite in the gapless regime, or the singlet and triplet gap are equal in the dimerized regime.

Impact of disorder on dynamics and ordering in the honeycomb-lattice iridate Na2IrO3

Kitaev's honeycomb spin-liquid model and its proposed realization in materials such as $\alpha$-RuCl$_3$, Li$_2$IrO$_3$ and Na$_2$IrO$_3$ continue to present open questions about how the dynamics of a spin-liquid are modified in the presence of non-Kitaev interactions as well as the presence of inhomogeneities. Here we use $^{23}$Na nuclear magnetic resonance to probe both static and dynamical magnetic properties in single crystal Na$_2$IrO$_3$. We find that the NMR shift follows the bulk susceptibility above 30 K but deviates from it below; moreover below $T_N$ the spectra show a broad distribution of internal magnetic fields. Both of these results provide evidence for inequivalent magnetic sites at low temperature, suggesting inhomogeneities are important for the magnetism. The spin-lattice relaxation rate is isotropic and diverges at $T_N$, suggesting that the Kitaev cubic axes may control the critical quantum spin fluctuations. In the ordered state, we observe gapless excitations, which may arise from site substitution, emergent defects from milder disorder, or possibly be associated with nearby quantum paramagnetic states distinct from the Kitaev spin liquid.

Spin liquids from Majorana zero modes in a Cooper-pair box

We propose a path for constructing diverse interacting spin systems from topological nanowires in Cooper Boxes. The wires are grouped into a three-wire building block called an 'hexon', consisting of six Majorana zero modes. In the presence of a strong charging energy, the hexon becomes a Cooper box equivalent to two spin-$1/2$ degrees of freedom. By considering arrays of hexons and controlling the distances between the various wires, one can tune the Hamiltonian governing the low-energy spins, thus providing a route for controllably constructing interacting spin systems in one- and two-dimensions. We explicitly present realizations of the one-dimensional spin-$1/2$ XXZ chain, as well as the transverse field Ising model. We propose an experiment capable of revealing the nature of critical points in such effective spin systems by applying a local gate voltage and measuring the induced charge at a distance. To demonstrate the applicability of this approach to two-dimensions, we provide a scheme for realizing the topologically ordered Yao-Kivelson spin-liquid model, which has a collective Majorana edge mode, similar to the B-phase of Kitaev's honeycomb model.

Gapped SU(3) spin liquid with Z3 topological order

We construct a topological spin liquid (TSL) model on the kagome lattice, with SU(3) symmetry with the fundamental representation at each lattice site, based on Projected Entangled Pair States (PEPS). Using the PEPS framework, we can adiabatically connect the model to a fixed point model (analogous to the dimer model for Resonating Valence Bond states) which we prove to be locally equivalent to a $Z_3$ quantum double model. Numerical study of the interpolation reveals no sign of a phase transition or long-range order, characterizing the model conclusively as a gapped TSL. We further study the entanglement spectrum of the model and find that while it is gapped, it exhibits branches with vastly different velocities, with the slow branch matching the counting of a chiral $SU(3)_1$ CFT, suggesting that it can be deformed to a model with chiral $SU(3)_1$ entanglement spectrum.

Pinch point singularities of tensor spin liquids

Recently, a new class of three-dimensional spin liquid models have been theoretically discovered, which feature generalized Coulomb phases of emergent symmetric tensor $U(1)$ gauge theories. These ``higher rank'' tensor models are particularly intriguing due to the presence of quasiparticles with restricted mobility, such as fractons. We investigate universal experimental signatures of tensor Coulomb phases. Most notably, we show that tensor Coulomb spin liquids (both quantum and classical) feature characteristic pinch point singularities in their spin-spin correlation functions, accessible via neutron scattering, which can be readily distinguished from pinch points in conventional $U(1)$ spin liquids. These pinch points can thus serve as a crisp experimental diagnostic for such phases. We also tabulate the low-temperature heat capacity of various tensor Coulomb phases, which serves as a useful additional diagnostic in certain cases.

Edge magnetization and spin transport in an SU(2)-symmetric Kitaev spin liquid

We investigate the edge magnetism and the spin transport properties of an SU(2)-symmetric Kitaev spin liquid (KSL) model put forward by Yao and Lee [Phys. Rev. Lett. 107, 087205 (2011)] on the honeycomb lattice. In this model, the spin degrees of freedom fractionalize into a ${\mathbb{Z}}_{2}$ static gauge field and three species of either gapless (Dirac) or gapped (chiral) Majorana fermionic excitations. We find that, when a magnetic field is applied on a zigzag edge, the Dirac KSL exhibits a nonlocal magnetization associated with the existence of zero-energy edge modes. The application of a spin bias $V={\ensuremath{\mu}}_{\ensuremath{\uparrow}}\ensuremath{-}{\ensuremath{\mu}}_{\ensuremath{\downarrow}}$ at the interface of the spin system with a normal metal produces a spin current into the KSL, which depends as a power law on $V$, in the zero-temperature limit, for both Dirac and chiral KSLs, but with different exponents. Lastly, we study the longitudinal spin Seebeck effect, in which a spin current is driven by the combined action of a magnetic field perpendicular to the plane of the honeycomb lattice and a thermal gradient at the interface of the KSL with a metal. Our results suggest that edge magnetization and spin transport can be used to probe the existence of charge-neutral edge states in quantum spin liquids.