Detailed first-order post-Newtonian analysis of massive Brans-Dicke theories: Numerical constraints and the meaning of the $\beta$ parameter

Massive Brans-Dicke (BD) theory is among the simplest general relativity extensions. It is commonly found as the weak-field limit of other gravitational theories. Here we do a detailed post-Newtonian analysis of massive BD theories. We start by expanding the massive BD field equations following the Will-Nordtvedt parametrized post-Newtonian (PPN) formalism, without point-particle approximations. A single potential that is not present in the standard PPN formalism is found. This new potential hinders immediate PPN conclusions. To proceed, we do a complete first-order post-Newtonian analysis and explicitly derive all the conserved quantities. After demanding that there exists a Newtonian limit by requiring the BD mass to be sufficiently large, we find, as expected, that $\gamma =1$; but there is no effective $\beta$ parameter that can have the same physical role of the standard $\beta$ in PPN formalism. All the others standard PPN parameters can be extended to the massive BD case without issues and are shown to have the same values of general relativity. At last, we consider numerical relations on the periastron advance and the BD mass in two different physical contexts, the orbit of Mercury about the Sun and the orbit of the star S2 about the expected supermassive black hole in the Milky Way.