Interacting dark energy after DESI DR2: a challenge for $\Lambda$CDM paradigm?

We investigate the scenario of interacting dark energy through a detailed confrontation with various observational datasets. We quantify the interaction in a general way, through the deviation from the standard scaling of the dark matter energy density. We use observational Hubble Data from Cosmic Chronometers (CC), data from Baryon Acoustic Oscillations (BAO) from the recently released DESI DR2, data from Supernova Type Ia (SNIa), both with and without calibration by SH0ES, and finally the CMB data from Planck 2018. For the basic and simplest interacting model, we find that the data favor a non-zero interaction up to $2\sigma$. Notably, when we employ SH0ES-calibrated SNIa data the preference of non-zero interaction becomes more significant. However, comparison with $\Lambda$CDM scenario through AIC, BIC and Bayesian analysis, reveals a mixed picture.