Résumé
We analyse pre-DESI clustering data using a dark energy equation of state$w(z)$parametrised by$(w_0, w_a)$ , finding a$2.8-3.9\sigma$preference for evolving dark energy over the cosmological constant$\Lambda$when combined with cosmic microwave background data from Planck and supernova data from Pantheon+, Union3, or DESY5. Our constraints, consistent with DESI Y1 results, are derived from the power spectrum and bispectrum of SDSS/BOSS galaxies using the Effective Field Theory of Large Scale Structure (EFTofLSS) at one loop. The evidence remains robust across analysis variations but disappears without the one-loop bispectrum. When combining DESI baryon acoustic oscillations with BOSS full-shape data, while marginalising over the sound horizon in the latter to prevent unaccounted correlations, the significance increases to$3.7-4.4\sigma$ , depending on the supernova dataset. Using a data-driven reconstruction of$w(z)$ , we show that the evidence arises from deviations from$\Lambda$at multiple redshifts. In addition, our findings are interpreted within the Effective Field Theory of Dark Energy (EFTofDE), from which we explicitly track the non-standard time evolution in EFTofLSS predictions. For perturbatively stable theories in the$w < -1$regime, the evidence persists in the clustering limit$(c_s^2 \rightarrow 0)$when higher-derivative corrections are present, as well as in the quasi-static limit$(c_s^2 \rightarrow 1)$when additional EFTofDE parameters are considered.