ZhaoProfile#

class gammapy.astro.darkmatter.ZhaoProfile[source]#

Bases: DMProfile

Zhao Profile.

It is a generalization of the NFW profile. The volume density is parametrized with a double power-law. Scale radii smaller than the scale radius are described with a slope of \(-\gamma\) and scale radii larger than the scale radius are described with a slope of \(-\beta\). \(\alpha\) is a measure for the width of the transition region.

\[\rho(r) = \rho_s \left(\frac{r_s}{r}\right)^\gamma \left(1 + \left(\frac{r}{r_s}\right)^\frac{1}{\alpha} \right)^{(\gamma - \beta) \alpha}\]

Equation (1) from [1].

Parameters:
r_sQuantity

Scale radius, \(r_s\).

alphaQuantity

\(\alpha\).

beta: `~astropy.units.Quantity`

\(\beta\).

gammaQuantity

\(\gamma\).

rho_sQuantity

Characteristic density, \(\rho_s\).

References

Attributes Summary

DEFAULT_ALPHA

DEFAULT_BETA

DEFAULT_GAMMA

(alpha, beta, gamma) = (1,3,1) is NFW profile.

DEFAULT_SCALE_RADIUS

Methods Summary

evaluate(radius, r_s, alpha, beta, gamma, rho_s)

Evaluate the profile.

Attributes Documentation

DEFAULT_ALPHA = 1#
DEFAULT_BETA = 3#
DEFAULT_GAMMA = 1#

(alpha, beta, gamma) = (1,3,1) is NFW profile. Default scale radius as given in reference 2 (same as for NFW profile)

DEFAULT_SCALE_RADIUS = <Quantity 24.42 kpc>#

Methods Documentation

static evaluate(radius, r_s, alpha, beta, gamma, rho_s)[source]#

Evaluate the profile.

__init__(r_s=None, alpha=None, beta=None, gamma=None, rho_s=<Quantity 1. GeV / cm3>)[source]#
classmethod __new__(*args, **kwargs)#