Search

Twisted bilayer graphene is an excellent example of highly correlated system demonstrating a nearly flat electron band, the Mott transition and probably a spin liquid state. Besides the one-electron picture, analysis of Dirac points is performed in terms of spinon Fermi surface in the limit of strong correlations. Application of gauge field theory to describe deconfined spin liquid phase is treated. Topological quantum transitions, including those from small to large Fermi surface in the presence of van Hove singularities, are discussed., Comment: 4 pages, JETP Letters, accepted

We propose an exactly solvable quantum spin-1/2 model with time reversal invariance on a two dimensional brick-wall lattice, where each unit cell consists of three sites. We find that the ground states are algebraic quantum spin liquid states. The spinon excitations are gapless and the energy dispersion is linear around two Dirac points. The ground states are of three-fold topological degeneracy on a torus. Breaking the time reversal symmetry opens a bulk energy gap and the $Z_2$ vortices obey non-Abelian statistics., Comment: 5+2 pages, 4 figures

The interplay of symmetry and topology has been at the forefront of recent progress in quantum matter. Here we uncover an unexpected connection between band topology and the description of competing orders in a quantum magnet. Specifically we show that aspects of band topology protected by crystalline symmetries determine key properties of the Dirac spin liquid (DSL) which can be defined on the honeycomb, square, triangular and kagom\'e lattices. At low energies, the DSL on all these lattices is described by an emergent Quantum Electrodynamics (QED$_3$) with $N_f=4$ flavors of Dirac fermions coupled to a $U(1)$ gauge field. However the symmetry properties of the magnetic monopoles, an important class of critical degrees of freedom, behave very differently on different lattices. In particular, we show that the lattice momentum and angular momentum of monopoles can be determined from the charge (or Wannier) centers of the corresponding spinon insulator. We also show that for DSLs on bipartite lattices, there always exists a monopole that transforms trivially under all microscopic symmetries owing to the existence of a parent SU(2) gauge theory. We connect our results to generalized Lieb-Schultz-Mattis theorems and also derive the time-reversal and reflection properties of monopoles. Our results indicate that recent insights into free fermion band topology can also guide the description of strongly correlated quantum matter., Comment: 26 pages, 7 figures

We consider the Abelian-Higgs model in 2+1 dimensions with instanton-monopole defects. This model is closely related to the phases of quantum anti-ferromagnets. In the presence of $\mathbb{Z}_2$ preserving monopole operators, there are two confining ground states in the monopole phase, corresponding to the Valence Bond Solid (VBS) phase of quantum magnets. We show that the domain-wall carries a 't Hooft anomaly in this case. The anomaly can be saturated by, e.g., charge-conjugation breaking on the wall or by the domain wall theory becoming gapless (a gapless model that saturates the anomaly is $SU(2)_1$ WZW). Either way the fundamental scalar particles (i.e. spinons) which are confined in the bulk are deconfined on the domain-wall. This $\mathbb{Z}_2$ phase can be realized either with spin-1/2 on a rectangular lattice, or spin-1 on a square lattice. In both cases the domain wall contains spin-1/2 particles (which are absent in the bulk). We discuss the possible relation to recent lattice simulations of domain walls in VBS. We further generalize the discussion to Abrikosov-Nielsen-Olsen (ANO) vortices in a dual superconductor of the Abelian-Higgs model in 3+1 dimensions, and to the easy-plane limit of anti-ferromagnets. In the latter case the wall can undergo a variant of the BKT transition (consistent with the anomalies) while the bulk is still gapped. The same is true for the easy-axis limit of anti-ferromagnets. We also touch upon some analogies to Yang-Mills theory., Comment: 13 pages, 8 figures

Surprising properties of doped Mott insulators are at the heart of many quantum materials, including transition metal oxides and organic materials. The key to unraveling complex phenomena observed in these systems lies in understanding the interplay of spin and charge degrees of freedom. One of the most debated questions concerns the nature of charge carriers in a background of fluctuating spins. To shed new light on this problem, we suggest a simplified model with mixed dimensionality, where holes move through a Mott insulator unidirectionally while spin exchange interactions are two dimensional. By studying individual holes in this system, we find direct evidence for the formation of mesonic bound states of holons and spinons, connected by a string of displaced spins -- a precursor of the spin-charge separation obtained in the 1D limit of the model. Our predictions can be tested using ultracold atoms in a quantum gas microscope, allowing to directly image spinons and holons, and reveal the short-range hidden string order which we predict in this model., Comment: 7 pages, 6 figures, 4 pages methods section

Resonant inelastic x-ray scattering (RIXS) at the oxygen $K$-edge has recently accessed multi-spinon excitations in the one-dimensional antiferromagnet (1D-AFM) \sco, where four-spinon excitations are resolved separately from the two-spinon continuum. This technique, therefore, provides new opportunities to study fractionalized quasiparticle excitations in doped 1D-AFMs. To this end, we carried out exact diagonalization studies of the doped $t$-$J$ model and provided predictions for oxygen $K$-edge RIXS experiments on doped 1D-AFMs. We show that the RIXS spectra are rich, containing distinct two- and four-spinon excitations, dispersive (anti)holon excitations, and combinations thereof. Our results highlight how RIXS complements inelastic neutron scattering experiments by accessing additional charge and spin components of fractionalized quasiparticles., Comment: 5 pages, 4 figures

The Kitaev model with an applied magnetic field in the [Formula: see text] direction shows two transitions: from a nonabelian gapped quantum spin liquid (QSL) to a gapless QSL at [Formula: see text] and a second transition at a higher field [Formula: see text] to a gapped partially polarized phase, where [Formula: see text] is the strength of the Kitaev exchange interaction. We identify the intermediate phase to be a gapless U(1) QSL and determine the spin structure function [Formula: see text] and the Fermi surface [Formula: see text] of the gapless spinons using the density matrix renormalization group (DMRG) method for large honeycomb clusters. Further calculations of static spin-spin correlations, magnetization, spin susceptibility, and finite temperature-specific heat and entropy corroborate the gapped and gapless nature of the different field-dependent phases. In the intermediate phase, the spin-spin correlations decay as a power law with distance, indicative of a gapless phase.

