Quantum Monte Carlo studies of spinons in one-dimensional spin systems

Abstract

Observing constituent particles with fractional quantum numbers in confined and deconfined states is an interesting and challenging problem in quantum many-body physics. Here, we further explore a computational scheme [Y. Tang and A. W. Sandvik, Phys. Rev. Lett. 107, 157201 (2011)] based on valence-bond quantum Monte Carlo simulations of quantum spin systems. Using several different one-dimensional models, we characterize $S=\frac{1}{2}$ spinon excitations using the intrinsic spinon size $\ensuremath{\lambda}$ and confinement length $\mathrm{\ensuremath{\Lambda}}$ (the size of a bound state). The spinons have finite size in valence-bond-solid states, infinite size in the critical region (with overlaps characterized by power laws), and become ill defined (completely unlocalizable) in the N\'eel state (which we stabilize in one dimension by introducing long-range interactions). We also verify that pairs of spinons are deconfined in uniform spin chains but become confined upon introducing a pattern of alternating coupling strengths (dimerization) or coupling two chains (forming a ladder). In the dimerized system, an individual spinon can be small when the confinement length is large; this is the case when the imposed dimerization is weak but the ground state of the corresponding uniform chain is a spontaneously formed valence-bond solid (where the spinons are deconfined). Based on our numerical results, we argue that a system with $\ensuremath{\lambda}\ensuremath{\ll}\mathrm{\ensuremath{\Lambda}}$ is associated with weak repulsive short-range spinon-spinon interactions. In principle, both the length scales $\ensuremath{\lambda}$ and $\mathrm{\ensuremath{\Lambda}}$ can still be individually tuned from small to infinite (with $\ensuremath{\lambda}\ensuremath{\le}\mathrm{\ensuremath{\Lambda}}$) by varying model parameters. In contrast, in the ladder system the two lengths are always similar, and this is the case also in the weakly dimerized systems when the corresponding uniform chain is in the critical phase. In these systems, the effective spinon-spinon interactions are purely attractive and there is only a single large length scale close to criticality, which is reflected in the standard spin correlations as well as in the spinon characteristics.