Phantom of RAMSES (POR): A new Milgromian dynamics N-body code

Since its first formulation in 1983, Milgromian dynamics (MOND) has been very successful in predicting the gravitational potential of galaxies from the distribution of baryons alone, including general scaling relations and detailed rotation curves of large statistical samples of individual galaxies covering a large range of masses and sizes. Most predictions, however, rely on static models, and only a handful of N-body codes have been developed over the years to investigate the consequences of the Milgromian framework for the dynamics of complex evolving dynamical systems. In this work, we present a new Milgromian N-body code, which is a customized version of the RAMSES code (R. Teyssier. Astron. Astrophys. 385, 337 (2002). doi: 10.1051/0004-6361:20011817) and thus comes with all its features: it includes particles and gas dynamics, and importantly allows for high spatial resolution of complex systems owing to the adaptive mesh refinement technique. It further allows the direct comparison between Milgromian simulations and standard Newtonian simulations with dark matter particles. We provide basic tests of this customized code and demonstrate its performance by presenting N-body computations of dark-matter-free spherical equilibrium models as well as dark-matter-free disk galaxies in Milgromian dynamics.