Synthetic O – C diagrams obtained after performing forward numerical N-body simulations with TRADES (integration of 5.2 years). The C represents the timings calculated from the linear ephemerides in Table 2. The MAP model is plotted with a black line, while the confidence intervals at 1, 2, and 3σ are marked as shaded grey regions and come from 200 random samples drawn from the merged posterior distributions derived from the joint RV and TTV dynamical analysis.