High-Order Coupled Cluster Method (CCM) Calculations for Quantum Magnets with Valence-Bond Ground States

In this article, we prove that exact representations of dimer and plaquette valence-bond ket ground states for quantum Heisenberg antiferromagnets may be formed via the usual coupled cluster method (CCM) from independent-spin product (e.g. Néel) model states. We show that we are able to provide good results for both the ground-state energy and the sublattice magnetization for dimer and plaquette valence-bond phases within the CCM. As a first example, we investigate the spin-half J ₁–J ₂ model for the linear chain, and we show that we are able to reproduce exactly the dimerized ground (ket) state at J ₂/J ₁=0.5. The dimerized phase is stable over a range of values for J ₂/J ₁ around 0.5, and results for the ground-state energies are in good agreement with the results of exact diagonalizations of finite-length chains in this regime. We present evidence of symmetry breaking by considering the ket- and bra-state correlation coefficients as a function of J ₂/J ₁. A radical change is also observed in the behavior of the CCM sublattice magnetization as we enter the dimerized phase. We then consider the Shastry-Sutherland model and demonstrate that the CCM can span the correct ground states in both the Néel and the dimerized phases. Once again, very good results for the ground-state energies are obtained. We find CCM critical points of the bra-state equations that are in agreement with the known phase transition point for this model. The results for the sublattice magnetization remain near to the “true” value of zero over much of the dimerized regime, although they diverge exactly at the critical point. Finally, we consider a spin-half system with nearest-neighbor bonds for an underlying lattice corresponding to the magnetic material CaV₄O₉ (CAVO). We show that we are able to provide excellent results for the ground-state energy in each of the plaquette-ordered, Néel-ordered, and dimerized regimes of this model. The exact plaquette and dimer ground states are reproduced by the CCM ket state in their relevant limits. Furthermore, we estimate the range over which the Néel order is stable, and we find the CCM result is in reasonable agreement with the results obtained by other methods. Our new approach has the dual advantages that it is simple to implement and that existing CCM codes for independent-spin product model states may be used from the outset. Furthermore, it also greatly extends the range of applicability to which the CCM may be applied. We believe that the CCM now provides an excellent choice of method for the study of systems with valence-bond quantum ground states.