Virial radius

 Posts: 183
 Joined: September 24 2004
 Affiliation: Brookhaven National Laboratory
 Contact:
Virial radius
What should be the correct convention?
It makes about 50% difference for the two most extreme cases...
It makes about 50% difference for the two most extreme cases...

 Posts: 27
 Joined: September 25 2004
 Affiliation: McGill University
Virial radius
What happened to the old virial radius, defined through the spherical collapse model? That was the basis of the argument for using r_{180} for omega_m=1, but you end up with something close to 100 times the critical density for omega_m=0.3. It doesn't seem to have had much success in defining mass functions, but it always seemed like the one with the best physical foundation.

 Posts: 8
 Joined: December 29 2004
 Affiliation: University of Chicago
Virial radius
As long as people agree on a consistent definition, the exact choice seems arbitrary  until you want to compare with observations, at which point you want a quantity whose definition doesn't depend on other parameters (or depends on as few as possible). And if we're not going to tie the quantity directly to the spherical collapse overdensity (because that does vary with cosmology), then we might as well use a nice round number (rather than 178 or 180). Thus my vote for [tex]200*\rho_c[/tex].

 Posts: 8
 Joined: August 17 2006
 Affiliation: University of Sussex
 Contact:
Virial radius
Actually, the correct answer is probably "none of the above". The physicallymeaningful mean density for halos is (roughly) proportional to the mean density at its time of formation. We have written some time ago a paper where we generalized the simple tophat model to something more physical, including centrallyconcentrated density profile and dependence on Lambda (see http://adsabs.harvard.edu/cgibin/nphb ... 4c09000496). This gave mean halo overdensity of ~130, weakly varying with Lambda. Numerical simulations also give more extended halos and lower mean density (see e.g. recent papers by Anatoly Klypin).

 Posts: 144
 Joined: September 24 2004
 Affiliation: University College London (UCL)
 Contact:
Virial radius
It is a bit confusing that these definitions are referred to as the "virial radius", when indeed Ilian and Tom point out the two opposing uses: "something related to virialisation" vs "a convenient definition of the mass of a halo".
Maybe we should use the term e.g. "200critical radius" for the latter instead of "virial radius"?
I agree with Tom that it is handy to use a round number, and good if it doesn't change too much with cosmology.
Maybe we should use the term e.g. "200critical radius" for the latter instead of "virial radius"?
I agree with Tom that it is handy to use a round number, and good if it doesn't change too much with cosmology.

 Posts: 8
 Joined: August 17 2006
 Affiliation: University of Sussex
 Contact:
Virial radius
The appropriate definition apparently depends on the problem. If one needs merely to locate a halo and give it some rough mass, then almost any resonable definition would do. For things that depend on the inner halo parts (e.g. Xray emission) even higher overdensities would be approporiate. But if one wants to know what mass is virialized and the total mass of a halo, then the boundary moves quite a bit out and the overdensity is lower (maybe ~100 or even less, this is still not settled).

 Posts: 144
 Joined: September 24 2004
 Affiliation: University College London (UCL)
 Contact:
Virial radius
Sheth & Tormen give their mass function (astroph/9901122), but I can't see from their paper how they define the mass of a halo. Am I missing it?
Jenkins et al. mention "The mean overdensity was set to ... 324 in the ... LambdaCDM" which sounds promising, but still doesn't seem to contain all the information.
In particular, is this number relative to the critical density or to the mean density? (Should it be obvious? 324 seems v different to the numbers Ilian is giving..)
And do Sheth & Tormen define halo mass the same as Jenkins et al anyway? Thanks very much for any clues as I would like to be able to use the Sheth & Tormen mass function correctly!
Jenkins et al. mention "The mean overdensity was set to ... 324 in the ... LambdaCDM" which sounds promising, but still doesn't seem to contain all the information.
In particular, is this number relative to the critical density or to the mean density? (Should it be obvious? 324 seems v different to the numbers Ilian is giving..)
And do Sheth & Tormen define halo mass the same as Jenkins et al anyway? Thanks very much for any clues as I would like to be able to use the Sheth & Tormen mass function correctly!

 Posts: 5
 Joined: December 17 2004
 Affiliation: UCL
Virial radius
Sheth and Tormen (1999) paper defines the critical value of the initial overdensity as [tex]\delta_c=1.68647[/tex] (page 3, equations 9 and 10). Thus, their [tex]\delta_{vir}[/tex] must be [tex]\approx[/tex]180200 (mean or critical density, it does not matter since this value is for an Einsteinde Sitter universe). Or?

 Posts: 8
 Joined: August 17 2006
 Affiliation: University of Sussex
 Contact:
Virial radius
This is a rather complex question, which is often glossed over and not really discussed much. The definition of what is a halo in Nbody simulations depends on the halo finder used. Two popular ones are the friendsoffriends (FOF) and spherical overdensity, but there are also others. Each of these has free parameters, the linking length in the first case and the overdensity in the second. The overdensity is of course what we are discussing here, while the linking length (usually chosen to be 0.2) is only loosely connected to a particular overdensity. The value of the free parameter is chosen rather empirically so as to give reasonable results for the halos. The halo mass depends on the value chosen, so there is obviously a certain freedom here. It should also be noted that while these methods agree fairly well neither method finds the halos perfectly and each has its difficulties in some situations.
If I am not mistaken, the Sheth and Tormen model does not point to a particular overdensity, but is just calibrated against numerical simulations, so has much the same freedoms the Nbody simulations have.
A note regarding the previous post  the critical density δc is not an initial overdensity, but is the value of the lineartheory overdensity reached when the nonlinear tophat pertubation collapses to infinite density. The value δc=1.68 is correct for Einsteinde Sitter model, while for LambdaCDM it is slightly lower.
If I am not mistaken, the Sheth and Tormen model does not point to a particular overdensity, but is just calibrated against numerical simulations, so has much the same freedoms the Nbody simulations have.
A note regarding the previous post  the critical density δc is not an initial overdensity, but is the value of the lineartheory overdensity reached when the nonlinear tophat pertubation collapses to infinite density. The value δc=1.68 is correct for Einsteinde Sitter model, while for LambdaCDM it is slightly lower.

 Posts: 5
 Joined: December 17 2004
 Affiliation: UCL
Virial radius
Sorry for being a bit unclear in my previous post.
Ilian's post reads:
where (I quote from their paper) [tex]\delta_c(z)[/tex] is the critical value of the initial overdensity which is required to collapse at [tex]z[/tex], computed using the spherical collapse model and [tex]\sigma(m)[/tex] is the value of the rms fluctuation in spheres which on average contain mass [tex]m[/tex] at inital time, extrapolated using linear theory to [tex]z[/tex].
In section 3 of ST99, it reads:
Ilian's post reads:
Sheth and Tormen (ST99, 1999) find a fit to [tex]\nu f(\nu)[/tex] and [tex]\nu\equiv[\delta_c(z)/\sigma(m)]^2[/tex]If I am not mistaken, the Sheth and Tormen model does not point to a particular overdensity, but is just calibrated against numerical simulations, so has much the same freedoms the Nbody simulations have.
where (I quote from their paper) [tex]\delta_c(z)[/tex] is the critical value of the initial overdensity which is required to collapse at [tex]z[/tex], computed using the spherical collapse model and [tex]\sigma(m)[/tex] is the value of the rms fluctuation in spheres which on average contain mass [tex]m[/tex] at inital time, extrapolated using linear theory to [tex]z[/tex].
In section 3 of ST99, it reads:
And in Section 2.2 of Tormen (1998, MNRAS, 297, 648), it reads:We measure mass functions in the simulations in the usual way, using spherical overdensity group finder (see Tormen 1998 for details).
It seems to me that they have a particular virial radius for the simulations ([tex]\delta_v=178[/tex]). Any thoughts on this would be much appreciated.We defined [dark matter] lumps by an overdensity criterion and included all particles within a sphere of mean overdensity [tex]\delta_v=178[/tex], centred on the particle with lowest potential energy. The value [tex]\delta_v[/tex] corresponds to the viral overdensity in the model of a spherical tophat collapse in an Einstein de Sitter universe.

 Posts: 8
 Joined: August 17 2006
 Affiliation: University of Sussex
 Contact:
Virial radius
The value 178 (or more precisely 18[tex]\pi^2[/tex]) is the standard value coming from the classic tophat model. As far as I am aware, simulations do not point specifically to this value. In the spherical overdensity halo finder it is a free parameter, which could be varied (within some limits). If one picks 200 or 130, this would be just as reasonable. As I mentioned before, some newer simulations that were specifically looking for the halo boundaries (virialization shocks in gas/caustics in dark matter) were finding rather smaller overdensities, maybe ~100 or less.

 Posts: 144
 Joined: September 24 2004
 Affiliation: University College London (UCL)
 Contact:
Virial radius
Thanks very much indeed for your replies. To be honest, at the moment I'm not too bothered about the definition of the virial radius (as defined by turnaround), but just want to know what "mass" is in the mass function I get from the ShethTormen equations (i.e. once they had selected a halo in their simulations, out to what radius did they sum up their mass?).
Thanks for the Tormen (astroph/9802290) ref, which suggests use of 178 in the EdS simulations in that paper. Does anyone know if 178 was also used in Sheth & Tormen?
If indeed Sheth & Tormen use 178 then they must have had to make a decision when applying this to the nonEdS simulations they use in their paper: is this 178 relative to the mean density at the halo redshift, or the critical density at the halo redshift?
Thanks very much for any more comments.
Thanks for the Tormen (astroph/9802290) ref, which suggests use of 178 in the EdS simulations in that paper. Does anyone know if 178 was also used in Sheth & Tormen?
If indeed Sheth & Tormen use 178 then they must have had to make a decision when applying this to the nonEdS simulations they use in their paper: is this 178 relative to the mean density at the halo redshift, or the critical density at the halo redshift?
Thanks very much for any more comments.

 Posts: 17
 Joined: September 26 2004
 Affiliation: University Of Nottingham
 Contact:
Virial radius
I'm not an expert, but I think that mass is usually taken to be the mass within the virial radius, so it does boil down to what definition is used for the virial radius.
Bryan and Norman (astroph/9710107) have fitting functions (which are often cited) for the dependence of the virial over density, defined in terms of the collapse of a spherical perturbation, on cosmology (eq. 6).
A related problem that's been bothering me lately is the time dependence of a halo's properties (radius and mass) if you use the "standard" overdensity of 200; as the Universe evolves the background density
(total or matter) decreases so that the virial radius, and hence the nominal mass, of a given halo increases even if it remains physically unchanged. I'd be interested to know if anyone has any thoughts on this problem.
Bryan and Norman (astroph/9710107) have fitting functions (which are often cited) for the dependence of the virial over density, defined in terms of the collapse of a spherical perturbation, on cosmology (eq. 6).
A related problem that's been bothering me lately is the time dependence of a halo's properties (radius and mass) if you use the "standard" overdensity of 200; as the Universe evolves the background density
(total or matter) decreases so that the virial radius, and hence the nominal mass, of a given halo increases even if it remains physically unchanged. I'd be interested to know if anyone has any thoughts on this problem.

 Posts: 49
 Joined: December 17 2004
 Affiliation: Perimeter Institute/ University of Waterloo
 Contact:
Virial radius
A related problem that's been bothering me lately is the time dependence of a halo's properties (radius and mass) if you use the "standard" overdensity of 200; as the Universe evolves the background density
(total or matter) decreases so that the virial radius, and hence the nominal mass, of a given halo increases even if it remains physically unchanged. I'd be interested to know if anyone has any thoughts on this problem.
This is correct. At the same time the NFW concentration parameter increases, so that the profile remains the same, although the mass is increasing. This is also why smaller haloes that form earlier have a larger concentration parameter.

 Posts: 17
 Joined: September 26 2004
 Affiliation: University Of Nottingham
 Contact:
Virial radius
I know/understand that. It doesn't solve the 'problem' though!At the same time the NFW concentration parameter increases, so that the profile remains the same, although the mass is increasing. This is also why smaller haloes that form earlier have a larger concentration parameter.
Putting it another way: is there really a sharp transition at the virial radius between bound particles and particles which are being accreted? I think I need to study astroph/0506432 and astroph/0509897 carefully. I'd still be interested to hear what simulation experts think though.