We investigate the nature of the ground state of the spin-12 Heisenberg antiferromagnet on the shuriken lattice by complementary state-of-the-art numerical techniques, such as variational Monte Carlo (VMC) with versatile Gutzwiller-projected Jastrow wave functions, unconstrained multivariable variational Monte Carlo (mVMC), and pseudofermion/pseudo-Majorana functional renormalization group (PFFRG/PMFRG) methods. We establish the presence of a quantum paramagnetic ground state and investigate its nature, by classifying symmetric and chiral quantum spin liquids, and inspecting their instabilities towards competing valence bond crystal (VBC) orders. Our VMC analysis reveals that a VBC with a pinwheel structure emerges as the lowest-energy variational ground state, and it is obtained as an instability of the U(1) Dirac spin liquid. Analogous conclusions are drawn from mVMC calculations employing accurate BCS pairing states supplemented by symmetry projectors, which confirm the presence of pinwheel VBC order by a thorough analysis of dimer-dimer correlation functions. Our work highlights the nontrivial role of quantum fluctuations via the Gutzwiller projector in resolving the subtle interplay between competing orders.

We assess the ground-state phase diagram of the J1-J2 Heisenberg model on the kagome lattice by employing Gutzwiller-projected fermionic wave functions. Within this framework, different states can be represented, defined by distinct unprojected fermionic Hamiltonians that include hopping and pairing terms, as well as a coupling to local Zeeman fields to generate magnetic order. For J2 = 0, the so-called U(1) Dirac state, in which only hopping is present (such as to generate a π-flux in the hexagons), has been shown to accurately describe the exact ground state [Y.

No abstract