Phase-folded plots. Top panel: global fit result for TESS and the ground-based transits. Both LCOGT light curves are shifted on the y-axis for clarity, and their respective filter band are indicated in the legend, while the superimposed points represent ~30-minute bins. Bottom panels: phase-folded HARPS-N RVs to the period of planet b and candidate c, along with their residuals. The red circles represent the average value of ~10 phased RV data points at a time.