Numerical contractor renormalization method for quantum spin models -: art. no. 104424

We demonstrate the utility of the numerical contractor renormalization (CORE) method for quantum spin systems by studying one- and two-dimensional model cases. Our approach consists of two steps: (i) building an effective Hamiltonian with longer ranged interactions up to a certain cutoff using the CORE algorithm and (ii) solving this new model numerically on finite clusters by exact diagonalization and performing finite-size extrapolations to obtain results in the thermodynamic limit. This approach, giving complementary information to analytical treatments of the CORE Hamiltonian, can be used as a semiquantitative numerical method. For ladder-type geometries, we explicitely check the accuracy of the effective models by increasing the range of the effective interactions until reaching convergence. Our results in the perturbative regime and also away from it are in good agreement with previously established results. In two dimensions we consider the plaquette lattice and the kagome lattice as nontrivial test cases for the numerical CORE method. As it becomes more difficult to extend the range of the effective interactions in two dimensions, we propose diagnostic tools (such as the density matrix of the local building block) to ascertain the validity of the basis truncation. On the plaquette lattice we have an excellent description of the system in both the disordered and the ordered phases, thereby showing that the CORE method is able to resolve quantum phase transitions. On the kagome lattice we find that the previously proposed twofold degenerate S=1/2 basis can account for a large number of phenomena of the spin 1/2 kagome system. For spin 3/2, however, this basis does not seem to be sufficient. In general we are able to simulate system sizes which correspond to an 8x8 lattice for the plaquette lattice or a 48-site kagome lattice, which are beyond the possibilities of a standard exact diagonalization approach.