the profiles module¶

pNbody.profiles.burkert_mr(r, rs, rho0=1.0)¶

Burkert profile rhob = rho0 / ( ( 1 + r/rs ) * ( 1 + (r/rs)**2 ) )

Burkert, Astrophys. J. 447 (1995) L25.

pNbody.profiles.burkert_profile(r, rs, rho0=1.0)¶

Burkert profile rhob = rho0 / ( ( 1 + r/rs ) * ( 1 + (r/rs)**2 ) )

Burkert, Astrophys. J. 447 (1995) L25.

pNbody.profiles.generic2c_mr(r, rs, a, b, rho0=1.0)¶: Mass in the radius r for the distribution rho = 1/( (r/rs)**a * (1+r/rs)**(b-a) )

pNbody.profiles.generic2c_profile(r, rs, a, b, rho0=1.0)¶: generic2c profile rho = 1/( (r/rs)**a * (1+r/rs)**(b-a) )

pNbody.profiles.hernquist_mR(R, rs, rho0=1)¶

Mass in the projected radius R for the distribution rho = 1/( (r/rs) * (1+r/rs)**3 )

(Hernquist 90, Eq. 37)

Warning : the function diverges in r=0 and r/rs=1. Warning : it is badly implemented for arrays

pNbody.profiles.hernquist_mr(r, rs, rho0=1.0)¶: Mass in the radius r for the distribution rho = 1/( (r/rs) * (1+r/rs)**3 )

pNbody.profiles.hernquist_profile(r, rs, rho0=1.0)¶: hernquist profile rho = 1/( (r/rs) * (1+r/rs)**3 )

pNbody.profiles.jaffe_mr(r, rs, rho0=1.0)¶: Mass in the radius r for the distribution rho = 1/( (r/rs)**2 * (1+r/rs)**2 )

pNbody.profiles.jaffe_profile(r, rs, rho0=1.0)¶: jaffe profile rho = 1/( (r/rs)**2 * (1+r/rs)**2 )

pNbody.profiles.king_Rc(rs, rt)¶

Core radius Find R such that

Sigma(Rc) = Sigma(0)/2.

pNbody.profiles.king_profile(r, rs, rt)¶: King profile (see King 62)

pNbody.profiles.king_profile_Rz(R, z, rs, rt)¶: King profile in cylindrical coord (needed for surface density computation) (see King 62)

pNbody.profiles.king_surface_density(R, rs, rt)¶: Surface density of King profile (see King 62)

pNbody.profiles.king_surface_density_old(R, rs, rt)¶: Obsolete implementation

pNbody.profiles.nfw_mr(r, rs, rho0=1.0)¶: Mass in the radius r for the distribution rho = rho0/((r/rs)*(1+r/rs)**2)

pNbody.profiles.nfw_profile(r, rs, rho0=1.0)¶: NFW profile rho = rho0/((r/rs)*(1+r/rs)**2)

pNbody.profiles.nfwg_mr(r, rs, gamma, rho0=1.0)¶: Mass in the radius r for the distribution rho = rho0/((r/rs)**(gamma)*(1+(r/rs)**2)**(0.5*(3.-gamma)))

pNbody.profiles.nfwg_profile(r, rs, gamma, rho0=1.0)¶: NFW modified profile rho = rho0/((r/rs)**(gamma)*(1+(r/rs)**2)**(0.5*(3.-gamma)))

pNbody.profiles.nfws_mr(r, rhos, rs, r0)¶: Mass in the radius r for the distribution rho = rhos/((r/rs)*(1+r/rs)**2)

pNbody.profiles.nfws_profile(r, rhos, rs, r0)¶: NFW softened profile rho = rhos/(((r+r0)/rs)*(1+r/rs)**2)

pNbody.profiles.pisothm_mr(r, rs, rho0=1.0)¶: Mass in the radius r for the distribution rho = 1/(1+(r/rs)**2)

pNbody.profiles.pisothm_profile(r, rs, rho0=1.0)¶: Pseudo-isothermal profile rho = 1/(1+(r/rs)**2)

pNbody.profiles.plummer_mr(r, rc, rho0=1.0)¶: Mass in the radius r for the distribution rho = 1/(1+(r/rc)**2)**(5/2)

pNbody.profiles.plummer_profile(r, rc, rho0=1.0)¶: Plummer profile rho = 1/(1+(r/rc)**2)**(5/2)