Testing cosmic anisotropy with cluster scaling relations

We test claims of large-scale anisotropy in the local expansion rate using cluster scaling relations as distance indicators. Using a Bayesian forward model, we jointly fit the X-ray luminosity--temperature (LT) and thermal Sunyaev-Zel'dovich--temperature (YT) relations, marginalising over the latent cluster distances and modelling selection effects as well as peculiar velocities. The latter are modelled using reconstructions of the local peculiar velocity field where we self-consistently account for possible anisotropic redshift--distance relations via an approximate scheme. This treatment proves crucial to the inferred anisotropy and breaks the degeneracy between anisotropy in scaling relation normalisations and underlying cosmological anisotropy. We apply our method to 312 clusters at $z \lesssim 0.2$, testing dipolar, quadrupolar and general (pixelised) anisotropy models. Bayesian model selection finds no more than weak evidence for any anisotropic model. For dipole models, we obtain upper limits of $δH_0 / H_0 < 3.2\%$ and bulk flow magnitude $< 1300\,\mathrm{km\,s^{-1}}$. Our results contrast with previous claims of statistically significant anisotropy from the same data, which we attribute to our principled forward modelling of both redshifts and scaling relation observables through latent distances and our treatment of the impact of anisotropic redshift--distance relations when modelling the local peculiar velocity field. Our work highlights the importance of accurately modelling peculiar velocities when testing isotropy with distance indicators, and motivates the further development of reconstructions that self-consistently treat large-scale deviations from the Hubble flow.