This should be a straightforward exercise, but I'm struggling. Help is much appreciated

I would like to parameterize a Navarro-Frenk-White profile with one parameter, using the correlation between total mass and the concentration.

As a starting point, I have the density in terms of two parameters: the total dark matter mass inside the virial radius $r_{vir}$, and the concentration parameter $c=r_{vir}/r_{s}$:

$\rho_{DM}(r) = \frac{M_{DM}}{4\pi A(c)}\frac{1}{r(r_{s}+r)^{2}}$

where the function $A(c)$ is given by

$A(c) = \ln{(1+c)} - \frac{c}{1+c} $

Using the correlation between $M_{DM}$ and $c$ (from e.g. Napolitano et al. 2005)

$c(M_{DM}) \approx 9.195 \left(\frac{M_{DM}}{h^{-1}10^{12} M_{\odot}} \right)^{-0.094}$

it should then in principle be possible to parameterize $\rho_{DM}$ as a function of $M_{DM}$ only. This is how they do it in e.g. Williams et al. 2009.

However, I'm not managing to reproduce that. Sure, given the correlation above, $M_{DM}$ determines $A(c)$, but won't $\rho_{DM}$ still depend on $r_{s}$?

  • $\begingroup$ Can you provide the equation you end up getting and the equation from Williams et al. 2009 you're trying to match? $\endgroup$ – zephyr May 26 '17 at 14:38
  • $\begingroup$ What is $r_s$? Is it small enough to be neglected? $\endgroup$ – probably_someone May 26 '17 at 22:55
  • $\begingroup$ @zephyr, Williams et al. mention the equations that I state above allow them to rewrite it, but they don't give the rewritten expression. I've solved it by now actually, doing essentially what Walter below also says. Thanks! $\endgroup$ – user1991 May 27 '17 at 16:27
  • $\begingroup$ @probably_someone $r_{s}$ is the scale radius of the NFW profile, the radius at which the functional behaviour changes significantly. It is not negligible, no. Thanks! $\endgroup$ – user1991 May 27 '17 at 16:28

You are missing one last ingredient: the correlation between the virial radius $r_{\mathrm{vir}}$ and the virial mass $M_{\mathrm{vir}}$. The idea is that the orbital time scale of a particle orbiting in the halo becomes roughly equal to the age of the universe at $r_{\mathrm{vir}}$. This implies a relation between the (mean) density of the halo and of the universe: $$ \bar{\rho}_{\mathrm{vir}} \equiv \frac{3M_{\mathrm{vir}}}{4\pi r_{\mathrm{vir}}^3} = \Delta \rho_{\mathrm{crit}}(z) $$ where $\bar{\rho}_{\mathrm{vir}}$ is the halo's mean density interior of $r_{\mathrm{vir}}$, $\rho_{\mathrm{crit}}$ is the critical density of the universe (as function of redshift $z$), and $\Delta$ a parameter that depends on the cosmology. Often $\Delta=200$ is used and the corresponding quantities labeled $r_{200}$ and $M_{200}$ instead. See also this Wikipedia page.

So the way to construct a halo is as follows:

  1. Pick $M_{\mathrm{vir}}$ (or $M_{200}$) and redshift $z$ (and your favorite cosmology)
  2. Use above relation to find $r_{\mathrm{vir}}$ (or $r_{200}$)
  3. Use the relation from Napolitano et al. (2005), or equivalent, to find $c$ and hence $r_s$
  4. Obtain $\rho_{\mathrm{DM}}$ (if you still need it)
|improve this answer|||||

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.