The evolution of disk galaxies in modified gravity is studied by using high-resolution N-body simulations. More specifically, we use the weak field limit of two modified gravity theories, i.e., nonlocal gravity (NLG) and scalar-tensor-vector gravity known as MOG, and ignore the existence of dark matter halo. On the other hand, we construct the same models in the standard dark matter model and compare their dynamics with the galactic models in modified gravity. It turns out that there are serious differences between galactic models in these different viewpoints. For example, we explicitly show that the galactic models in modified gravity, host faster bars compared to the dark matter case. On the other hand, final stellar bars are weaker in modified gravity. These facts are not new and have already been reported in our previous simulations for exponential galactic models. Therefore, our main purpose in this study is to show that the above-mentioned differences, with emphasis on the speed of the bars, are independent of the initial density profile of the adopted disk/halo. To do so, we employ different profiles for the disk and halo and show that the results remain qualitatively independent of the initial galactic models. Moreover, a more accurate method has been used to quantify the kinematic properties of the stellar bar. Our results imply that contrary to the dark matter models, bars in modified gravity are fast rotators which never leave the fast-bar region until the end of the simulation.

The local stability of stellar and fluid discs, under a new modified dynamical model, is surveyed by using WKB approximation. The exact form of the modified Toomre criterion is derived for both types of systems and it is shown that the new model is, in all situations, more locally stable than Newtonian model. In addition, it has been proved that the central surface density of the galaxies plays an important role in the local stability in the sense that LSB galaxies are more stable than HSBs. Furthermore, the growth rate in the new model is found to be lower than the Newtonian one. We found that, according to this model, the local instability is related to the ratio of surface density of the disc to a critical surface density $\Sigma^{crit}$. We provide observational evidence to support this result based on star formation rate in HSBs and LSBs.

We study the global stability of a self-gravitating disk in the context of Modified Gravity (MOG) using N-body simulations. This theory is a relativistic scalar-tensor-vector theory of gravity and presented to address the dark matter problem. In the weak field limit MOG possesses two free parameters $\alpha$ and $\mu_0$ which have been already determined using rotation curve data of spiral galaxies. The evolution of a stellar self-gravitating disk and more specifically the bar instability in MOG is investigated and compared to a Newtonian case. Our models have exponential and Mestel-like surface densities as $\Sigma\propto \exp(-r/h)$ and $\Sigma\propto 1/r$. It is found out that, surprisingly, the disks are more stable against the bar mode in MOG than in Newtonian gravity. In other words, the bar growth rate is effectively slower than the Newtonian disks. Also we show that both free parameters, i.e. $\alpha$ and $\mu_0$, have stabilizing effects. In other words, increase in these parameters will decrease the bar growth rate.