Difference between the stellar radial velocity (computed with the n-body integration) and a two-Keplerian fit to the model prediction, computed over 100 random MCMC samples from the full photodynamic model. The black line and the grey shaded region represent the median and the 68.3% confidence interval, respectively. The smaller panel is an enlargementaround the maximum difference, marked in red in the main panel.