We analyze the crossover from Kondo to weak-link regime by means of a model of tunable bond impurities in the middle of a spin-1/2 XXZ Heisenberg chain. We study the Kondo screening cloud and estimate the Kondo length by combining perturbative renormalization group approach with the exact numerical calculation of the integrated real-space spin-spin correlation functions. We show that, when the spin impurity is symmetrically coupled to the two parts of the chain with realistic values of the Kondo coupling strengths and spin-parity symmetry is preserved, the Kondo length takes values within the reach of nowadays experimental technology in ultracold-atom setups. In the case of non-symmetric Kondo couplings and/or spin parity broken by a nonzero magnetic field applied to the impurity, we discuss how Kondo screening redistributes among the chain as a function of the asymmetry in the couplings and map out the shrinking of the Kondo length when the magnetic field induces a crossover from Kondo impurity to weak-link physics., 27 pages, 12 .eps figures, version published in PRB

Different scenarios of the implementation of the two-band model in strongly correlated electrons systems, including frustrated magnets, high-temperature superconductors, and Kondo lattices, are considered. The interaction of current carriers with magnetic moments in the representations of pseudofermions or Schwinger bosons describing the spinon excitations is studied on the basis of the derived Hamiltonians of the $s-d$ exchange and $t-J$ models within the formalism of many-electron Hubbard operators., Comment: 5 pages

It has recently been found that bosonic excitations of ordered media, such as phonons or spinons, can exhibit topologically nontrivial band structures. Of particular interest are magnon and triplon excitations in quantum magnets, as they can easily be manipulated by an applied field. Here we study triplon excitations in an S=1/2 quantum spin ladder and show that they exhibit nontrivial topology, even in the quantum-disordered paramagnetic phase. Our analysis reveals that the paramagnetic phase actually consists of two separate regions with topologically distinct triplon excitations. We demonstrate that the topological transition between these two regions can be tuned by an external magnetic field. The winding number that characterizes the topology of the triplons is derived and evaluated. By the bulk-boundary correspondence, we find that the non-zero winding number implies the presence of localized triplon end states. Experimental signatures and possible physical realizations of the topological paramagnetic phase are discussed., 6+4 pages; References, footnotes, small clarification added in conclusions and suppl. mat (v2); Minor modifications, close to published version (v3)

This manuscript deals with the study of low dimensional cuprates, namely, the spin chains systems SrCuO₂ and Sr₂CuO₃. These two compounds exhibit highly anisotropic thermal conduction properties along the spin-chains direction, where magnetic thermal conduction contributes to the heat transport process via spinon excitations. Our study aims to highlight the mechanisms that govern the heat transport properties, particularly through the study of the scattering channels involving spinon, phonon and defects. The spinon (phonon)–defect scattering was probed by the intentional introduction of nonmagnetic dopants (1-2%) on the copper site, by: Mg²⁺, Zn²⁺, Pd²⁺ or Ni²⁺, or by the introduction of elements carrying different oxidation level on the strontium site, by: La³⁺ or K⁺. Single crystals of the pure and doped materials have been grown by the travelling solvent floating zone method. The structural, magnetic and thermal characterizations of the pure and doped compounds were performed. The magnetic excitation spectra of the compounds were determined by inelastic neutron scattering, NMR spectroscopy, and angle resolved photoemission spectroscopy to reveal the impact of the substitution on the spin dynamics of the doped compounds. The study of phonon spectra has also been performed by inelastic neutron scattering. Results from inelastic neutron scattering have been correlated to the heat transport properties of the pristine and substituted materials., Ce manuscrit porte sur l’étude de cuprates de basse dimensionnalité, les systèmes à chaînes de spins SrCuO₂ et Sr₂CuO₃. Un des intérêts de ces deux composés est qu’ils présentent des conductions thermiques fortement anisotropes. Celles-ci comportent une contribution magnétique due au transport de la chaleur via les excitations de spinons qui se manifeste uniquement dans la direction des chaînes de spins. Notre étude a pour objectif la mise en évidence des mécanismes qui gouvernent ces propriétés de transport, notamment à travers l’étude des interactions entre les spinons, les phonons et les défauts. Les interactions spinons (phonons)-défauts ont été sondées par l’introduction intentionnelle de dopants (1-2%) non-magnétiques sur le site du cuivre : Mg²⁺, Zn²⁺, Pd²⁺ ou Ni²⁺, ou encore par l’introduction d’éléments possédant des degrés d’oxydation différents sur le site du strontium : La³⁺ ou K⁺. Les composés ont été synthétisés sous leur forme monocristalline par la méthode de fusion de la zone solvante. Des caractérisations structurales, magnétiques et thermiques des composés purs et dopés ont été réalisées. Les spectres d’excitations magnétiques de ces cuprates ont été déterminés par diffusion inélastique de neutrons, spectroscopie RMN et spectroscopie de photoémission résolue en angle afin de révéler l’impact de la substitution. L’étude des spectres de phonons a également été réalisée par diffusion inélastique de neutrons. Les résultats de ces mesures sont corrélés aux propriétés de conduction thermique des composés purs et dopés..

We study the dynamic spin structure factor of the spin-$1/2$ square-lattice Heisenberg antiferromagnet and of the $J$-$Q$ model (with 4-spin interactions $Q$ and Heisenberg exchange $J$). Using an improved method for stochastic analytic continuation of imaginary-time correlation functions computed with QMC simulations, we can treat the sharp ($\delta$-function) contribution from spinwaves (magnons) and a continuum at higher energy. The results for the Heisenberg model agree with neutron scattering experiments on Cu(DCOO)$_2$$\cdot$4D$_2$O, where a broad spectral-weight continuum at $q=(\pi,0)$ was interpreted as deconfined spinons. Our results at $(\pi,0)$ show a similar reduction of the magnon weight and a large continuum, while the continuum is much smaller at $q=(\pi/2,\pi/2)$ (as also seen experimentally). Turning on $Q$, we observe a rapid reduction of the $(\pi,0)$ magnon weight to zero, well before the deconfined quantum phase transition into a spontaneously dimerized state. We re-interpret the picture of deconfined spinons at $(\pi,0)$ in the experiments as nearly deconfined spinons---a precursor to deconfined quantum criticality. To further elucidate the picture of a fragile $(\pi,0)$-magnon in the Heisenberg model and its depletion in the $J$-$Q$ model, we introduce an effective model in which a magnon can split into two spinons that do not separate but fluctuate in and out of the magnon space (in analogy with the resonance between a photon and a particle-hole pair in the exciton-polariton problem). The model reproduces the $(\pi,0)$ and $(\pi/2,\pi/2)$ features of the Heisenberg model. It can also account for the rapid loss of the $(\pi,0)$ magnon with increasing $Q$ and a remarkable persistence of a large magnon pole at $q=(\pi/2,\pi/2)$ even at the deconfined critical point., Comment: 25 pages, 24 figures, some additional discussion in version 3, to be appeared in PRX

