Radial velocity residuals time series (black dots), after subtracting our best-fit orbital solution for Proxima b. The blue line with grey shaded 1σ regions represents our best-fit GP quasi-periodic model for the correlated stellar noise. The upper plot shows the complete dataset, while the two plots in the second row show selected epochs, to easier visualize the agreement between the model and the data.