We study the plaquette valence-bond solid phase of the spin-$1/2$${J}_{1}$-${J}_{2}$ antiferromagnet Heisenberg model on the square lattice within the bond-operator theory. We start by considering four $S=1/2$ spins on a single plaquette and determine the bond operator representation for the spin operators in terms of singlet, triplet, and quintet boson operators. The formalism is then applied to the ${J}_{1}$-${J}_{2}$ model and an effective interacting boson model in terms of singlets and triplets is derived. The effective model is analyzed within the harmonic approximation and the previous results of Zhitomirsky and Ueda [Phys. Rev. B 54, 9007 (1996)] are recovered. By perturbatively including cubic (triplet-triplet-triplet and singlet-triplet-triplet) and quartic interactions, we find that the plaquette valence-bond solid phase is stable within the parameter region $0.34<{J}_{2}/{J}_{1}<0.59$, which is narrower than the harmonic one. Differently from the harmonic approximation, the excitation gap vanishes at both critical couplings ${J}_{2}=0.34{J}_{1}$ and ${J}_{2}=0.59{J}_{1}$. Interestingly, for ${J}_{2}<0.48{J}_{1}$, the excitation gap corresponds to a singlet-triplet excitation at the $\ensuremath{\Gamma}$ point while, for ${J}_{2}>0.48{J}_{1}$, it is related to a singlet-singlet excitation at the $\mathbf{X}=(\ensuremath{\pi}/2,0)$ point of the tetramerized Brillouin zone.
Read full abstract