Hi TNG team,
I would like to know the stellar gravitational potential in a galaxy at any coordinate, which is potential profile. Is there any easy way to get it? I can calculate with all stellar particles by -GM/r, but it is very time-consuming for computer. Besides, we have to set a zero point of the potential. How to set zero point in TNG?
A value of the gravitational potential is available in the Potential field for snapshots, but this (i) includes all matter, not just stars, and (ii) is a global potential, not a galaxy or halo local value.
If you want just stars, you need to calculate this yourself, as you say.
As for zero point, I would think that most likely you are interested in relative potential differences?
Actually I am figuring out how to calculate j_circ(E). I need to find a radius, where \Phi(r)+0.5*v_c^2(r)=E. For v_c, I can calculate it from \sqrt(GM(<r)/r), like what we did to get rotation curves. For \Phi, maybe we should estimate the potential profile. Or is there more easier and more standard way to calculate it?
Yes Phi is the gravitational potential. I think you can ask Ling, or Shy, for details about how to calculate j_circ(E). You could compare your results to existing results under the "stellar circularities" catalog.