Strongly correlated low-dimensional quantum spin models provide a well-established frame- \nwork to study magnetic properties of insulators, and are of great theoretical interest and experimental relevance in condensed-matter physics. In this thesis, I use quantum Monte Carlo methods to numerically study quantum critical behavior in low-dimensional quantum spin models and wavefunctions. \n \nFirst, I study spinons &ndash emergent spin-1/2 bosonic excitations &ndash at certain one- and two-dimensional quantum phase transitions (QPTs) in spin models, by characterizing their size and confinement length quantitatively. In particular, I focus on the QPT from an antiferromagnetic (AFM) phase into a valence-bond solid (VBS) phase, which is an example of a violation of the standard Landau-Ginzburg-Wilson paradigm for phase transitions. This transition in two dimensions (2D) is instead likely described by a novel theory called "deconfined quantum criticality" (DQC). According to the theory, spinons should be deconfined. The degree of deconfinement is quantified in my calculations. \n \nSecond, I present a comprehensive study of so-called short-bond resonating-valence-bond (RVB) spin liquids in 2D, which have been suggested as a good starting point for understanding the spin physics of high-temperature cuprates. I find that these RVB states can also be classified as quantum-critical VBS states, which indicates that RVB is less disordered than expected. This work suggests a possible mapping from the quantum RVB states to classical dimer models via a classical continuum field theory--the height model. This map explicitly bridges well-established classical results to future quantum studies. \n \nThird, I consider 1D amplitude product (AP) states, which are generalized versions of RVB states, with different wavefunction weightings of bonds according to their lengths. AP states constitute a good ansatz for certain Hamiltonians and are of broad interest in quantum magnetism. I study phase transitions from AFM-VBS phases in AP states by tuning their amplitudes, and obtain continuously varying critical exponents. In addition, I classify the 1D AP states through entanglement entropy calculations of the central charge in (1+1)D conformal field theory. This new classification could serve as guide for AP states as trial wavefunctions to search for ground states of corresponding quantum spin models.
We consider the
S=1/2 Heisenberg model with nearest-neighbor interaction
J and an additional multispin interaction
Q3 on the square lattice. The
Q3 term consists of three bond-singlet projectors and is chosen to favor the formation of a valence-bond solid (VBS) where the valence bonds (singlet pairs) form a staggered pattern. The model exhibits a quantum phase transition from the N\'eel state to the VBS as a function of
Q3/J. We study the model using quantum Monte Carlo (stochastic series expansion) simulations. The N\'eel-VBS transition in this case is strongly first order in nature, in contrast to similar previously studied models with continuous transitions into columnar VBS states. The qualitatively different transitions illustrate the important role of an emerging
U(1) symmetry in the latter case, which is not possible in the present model due to the staggered VBS pattern (which does not allow local fluctuations necessary to rotate the local VBS order parameter).
We study the
S=1/2 Heisenberg (
J) model on the two-dimensional (2D) square lattice in the presence of additional higher-order spin interactions (
Q) which lead to a valence-bond-solid (VBS) ground state. Using quantum Monte Carlo simulations, we analyze the thermal VBS transition. We find continuously varying exponents, with the correlation-length exponent $\ensuremath{\nu}$ close to the Ising value for large
Q/J and diverging when
Q/J approaches the quantum-critical point (the critical temperature ${T}_{c}\ensuremath{\rightarrow}0$). We identify the transition with a class of conformal field theories with charge
c=1 and critical exponents varying between those of the 2D Ising model and the Kosterlitz-Thouless (KT) fixed point. We find explicit evidence for KT physics by studying the emergence of
U(1) symmetry of the order parameter at
T=Tc when ${T}_{c}\ensuremath{\rightarrow}0$.
We study two-dimensional Heisenberg antiferromagnets with additional multispin interactions which can drive the system into a valence-bond-solid state. For standard SU(2) spins, we consider both four- and six-spin interactions. We find continuous quantum phase transitions with the same critical exponents. Extending the symmetry to
SU(N), we also find continuous transitions for
N=3 and 4. In addition, we also study quantitatively the crossover of the order-parameter symmetry from
Z4 deep inside the valence-bond-solid phase to
U(1) as the phase transition is approached.
Various lattice geometries and boundary conditions are used to investigate valence-bond-solid (VBS) ordering in the ground state of an
S=1/2 square-lattice quantum spin model---the
J-
Q model, in which four- or six-spin interactions
Q are added to the standard Heisenberg exchange
J. Ground-state results for finite systems (with up to thousands of spins) are obtained using an unbiased projector quantum Monte Carlo method. It is found that great care has to be taken when extrapolating the order parameter to infinite lattice size, in particular, in cylinder geometry. Even though strong VBS order exists in two dimensions, and is established clearly with increasing system size on $L\ifmmode\times\else\texttimes\fi{}L$ lattices (or ${L}_{x}\ifmmode\times\else\texttimes\fi{}{L}_{y}$ lattices with a fixed aspect ratio
Lx/Ly of order 1), only short-range VBS correlations are observed on long cylinders (when ${L}_{x}\ensuremath{\rightarrow}\ensuremath{\infty}$ at fixed
Ly). The correlation length increases with the cylinder width, until long-range order sets in at a ``critical'' width. This width is very large even when the 2D order is relatively strong. For example, for a system in which the order parameter is
70 of the largest possible value,
Ly=8 is required for ordering. Extrapolations of the VBS order parameter based on correlation functions (the square of the order parameter) for small $L\ifmmode\times\else\texttimes\fi{}L$ lattices can also be misleading. For a
20-ordered system, results for
L up to $\ensuremath{\approx}$20 appear to extrapolate clearly to a vanishing order parameter, while for larger lattices the scaling behavior crosses over and extrapolates to a nonzero value (with exponentially small finite-size corrections). The VBS order parameter also exhibits interesting edge effects related the known emergent U(1) symmetry close to a ``deconfined'' critical point, which, if not considered properly, can lead to wrong conclusions for the thermodynamic limit. The observed finite-size behavior for small $L\ifmmode\times\else\texttimes\fi{}L$ lattices and long cylinders is very similar to that predicted for a
Z2 spin liquid. The results therefore raise concerns about recent numerical work claiming
Z2 spin-liquid ground states in 2D frustrated quantum spin systems, in particular, the Heisenberg model with nearest and next-nearest-neighbor couplings. Based on the results presented so far, a VBS state in this system cannot be ruled out.
We calculate the bipartite von Neumann and second R\'enyi entanglement entropies of the ground states of spin-1/2 dimerized Heisenberg antiferromagnets on a square lattice. Two distinct dimerization patterns are considered: columnar and staggered. In both cases, we concentrate on the valence bond solid (VBS) phase and describe such a phase with the bond-operator representation. Within this formalism, the original spin Hamiltonian is mapped into an effective interacting boson model for the triplet excitations. We study the effective Hamiltonian at the harmonic approximation and determine the spectrum of the elementary triplet excitations. We then follow an analytical procedure, which is based on a modified spin-wave theory for finite systems and was originally employed to calculate the entanglement entropies of magnetic ordered phases, and calculate the entanglement entropies of the VBS ground states. In particular, we consider one-dimensional (line) subsystems within the square lattice, a choice that allows us to consider line subsystems with sizes up to ${L}^{\ensuremath{'}}\phantom{\rule{0.16em}{0ex}}=\phantom{\rule{0.16em}{0ex}}1000$. We combine such a procedure with the results of the bond-operator formalism at the harmonic level and show that, for both dimerized Heisenberg models, the entanglement entropies of the corresponding VBS ground states obey an area law as expected for gapped phases. For both columnar-dimer and staggered-dimer models, we also show that the entanglement entropies increase but they seem to not diverge as the dimerization decreases and the system approaches the N\'eel-VBS quantum phase transition. Finally, the entanglement spectra associated with the VBS ground states are presented.