Simulating disk galaxies and interactions in Milgromian dynamics

Since its publication 1983, Milgromian dynamics (aka MOND) has been very successful in modeling the gravitational potential of galaxies from baryonic matter alone. However, the dynamical modeling has long been an unsolved issue. In particular, the setup of a stable galaxy for Milgromian N-body calculations has been a major challenge. Here, we show a way to set up disc galaxies in MOND for calculations in the PHANTOM OF RAMSES (PoR) code by L\"ughausen (2015) and Teyssier (2002). The method is done by solving the QUMOND Poisson equations based on a baryonic and a phantom dark matter component. The resulting galaxy models are stable after a brief settling period for a large mass and size range. Simulations of single galaxies as well as colliding galaxies are shown.