Supersolidity in the triangular lattice spin-1/2 XXZ model: a variational perspective.

We study the spin-1/2 XXZ model on the triangular lattice with a nearest neighbor antiferromagnetic Ising coupling J(z) > 0 and unfrustrated (J(perpendicular) < 0) or frustrated (J(perpendicular) >0) kinetic terms in a zero magnetic field. Incorporating long-range Jastrow correlations over a mean-field spin state, we obtain the variational phase diagram of this model on large lattices for arbitrary J(z) and either sign of J(perpendicular). For J(perpendicular) < 0, we find a square root(3) x square root(3) supersolid for J(z)/J(perpendicular)| approximately > 4.7, in excellent agreement with quantum Monte Carlo data. For J(perpendicular) > 0, a distinct square root(3) x square root(3) supersolid is found to emerge for J(z)/J(perpendicular) > or = 1. Both supersolids exhibit a spontaneous density deviation from half-filling. At J(z)/J(perpendicular) = infinity, the crystalline order parameters of these two supersolids are nearly identical, consistent with exact results.