ISCOs and OSCOs in the presence of a positive cosmological constant in massive gravity

We study the impact of a non-vanishing (positive) cosmological constant on the innermost and outermost stable circular orbits (ISCOs and OSCOs, respectively) within massive gravity in four dimensions. The gravitational field generated by a point-like object within this theory is known, generalizing the usual Schwarzschild--de Sitter geometry of General Relativity. In the non-relativistic limit, the gravitational potential differs by the one corresponding to the Schwarzschild--de Sitter geometry by a term that is linear in the radial coordinate with some prefactor $gamma$, which is the only free parameter. Starting from the geodesic equations for massive test particles and the corresponding effective potential, we obtain a polynomial of fifth order that allows us to compute the innermost and outermost stable circular orbits. Next, we numerically compute the real and positive roots of the polynomial for several different structures (from the hydrogen atom to stars and globular clusters to galaxies and galaxy clusters) considering three distinct values of the parameter $gamma$, determined using physical considerations, such as galaxy rotation curves and orbital precession. Similarly to the Kottler spacetime, both ISCOs and OSCOs appear. Their astrophysical relevance as well as the comparison with the Kottler spacetime are briefly discussed.