In a recent paper we applied a rigorous perturbed matching framework to show the amendment of the mass of rotating stars in Hartle's model. Here, we apply this framework to the tidal problem in binary systems. Our approach fully accounts for the correction to the Love numbers needed to obtain the universal $I$-Love-$Q$ relations. We compute the corrected mass vs radius configurations of rotating quark stars, revisiting a classical paper on the subject. These corrections allow us to find a universal relation involving the second-order contribution to the mass $\delta M$. We thus complete the set of universal relations for the tidal problem in binary systems, involving four perturbation parameters, namely $I$, Love, $Q$, and $\delta M$. These relations can be used to obtain the perturbation parameters directly from observational data.

We show how the correction to the calculation of the mass in the original relativistic model of a rotating star by Hartle [6], found recently [10], appears in the Newtonian limit, and that the correcting term is indeed present, albeit hidden, in the original Newtonian approach by Chandrasekhar [2].