CosmoCoffee

 FAQ   Search  SmartFeed   Memberlist    Register Profile   Log in Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

 Virial radius Goto page 1, 2  Next

Virial radius is defined. as the radius of a sphere that has a mean enclosed density of
180 times the mean density  23%  [ 3 ]
200 times the mean density  30%  [ 4 ]
180 times the critical density  15%  [ 2 ]
200 times the critical density  30%  [ 4 ]

Author Message
Anze Slosar

Joined: 24 Sep 2004
Posts: 198
Affiliation: Brookhaven National Laboratory

 Posted: October 01 2005 What should be the correct convention? It makes about 50% difference for the two most extreme cases...
Gil Holder

Joined: 25 Sep 2004
Posts: 27
Affiliation: McGill University

 Posted: October 02 2005 What happened to the old virial radius, defined through the spherical collapse model? That was the basis of the argument for using r180 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.
Tom Crawford

Joined: 29 Dec 2004
Posts: 9
Affiliation: University of Chicago

 Posted: October 03 2005 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 200 * ρc.
Ilian Iliev

Joined: 17 Aug 2006
Posts: 8
Affiliation: University of Sussex

 Posted: August 17 2006 Actually, the correct answer is probably "none of the above". The physically-meaningful 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 centrally-concentrated density profile and dependence on Lambda (see http://adsabs.harvard.edu/cgi-bin/nph-bib_query?bibcode=2001MNRAS.325..468I&db_key=AST&data_type=HTML&format=&high=3ecea4c09000496). 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).
Sarah Bridle

Joined: 24 Sep 2004
Posts: 147
Affiliation: University College London (UCL)

 Posted: August 18 2006 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.
Ilian Iliev

Joined: 17 Aug 2006
Posts: 8
Affiliation: University of Sussex

 Posted: August 19 2006 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. X-ray 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).
Sarah Bridle

Joined: 24 Sep 2004
Posts: 147
Affiliation: University College London (UCL)

 Posted: August 26 2006 Sheth & Tormen give their mass function (astro-ph/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!
Pirin Erdogdu

Joined: 17 Dec 2004
Posts: 7
Affiliation: UCL

 Posted: August 26 2006 Sheth and Tormen (1999) paper defines the critical value of the initial overdensity as δc = 1.68647 (page 3, equations 9 and 10). Thus, their δvir must be $\approx$180−200 (mean or critical density, it does not matter since this value is for an Einstein-de Sitter universe). Or?
Ilian Iliev

Joined: 17 Aug 2006
Posts: 8
Affiliation: University of Sussex

 Posted: August 27 2006 This is a rather complex question, which is often glossed over and not really discussed much. The definition of what is a halo in N-body simulations depends on the halo finder used. Two popular ones are the friends-of-friends (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 N-body simulations have. A note regarding the previous post - the critical density δc is not an initial overdensity, but is the value of the linear-theory overdensity reached when the nonlinear top-hat pertubation collapses to infinite density. The value δc=1.68 is correct for Einstein-de Sitter model, while for LambdaCDM it is slightly lower.
Pirin Erdogdu

Joined: 17 Dec 2004
Posts: 7
Affiliation: UCL

Posted: August 27 2006

Sorry for being a bit unclear in my previous post.

 Quote: 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 N-body simulations have.

Sheth and Tormen (ST99, 1999) find a fit to νf(ν) and ν≡[δc(z) / σ(m)]2
where (I quote from their paper) δc(z) is the critical value of the initial overdensity which is required to collapse at z, computed using the spherical collapse model and σ(m) is the value of the rms fluctuation in spheres which on average contain mass m at inital time, extrapolated using linear theory to z.

In section 3 of ST99, it reads:
 Quote: We measure mass functions in the simulations in the usual way, using spherical overdensity group finder (see Tormen 1998 for details).

And in Section 2.2 of Tormen (1998, MNRAS, 297, 648), it reads:
 Quote: We defined [dark matter] lumps by an overdensity criterion and included all particles within a sphere of mean overdensity δv = 178, centred on the particle with lowest potential energy. The value δv corresponds to the viral overdensity in the model of a spherical top-hat collapse in an Einstein- de Sitter universe.

It seems to me that they have a particular virial radius for the simulations (δv = 178). Any thoughts on this would be much appreciated.
Ilian Iliev

Joined: 17 Aug 2006
Posts: 8
Affiliation: University of Sussex

 Posted: August 28 2006 The value 178 (or more precisely 18π2) is the standard value coming from the classic top-hat 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.
Sarah Bridle

Joined: 24 Sep 2004
Posts: 147
Affiliation: University College London (UCL)

 Posted: August 28 2006 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 turn-around), but just want to know what "mass" is in the mass function I get from the Sheth-Tormen 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 (astro-ph/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 non-EdS 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.
Anne Green

Joined: 26 Sep 2004
Posts: 10
Affiliation: University Of Nottingham

 Posted: August 29 2006 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 (astro-ph/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.
Niayesh Afshordi

Joined: 17 Dec 2004
Posts: 62
Affiliation: Perimeter Institute/ University of Waterloo

Posted: August 29 2006

 Quote: 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.
Anne Green

Joined: 26 Sep 2004
Posts: 10
Affiliation: University Of Nottingham

Posted: August 29 2006

 Quote: 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.

I know/understand that. It doesn't solve the 'problem' though!

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 astro-ph/0506432 and astro-ph/0509897 carefully. I'd still be interested to hear what simulation experts think though.
 Display posts from previous: All Posts1 Day7 Days2 Weeks1 Month3 Months6 Months1 Year Oldest FirstNewest First
 All times are GMT + 5 HoursGoto page 1, 2  Next Page 1 of 2

 Jump to: Select a forum Arxiv paper discussion----------------arXiv papers Topic discussion----------------Early UniverseCosmological ModelCosmological Observations  Projects and Resources----------------Computers and softwareTeaching, Papers and PresentationsResearch projectsiCosmo Coming up----------------Job vacanciesConferences and meetings Management----------------CosmoCoffee
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum