The coupled cluster method (CCM) is a method of quantum many-body theory that may provide accurate results for the ground-state properties of lattice quantum spin systems even in the presence of strong frustration and for lattices of arbitrary spatial dimensionality. Here we present a significant extension of the method by introducing a new approach that allows an efficient parallelization of computer codes that carry out "high-order" CCM calculations. We find that we are able to extend such CCM calculations by an order of magnitude higher than ever before utilized in a high-order CCM calculation for an antiferromagnet. Furthermore, we use only a relatively modest number of processors, namely, eight. Such very high-order CCM calculations are possible only by using such a parallelized approach. An illustration of the new approach is presented for the ground-state properties of a highly frustrated two-dimensional magnetic material, CaV4O9. Our best results for the ground-state energy and sublattice magnetization for the pure nearest-neighbor model are given by Eg/N = −0.5534 and M = 0.19, respectively, and we predict that there is no Néel ordering in the region 0.2 ≤ J2/J1 ≤ 0.7. These results are shown to be in excellent agreement with the best results of other approximate methods.A new procedure for carrying out high-order coupled cluster method (CCM) [1,2,3,4,5,6,7,8] calculations via parallel processing is presented here. The CCM may be applied to systems demonstrating strong frustration and for arbitrary spatial dimension of the lattice. We illustrate our new approach by applying it to CaV 4 O 9 [9,10,11,12,13,14,15] at zero temperature. An increasing number of insulating quantum magnetic systems for lattices of low spatial dimensionality are being studied experimentally. Indeed, the calcium vanadium oxide (CAVO) materials are one particularly useful example. They exhibit strong frustration for a number of different crystallographic lattices with respect to varying chemical composition, and they may demonstrate "novel" groundstate ordering. Indeed, theoretical evidence suggests that CAVO systems may contain a number of differing quantum ground states (see, e.g., Refs. [9,11]) at zero temperature as a function of varying nearest-neighbor nextnearest-neighbor bond strengths. The relevant Hamiltonian is given bywhere i runs over all lattice sites, and j and k run over all nearest-neighbor and next-nearest-neighbor sites to i, respectively, counting each bond once and once only. The bond strengths are given by J 1 and J 2 for nearestneighbor and next-nearest-neighbor terms, respectively. The lattice and exchange "bonds" are illustrated graphically in Fig. 1. For this model, collinear Néel ordering is observed with respect to nearest neighbors at J 1 > 0 and J 2 = 0 and with respect to next-nearest neighbors at J 2 > 0 and diate phase (or phases) for 0.2 < J 2 /J 1 < 0.7, which is also in good agreement with results of a "spin-wave theory-like" treatment [9] that predict such an intermediate regime for 0.25 < J 2 ...