“…In this paper, we study the observational viability of different classes of interacting scenarios, including both linear and non-linear interaction terms. In order to observationally distinguish between the different models we perform a Bayesian model selection analysis using current data of type Ia supernova (SN Ia) [17,18], measurements of the baryon acoustic oscillation (BAO) from the 6dFGS [19], SDSS-MGS [20], BOSS-LOWZ [21], BOSS-CMASS [21], BOSS-DR12 [22], eBOSS [23,24] and BOSS-Lyα [25] surveys, along with BAO measurements obtained from SDSS-DR7 [26], SDSS-DR10 [27], SDSS-DR11 [28], SDSS-DR12Q [29]; measurements of the expansion rate from cosmic chronometers [30] and the estimate of the sound horizon scale at the last scattering reported by the Planck collaboration [31]. We find weak and moderate evidence against some of the interacting scenarios studied with respect to ΛCDM when the full joint analysis is considered.…”