Benchmarks and explanations for deep learning estimates of X-ray galaxy cluster masses

We evaluate the effectiveness of deep learning (DL) models for reconstructing the masses of galaxy clusters using X-ray photometry data from next-generation surveys. We establish these constraints using a catalogue of realistic mock eROSITA X-ray observations which use hydrodynamical simulations to model realistic cluster morphology, background emission, telescope response, and active galactic nucleus (AGN) sources. Using bolometric X-ray photon maps as input, DL models achieve a predictive mass scatter of $\sigma _{\ln M_\mathrm{500c}} = 17.8~{{\ \rm per\ cent}}$, a factor of two improvements on scalar observables such as richness N_gal, 1D velocity dispersion σ_v,1D, and photon count N_phot as well as a 32 per cent improvement upon idealized, volume-integrated measurements of the bolometric X-ray luminosity L_X. We then show that extending this model to handle multichannel X-ray photon maps, separated in low, medium, and high energy bands, further reduces the mass scatter to 16.2 per cent. We also tested a multimodal DL model incorporating both dynamical and X-ray cluster probes and achieved marginal gains at a mass scatter of 15.9 per cent. Finally, we conduct a quantitative interpretability study of our DL models and find that they greatly down-weight the importance of pixels in the centres of clusters and at the location of AGN sources, validating previous claims of DL modelling improvements and suggesting practical and theoretical benefits for using DL in X-ray mass inference.