K3Cu3AlO2(SO4)4 is a highly one-dimensional spin-1/2 inequilateral diamond-chain antiferromagnet. Spinon continuum and spin-singlet dimer excitations are observed in the inelastic neutron scattering spectra, which is in excellent agreement with a theoretical prediction: a dimer-monomer composite structure, where the dimer is caused by strong antiferromagnetic (AFM) coupling and the monomer forms an almost isolated quantum AFM chain controlling low-energy excitations. Moreover, muon spin rotation/relaxation spectroscopy shows no long-range ordering down to 90 mK, which is roughly three orders of magnitude lower than the exchange interaction of the quantum AFM chain. K3Cu3AlO2(SO4)4 is, thus, regarded as a compound that exhibits a Tomonaga-Luttinger spin liquid behavior at low temperatures close to the ground state.

We analyze recent heat capacity measurements in herbertsmithite $\rm ZnCu_3(OH)_6Cl_2$ single crystal samples subjected to strong magnetic fields. We show that the temperature dependence of specific heat $C_{mag}$ formed by quantum spin liquid at different magnetic fields $B$ resembles the electronic heat capacity $C_{el}$ of the HF metal $\rm YbRh_2Si_2$. We demonstrate that the spinon effective mass $M^*_{mag}\propto C_{mag}/T$ exhibits a scaling behavior like that of $C_{el}/T$. We also show that the recent measurements of $C_{mag}$ are compatible with those obtained on powder samples. These observations allow us to conclude that $\rm ZnCu_3(OH)_6Cl_2$ holds a stable strongly correlated quantum spin liquid, and a possible gap in the spectra of spinon excitations is absent even under the application of very high magnetic fields of 18 T., 5 pages, 5 figures

Two-dimensional density-matrix renormalization group method is employed to examine the ground state phase diagram of the Hubbard model on the triangular lattice at half filling. The calculation reveals two discontinuities in the double occupancy with increasing the repulsive Hubbard interaction U at Uc1 = 7.55 t and Uc2 = 9.65 t (t being the hopping integral), indicating that there are three phases separated by first order transitions. The absence of any singularity in physical quantities for 0 < U < Uc1 implies that this phase corresponds to a metallic phase. The local spin density induced by an applied pinning magnetic field for U > Uc2 exhibits a three sublattice feature, which is compatible with the Neel ordered state realized in the strong coupling limit. For Uc1 < U < Uc2, a response to the applied pinning magnetic field is comparable to that in the metallic phase but a relatively large spin correlation length is found with neither valence bond nor chiral magnetic order, suggesting a paramagnetic nature which resembles gapless spin liquid. The calculation also finds that the pair- ing correlation function monotonically decreases with increasing U and thus the superconductivity is unlikely in the intermediate phase., 11 pages, 13 figures

We use Schwinger-boson mean-field theory (SBMFT) to study the ground state of the spin-$S$ triangular-lattice Heisenberg model with nearest (${J}_{1}$) and next-nearest (${J}_{2}$) neighbor antiferromagnetic interactions. Previous work on the $S=1/2$ model leads us to consider two spin-liquid Ans\"atze, one symmetric and one nematic, which upon spinon condensation give magnetically ordered states with ${120}^{\ensuremath{\circ}}$ order and collinear stripe order, respectively. The SBMFT contains the parameter $\ensuremath{\kappa}$, the expectation value of the number of bosons per site, which in the exact theory equals $2S$. For $\ensuremath{\kappa}=1$ there is a direct, first-order transition between the ordered states as ${J}_{2}/{J}_{1}$ increases. Motivated by arguments that in SBMFT, smaller $\ensuremath{\kappa}$ may be more appropriate for describing the $S=1/2$ case qualitatively, we find that in a $\ensuremath{\kappa}$ window around 0.6, a region with the (gapped ${Z}_{2}$) symmetric spin liquid opens up between the ordered states. As a consequence, the static structure factor has the same peak locations in the spin liquid as in the ${120}^{\ensuremath{\circ}}$ ordered state, and the phase transitions into the ${120}^{\ensuremath{\circ}}$ and collinear stripe-ordered states are continuous and first order, respectively.

We describe a theoretical approach for finding spontaneously symmetry-broken electronic phases due to strong electronic interactions when using recently developed slave-particle (slave-boson) approaches based on occupation numbers. We describe why, to date, spontaneous symmetry breaking has proven difficult to achieve in such approaches. We then provide a total-energy based approach for introducing auxiliary symmetry breaking fields into the solution of the slave-particle problem that leads to lowered total energies for symmetry broken phases. We point out that not all slave-particle approaches yield to energy lowering: the slave-particle model being used must explicitly describe the degrees of freedom that break symmetry. Finally, our total energy approach permits us to greatly simplify the formalism used to achieve a self-consistent solution between spinon and slave modes while increasing numerical stability and greatly speeding up the calculations.

Motivated by the recent developments on cluster Mott insulating materials such as the cluster magnet LiZn$_2$Mo$_3$O$_8$, we consider the strong plaquette charge ordered regime of the extended Hubbard model on a breathing Kagome lattice and reveal the properties of the cluster Mottness. The plaquette charge order arises from the inter-site charge interaction and the collective motion of three localized electrons on the hexagon plaquettes. This model leads naturally to a reduction of the local moments by 2/3 as observed in LiZn$_2$Mo$_3$O$_8$. Furthermore, at low temperatures each hexagon plaquette contains an extra orbital-like degree of freedom in addition to the remaining spin 1/2. We explore the consequence of this emergent orbital degree of freedom. We point out the interaction between the local moments is naturally described by a Kugel-Khomskii spin-orbital model. We develop a parton approach and suggest a spin liquid ground state with spinon Fermi surfaces for this model. We further predict an emergent orbital order when the system is under a strong magnetic field. Various experimental consequences for LiZn$_2$Mo$_3$O$_8$ are discussed, including an argument that the charge ordering much be short ranged if the charge per Mo is slightly off stoichiometry., 12 pages, 13 figures

