The robustness of cluster clustering covariance calibration

Ongoing and upcoming wide-field surveys at different wavelengths will measure the distribution of galaxy clusters with unprecedented precision, demanding accurate models for the two-point correlation function (2PCF) covariance. In this work we assessed a semi-analytical framework for the cluster 2PCF covariance that employs three nuisance parameters to account for non-Poissonian shot noise, residual uncertainties in the halo bias model, and sub-leading noise terms. We calibrated these parameters on a suite of fast approximate simulations generated by PINOCCHIO as well as full N-body simulations from OpenGADGET3. We demonstrate that PINOCCHIO can reproduce the 2PCF covariance measured in OpenGADGET3 at the few percent level, provided the mass functions are carefully rescaled. Resolution tests confirm that high particle counts are necessary to capture shot-noise corrections, especially at high redshifts. We performed the parameter calibration across multiple cosmological models, showing that one of the nuisance parameters, the non-Poissonian shot-noise correction (α), depends mildly on the amplitude of matter fluctuations (σ8). In contrast, the remaining two parameters, β, which controls the bias correction, and γ, which controls the secondary shot-noise correction, exhibit more significant variation with redshift and halo mass. Overall, our results underscore the importance of calibrating covariance models on realistic mock catalogs that replicate the selection function of forthcoming surveys and highlight that approximate methods, when properly tuned, can effectively complement full N-body simulations for precision cluster cosmology.