RVs and photometry light curves of TOI-5624

By using TESS and CHEOPS photometry along with HARPS-N and SOPHIE high-resolution spectra, we determined the planet properties and performed a dynamical analysis of the system.

After analysing the photometric data, we extracted and modelled the radial velocity (RV) time series using two independent methodologies both within a Markov chain Monte Carlo framework. We further integrated the N-body equations of motion, while simultaneously fitting the transit times and the detrended RVs, to dynamically characterise the system. We present the discovery of four transiting sub-Neptunes with radii of R_b_=2.314+/-0.035R_{Earth}_, R_c_=2.474+/-0.042R_{Earth}_, R_d_=3.584_-0.050_^+0.051^ R_{Earth}_, and R_e_=3.247_-0.043_^+0.042^ R_{Earth}_ and masses of M_b_=9.4+/-1.4M_{Earth}_, M_c_=4.8+/-1.9M_{Earth}, M_d_=4.9+/-2.2M_{Earth}_, and M_e_=8.9_-3.0_^+2.9^M_{Earth}_.