Ferrers potential¶

class
galpy.potential.
FerrersPotential
(amp=1.0, a=1.0, n=2, b=0.35, c=0.2375, omegab=0.0, pa=0.0, normalize=False, ro=None, vo=None)[source]¶ Class that implements triaxial Ferrers potential for the ellipsoidal density profile with the short axis along the zdirection
\[\rho(x,y,z) = \frac{\mathrm{amp}}{\pi^{1.5} a^3 b c} \frac{\Gamma(n+\frac{5}{2})}{\Gamma(n+1)}\,(1(m/a)^2)^n\]with
\[m^2 = x'^2 + \frac{y'^2}{b^2}+\frac{z'^2}{c^2}\]and \((x',y',z')\) is a rotated frame wrt \((x,y,z)\) so that the major axis is aligned with \(x'\).
Note that this potential has not yet been optimized for speed and has no C implementation, so orbit integration is currently slow.

__init__
(amp=1.0, a=1.0, n=2, b=0.35, c=0.2375, omegab=0.0, pa=0.0, normalize=False, ro=None, vo=None)[source]¶ NAME:
__init__PURPOSE:
initialize a Ferrers potentialINPUT:
amp  total mass of the ellipsoid determines the amplitude of the potential; can be a Quantity with units of mass or Gxmass
a  scale radius (can be Quantity)
n  power of Ferrers density (n > 0)
b  ytox axis ratio of the density
c  ztox axis ratio of the density
omegab  rotation speed of the ellipsoid (can be Quantity)
pa= (None) If set, the position angle of the x axis (rad or Quantity)
normalize  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=, vo= distance and velocity scales for translation into internal units (default from configuration file)
OUTPUT:
(none)