The cluster slave-spin method is introduced to systematically investigate the solutions of the Hubbard model including the symmetry-broken phases. In this method, the electron operator is factorized into a fermionic spinon describing the physical spin and a slave-spin describing the charge fluctuations. Following the $U(1)$ formalism derived by Yu and Si [Phys. Rev. B 86, 085104 (2012)], it is shown that the self-consistent equations to explore various symmetry-broken density wave states can be constructed in general with a cluster of multiple slave-spin sites. We employ this method to study the antiferromagnetic (AFM) state in the single band Hubbard model with the two- and four-site clusters of slave spins. While the Hubbard gap, the charge gap due to the doubly occupied states, scales with the Hubbard interaction $U$ as expected, the AFM gap $\mathrm{\ensuremath{\Delta}}$, the gap in the spinon dispersion in the AFM state, exhibits a crossover from the weak- to strong-coupling behaviors as $U$ increases. Our cluster slave-spin method reproduces not only the traditional mean-field behavior of $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}U$ in the weak-coupling limit, but also the behavior of $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}{t}^{2}/U$ predicted by the superexchange mechanism in the strong-coupling limit. In addition, the holon-doublon correlator as functions of $U$ and doping $x$ is also computed, which exhibits a strong tendency toward the holon-doublon binding in the strong coupling regime. We further show that the quasiparticle weight obtained by the cluster slave-spin method is in a good agreement with the generalized Gutzwiller approximation in both AFM and paramagnetic states, and the results can be improved beyond the generalized Gutzwiller approximation as the cluster is enlarged from a single site to four sites. Our results demonstrate that the cluster slave-spin method can be a powerful tool to systematically investigate the strongly correlated system.

We propose that resonant inelastic X-ray scattering (RIXS) is an effective probe of the fractionalized excitations in three-dimensional (3D) Kitaev spin liquids. While the non-spin-conserving RIXS responses are dominated by the gauge-flux excitations and reproduce the inelastic-neutron-scattering response, the spin-conserving (SC) RIXS response picks up the Majorana-fermion excitations and detects whether they are gapless at Weyl points, nodal lines, or Fermi surfaces. As a signature of symmetry fractionalization, the SC RIXS response is suppressed around the $\Gamma$ point. On a technical level, we calculate the exact SC RIXS responses of the Kitaev models on the hyperhoneycomb, stripyhoneycomb, hyperhexagon, and hyperoctagon lattices, arguing that our main results also apply to generic 3D Kitaev spin liquids beyond these exactly solvable models., Comment: 5+5 pages, 3+4 figures

Motivated by the recent progress in the spin-orbit-coupled triangular lattice spin liquid candidate ${\mathrm{YbMgGaO}}_{4}$, we carry out a systematic projective symmetry group analysis and mean-field study of candidate $U(1)$ spin-liquid ground states. Due to the spin-orbital entanglement of the Yb moments, the space-group symmetry operation transforms both the position and the orientation of the local moments, and hence it brings different features for the projective realization of the lattice symmetries from the cases with spin-only moments. Among the eight $U(1)$ spin liquids that we find with the fermionic parton construction, only one spin-liquid state, which was proposed and analyzed by Yao Shen et al. [Nature (London) 540, 559 (2016)] and labeled as U1A00 in the present work, stands out and gives a large spinon Fermi surface and provides a consistent explanation for the spectroscopic results in ${\mathrm{YbMgGaO}}_{4}$. Further connection of this spinon Fermi surface $U(1)$ spin liquid with ${\mathrm{YbMgGaO}}_{4}$ and the future directions are discussed. Finally, our results may apply to other spin-orbit-coupled triangular lattice spin-liquid candidates, and more broadly, our general approach can be well extended to spin-orbit-coupled spin-liquid candidate materials.

Motivated by the rapid experimental progress of quantum spin ice materials, we study the dynamical properties of pyrochlore spin ice in the U(1) spin liquid phases. In particular, we focus on the spinon excitations that appear at high energies and show up as an excitation continuum in the dynamic spin structure factor. The keen connection between the crystal symmetry fractionalization of the spinons and the spectral periodicity of the spinon continuum is emphasized and explicitly demonstrated. When the spinon experiences a background $\ensuremath{\pi}$ flux and the spinon continuum exhibits an enhanced spectral periodicity with a folded Brillouin zone, this spectral property can then be used to detect the spin quantum number fractionalization and U(1) spin liquid. Our prediction can be immediately examined by inelastic neutron-scattering experiments among quantum spin ice materials with Kramers' doublets. Further application to the non-Kramers' doublets is discussed.

Quantum fluctuations of the N\'eel state of the square lattice antiferromagnet are usually described by a $\mathbb{CP}^1$ theory of bosonic spinons coupled to a U(1) gauge field, and with a global SU(2) spin rotation symmetry. Such a theory also has a confining phase with valence bond solid (VBS) order, and upon including spin-singlet charge 2 Higgs fields, deconfined phases with $\mathbb{Z}_2$ topological order possibly intertwined with discrete broken global symmetries. We present dual theories of the same phases starting from a mean-field theory of fermionic spinons moving in $\pi$-flux in each square lattice plaquette. Fluctuations about this $\pi$-flux state are described by 2+1 dimensional quantum chromodynamics (QCD$_3$) with a SU(2) gauge group and $N_f=2$ flavors of massless Dirac fermions. It has recently been argued by Wang et al. (arXiv:1703.02426) that this QCD$_3$ theory describes the N\'eel-VBS quantum phase transition. We introduce adjoint Higgs fields in QCD$_3$, and obtain fermionic dual descriptions of the phases with $\mathbb{Z}_2$ topological order obtained earlier using the bosonic $\mathbb{CP}^1$ theory. We also present a fermionic spinon derivation of the monopole Berry phases in the U(1) gauge theory of the VBS state. The global phase diagram of these phases contains multi-critical points, and our results imply new boson-fermion dualities between critical gauge theories of these points., Comment: Version 2: 32 pages, 3 figures, 12 tables; fixed typos, merged figures

