We use very large cosmological N-body simulations to obtain accurate predictions for the two-point correlations and power spectra of mass-limited samples of galaxy clusters. We consider two currently popular cold dark matter (CDM) cosmogonies, a critical density model (tCDM) and a flat low density model with a cosmological constant (LCDM). Our simulations each use 10 particles to follow the mass distribution within cubes of side 2 h Gpc (tCDM) and 3 h Gpc (LCDM) with a force resolution better than 10 of the cube side. We investigate how the predicted cluster correlations increase for samples of increasing mass and decreasing abundance. Very similar behaviour is found in the two cases. The correlation length increases from r0 12±13 h Mpc for samples with mean separation dc 30 h Mpc to r0 22±27 h Mpc for samples with dc 100 h Mpc: The lower value here corresponds to tCDM and the upper to LCDM. The power spectra of these cluster samples are accurately parallel to those of the mass over more than a decade in scale. Both correlation lengths and power spectrum biases can be predicted to better than 10 per cent using the simple model of Sheth, Mo & Tormen. This prediction requires only the linear mass power spectrum and has no adjustable parameters. We compare our predictions with published results for the automated plate measurement (APM) cluster sample. The observed variation of correlation length with richness agrees well with the models, particularly for LCDM. The observed power spectrum (for a cluster sample of mean separation dc 31 h Mpc) lies significantly above the predictions of both models.