Hello ! I have some problem in understanding how you calculate your virial Radius.
I'm currently working with Illustris 3 at redshift 0.
I retrieve the virial mass (Tophat200) and i calculate the virial radius with the following code :
Mvir = 4/3piRvir*3 deltac * pc with deltac calculated with the formula from Bryan and Norman (1998)).
Om = WMAP9.Om(z)
E = WMAP9.efunc(z)
rhoc = WMAP9.critical_density(z)
deltac = 18pi**2+82(1-Om)-(39(1-Om)*2)
radius = (M/(deltacrhoc(4/3)pi))*(1./3)
for i in range(M.size):
For example for M(TOPHAT200) = 3.63 * 1e14 solar mass i found a virial radius of approx 4.5 Mpc
for the same halo however, R(TOPHAT200) = 1310 ckpc/h = 1.8 Mpc at z = 0
Even if it not a fiducial cosmology it shouldn't lead to such difference right ? so i was wondering where i was wrong in my calculation
And a last question : The actual total mass is always greater than the tophat mass. Does it means that our actual mean density is greater than deltac ?
So does it mean so the group is already virialized since a given time t ?
Have a great day and thanks you for reading this !!
my bad i have managed to calculate it afterward