We study the gravitational potential arising from static, spherically symmetric sources in a modified gravity theory described by a quadratic f (R) model. The linearized field equations lead to a fourth-order differential equation, whose solution includes Newtonian and Yukawa-like terms. Imposing regularity at the origin and asymptotic flatness fixes all integration constants, making the potential fully determined by the mass distribution. We analyze classical profiles such as Plummer, Hernquist, and Navarro-Frenk-White (NFW), along with some new analytic models. For each case, we derive explicit expressions for the modified potential and analyze its dependence on the f (R) parameter α. The solutions demonstrate that while f (R) corrections generically introduce r^-1 divergences at finite α, these are cancelled in the GR limit, that is, when the scalar field mass becomes infinite. The new profiles offer alternatives for modeling astrophysical systems, particularly in regimes where deviations from Newtonian gravity may play a role. These results may serve as a theoretical groundwork for future studies of modified gravity in astrophysical contexts.
Roger A. Hurtado (Wed,) studied this question.