In spin lattices, frustration happens when there exists a competition between different terms of the Hamiltonian, which cannot be simultaneously minimized in all local physical bases. Frustration can destabilize the bipartite Neel order, which is the conventional groundstate of antiferromagnetic Heisenberg-type models in two dimensions. Hence, the focus of two-dimensional quantum magnetism studies has been directed toward elucidating the zero-temperature properties of Heisenberg-type models on the Archimedean lattices that exhibit geometrical frustration. This type of frustration is of great interest since it can lead to exotic forms of the quantum matter, such as topological spin-liquid (SL) and many-sublattice long-range-ordered phases. The spin-1/2 triangular Heisenberg model (THM) is a prototypical model with geometrical frustration, which has the highest coordination number of all Archimedean lattices. However, there exist two major, distinct difficulties for numerical methods to efficiently determine the groundstate properties of the THM: the immense dimension of the Hilbert space and the geometrical frustration phenomenon. In this thesis, to deal with such difficulties, we employ the state-of-the- art non-Abelian matrix product states (MPS) and density matrix renormalization group (DMRG) methods. We exploit the symmetries of the Hamiltonian to reduce the working dimension of the Hilbert space and, to lessen the effects of the geometrical frustration, we depend on the local entanglement (area-law) nature of the MPS as it scales with the size of a lattice subspace with reduced spatial dimensions. In such a way, we thoroughly examine the phase diagram of the THM with nearest and next-nearest neighbor (NN and NNN) exchange couplings on finite- and infinite-length cylinders of widths up to 12 sites. In addition, we develop new numerical tools to analyze topological and symmetry-broken phases in the context of the SU(2)-symmetric translation-invariant MPS algorithm. We establish that on infinite cylinders, the anyonic sectors of a topological order can be classified through the fractionalization of global symmetries and degeneracy patterns of the entanglement spectrum (ES). On the other hand, we detect symmetry-breaking phase transitions by a combination of the correlation lengths and second and fourth cumulants of the magnetic order parameters, even though symmetry implies that the order parameter itself is strictly zero. Furthermore, the appearance of Nambu-Goldstone modes in the excitation spectrum can be detected using “tower of states” (TOS) levels in the momentum-resolved ES. By applying the new tools to the THM, we discover a variety of exotic groundstates such as a Z2-gauge toric-code-type topological order (with vison and spinon excitations), a NNN Majumdar-Gosh, and a long-range 120◦-ordered state.

Metal organic magnets have enormous potential to host a variety of electronic and magnetic phases that originate from a strong interplay between the spin, orbital, and lattice degrees of freedom. We control this interplay in the quantum magnet ${\mathrm{CuF}}_{2}$(${\mathrm{D}}_{2}\mathrm{O}$)${}_{2}(\mathrm{pyz})$ by using high pressure to drive the system through structural and magnetic phase transitions. Using neutron scattering, we show that the low pressure state, which hosts a two-dimensional square lattice with spin-wave excitations and a dominant exchange coupling of 0.89 meV, transforms at high pressure into a one-dimensional spin chain hallmarked by a spinon continuum and a reduced exchange interaction of 0.43 meV. This direct microscopic observation of a magnetic dimensional crossover as a function of pressure opens up new possibilities for studying the evolution of fractionalised excitations in low-dimensional quantum magnets and eventually pressure-controlled metal--insulator transitions.

Integrable models support elementary excitations with infinite lifetimes. In the spin-1/2 Heisenberg chain these are known as spinons. We consider the stability of spinons when a weak integrability breaking perturbation is added to the Heisenberg chain in a magnetic field. We focus on the case where the perturbation is a next nearest neighbour exchange interaction. We calculate the spinon decay rate in leading order in perturbation theory using methods of integrability and identify the dominant decay channels. The decay rate is found to be small, which indicates that spinons remain well-defined excitations even though integrability is broken., 29 pages, 8 figures

The Kagome-lattice-based material, Volborthite, $\mathrm{Cu_3 V_2 O_7 (OH)_2 \cdot 2 H_2 O}$, has been considered as a promising platform for discovery of unusual quantum ground states due to the frustrated nature of spin interactions. Here we explore possible quantum spin liquid and magnetically ordered phases in a two-dimensional non-symmorphic lattice described by $p2gg$ layer space group, which is consistent with the spatial anisotropy of the spin model derived from density functional theory (DFT) for Volborthite. Using the projective symmetry group (PSG) analysis and Schwinger boson mean field theory, we classify possible spin liquid phases with bosonic spinons and investigate magnetically ordered phases connected to such states. It is shown, in general, that only translationally invariant mean field states are allowed in two-dimensional non-symmorphic lattices, which simplifies the classification considerably. The mean field phase diagram of the DFT-derived spin model is studied and it is found that possible quantum spin liquid phases are connected to two types of magnetically ordered phases, a coplanar incommensurate $(q,0)$ spiral order as the ground state and a closely competing coplanar commensurate $(\pi,\pi)$ spin density wave order. In addition, periodicity enhancement of the two-spinon continuum, a signature of symmetry fractionalization, is found in the spin liquid phases connected to the $(\pi,\pi)$ spin density wave order. We discuss relevance of these results to recent and future experiments on Volborthite., Comment: 16 pages, 6 figures, 2 tables; typos corrected

We study the effect of nonmagnetic ${\mathrm{Zn}}^{2+}$ (spin-0) and magnetic ${\mathrm{Ni}}^{2+}$ (spin-1) impurities on the ground state and low-lying excitations of the quasi-one-dimensional spin-$1/2$ Heisenberg antiferromagnet ${\mathrm{Sr}}_{2}{\mathrm{CuO}}_{3}$ using inelastic neutron scattering, specific heat, and bulk magnetization measurements. We show that 1% ${\mathrm{Ni}}^{2+}$ doping in ${\mathrm{Sr}}_{2}{\mathrm{CuO}}_{3}$ results in a sizable spin gap in the spinon excitations, analogous to the case of Ni doped ${\mathrm{SrCuO}}_{2}$ previously reported [Simutis et al., Phys. Rev. Lett. 111, 067204 (2013)]. However, a similar level of ${\mathrm{Zn}}^{2+}$ doping in ${\mathrm{SrCuO}}_{2}$, investigated here for comparison, did not reveal any signs of a spin gap. Magnetic ordering temperature was found to be suppressed in the presence of both ${\mathrm{Zn}}^{2+}$ and ${\mathrm{Ni}}^{2+}$ impurities; however, the rate of suppression due to ${\mathrm{Ni}}^{2+}$ was found to be much more pronounced than for ${\mathrm{Zn}}^{2+}$. Effect of magnetic field on the ordering temperature is investigated. We found that with increasing magnetic field, not only the magnetic ordering temperature gradually increases but the size of specific heat anomaly associated with the magnetic ordering also progressively enhances, which can be qualitatively understood as due to the field induced suppression of quantum fluctuations.

