Fornax dwarf spheroidal in MOND: Its formation and the survival of its globular clusters

The Fornax dwarf spheroidal galaxy has five massive globular clusters (GCs). They are often used for testing different dark matter and modified gravity theories, because it is difficult to reconcile their old stellar ages with the short time they need to settle in the center of the galaxy due to dynamical friction. Using high resolution N-body simulations with the Phantom of Ramses code, we investigate whether the GCs of Fornax can be reconciled with the modified Newtonian dynamics (MOND), namely its QUMOND formulation. Observational data interpreted in MOND indicate that Fornax is a tidal dwarf galaxy formed at redshift z = 0.9 in a flyby of the Milky Way (MW) and Andromeda galaxies, and that its GCs were initially massive star clusters in the disk of the MW. This helps us to set up and interpret the simulations. In the simulations, a point-mass GC orbits Fornax, and they both orbit the MW. When we ran multiple simulations with varying initial conditions for the GC, we found a 20% probability of Fornax being observed with five unsunk GCs. The unsunk GCs have the observed radial distribution. Moreover, we found: 1) In MOND, Fornax has an orbit around the MW such that the pericenters coincide with the observed peaks in the star formation history of Fornax; 2) The simulations reproduce the observed “diffuse stellar halo” of Fornax; 3) The simulations predict that Fornax has a stellar stream, which could be detectable in the existing data. 4) An extra simulation shows that if Fornax was initially a rotating disky tidal dwarf galaxy, the gravitational influence of the MW would be able to transform it into a nonrotating spheroidal. 5) Sometimes Phantom of Ramses does not conserve angular momentum. This makes the GC sink too fast if it is simulated as an N-body object.