MOG in PySCo

We present the first computational simulations of John Moffat’s Modified Gravity (MOG) theory. MOG is a promising post-Einsteinean gravity theory, already having passed validation checks in several regimes where other models have failed. We compute the N-body power spectrum for a popular manifestation of MOG using PySCo, a high performance Python-based multigrid solver, and we discuss the implications for the Hubble tension. The addition of MOG to PySCo opens the door to computationally comparing scalar-vector-tensor theories against simulated data from Newtonian gravity, f(R) gravity, and other approaches.