We study the U(1) quantum spin liquid on the pyrochlore spin ice systems. For the non-Kramers doublets such as Pr$^{3+}$ and Tb$^{3+}$, we point out that the inelastic neutron scattering not only detects the low-energy gauge photon, but also contains the continuum of the "magnetic monopole" excitations. Unlike the spinons, these "magnetic monopoles" are purely of quantum origin and have no classical analogue. We further point out that the "magnetic monopole" experiences a background dual "$\pi$" flux due to the spin-1/2 nature of the local moment when the "monopole" hops on the dual diamond lattice. We then predict that the "monopole" continuum has an enhanced spectral periodicity. This prediction can be examined among the existing data on the non-Kramers doublet spin liquid candidate materials like Pr$_2$TM$_2$O$_7$ and Tb$_2$TM$_2$O$_7$ (with TM = "transition metal"). The application to the Kramers doublet systems and numerical simulation is further discussed. Finally, we present a general classification of distinct symmetry enriched U(1) quantum spin liquids based on the translation symmetry fractionalization patterns of "monopoles" and "spinons"., Comment: 8 pages, 2 figures, 2 tables, update title, small correction

We describe electron spin resonance in a quantum spin liquid with significant spin-orbit coupling. We find that the resonance directly probes spinon continuum which makes it an efficient and informative probe of exotic excitations of the spin liquid. Specifically, we consider spinon resonance of three different spinon mean-field Hamiltonians, obtained with the help of projective symmetry group analysis, which model a putative quantum spin liquid state of the triangular rare-earth antiferromagnet YbMgGaO4. The band of absorption is found to be very broad and exhibit strong van Hove singularities of single spinon spectrum as well as pronounced polarization dependence., 14 pages, 9 figures (including supplement)

The ground state of the quantum kagome antiferromagnet Zn-brochantite, ZnCu$_3$(OH)$_6$SO$_4$, which is one of only a few known spin-liquid (SL) realizations in two or three dimensions, has been described as a gapless SL with a spinon Fermi surface. Employing nuclear magnetic resonance in a broad magnetic-field range down to millikelvin temperatures, we show that in applied magnetic fields this enigmatic state is intrinsically unstable against a SL with a full or a partial gap. A similar instability of the gapless Fermi-surface SL was previously encountered in an organic triangular-lattice antiferromagnet, suggesting a common destabilization mechanism that most likely arises from spinon pairing. A salient property of this instability is that an infinitesimal field suffices to induce it, as predicted theoretically for some other types of gapless SL's., minor corrections

We explore the thermodynamic properties of the layered copper(II) carbodiimide CuNCN by heat-capacity measurements and investigate the corresponding thermal atomic motions by means of neutron powder diffraction as well as inelastic neutron scattering. The experiments are complemented by a combination of density-functional calculations, phonon analysis and analytic theory. The existence of a soft flexural mode-bending of the layers, characteristic for the material structure-is established in the phonon spectrum of CuNCN by giving characteristic temperature-dependent contributions to the heat capacity and atomic displacement parameters. The agreement with the neutron data allows us to extract a residual-on top of the lattice-presumably spinon contribution to the heat capacity [Formula: see text], speaking in favor of the spin-liquid picture of the electronic phases of CuNCN.

Quantum spin ice, modeled for magnetic rare-earth pyrochlores, has attracted great interest for hosting a U(1) quantum spin liquid, which involves spin-ice monopoles as gapped deconfined spinons, as well as gapless excitations analogous to photons. However, the global phase diagram under a [111] magnetic field remains open. Here we uncover by means of unbiased quantum Monte-Carlo simulations that a supersolid of monopoles, showing both a superfluidity and a partial ionization, intervenes the kagom\'e spin ice and a fully ionized monopole insulator, in contrast to classical spin ice where a direct discontinuous phase transition takes place. We also show that on cooling, kagom\'e spin ice evolves towards a valence bond solid similar to what appears in the associated kagom\'e lattice model [S. V. Isakov et al., Phys. Rev. Lett. 97, 147202 (2006)]. Possible relevance to experiments is discussed., Comment: 5 pages, 4 figures; accepted for publication in PRL

Motivated by the recent proposal of the spinon Fermi surface spin liquids for several candidate materials such as YbMgGaO4, we explore the experimental consequences of the external magnetic fields on this exotic state. Specifically, we focus on the weak field regime where the spin liquid state is well preserved and the spinon remain to be a good description of the magnetic excitations. From the spin-1/2 nature of the spinon excitation, we predict the unique features of the spinon continuum when the weak magnetic field is applied to the system. Due to the small energy scale of the exchange interactions between the local moments in the spin liquid candidate like YbMgGaO4, our proposal for the spectral weight shifts and spectral crossing in the magnetic fields can be immediately tested by inelastic neutron scattering experiments. Several other experimental aspects about the spinon Fermi surface and the spinon excitations are discussed and proposed. Our work provides an experimental scheme to examine the fractionalized spinon excitation and the candidate spin liquid states in YbMgGaO4, the 6H-B phase of Ba3NiSb2O9 and other relevant materials., 9 pages, 5 figures, modified title and discussion

We report a new kagome quantum spin liquid candidate Cu$_3$Zn(OH)$_6$FBr, which does not experience any phase transition down to 50 mK, more than three orders lower than the antiferromagnetic Curie-Weiss temperature ($\sim$ 200 K). A clear gap opening at low temperature is observed in the uniform spin susceptibility obtained from $^{19}$F nuclear magnetic resonance measurements. We observe the characteristic magnetic field dependence of the gap as expected for fractionalized spin-1/2 spinon excitations. Our experimental results provide firm evidence for spin fractionalization in a topologically ordered spin system, resembling charge fractionalization in the fractional quantum Hall state., 5 pages and 3 figures. Supplementary Materials file is in the source file and also can be found in CPL published online version

