Hello, I would appreciate very much some help on this problem.(adsbygoogle = window.adsbygoogle || []).push({});

I have a Hernquist model (distribution function) and I am trying to generate equal mass initial conditions. The way I do it is this. I have the mass profile M(r) and I pick say 10^6 random numbers between 0 and 1=M(r->inf). Then I solve the inverse equation r=r(m_random) to find the radius. The next step is to draw the velocities. For each radius I calculate f(v) = f(ri,v)=f(0.5*v^2+V(r))=f(E) and then I use a rejection method to draw the velocity from f(v). Then I divide the space r (radius) in a logarithmic scale and within each bin I calculate the <v^2> from all the v's within this bin therefore I have the function <v^2>(r) and draw it and compare it with <vr^2>(r) form analytic from of the paper.

THE PROBLEM IS that the two curves are THE SAME but they shouldn't because v^2=vr^2+vphi^2+vth^2 => <v^2>=3<vr^2> because the model is isotropic i.e. <vr^2>=<vphi^2>=<vth^2>. Therefore there should be a difference by a factor of 3 in the 2 curves but there isn't any difference.

Any ideas?

**Physics Forums - The Fusion of Science and Community**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Velocity dispersion problem

Loading...

Similar Threads - Velocity dispersion problem | Date |
---|---|

I Why use velocity dispersion in Faber Jackson relation | Jun 5, 2017 |

B Black hole mass and sigma (velocity dispersion) | Jan 11, 2017 |

A Velocity Dispersion | Nov 4, 2016 |

Velocity Dispersion | Dec 22, 2012 |

Velocity dispersion - confused. | May 14, 2012 |

**Physics Forums - The Fusion of Science and Community**