| 110 | Now, this is where things get hard for me to understand. First, I see the range of phi, and magnitude of phi is drastically different for this matched ambient case. This may make sense, based on more mass in the ambient meaning stronger gravitational potential due to the ambient. However, I would expect that at very small radii, just above Rbe, the approximation of the ambient as a uniform sphere should be very good. This is because there is SOOOO much more mass in the ambient than in the sphere (see above calculation). Thus, very early on, I would expect the potential to be growing more and more negative as more and more mass is enclosed in this uniform sphere. Indeed, we see this behavior from the function derived above, but NOT in the simulation. Can someone please explain this discrepancy? |