Material line defects are one-dimensional structures but the search and proof of electron behaviour consistent with the reduced dimension of such defects has been so far unsuccessful. Here we show using angle resolved photoemission spectroscopy that twin-grain boundaries in the layered semiconductor MoSe2 exhibit parabolic metallic bands. The one-dimensional nature is evident from a charge density wave transition, whose periodicity is given by kF/p, consistent with scanning tunnelling microscopy and angle resolved photoemission measurements. Most importantly, we provide evidence for spin- and charge-separation, the hallmark of one-dimensional quantum liquids. Our studies show that the spectral line splits into distinctive spinon and holon excitations whose dispersions exactly follow the energy-momentum dependence calculated by a Hubbard model with suitable finite-range interactions. Our results also imply that quantum wires and junctions can be isolated in line defects of other transition metal dichalcogenides, which may enable quantum transport measurements and devices., The USF group acknowledges support from the National Science Foundation (DMR-1204924). V.K., R.D. and M.-H. P. acknowledges support from the Army Research Office (W911NF-15-1-0626) and thank Prof. Hari Srikanth for resistance measurements in his laboratory. M.C.A., J.A. and C.C. thank enlightening exchanges with Gabriel Kotliar and Zhi-Xun Shen. The Synchrotron SOLEIL is supported by the Centre National de la Recherche Scientifique (CNRS) and the Commissariat a` l’Energie Atomique et aux Energies Alternatives (CEA), France. T.Cˇ. and J.M.P.C. thank Eduardo Castro, Hai-Qing Lin and Pedro D. Sacramento for illuminating discussions. The theory group acknowledges the support from NSAF U1530401 and computational resources from CSRC (Beijing), the Portuguese FCT through the Grant UID/FIS/04650/2013 and the NSFC Grant 11650110443., info:eu-repo/semantics/publishedVersion

${\mathrm{Na}}_{4}{\mathrm{Ir}}_{3}{\mathrm{O}}_{8}$ provides a material platform to study three-dimensional quantum spin liquids in the geometrically frustrated hyperkagome lattice of ${\mathrm{Ir}}^{4+}$ ions. In this work, we consider quantum spin liquids on a hyperkagome lattice for generic spin models, focusing on the effects of anisotropic spin interactions. In particular, we classify possible ${\mathbb{Z}}_{2}$ and $U(1)$ spin liquid states, following the projective symmetry group analysis in the slave-fermion representation. There are only three distinct ${\mathbb{Z}}_{2}$ spin liquids, together with 2 different $U(1)$ spin liquids. The nonsymmorphic space group symmetry of the hyperkagome lattice plays a vital role in simplifying the classification, forbidding ``$\ensuremath{\pi}$-flux'' or ``staggered-flux'' phases in contrast to symmorphic space groups. We further prove that both $U(1)$ states and one ${\mathbb{Z}}_{2}$ state among all 3 are symmetry-protected gapless spin liquids, robust against any symmetry-preserving perturbations. Motivated by the ``spin-freezing'' behavior recently observed in ${\mathrm{Na}}_{4}{\mathrm{Ir}}_{3}{\mathrm{O}}_{8}$ at low temperatures, we further investigate the nearest-neighbor spin model with the dominant Heisenberg interaction subject to all possible anisotropic perturbations from spin-orbit couplings. We find that a $U(1)$ spin liquid ground state with spinon Fermi surfaces is energetically favored over ${\mathbb{Z}}_{2}$ states. Among all spin-orbit coupling terms, we show that only the Dzyaloshinskii-Moriya interaction can induce spin anisotropy in the ground state when perturbing from the isotropic Heisenberg limit. Our work paves the way for a systematic study of quantum spin liquids in various materials with a hyperkagome crystal structure.

Inelastic neutron scattering is used to study the low-energy magnetic excitations in the spin-1 triangular lattice of the 6H-B phase of Ba$_3$NiSb$_2$O$_9$. We study two powder samples: Ba$_3$NiSb$_2$O$_9$ synthesized under high pressure and Ba$_{2.5}$Sr$_{0.5}$NiSb$_2$O$_9$ in which chemical pressure stabilizes the 6H-B structure. The measured excitation spectra show broad gapless and nondispersive continua at characteristic wave vectors. Our data rules out most theoretical scenarios that have previously been proposed for this phase, and we find that it is well described by an exotic quantum spin liquid with three flavors of unpaired fermionic spinons, forming a large spinon Fermi surface.

We study quantum disordered ground states of the two dimensional Heisenberg-Kitaev model on the triangular lattice using a Schwinger boson approach. Our aim is to identify and characterize potential gapped quantum spin liquid phases that are stabilized by anisotropic Kitaev interactions. For antiferromagnetic Heisenberg- and Kitaev couplings and sufficiently small spin $S$ we find three different symmetric $Z_2$ spin liquid phases, separated by two continuous quantum phase transitions. Interestingly, the gap of elementary excitations remains finite throughout the transitions. The first spin liquid phase corresponds to the well known zero-flux state in the Heisenberg limit, which is stable with respect to small Kitaev couplings and develops $120^\circ$ order in the semi-classical limit at large $S$. In the opposite Kitaev limit we find a different spin liquid ground-state, which is a quantum disordered version of a magnetically ordered state with antiferromagnetic chains, in accordance with results in the classical limit. Finally, at intermediate couplings we find a spin liquid state with unconventional spin correlations. Upon spinon condensation this state develops Bragg peaks at incommensurate momenta in close analogy to the magnetically ordered Z2 vortex crystal phase, which has been analyzed in recent theoretical works., Comment: 14 pages, 13 figures

Spin liquids are exotic quantum states characterized by the existence of fractional and deconfined quasiparticle excitations, referred to as spinons and visons. Their fractional nature establishes topological properties such as a protected ground-state degeneracy. This work investigates spin-orbit coupled spin liquids where, additionally, topology enters via non-trivial band structures of the spinons. We revisit the $Z_2$ spin-liquid phases that have recently been identified in a projective symmetry-group analysis on the square lattice when spin-rotation symmetry is maximally lifted [Phys. Rev. B 90, 174417 (2014)]. We find that in the case of nearest neighbor couplings only, $Z_2$ spin liquids on the square lattice always exhibit trivial spinon bands. Adding second neighbor terms, the simplest projective symmetry-group solution closely resembles the Bernevig-Hughes-Zhang model for topological insulators. Assuming that the emergent gauge fields are static we investigate vison excitations, which we confirm to be deconfined in all investigated spin phases. Particularly, if the spinon bands are topological, the spinons and visons form bound states consisting of several spinon-Majorana zero modes coupling to one vison. The existence of such zero modes follows from an exact mapping between these spin phases and topological $p+ip$ superconductors with vortices. We propose experimental probes to detect such states in real materials., Comment: 16 pages, 7 figures

The search for exotic quantum spin liquid states in simple yet realistic spin models remains a central challenge in the field of frustrated quantum magnetism. Here we consider the canonical nearest-neighbor kagome Heisenberg antiferromagnet restricted to a quasi-1D strip consisting entirely of corner-sharing triangles. Using large-scale density matrix renormalization group calculations, we identify in this model an extended gapless quantum phase characterized by central charge $c=2$ and power-law decaying spin and bond-energy correlations which oscillate at tunably incommensurate wave vectors. We argue that this intriguing spin liquid phase can be understood as a marginal instability of a two-band spinon Fermi surface coupled to an emergent U(1) gauge field, an interpretation which we substantiate via bosonization analysis and Monte Carlo calculations on model Gutzwiller variational wave functions. Our results represent one of the first numerical demonstrations of emergent fermionic spinons in a simple SU(2) invariant nearest-neighbor Heisenberg model beyond the strictly 1D (Bethe chain) limit., Comment: 14 pages, 12 figures

We develop a theory of the dynamical response of a minimal model of quantum spin ice (QSI) by means of inelastic light scattering. In particular, we are interested in the Raman response of the fractionalized U(1) spin liquid realized in the XXZ QSI. We show that the low-energy Raman intensity is dominated by spinon and gauge fluctuations. We find that the Raman response in the QSI state of that model appears only in the $T_{2g}$ polarization channel. We show that the Raman intensity profile displays a broad continuum from the spinons and coupled spinon and gauge fluctuations, and a low-energy peak arising entirely from gauge fluctuations. Both features originate from the exotic interaction between photon and the fractionalized excitations of QSI. Our theoretical results suggest that inelastic Raman scattering can in principle serve as a promising experimental probe of the nature of a U(1) spin liquid in QSI., Comment: 16 pages, 4 figures

We present an exact ground state solution of a quantum dimer model introduced in Ref.[1], which features ordinary bosonic spin-singlet dimers as well as fermionic dimers that can be viewed as bound states of spinons and holons in a hole-doped resonating valence bond liquid. Interestingly, this model captures several essential properties of the metallic pseudogap phase in high-$T_c$ cuprate superconductors. We identify a line in parameter space where the exact ground state wave functions can be constructed at an arbitrary density of fermionic dimers. At this exactly solvable line the ground state has a huge degeneracy, which can be interpreted as a flat band of fermionic excitations. Perturbing around the exactly solvable line, this degeneracy is lifted and the ground state is a fractionalized Fermi liquid with a small pocket Fermi surface in the low doping limit., Comment: Revised version, 8 pages

The free fermion nature of interacting spins in one dimensional (1D) spin chains still lacks a rigorous study. In this letter we show that the length-$1$ spin strings significantly dominate critical properties of spinons, magnons and free fermions in the 1D antiferromagnetic spin-1/2 chain. Using the Bethe ansatz solution we analytically calculate exact scaling functions of thermal and magnetic properties of the model, providing a rigorous understanding of the quantum criticality of spinons. It turns out that the double peaks in specific heat elegantly mark two crossover temperatures fanning out from the critical point, indicating three quantum phases: the Tomonaga-Luttinger liquid (TLL), quantum critical and fully polarized ferromagnetic phases. For the TLL phase, the Wilson ratio $R_W=4K_s$ remains almost temperature-independent, here $K_s$ is the Luttinger parameter. Furthermore, applying our results we precisely determine the quantum scalings and critical exponents of all magnetic properties in the ideal 1D spin-1/2 antiferromagnet Cu(C${}_4$H${}_4$N${}_2$)(NO${}_3$)${}_2$ recently studied in Phys. Rev. Lett. {\bf 114}, 037202 (2015)]. We further find that the magnetization peak used in experiments is not a good quantity to map out the finite temperature TLL phase boundary., Comment: 6 pages + 10 pages, 4 figures + 4 figures

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., Comment: 10 pages, 5 figures

Throughout this dissertation, we have made the claim that the low-energy degrees of freedom in quantum spin chains are effective particles that move against a strongly-correlated background. These quasiparticles cannot be understood as dressed versions of the microscopic degrees of freedom – as is the case in e.g. Fermi-liquid theory – but they are collective modes that often have exotic emergent properties.

Quantum spin liquids (QSLs) are exotic phases of matter exhibiting long-range entanglement and supporting emergent gauge fields. A vigorous search for experimental realizations of these states has identified several materials with properties hinting at QSL physics. A key issue in understanding these QSL candidates is often the interplay of weak disorder of the crystal structure with the spin liquid state. It has recently been pointed out that in at least one important class of candidate QSLs - pyrochlore magnets based on non-Kramers ions such as Pr$^{3+}$ or Tb$^{3+}$- structural disorder can actually promote a $U(1)$ QSL ground state. Here we set this proposal on a quantitative footing by analyzing the stability of the QSL state in the minimal model for these systems: a random transverse field Ising model. We consider two kinds of instability, which are relevant in different limits of the phase diagram: condensation of spinons and confinement of the $U(1)$ gauge fields. Having obtained stability bounds on the QSL state we apply our results directly to the disordered candidate QSL Pr$_2$Zr$_2$O$_7$. We find that the available data for currently studied samples of Pr$_2$Zr$_2$O$_7$ is most consistent with it a ground state outside the spin liquid regime, in a paramagnetic phase with quadrupole moments near saturation due to the influence of structural disorder., Comment: Significant revisions including new calculation of confinement instability of the QSL. Accepted for publication in Physical Review Letters

Historically, the Heisenberg chain has been the first exactly solved (interacting) model. It is a fairly realistic model for a one-dimensional quantum magnet and it has been instrumental in moving beyond the classical Ising-type models. We will solve this chain following the same steps we introduced in the previous chapter for the Lieb-Liniger model: in Sect. 3.2 we introduce the chain’s fundamental excitation, the magnon, in Sect. 3.3 we study the two-body problem and in Sect. 3.4 we write the coordinate Bethe Ansatz solution and the Bethe equations. We will see that the latter admit a mixture of complex and real solutions and discuss the challenges this implies. For infinitely long chains, it is believe that the string hypothesis allows to organize the different bound solutions: this is the topic of Sect. 3.4.1. Having organized the Hilbert space of the model, we discuss the ground state and the low energy excitations for the ferromagnetic and antiferromagnetic order in Sects. 3.5 and 3.6, respectively. For the latter, we introduce and discuss the role of spinons as emergent quasi-particles. Finally, in Sect. 3.7 we study the effect of switching on an external magnetic field.

In a spontaneously dimerized quantum antiferromagnet, spin-1/2 excitations (spinons) are confined in pairs by strings akin to those confining quarks in non-abelian gauge theories. The system has multiple degenerate ground states (vacua) and domain walls between regions of different vacua. For two vacua, we demonstrate that spinons on a domain wall are liberated, in a mechanism strikingly similar to domain-wall deconfinement of quarks in variants of quantum chromodynamics. This observation not only establishes a novel phenomenon in quantum magnetism, but also provides a new direct link between particle physics and condensed-matter physics. The analogy opens doors to improving our understanding of particle confinement and deconfinement by computational and experimental studies in quantum magnetism., Comment: Title changed, and overall condensation of the text. 5 pages, 4 figures