Twenty-five years after the first proposal, the question whether the ground state of a frustrated spin-half system is well described by a spin-liquid Resonating Valence Bond (RVB) wave function is still controversial. A physically transparent representation of a RVB state can be obtained in fermionic representation with a standard BCS-type pairing wave function, working in the subspace with fixed number of electrons and no double occupancies. In this work, we show that, using this variational wave function with a careful parameterization of the pairing function, it is possible to obtain an extremely accurate {\em ansatz} for the ground state of the Heisenberg antiferromagnet with next-nearest neighbors interactions (
J1−J2 model) in the regime of strong frustration. Indeed, in the spin-half realization of this model, it is known that the combined effect of frustration and zero-point motion interferes with the mechanism of spontaneously broken symmetry, giving rise to a non-magnetic phase of purely quantum-mechanical nature (
J2/J1≃0.5). This wave function is proposed to represent the generic spin-half RVB ground state in spin liquids.
We investigate the ground-state properties of two strongly-correlated systems: the two-dimensional
t−J model that can be used to study the high-temperature superconducting phase in the doped antiferromagnet, and the frustrated N\'eel antiferromagnet described by the
J1−J2 Heisenberg model on the square lattice, which is widely considered as the prototype model for spin frustration. In this thesis, we apply to these two systems state-of-the-art quantum Monte Carlo techniques, including the variational and the Green's function Monte Carlo with the fixed-node approximation. Few Lanczos steps are used to systematically improve the accuracy of the trial wave functions. By introducing a suitable regularization scheme for the variational calculations with few Lanczos steps, stable and controllable simulations can be performed up to very large cluster sizes with very good accuracy. In the two-dimensional
t−J model at
J/t=0.4, we show that the accuracy of the Gutzwiller-projected variational state (containing
dx2−y2 pairing) can be improved much by few Lanczos steps; in addition, the fixed-node Monte Carlo with these systematically improvable trial wave functions gives results that are comparable with the best accurate DMRG ones. Our main outcome is that the ground state is homogeneous and no evidence of stripes is detected around the doping
δ=1/8. Indeed, our best approximation to the ground state does not show any tendency towards charge inhomogeneity. Furthermore, our results show that a uniform state containing superconductivity and antiferromagnetism is stabilized at low hole doping, i.e.,
δ≲0.1. In the
J1−J2 Heisenberg model on the square lattice, we use the projected mean-field state that is built from Abrikosov fermions having a
Z2 gauge structure and four Dirac points in the spinon spectrum. No spin or dimer order is found in the strongly frustrated regime and our calculations imply that a spin liquid phase may faithfully represent the exact ground state around
J2/J1=0.5. The few Lanczos step technique is used to systematically improve the accuracy of the variational states both for the ground state and for few relevant low-energy excitations. This procedure allows us to estimate, in a valuable and systemic way, the spin gaps within thermodynamical limit and to show a solid evidence of an unconventional gapless excitation spectrum in the strongly frustrated regime, i.e.,
J2/J1≃0.5. In particular, we found gapless triplet excitations at momenta
(π,0) and
(0,π), which are compatible with the presence of four Dirac points at momenta
(±2π,±2π) in the spinon spectrum.
The strong interest in strongly correlated systems in condensed matter physics has continued unabated for the past few decades. In recent years, the number of novel, exotic quantum phases found in theoretical studies has seen a phenomenal rise. Among those interesting quantum states are bose liquids and spin liquids, where strong quantum fluctuations have prevented the systems from developing a long range order. Our work in this thesis seeks to further the understanding of frustrated systems. In the study of a hard-core boson model with ring-only exchange interactions on a square lattice, we obtain concrete numerical realization of the unconventional Exciton Bose Liquid (EBL) phase, which possesses interesting properties such as a "Bose surface" which resembles the Fermi surface in a metal, as well as unusual thermodynamic properties such as a T Log T dependence for specific heat. An equally important result from this work is the demonstration that the widely used Gutzwiller projection on slave-particle wave functions may generally fail to capture the correct long wavelength physics in the respective systems. For the Heisenberg antiferromagnet on the kagome lattice, which is a promising candidate for realizing a spin-disordered ground state, our variational study shows that the projected Schwinger boson wave function is energetically better than the Dirac spin liquid wave function when a small antiferromagnetic second-neighbor spin coupling is added to the nearest-neighbor model. We also study the anisotropic triangular Heisenberg antiferromagnetic in magnetic field, and find simple, yet accurate wave functions for various regions of the surprisingly rich phase diagram, thus providing insights into the energetics of the competing phases in this interesting model. Finally, our work also highlights permanent-type wave functions as potentially useful constructions in variational studies of systems with short-ranged correlations, e.g., a Mott insulator and a gapped spin liquid.
This thesis is devoted to the numerical investigation of quantum spin models which describe the low-energy physics of frustrated magnets. At extremely low temperatures, these systems can host the so-called spin liquid phase, an unconventional state of matter characterized by a high degree of quantum entanglement and the absence of magnetic order. The experimental identification of the spin liquid phase in actual materials relies on the detection of its distinctive excitations (named spinons), which possess fractional quantum numbers and can be probed by inelastic neutron scattering experiments. From the theoretical point of view, variational methods have been largely employed to tackle ground state properties of frustrated spin models. In particular, variational Monte Carlo techniques based on Gutzwiller-projected fermionic wave functions have been shown to provide accurate results for several frustrated systems. In this thesis, we pursue an extension of this variational scheme to target dynamical spectra, which are directly measured by inelastic neutron scattering experiments. Specifically, we compute the dynamical structure factor by constructing approximate excited states, which are obtained by applying two-spinon operators to the ground state wave function. Our results prove that this variational method can accurately describe the spectral features of different spin systems. Focusing on prototypical frustrated models on the square and triangular lattices, we observe how the dynamical structure factor reflects the phase transition between a magnetically ordered phase with spin wave excitations, to a spin liquid state with fractional degrees of freedom. In addition to spectral properties, we also explore new directions to improve the accuracy of Gutzwiller-projected states by the application of a neural network correlator, in the form of a restricted Boltzmann machine. While this hybrid variational scheme provides a considerable improvement of the variational energy in the case of unfrustrated spin models, less satisfactory results are obtained for frustrated systems, which call for further refinements of the neural network variational Ansätze.
Variational wave functions have enabled exceptional scientific breakthroughs related to the understanding of novel phases of matter. Examples include the Bardeen-Cooper-Schrieffer theory of superconductivity, the description of the fractional quantum Hall effect through the Laughlin state, and Feynman's variational understanding of large-scale quantum effects in liquid helium. More recently, Gutzwiller-projected wave functions, typically constructed from fermionic degrees of freedom, have been employed to examine quantum spin models in the presence of competing interactions, where exotic phases with no spontaneous symmetry breaking and fractional excitations may exist. In this work, we investigate the aforementioned fermionic wave functions supplemented with neural networks, specifically with the so-called restricted Boltzmann machine (RBM), to boost their accuracy and obtain reliable approximations to the ground state of generic spin models. In particular, we apply our neural augmented fermionic construction to the description of both magnetically ordered and disordered phases of increasing complexity, including cases where the ground state displays a nontrivial sign structure. Even though the RBM state is by far more effective for N\'eel states endowed with a particularly simple sign structure, it provides a significant improvement over the original fermionic state in highly frustrated regimes where a complex sign structure is anticipated, thus marking the path to an understanding of strongly correlated spin models on the lattice via neural Gutzwiller-projected variational wave functions.