Sunday 25 March 2012

galaxy - Mass resolution - Astronomy

Usually, in an N-body or SPH simulation the term "mass resolution" refers to the mass of a single particle, which usually all have the same mass.



A single particle can always be "detected" in the simulation, since we have control of the coordinates of all particles, but a structure of several particles becomes ill-defined if the number of particles is too small



The mass of the smallest resolved structure depends on your definition of "resolved". That could be, say, 10 or 100 particles, depending on what quantity you are interested in measuring, and how accurately you want it. For instance, to define the mass of the structure, you need to be able to count the number of particles in the structure. But how does one know where to stop counting in a more or less continuous field of particles? One way is to determine an approximate center of mass $x_mathrm{CM}$ (approximate, since it can only be exact once all associated particles are defined), calculate the average density $langlerhorangle$ inside a sphere centered on $x_mathrm{CM}$ (which will be higher than the global average, since you started out on an overdensity), and increase the radius until $langlerhorangle$ falls below some threshold (e.g. 200 times the global average). If the number of particles is too small, $langlerhorangle$ will change a lot between each iteration, making your mass inaccurate. I think for this purpose, at least ten particles are needed.



If you're interested in the structure of the interstellar medium in a galaxy in a hydrodynamic (i.e. SPH) simulation, you probably need more particles than this. But I thinks it's fair to say that people disagree on how many particles are needed to resolve a galaxy.

No comments:

Post a Comment