Power-law density spherical potential with an exponential cut-off

class galpy.potential.PowerSphericalPotentialwCutoff(amp=1.0, alpha=1.0, rc=1.0, normalize=False, r1=1.0, ro=None, vo=None)[source]

Class that implements spherical potentials that are derived from power-law density models

\[\rho(r) = \mathrm{amp}\,\left(\frac{r_1}{r}\right)^\alpha\,\exp\left(-(r/rc)^2\right)\]
__init__(amp=1.0, alpha=1.0, rc=1.0, normalize=False, r1=1.0, ro=None, vo=None)[source]

Initialize a power-law-density potential.

Parameters:

  • amp (float or Quantity, optional) – Amplitude to be applied to the potential. Can be a Quantity with units of mass density or Gxmass density.

  • alpha (float, optional) – Inner power.

  • rc (float or Quantity, optional) – Cut-off radius.

  • r1 (float or Quantity, optional) – Reference radius for amplitude. Default is 1.0. Can be Quantity.

  • normalize (bool or float, optional) – If True, normalize such that vc(1.,0.)=1., or, if given as a number, such that the force is this fraction of the force necessary to make vc(1.,0.)=1.

  • ro (float, optional) – Distance scale for translation into internal units (default from configuration file).

  • vo (float, optional) – Velocity scale for translation into internal units (default from configuration file).

Notes

  • 2013-06-28 - Written - Bovy (IAS)