Potential-density pairs for axisymmetric galaxies: the influence of   scalar fields
  We present a formulation for potential-density pairs to describe axisymmetric galaxies in the Newtonian limit of scalar-tensor theories of gravity. The scalar field is described by a modified Helmholtz equation with a source that is coupled to the standard Poisson equation of Newtonian gravity. The net gravitational force is given by two contributions: the standard Newtonian potential plus a term stemming from massive scalar fields. General solutions have been found for axisymmetric systems and the multipole expansion of the Yukawa potential is given. In particular, we have computed potential-density pairs of galactic disks for an exponential profile and their rotation curves. 
M. A. Rodriguez-Meza, Jorge L. Cervantes-Cota, M. I. Pedraza, J. F.   Tlapanco, and E. M. De la Calleja