We study the energy and the static spin structure factor of the ground state of the spin-1/2 quantum Heisenberg antiferromagnetic model on the kagome lattice. By the iterative application of a few Lanczos steps on accurate projected fermionic wave functions and the Green's function Monte Carlo technique, we find that a gapless (algebraic) U (1) Dirac spin liquid is competitive with previously proposed gapped (topological) Z2 spin liquids. By performing a finite-size extrapolation of the ground-state energy, we obtain an energy per site E/J = −0.4365(2), which is equal, within three error bars, to the estimates given by the density-matrix renormalization group (DMRG). Our estimate is obtained for a translationally invariant system, and, therefore, does not suffer from boundary effects, like in DMRG. Moreover, on finite toric clusters at the pure variational level, our energies are lower compared to those from DMRG calculations.PACS numbers: 75.10. Jm, 75.10.Kt, 75.40.Mg, 75.50.Ee Introduction. The spin-1/2 quantum Heisenberg antiferromagnet (QHAF) on the kagome lattice provides a conducive environment to stabilize a quantum paramagnetic phase of matter down to zero temperature, [1][2][3] a fact that has been convincingly established theoretically from several studies, including exact diagonalization, 4-8 series expansion, 9,10 quantum Monte Carlo,
11and analytical techniques. 12 The question of the precise nature of the spin-liquid state of the kagome spin-1/2 QHAF has been intensely debated on the theoretical front, albeit without any definitive conclusions. Different approximate numerical techniques have claimed a variety of ground states. On the one hand, densitymatrix renormalization group (DMRG) calculations have been claimed for a fully gapped (nonchiral) Z 2 topological spin-liquid ground state that does not break any point group symmetry.13,14 On the other hand, an algebraic and fully symmetric U (1) Dirac spin liquid has been proposed as the ground state, by using projected fermionic wave functions and the variational Monte Carlo (VMC) approach. [15][16][17][18][19][20] In addition, valence bond crystals have been suggested from many other techniques. In particular, a 36-site unit cell valence-bond crystal [21][22][23] was proposed using quantum dimer models, 24-28 series expansion 29,30 and multiscale entanglement renormalization ansatz (MERA) 31 techniques. Finally, a recent coupled cluster method (CCM) suggested a q = 0 (uniform) state.
32On general theoretical grounds, the Z 2 spin liquids in two spatial dimensions are known to be stable phases, 33-35 as compared to algebraic U (1) spin liquids, which are known to be only marginally stable.36 However, explicit numerical calculations using projected wave functions have shown the U (1) Dirac spin liquid to be stable (locally and globally) with respect to dimerizing into all known valence-bond crystal phases. 15,17,18,20 Furthermore, it was shown that, within this class of Gutzwiller projected wave functions, all the fully symmetric, gapped Z 2 spin ...