RV analysis performed with pyaneti. Top-left: Time series of the HARPS-N TERRA RV measurements (blue points). The inferred models for the five-planet and GP signals are shown with red and blue thick lines, respectively, whereas the combined signal is displayed in black. Topright: Time series of the FWHM of the DRS cross-correlation functions (blue points). The black thick line marks the inferred GP model, whereas the grey shaded areas show its 1σ and 2σ credible intervals. Second and third rows, from top to bottom and left to right: Phase-folded RV curves of TOI-5624 b, c, d, e, and f, and best-fitting Keplerian models (thick black lines). In each panel, the nominal uncertainties of the data points and the jitter contribution are shown with thin and thick blue lines, respectively.