We introduce Szβ spin-projection based on cluster mean-field theory and
apply it to the ground state of strongly-correlated spin systems. In cluster
mean-field, the ground state wavefunction is written as a factorized tensor
product of optimized cluster states. In previous work, we have focused on
unrestricted cluster mean-field, where each cluster is Szβ symmetry adapted.
We here remove this restriction by introducing a generalized cluster mean-field
(GcMF) theory, where each cluster is allowed to access all Szβ sectors,
breaking Szβ symmetry. In addition, a projection scheme is used to restore
global Szβ, which gives rise to Szβ spin-projected generalized cluster
mean-field (SzβGcMF). Both of these extensions contribute to accounting for
inter-cluster correlations. We benchmark these methods on the 1D, quasi-2D, and
2D J1ββJ2β and XXZ Heisenberg models. Our results indicate that the new
methods (GcMF and SzβGcMF) provide a qualitative and semi-quantitative
description of the Heisenberg lattices in the regimes considered, suggesting
them as useful references for further inter-cluster correlations, which are
discussed in this work