A way to probe alternative theories of gravitation is to study if they could account for the structures of the universe. We then modified the well-known Gadget-2 code to probe alternative theories of gravitation through galactic dynamics. As an application, we simulate the evolution of spiral galaxies to probe alternative theories of gravitation whose weak field limits have a Yukawa-like gravitational potential. These simulations show that galactic dynamics can be used to constrain the parameters associated with alternative theories of gravitation. It is worth stressing that the recipe given in the present study can be applied to any other alternative theory of gravitation in which the superposition principle is valid.

In the standard scenario of the Newtonian gravity, a late-type galaxy (i.e., a spiral galaxy) is well described by a disk and a bulge embedded in a halo mainly composed by dark matter. In Brownstein-Moffat gravity, there is a claim that late-type galaxy systems would not need to have halos, avoiding as a result the dark matter problem, i.e., a modified gravity (non-Newtonian) would account for the galactic structure with no need of dark matter. In the present paper, we probe this claim via numerical simulations. Instead of using a "static galaxy," where the centrifugal equilibrium is usually adopted, we probe the Brownstein-Moffat gravity dynamically via numerical $N$-body simulations.

Since the Newtonian gravitation is largely used to model with success the structures of the universe, such as galaxies and clusters of galaxies, for example, a way to probe and constrain alternative theories, in the weak field limit, is to apply them to model the structures of the universe. We then modified the well known Gadget-2 code to probe alternative theories of gravitation through galactic dynamics. In particular, we modified the Gadget-2 code to probe alternatives theories whose weak field limits have a Yukawa-like gravitational potential. As a first application of this modified Gadget-2 code we simulate the evolution of elliptical galaxies. These simulations show that galactic dynamics can be used to constrain the parameters associated with alternative theories of gravitation.

In the weak field limit general relativity reduces, as is well known, to the Newtonian gravitation. Alternative theories of gravity, however, do not necessarily reduce to Newtonian gravitation; some of them, for example, reduce to Yukawa-like potentials instead of the Newtonian potential. Since the Newtonian gravitation is largely used to model with success the structures of the universe, such as for example galaxies and clusters of galaxies, a way to probe and constrain alternative theories, in the weak field limit, is to apply them to model the structures of the universe. In the present study, we consider how to probe Yukawa-like potentials using N-body numerical simulations.